Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913909,4000CA71.2% 96.7 / 50.0 66M550I (ATG→ATT) poxBpyruvate dehydrogenase, thiamine triphosphate‑binding, FAD‑binding
Reads supporting (aligned to +/- strand):  ref base C (10/9);  new base A (20/27);  total (30/36)
Fisher's exact test for biased strand distribution p-value = 5.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GTCGGGTAACGGTATCACTGCGTAAATCAATCATGGCATGTCCTTATTATGACGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCGCCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGTTGCAGGGCTTCATCAACTTCAGACGCTTTT  >  NC_000913/909261‑909544
                                                                                                                                           |                                                                                                                                                
gTCGGGTAACGGTATCACTGCGTAAATCAATCATGGCATGTCCTTATTATGACGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCGCCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGAATATACAg                                                                                                                                          <  2:188776/148‑1 (MQ=255)
gTCGGGTAACGGTATCACTGCGTAAATCAATCATGGCATGTCCTTATTATGACGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCACCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAgg                                                                                                                                         >  7:34313/1‑149 (MQ=255)
     gTAACGGTATCACTGCGTAAATCAATCATGGCATGTCCTTATTATGACGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCGCCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGAATATACAGGCTGa                                                                                                                                     >  7:271811/1‑148 (MQ=255)
         cGGTATCACTGCGTAAATCAATCATGGCATGTCCTTATTATGACGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCGCCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGAATATACAGGCTGAAACCt                                                                                                                                <  7:217807/149‑1 (MQ=255)
           gTATCACTGCGTAAATCAATCATGGCATGTCCTTATTATGACGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCGCCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGAATATACAGGCTGAACCCttt                                                                                                                              >  2:408111/1‑149 (MQ=255)
           gTATCACTGCGTAAATCAATCATGGCATGTCCTTATTATGACGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCGCCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGAATATACAGGCTGAAACCttt                                                                                                                              <  2:315536/149‑1 (MQ=255)
            tATCACTGCGTAAATCAATCATGGCATGTCCTTATTCTGACGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCGCCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGAATATACAGGCTGAAACCTTTg                                                                                                                             <  5:159405/149‑1 (MQ=255)
            tATCACTGCGTAAATCAATCATGGCATGTCCTTATTATGACGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCGCCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGAATATACAGGCTGAAACCTTTg                                                                                                                             <  1:207298/149‑1 (MQ=255)
            tATCACTGCGTAAATCAATCATGGCATGTCCTTATTATGACGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCGCCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGAATATACAGGCTGAAACCTTTg                                                                                                                             <  1:213423/149‑1 (MQ=255)
            tATCACTGCGTAAATCAATCATGGCATGTCCTTATTATGACGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCACCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTg                                                                                                                             <  7:360313/149‑1 (MQ=255)
                  tGCGTAAATCAATCATGGCATGTCCTTAATATGACGGGAAATGCCACCATTTTTACCTTAGCCAGTATCTTTTCGCCAGATCGATCACTTCATCACCGCGACCGCAGAAGATAGCGCGCAGAATATACAGGCTGAAACCTTTGGCCTGt                                                                                                                       >  2:318408/1‑149 (MQ=255)
                      tAAATCAATCATGGCCTGTCCTTATTATGACGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCGCCAGTTCGATCACTTCATCACCGCGTCCTCTGATGAATGCGCGCAGAATATACAGGCTGAATCCTTTGGCCTGTTCGa                                                                                                                   <  6:10866/149‑1 (MQ=255)
                       aaaTCAATCATGGCATGTCCTTATTATGACGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCGCCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGAATATACAGGCTGAAACCTTTGGCGTGTTCGAg                                                                                                                  >  2:338303/1‑149 (MQ=255)
                       aaaTCAATCATGGCATGTCCTTATTATGACGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCGCCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGAATATACAGGCTGAAACCTTTGGCCTGTTCGAg                                                                                                                  >  4:160430/1‑149 (MQ=255)
                         atcaatcaTGGCATGTCCTTATTATGACGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCGCCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGAATATACAGGCTGAAACCTTTGGCCTGTTCGAGtt                                                                                                                >  1:357892/1‑149 (MQ=255)
                         atcaatcaTGGCATGTCCTTATTATGACGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCACCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGAGtt                                                                                                                >  2:346961/1‑149 (MQ=255)
                          tcaatcaTGGCATGTCCTTATTATGACGGGAAATGCCCCCCTTTTTACCTTAGCCAGTTTGTTTTCGCCAGTTCGATCACTTCATCACCGCGTCCGCTGATGAGTGCGCGCAGAATATACAGGCTTAAACCTTTGGCCTGTGCGAGttt                                                                                                               >  2:287285/1‑149 (MQ=255)
                             atcaTGGCATGTCCTTATTATGACGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCACCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGAGTTTGAt                                                                                                            <  1:346961/149‑1 (MQ=255)
                                     aTGTCCTTATTATGACGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCGCAAGTTAGATCACTTCATCACCGCTTCCGCTGATGATTGCGCGCAGAATATACAGGCTTAAACCTTTGGCCTGTTCGAGTTTGATCTGCGGTg                                                                                                    >  3:247993/1‑149 (MQ=255)
                                      tGTCCTTATTATGACGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCACCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGAGTTTGATCTGCGGTgg                                                                                                   <  7:297906/149‑1 (MQ=255)
                                        tCCTTATTATGACGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCACCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGAGTTTGATCTGCGGTGGaa                                                                                                 <  5:253790/149‑1 (MQ=255)
                                          cTTATTATGACTGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCGCCAGTTCGATCACTTCATCACCGCGTCCGCTGATGTTTGCGCGCAGAAAATACAGGCTGAAACCTTTGGCCTGTTCGAGTTTGATCTGCGGTGGAAt                                                                                                <  1:318408/148‑1 (MQ=255)
                                          cTTATTATGACGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCGCCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGAATATACAGGCTGAAACCTTTGGCCTGTTCGAGTTTGATCTGCGGTGGAAt                                                                                                <  8:271811/148‑1 (MQ=255)
                                          cTTATTATGACGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCGCCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGAATATACAGGCTGAAACCTTTGGCCTGTTCGAGTTTGATCTGCGGTGGAAt                                                                                                <  1:266929/148‑1 (MQ=255)
                                            tattatGACGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCGCCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGAATATACAGGCTGAAACCTTTGGCCTGTTCGAGTTTGATCTGCGGTGGAATGGc                                                                                             >  2:198152/1‑149 (MQ=255)
                                              ttatGACGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCGCCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGAATATACAGGCTGAAACCTTTGGCCTGTTCGAGTTTGATCTGCGGTGGAATGGCt                                                                                            >  1:393545/1‑148 (MQ=255)
                                                         caggCCACCCTTTTTCCCCTCCCCAGGATGTCTTCGCCCCGTCGAGCACTTCATCCCCGCGTCCGCTGATGATTGCGCGCAGAATATACCGGCTGAAACCTTTCGCCTGTTCGAGGTTGACCTGCGGTGGAATGGCTAACTCTTCTTTg                                                                                <  8:294204/146‑1 (MQ=255)
                                                         aaTGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCGCCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGAATATACAGGCTGAAACCTTTGGCCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTg                                                                                <  1:408111/149‑1 (MQ=255)
                                                          aTGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCGCCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGAATATACAGGCTGAAACCTTTGGCCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTgg                                                                               >  6:36042/1‑149 (MQ=255)
                                                             ccACCCTTTTTACCTTAGCCAGTTTGTTTTCGCCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGAATATACAGGCTGAAACCTTTGGCCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCg                                                                             <  7:367162/148‑1 (MQ=255)
                                                             ccACCCTTTTTACCTTAGCCAGTTTGTTTTCACCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGa                                                                            >  7:156211/1‑149 (MQ=255)
                                                                cccTTTTTACCTTAGCCAGTTTGTTTTCGCCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGAATATACAGGCTGAAACCTTTGGCCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAActc                                                                                       <  5:171401/135‑1 (MQ=255)
                                                                cccTTTTTACCTTAGCCAGTTTGTTTTCGCCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGAATATACAGGCTGAAACCTTTGGCCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAActc                                                                                       >  6:171401/1‑135 (MQ=255)
                                                                cccTTTTTACCTTAGCCAGTTTGTTTTCGCCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGAATATACAGGCTGAAACCTTTGGCCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGacca                                                                         <  2:357892/149‑1 (MQ=255)
                                                                        aCCTTAGCCAGTTTGTTTTCGCCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGAATATACAGGCTGAAACCTTTGGCCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATc                                                                  <  3:160430/148‑1 (MQ=255)
                                                                        aCCTTAGCCAGTTTGTTTTCACCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGAGTTTGATCTGCGGTGGaa                                                                                                 <  5:305649/117‑1 (MQ=255)
                                                                        aCCTTAGCCAGTTTGTTTTCACCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGAGTTTGATCTGCGGTGGaa                                                                                                 >  6:305649/1‑117 (MQ=255)
                                                                          cTTAGCCAGTTTGTTTTCGCCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGAATATACAGGCTGAAACCTTTGGCCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATccac                                                               <  2:393545/149‑1 (MQ=255)
                                                                          cTTAGCCAGTTTGGTTTCACCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATccac                                                               >  4:339537/1‑149 (MQ=255)
                                                                           ttAGCCAGTTTGTTTTCGCCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATccac                                                               >  7:204864/1‑148 (MQ=255)
                                                                           ttAGCCAGTTTGTTTTCGCCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGAATATACAGGCTGAAACCTTTGGCCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATccac                                                               >  4:275668/1‑148 (MQ=255)
                                                                               ccAGTTTGTTTTCACCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGCCCCCCCCATCCACCAaca                                                          >  6:114505/1‑149 (MQ=255)
                                                                                  gtttgtttTCGCCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGAATATACAGGCTGAAACCTTTGGCCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCg                                                       >  2:104556/1‑149 (MQ=255)
                                                                                       ttttCGCCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGAATATACAGGCTGAAACCTTTGGCCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCaa                                                            >  7:238706/1‑139 (MQ=255)
                                                                                       ttttCGCCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGAATATACAGGCTGAAACCTTTGGCCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCaa                                                            <  8:238706/139‑1 (MQ=255)
                                                                                           cGCCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGAATATACAGGCTGAAACCTTTGGCCTGTTCGAGTTTGATCTGCGGTGGAATGGCt                                                                                            >  4:57463/1‑103 (MQ=255)
                                                                                           cGCCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGAATATACAGGCTGAAACCTTTGGCCTGTTCGAGTTTGATCTGCGGTGGAATGGCt                                                                                            <  3:57463/103‑1 (MQ=255)
                                                                                             ccAGTTCGATCACTTCCTCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATg                                            <  5:114505/149‑1 (MQ=255)
                                                                                             ccAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGAATATACAGGCTGAAACCTTTGTCCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATg                                            <  1:198152/149‑1 (MQ=255)
                                                                                              cAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATgg                                           <  8:204864/149‑1 (MQ=255)
                                                                                              cAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATg                                            >  1:293335/1‑148 (MQ=255)
                                                                                              cAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGAATATACAGGCTGAAACCTTTGGCCTGTTCGAGTTTGATCTGCGGTGGAATGGCTaa                                                                                          >  6:259227/1‑102 (MQ=255)
                                                                                              cAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGAATATACAGGCTGAAACCTTTGGCCTGTTCGAGTTTGATCTGCGGTGGAATGGCTaa                                                                                          <  5:259227/102‑1 (MQ=255)
                                                                                              cAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGAATATACAGGCTGAAACCTTTGGCCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATgg                                           <  5:36042/149‑1 (MQ=255)
                                                                                                     aTCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGAATATACAGGCTGAAACCTTTGGCCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGgc                                    >  8:38327/1‑149 (MQ=255)
                                                                                                         cTTCATCACCGCGTCCGCTGATGATTGCGCGCAGAATATACAGGCTGAAACCTTTGGCCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGa                                              <  8:182392/135‑1 (MQ=255)
                                                                                                         cTTCATCACCGCGTCCGCTGATGATTGCGCGCAGAATATACAGGCTGAAACCTTTGGCCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGa                                              >  7:182392/1‑135 (MQ=255)
                                                                                                            catcaCCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGTTg                              <  8:156211/148‑1 (MQ=255)
                                                                                                               caCCGCGTCCGCTGATGATTGCGCGCAGAATATACAGGCTGAAACCTTTGGCCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGTTGCAgg                          >  1:353810/1‑149 (MQ=255)
                                                                                                                 ccGCGTCCGCTGATGATTGCGCGCAGAATATACAGGCTGAAACCTTTGGCCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGTTGCAGGGc                        <  2:353810/149‑1 (MQ=255)
                                                                                                                   gcgTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGTTGCAGGGCtt                      >  8:19035/1‑149 (MQ=255)
                                                                                                                   gcgTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGTTGCAGGGCtt                      <  2:293335/149‑1 (MQ=255)
                                                                                                                   gcgTCCGCTGATGATTGCGCGCAGAATATACAGGCTGAAACCTTTGGCCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGTTGCAGGGCtt                      <  3:275668/149‑1 (MQ=255)
                                                                                                                    gacgcggtgatgaTTGCGCGCAGAATATACAGGCTGAAACCTTTGGCCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGTTGCAGGGCTtc                     <  8:199851/142‑1 (MQ=255)
                                                                                                                    cgTCCGCTGATGATTGCGCGCAGAATATACAGGCTGAAACCTTTGGCCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACGACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGTTGCAGGGCTtc                     <  4:28289/149‑1 (MQ=255)
                                                                                                                           tgatgaTTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGTTGCAGGGCTTCATCAACt              >  8:97735/1‑149 (MQ=255)
                                                                                                                                       gcagaaTATACAGGCTGAAACCTTTGGCCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGTTGCAGGGCTTCATCAACTTCAGACGCtttt  >  6:205374/1‑149 (MQ=255)
                                                                                                                                           |                                                                                                                                                
GTCGGGTAACGGTATCACTGCGTAAATCAATCATGGCATGTCCTTATTATGACGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCGCCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGTTGCAGGGCTTCATCAACTTCAGACGCTTTT  >  NC_000913/909261‑909544

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑