Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,788,676 A→T 6.3% intergenic (+27/‑309) ygaQ → / → csiD uncharacterized protein; putative enzyme/carbon starvation protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,788,6760AT 6.3% 351.3 / 9.3 125intergenic (+27/‑309)ygaQ/csiDuncharacterized protein; putative enzyme/carbon starvation protein
Reads supporting (aligned to +/- strand):  ref base A (48/69);  new base T (4/4);  total (52/73)
Fisher's exact test for biased strand distribution p-value = 7.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.40e-01

Only 100 of 130 total aligned reads displayed.

AATTATCCCCCCTGGTAAAGTATTTACGCTTAGAGGCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTCCCCCTAACCTGTTGCTTTAGTTATTCATTTCCTGTCTCACTTTGCCTTAAT  >  NC_000913/2788529‑2788825
                                                                                                                                                   |                                                                                                                                                     
aaTTATCCCCCCTGGTAAAGTATTTACGCTTAGAGGCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACg                                                                                                                                                     >  2:852125/1‑150 (MQ=255)
aaTTATCCCCCCTGGTAAAGTATTTACGCTTAGAGGCGGAACTCTAAATATCAATATTCCTGGGCGAAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATTATTTAATGTTATATCTGCCCCGATAAAACg                                                                                                                                                     >  2:275856/1‑150 (MQ=255)
aaTTATCACCCCTGGTAAAGTATTTACGATTAGAGGCTGAACTCTAAATATCAATATTCCTGGGAGTAGCTCTCTTCTTTTATGGAAAACAGGAGAACCGCCTAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACg                                                                                                                                                     <  1:295951/150‑1 (MQ=255)
  ttATCCCCCCTGGTAAAGTATTTACGCTTAGAGGCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACggg                                                                                                                                                   <  2:1027100/150‑1 (MQ=255)
    aTCCCCCCTGGTAAAGTATTTACGCTTAGAGGCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGc                                                                                                                                                 <  1:1174567/150‑1 (MQ=255)
      ccccccTGGTAAAGTATTTACGCTTAGAGGCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAg                                                                                                                                               <  2:156108/150‑1 (MQ=255)
       cccccTGGTAAAGTATTTACGCTTAGAGTCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACAATCTCTATTTATAATTTAATGTTAAATCTGCCCAGATAAAACCGGGCAGa                                                                                                                                              >  2:255047/1‑150 (MQ=255)
       cccccTGGTAAAGTATTTACGCTTAGAGGCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGa                                                                                                                                              <  1:1382490/150‑1 (MQ=255)
       cccccTGGTAAAGTATTTACGCTTAGAGGCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTATTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATTTTATATCTGCCCCGATAAAACGGGGCAGa                                                                                                                                              <  1:1518476/150‑1 (MQ=255)
        ccccTGGTAAAGTATTTACGCTTAGAGGCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGat                                                                                                                                             >  1:954852/1‑150 (MQ=255)
        ccccTGGTAAAGTATTTACGCTTAGAGGCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGat                                                                                                                                             <  2:1225383/150‑1 (MQ=255)
        ccccTGGTAAAGTATTTACGCTTAGAGGCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGat                                                                                                                                             >  2:141709/1‑150 (MQ=255)
             ggTAAAGTATTTACGCTTAGAGGCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAatat                                                                                                                                        >  1:195089/1‑150 (MQ=255)
               tAAAGTATTTACGCTTAGAGGCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTAAAATTTATTGTTATATCTGCCCCGTTTTATCGGGGCAGatataacat                                                                                                                                      <  1:106167/150‑7 (MQ=255)
                aaaGTATTTACGCTTAGAGGCGGAACTCTAAATATCAATGTTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATgtt                                                                                                                                     >  2:887784/1‑150 (MQ=255)
                   gTATTTACGCTTAGAGGCGGTACTATAAATATCAATAATCCTGGGCGTAGCGCTCTTATTTTAGGGAAAACAGGAGAACCGACGAACTATCTCTATTTATAAATTTATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATgtttag                                                                                                                                  <  2:350921/150‑1 (MQ=255)
                   gTATTTACGCTTAGAGGCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATgtttag                                                                                                                                  <  2:1450844/150‑1 (MQ=255)
                     aTTTACGCTTAGAGGCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAACACAGGAGAACCGCCGAACTATCTCTATTGATAATTTAATGTTATATCTGCCCCGATAAAACGGGTCAGATAATATgtttagtt                                                                                                                                >  1:133109/1‑150 (MQ=255)
                     aTTTACGCTTAGAGGCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATgtttagtt                                                                                                                                <  1:986698/150‑1 (MQ=255)
                           gCTTAGAGGCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAATACAGGAGAAACGACGAACTATCTCTATTTATAATTTAATGTTATATCTGAACCGATAAAAAGGGGCAGATAATATGTTTAGTTTACTaa                                                                                                                          >  1:707625/1‑150 (MQ=255)
                               agagGCGGAACTCTAAATATGAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACATGAGAACCTCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGt                                                                                                                      <  1:1537379/150‑1 (MQ=255)
                               agagGCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTGTAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGt                                                                                                                      >  2:1596649/1‑150 (MQ=255)
                               agagGCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTCACTAACGGt                                                                                                                      <  2:1429755/150‑1 (MQ=255)
                                  ggCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCAt                                                                                                                   >  1:1331126/1‑150 (MQ=255)
                                  ggCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCAt                                                                                                                   >  1:84171/1‑150 (MQ=255)
                                     ggAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCAtttt                                                                                                                <  1:1617766/150‑1 (MQ=255)
                                       aaCTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTCGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGc                                                                                                              <  2:1518572/150‑1 (MQ=255)
                                       aaCTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTGTAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGc                                                                                                              <  1:1596649/150‑1 (MQ=255)
                                       aaCTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGc                                                                                                              <  2:804805/150‑1 (MQ=255)
                                        aCTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCa                                                                                                             <  1:673693/150‑1 (MQ=255)
                                        aCTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCa                                                                                                             <  2:782959/150‑1 (MQ=255)
                                          tctAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGt                                                                                                           >  1:761565/1‑150 (MQ=255)
                                            tAAATATCAATAATCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGa                                                                                                         <  2:1357246/150‑1 (MQ=255)
                                                tatCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGcc                                                                                                     >  1:342331/1‑150 (MQ=255)
                                                   cAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCAtt                                                                                                  >  2:1603740/1‑150 (MQ=255)
                                                   cAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGACGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCAtt                                                                                                  >  2:1672271/1‑150 (MQ=255)
                                                    aaTATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCAttt                                                                                                 <  1:758438/150‑1 (MQ=255)
                                                    aaTATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCAttt                                                                                                 <  2:1073967/150‑1 (MQ=255)
                                                      tatTCCTGGGCTTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTAc                                                                                               <  2:1375646/150‑1 (MQ=255)
                                                      tatTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTAc                                                                                               <  2:1702526/150‑1 (MQ=255)
                                                      tatTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTAc                                                                                               <  1:1629272/150‑1 (MQ=255)
                                                            tGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGttttt                                                                                        <  2:838394/150‑1 (MQ=255)
                                                             gggCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAAGCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGttttt                                                                                        >  1:1087566/1‑150 (MQ=255)
                                                             gggCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGttttt                                                                                        >  2:262532/1‑150 (MQ=255)
                                                             gggCGTAGCGCTCTTCATTTAGCGAAAACAGCATAACCGCCGAACTATCTCTATTTATAATTTAATGTTAAATCTGCCCCGATAAAACGCGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGttttt                                                                                        <  2:170608/150‑1 (MQ=255)
                                                                  tatctctctTCTTTTAGGGATAACTGGAGAACCGACGAACTATCTCTAATTATAATTTAATGTTAAATCAGCCCCGATAAAACGGGGCAGATATTATGTTTAGTTTACTAACTGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCg                                                                                   <  2:677663/145‑1 (MQ=255)
                                                                  actcgcTCTTCTTTTAGGGAAAACAGTAGAACCGCCGACCTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCg                                                                                   <  2:137119/147‑1 (MQ=255)
                                                                      gcTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCa                                                                               >  2:1656612/1‑150 (MQ=255)
                                                                       ctctTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAg                                                                              <  1:379198/150‑1 (MQ=255)
                                                                        tcttctTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGa                                                                             >  1:1624627/1‑150 (MQ=255)
                                                                            cttTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGTTAAt                                                                         >  2:841294/1‑150 (MQ=255)
                                                                            cttTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAAt                                                                         >  2:308670/1‑150 (MQ=255)
                                                                             ttttAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAAta                                                                        >  1:416550/1‑149 (MQ=255)
                                                                                 tgggAAAACAGGAGAACCGCCGAACTATCTATATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGtt                                                                    <  2:1450275/149‑1 (MQ=255)
                                                                                 aGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGtt                                                                    <  2:1729552/150‑1 (MQ=255)
                                                                                  gggAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTtc                                                                   <  2:615263/150‑1 (MQ=255)
                                                                                  gggAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTtc                                                                   <  1:685486/150‑1 (MQ=255)
                                                                                  gggAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTtc                                                                   <  2:682286/150‑1 (MQ=255)
                                                                                  gggAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTtc                                                                   <  2:1121205/150‑1 (MQ=255)
                                                                                   ggAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTtct                                                                  >  1:931748/1‑150 (MQ=255)
                                                                                    gAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTtcta                                                                 >  2:819908/1‑149 (MQ=255)
                                                                                       aaCAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATACAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTGTCTaa                                                              >  2:1472024/1‑150 (MQ=255)
                                                                                        aCAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAAt                                                             >  2:606891/1‑150 (MQ=255)
                                                                                          aGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGt                                                           <  2:719428/150‑1 (MQ=255)
                                                                                           ggAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACCGTTTTTTATCGACCAGATAATCTGTTCTCTAATGtt                                                          <  2:1483795/150‑1 (MQ=255)
                                                                                             agaACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTaa                                                        <  1:1611588/150‑1 (MQ=255)
                                                                                                  cgccgAACTATCTCTATTTATAATTTAAAGTTATATCTGCCCCGCTAAAACGGGGCATATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTccc                                                   <  2:761421/150‑1 (MQ=255)
                                                                                                          tATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTCCCCCTAACCt                                           >  1:743120/1‑150 (MQ=255)
                                                                                                          tATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTCCCACTAACCt                                           >  1:301400/1‑150 (MQ=255)
                                                                                                            tatctATATATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTAATGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTCCCCCTAACCTGt                                         <  2:68725/148‑1 (MQ=255)
                                                                                                             ctctATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTCCCCCTAACCTGtt                                        <  1:1689475/150‑1 (MQ=255)
                                                                                                              tctATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTAATAAAGTTAATTTTGCAGTGAAGACATTTACTGTTTTTTATCGACCAGATAATCTTTTATCTAATGTTAACTCACCCTAACCTGttt                                       >  1:72576/1‑149 (MQ=255)
                                                                                                               ctATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTCCCCCTAACCTGTTGc                                      >  2:10877/1‑150 (MQ=255)
                                                                                                               ctATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTCCCCCTAACCTGTTGc                                      >  2:465663/1‑150 (MQ=255)
                                                                                                               ctATTTATAATTTAATGGTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCATTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTCCCCCTAACCTGTTGc                                      <  2:1085723/150‑1 (MQ=255)
                                                                                                                  tttATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTCCCCCTAACCTGTTGCttt                                   <  2:342331/150‑1 (MQ=255)
                                                                                                                   ttATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTATCGGTCATTTTGCCGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTCCCCCTAACCTGTTGCTTTa                                  <  2:1104840/150‑1 (MQ=255)
                                                                                                                   ttATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTCCCCCTAACCTGTTGCTTTa                                  >  1:1397468/1‑150 (MQ=255)
                                                                                                                     ataATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTCCCCCTAACCTGTTGCTTTAGt                                <  2:587261/150‑1 (MQ=255)
                                                                                                                     ataATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCAGTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTCCCCCTAACCTGTTGCTTTAGt                                >  1:266414/1‑150 (MQ=255)
                                                                                                                      taATTTTATGTTTTATCTGGCCCGATAATAAGGGGCAGATAATATGTTTAGTTTATTTACGGTCATTTTGCAGTGTAGCGATTTACTGTTTTTTATCTACCAGATAATATGTTCTCTAATGTTAACTCCCCCTAACCTGTTGCTTTAGtt                               <  1:451864/150‑1 (MQ=255)
                                                                                                                      taATTTAATGTTATATCTGCCCTGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGGCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTCCCCCTAACCTGTTGCTTTAGtt                               <  1:346090/150‑1 (MQ=255)
                                                                                                                       aaTTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACTGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATGAGTTCTCTAATGTTAACTACCCCTAACCTGTTGCTTTAGTTa                              >  1:1127809/1‑150 (MQ=255)
                                                                                                                              tGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTCCCCCTAACCTGTTGCTTTAGTTATTCAttt                       >  1:1426684/1‑150 (MQ=255)
                                                                                                                                ttATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTATTTATCGACCAGATAATCTGTTCTCTAATGTTAACTCCACCTAACCTGTTGCTTTAGTTATTCATTTcc                     <  2:145343/150‑1 (MQ=255)
                                                                                                                                ttATATCTGCCCCGATAAAAAGGGTCAGATAATATGTTTAGTTTACTAACGTTAATTTTGCAGTGAAGCCATTTACTTTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTCCCCCTAACCTGTTGCTTTAGTTATTCATTTcc                     <  2:707625/150‑1 (MQ=255)
                                                                                                                                  ttatCTGCCCCGTTTTATCGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTCCCCCTAACCTGTTGCTTTAGTTATTCATTTCCTg                   >  1:966851/2‑150 (MQ=255)
                                                                                                                                  ttatCTGCCCCGTTTTATCGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTCCCCCTAACCTGTTGCTTTAGTTATTCATTTCCTg                   >  2:100827/2‑150 (MQ=255)
                                                                                                                                  atatCAGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTCCCCCTAACCTGTTGCTTTAGTTATTCATTTCCTg                   <  1:852125/150‑1 (MQ=255)
                                                                                                                                    atCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTCCCCCTAACCTGTTGCTTTAGTTATTCATTTCCTGtc                 <  2:1331126/150‑1 (MQ=255)
                                                                                                                                    atCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTCCCCCTAACCTGTTGCTTTAGTTATTCATTTCCTGtc                 >  2:1455664/1‑150 (MQ=255)
                                                                                                                                    atCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTAATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATTATCTGTTCTCTAATGTTAACTCCCCCTAACCTGTTGCTTTAGTTATTCATTTCCTGtc                 >  2:928121/1‑150 (MQ=255)
                                                                                                                                      cTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTCCACATAACCTGTGGCTTTAGTTATTCATTTCCTGtcta               >  2:835631/1‑149 (MQ=255)
                                                                                                                                             gATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTCCCCCTAACCTGTTGCTTTAGTTATTCATTTCCTGTCTCACTTTGc        <  2:133277/150‑1 (MQ=255)
                                                                                                                                              aTAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTCCCCCTAACCTGTTGCTTTAGTTATTCATTTCCTGTCTCACTTTGcc       <  2:1647681/150‑1 (MQ=255)
                                                                                                                                              aTAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTCCCCCTAACCTGTTGCTTTAGTTATTCATTTCCTGTCTCACTTTGcc       <  1:991184/150‑1 (MQ=255)
                                                                                                                                                 aaaCGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTCCCCCTAACCTGTTGCTTTATTTATTAATATCCTGTCTAACTTTGCCTTa    >  1:845868/1‑150 (MQ=255)
                                                                                                                                                 aaaCGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTCCCCCTAACCTGTTGCTTTAGTTATTCATTTCCTGTCTCACTTTGCCTTa    <  2:1693658/150‑1 (MQ=255)
                                                                                                                                                 aaaCGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTCCCCCTAACCTGTTGCTTTAGTTATTCATTTCCTGTCTCACTTTGCCTTa    <  1:1363659/150‑1 (MQ=255)
                                                                                                                                                   aCGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTCCCCCTAACCTGTTGCTTTAGTTATTCATTTCCTGTCTCACTTTGCCTTAAt  <  2:680848/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
AATTATCCCCCCTGGTAAAGTATTTACGCTTAGAGGCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTCCCCCTAACCTGTTGCTTTAGTTATTCATTTCCTGTCTCACTTTGCCTTAAT  >  NC_000913/2788529‑2788825

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