Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,788,674 A→T 9.7% intergenic (+25/‑311) ygaQ → / → csiD putative uncharacterized protein YgaQ/PF08943 family protein CsiD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,788,6740AT 9.7% 481.6 / 28.7 163intergenic (+25/‑311)ygaQ/csiDputative uncharacterized protein YgaQ/PF08943 family protein CsiD
Reads supporting (aligned to +/- strand):  ref base A (81/66);  new base T (7/9);  total (88/75)
Fisher's exact test for biased strand distribution p-value = 4.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 173 total aligned reads displayed.

GAAAAATTATCCCCCCTGGTAAAGTATTTACGCTTAGAGGCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTCCCCCTAACCTGTTGCTTTAGTTATTCATTTCCTGTCTCACTTTGCCTTA  >  NC_000913/2788525‑2788823
                                                                                                                                                     |                                                                                                                                                      
gAAAAATTATCCCCCCTGGTAAAGTATTTACGCTTAGAGGCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATaa                                                                                                                                                        <  1:1461739/150‑1 (MQ=255)
 aaaagTTATCCCCCCTGGTAAAGTATTTACGCTTAGAGGCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATaaa                                                                                                                                                       <  2:1709567/150‑1 (MQ=255)
 aaaaaTTATCCCCCCTGGTAAAGTATTTACGCTTAGAGGCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATaaa                                                                                                                                                       >  1:75879/1‑150 (MQ=255)
     aTTATCCCCCCTGGTAAAGTATTTACGCTTAGAGGCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACgg                                                                                                                                                   >  2:1578224/1‑150 (MQ=255)
      ttATCCCCCCTGGTAAAGTATTTACGCTTAGAGGCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACggg                                                                                                                                                  <  1:587991/150‑1 (MQ=255)
      ttATCCCCCCTGGTAAAGTATTTACGCTTAGAGGCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACggg                                                                                                                                                  <  2:384928/150‑1 (MQ=255)
      ttATCCCCCCTGGTAAAGTATTTACGCTTAGAGGCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACggg                                                                                                                                                  <  2:1383829/150‑1 (MQ=255)
        aTCCCCCCTGGTAAAGTATTTACGCTTAGAGGCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGc                                                                                                                                                <  1:478962/150‑1 (MQ=255)
        aTCCCCCCTGGTAAAGTATTTACGCTTAGAGGCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGc                                                                                                                                                <  2:653668/150‑1 (MQ=255)
        aTCCCCCCTGGTAAAGTATTTACGCTTAGAGGCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGc                                                                                                                                                <  1:662809/150‑1 (MQ=255)
         tCCCCCCTGGTAAAGTATTTACGCTTAGAGGCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGTTTTATCGGGGCa                                                                                                                                               <  2:1614571/150‑1 (MQ=255)
          ccccccTGGTATAGTTTTTACGCTTAGAGGCGGAACTCTAAATATCAATATTCCTGTGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATTTTATATCTGCCCCGATAAAACGGGGCAg                                                                                                                                              <  1:1337289/150‑1 (MQ=255)
          ccccccTGGTAAAGTATTTACGCTTAGAGGCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGTTTTATCGGGGCAg                                                                                                                                              <  1:1442393/150‑1 (MQ=255)
          ccccccTGGTAAAGTATTTACGCTTAGAGGCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAg                                                                                                                                              <  1:1677642/150‑1 (MQ=255)
           cccccTGGTAAAGTATTTACGCTTAGAGGCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGTTTTATCGGGGCAg                                                                                                                                              >  1:1424433/1‑148 (MQ=255)
            ccccTGGTAAAGTATTTACGCTTAGAGGCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTCTAGGGAAAACAGGAGAACCGCAGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGat                                                                                                                                            >  1:1226951/1‑150 (MQ=255)
                  gTAAAGTATTTACGCTTAGAGGCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATg                                                                                                                                      <  1:1116751/150‑1 (MQ=255)
                   tAAAGTATTTACGCTTAGAGGCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATgt                                                                                                                                     >  2:1877816/1‑150 (MQ=255)
                     aaGTATTTACGCTTAGAGGCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATgttt                                                                                                                                   <  2:67263/150‑1 (MQ=255)
                       gTATTTACGCTTAGAGGCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATgtttag                                                                                                                                 >  1:1255977/1‑150 (MQ=255)
                       gTATTTACGCTTAGAGGCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATgtttag                                                                                                                                 <  2:85548/150‑1 (MQ=255)
                        tATTTACGCTTAGAGGCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATgtttagt                                                                                                                                >  2:1465888/1‑150 (MQ=255)
                          tttACGCTTAGAGGCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCTCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATgtttagttt                                                                                                                              <  2:981838/150‑1 (MQ=255)
                          tttACGCTTAGAGGCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATgtttagttt                                                                                                                              >  2:851265/1‑150 (MQ=255)
                             aCGCTTAGAGGCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACt                                                                                                                           >  2:1409493/1‑150 (MQ=255)
                               gCTTAGAGGCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTaa                                                                                                                         >  2:1738633/1‑150 (MQ=255)
                               gCTTAGAGGCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCAGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTaa                                                                                                                         >  1:439492/1‑150 (MQ=255)
                                   agagGCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGt                                                                                                                     >  1:586738/1‑150 (MQ=255)
                                    gagGCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTc                                                                                                                    >  2:1438317/1‑150 (MQ=255)
                                      ggCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCAt                                                                                                                  >  1:29142/1‑150 (MQ=255)
                                         ggAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGACCTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCAtttt                                                                                                               <  2:1255977/150‑1 (MQ=255)
                                          gAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTg                                                                                                              >  2:182707/1‑150 (MQ=255)
                                               ctAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTg                                                                                                         <  2:1281451/150‑1 (MQ=255)
                                               ctAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTg                                                                                                         >  2:365202/1‑150 (MQ=255)
                                                tAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGa                                                                                                        <  1:866643/150‑1 (MQ=255)
                                                tAAATATCAATATTCATGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCTATAAAACGGGGCAGATAATATGTTTATTTTAATAACGGTCATTTTGCAGTGa                                                                                                        >  2:1820422/1‑150 (MQ=255)
                                                 aaaTATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGaa                                                                                                       >  2:213531/1‑150 (MQ=255)
                                                  aaTATCAATATTCCTGGGCGTATCGCTCTTCTTTTATGGAAAACAGGAGAACCGACGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGATTACTAACGGTCATTTTGCAGTGAAg                                                                                                      <  2:720333/150‑1 (MQ=255)
                                                     atCAATATTCCTGGGCGTAGCGCTCTTCTTTTATGGAATACAGCAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCa                                                                                                   <  2:1591722/150‑1 (MQ=255)
                                                      tCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCAt                                                                                                  >  2:204809/1‑150 (MQ=255)
                                                        aaTATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCAttt                                                                                                >  1:1573742/1‑150 (MQ=255)
                                                                tGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGtttt                                                                                       <  2:449439/150‑1 (MQ=255)
                                                                 gggCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGttttt                                                                                      <  1:1584576/150‑1 (MQ=255)
                                                                  ggCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACtttttttt                                                                                     >  2:1117722/1‑150 (MQ=255)
                                                                     gTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATc                                                                                  >  2:740691/1‑150 (MQ=255)
                                                                     gTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATc                                                                                  <  1:430503/150‑1 (MQ=255)
                                                                      tAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCg                                                                                 <  2:104913/150‑1 (MQ=255)
                                                                         cgcTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGAcc                                                                              <  2:74771/150‑1 (MQ=255)
                                                                          gcTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCa                                                                             >  1:1637500/1‑150 (MQ=255)
                                                                            tcttctTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCAttt                                                                                                                >  1:1701961/1‑114 (MQ=255)
                                                                               tcttTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATaa                                                                        >  1:1091309/1‑150 (MQ=255)
                                                                                  tttAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGtttt                                                                                       <  2:473472/132‑1 (MQ=255)
                                                                                    tAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGt                                                                   >  2:554536/1‑150 (MQ=255)
                                                                                      gggAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTtc                                                                 >  2:877036/1‑150 (MQ=255)
                                                                                      gggAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTtc                                                                 >  1:1854299/1‑150 (MQ=255)
                                                                                        taaaaCAGGAGAACCGCCTACCTATATCTATTTATAATTTAATGTTATATCTGCCCGGTTAAAACTGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTtctc                                                               <  2:1714170/149‑1 (MQ=255)
                                                                                         aaaaCAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTctct                                                              <  1:851265/150‑1 (MQ=255)
                                                                                             cAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTTATCGACCAGATAATCTGTTCTCTAAt                                                           >  1:1474458/1‑150 (MQ=255)
                                                                                                  gaACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCTATTAAACGGTGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAAc                                                     <  1:519610/150‑1 (MQ=255)
                                                                                                   aaCCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACt                                                    <  1:175910/150‑1 (MQ=255)
                                                                                                    aCCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTc                                                   >  2:1704650/1‑150 (MQ=255)
                                                                                                          caaCTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTCCCCCTa                                             <  1:1021578/149‑1 (MQ=255)
                                                                                                             cTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTCCCCCTAAcc                                          <  1:1607034/150‑1 (MQ=255)
                                                                                                             cTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTCCCCCTAAcc                                          <  2:1909609/150‑1 (MQ=255)
                                                                                                              tATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTCCCCCTAACCt                                         >  1:678275/1‑150 (MQ=255)
                                                                                                                 ctctATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTCCCCCTAACCTGtt                                      >  2:109412/1‑150 (MQ=255)
                                                                                                                      tttATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTCCCCCTAACCTGTTGCttt                                 <  2:1241390/150‑1 (MQ=255)
                                                                                                                      tttATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTCCCCCTAACCTGTTGCttt                                 >  2:1770487/1‑150 (MQ=255)
                                                                                                                        tataATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTCCCCCTAACCTGTTGCTTTAg                               >  2:249168/1‑150 (MQ=255)
                                                                                                                        tataATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTCCCCCTAACCTGTTGCTTTAg                               >  2:559554/1‑150 (MQ=255)
                                                                                                                        tataATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTCCCCCTAACCTGTTGCTTTAg                               >  2:1992041/1‑150 (MQ=255)
                                                                                                                           aaTTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTCCCCCTAACCTGTTGCTTTAGTTa                            >  2:1963324/1‑150 (MQ=255)
                                                                                                                            aTTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATTTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTCCCCCTAACCTGTTGCTTTAGTTAt                           <  1:1820422/150‑1 (MQ=255)
                                                                                                                            aTTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTCCCCCTAACCTGTTGCTTTAGTTAt                           <  2:75879/150‑1 (MQ=255)
                                                                                                                            aTTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTCCCCCTAACCTGTTGCTTTAGTTAt                           <  1:237589/150‑1 (MQ=255)
                                                                                                                               tAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTCCCCCTAACCTGTTGCTTTAGTTATTCa                        >  1:442019/1‑150 (MQ=255)
                                                                                                                                 aTGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTCCCCCTAACCTGTTGCTTTAGTTATTCAtt                      >  1:1020354/1‑150 (MQ=255)
                                                                                                                                  tGTTATATCTGCCCCGATAAAACGTGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTCCCCCTAACCTGTTGCTTTAGTTATTCAttt                     <  1:1769648/150‑1 (MQ=255)
                                                                                                                                  catattatCTGCCCCGTTTTATCGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTCCCCCTAACCTGTTGCTTTAGTTATTCAttt                     >  2:1285114/6‑150 (MQ=255)
                                                                                                                                   gTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTCCCCCTAACCTGTTGCTTTAGTTATTCATTTc                    >  2:384457/1‑150 (MQ=255)
                                                                                                                                   gTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTGCTCTAATGTTACCTCCCCCTAACCTGTTGCTTTAGTTATTCATTTc                    >  1:623819/1‑150 (MQ=255)
                                                                                                                                        atCTGCCCCGTTTTATCGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTCCCCCTAACCTGTTGCTTTAGTTATTCATTTCCTGtc               >  1:1123803/1‑150 (MQ=255)
                                                                                                                                        atCTGCCCCGTTTTATCGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTCCCCCTAACCTGTTGCTTTAGTTATTCATTTCCTGtc               <  1:819638/150‑1 (MQ=255)
                                                                                                                                        atCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTCCCCCTAACCTGTTGCTTTAGTTATTCATTTCCTGtc               >  2:655166/1‑150 (MQ=255)
                                                                                                                                        atCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTCCCCCTAACCTGTTGCTTTAGTTATTCATTTCCTGtc               >  2:1684089/1‑150 (MQ=255)
                                                                                                                                        atCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTCCCCCTAACCTGTTGCTTTAGTTATTCATTTCCTGtc               <  1:201229/150‑1 (MQ=255)
                                                                                                                                         tCTGCCCCGTTTTATCGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTCCCCCTAACCTGTTGCTTTAGTTATTCATTTCCTGtct              >  1:1882937/1‑150 (MQ=255)
                                                                                                                                         tCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTa                                                       <  2:272124/109‑1 (MQ=255)
                                                                                                                                         tCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTa                                                       >  1:272124/1‑109 (MQ=255)
                                                                                                                                         tCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTCCCCCTAACCTGTTGCTTTAGTTATTCATTTCCTGtct              >  2:610423/1‑150 (MQ=255)
                                                                                                                                         tCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCCGTTCTCTAATGTTAACTCCCCCTAACCTGTTGCTTTAGTTATTCATTTCCTGtct              >  2:1659980/1‑150 (MQ=255)
                                                                                                                                          cTGCCCCGTTTTATCGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTCCCCCTAACCTGTTGCTTTAGTTATTCATTTCCTGtctc             >  1:790851/1‑150 (MQ=255)
                                                                                                                                            gCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTCCCCCTAACCTGTTGCttt                                 <  2:164588/128‑1 (MQ=255)
                                                                                                                                             ccccGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTCCCCCTAACCTGTTGCTTTAGTTATTCATTTCCTGTCTCACt          >  2:898115/1‑150 (MQ=255)
                                                                                                                                                cGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTCCCCCTAACCTGTTGCTTTAGTTATTCATTTCCTGTCTCACTTTg       <  1:1963324/150‑1 (MQ=255)
                                                                                                                                                   tAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTCCCCCTAACCTGTTGCTTTAGTTATTCATTTCCTGTCTCACTTTGCCt    <  1:343995/150‑1 (MQ=255)
                                                                                                                                                    aaaaCGGGGCAGATAATATGTTTAGTTTACTAACTGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTCCCCCTAACCTGTTGCTTTAGTTATTCATTTCCTGTCTCACTTTGCCtt   <  1:1465888/150‑1 (MQ=255)
                                                                                                                                                    aaaaCGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTCCCCCTAACCTGTTGCTTTAGTTATTCATTTCCTGTCTCACTTTGCCtt   <  2:1020354/150‑1 (MQ=255)
                                                                                                                                                     aaaCGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTCCCCCTAACCTGTTGCTTTAGTTATTCATTTCCTGTCTCACTTTGCCTTa  >  1:1936961/1‑150 (MQ=255)
                                                                                                                                                     aaaCGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTCCCCCTAACCTGTTGCTTTAGTTATTCATTTCCTGTCTCACTTTGCCTTa  <  1:1734051/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
GAAAAATTATCCCCCCTGGTAAAGTATTTACGCTTAGAGGCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTCCCCCTAACCTGTTGCTTTAGTTATTCATTTCCTGTCTCACTTTGCCTTA  >  NC_000913/2788525‑2788823

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