Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,788,674 A→T 6.3% intergenic (+25/‑311) 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,6740AT 6.3% 357.1 / 9.3 125intergenic (+25/‑311)ygaQ/csiDuncharacterized protein; putative enzyme/carbon starvation protein
Reads supporting (aligned to +/- strand):  ref base A (47/69);  new base T (4/4);  total (52/73)
Fisher's exact test for biased strand distribution p-value = 7.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.91e-01

Only 100 of 131 total aligned reads displayed.

GAAAAATTATCCCCCCTGGTAAAGTATTTACGCTTAGAGGCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTCCCCCTAACCTGTTGCTTTAGTTATTCATTTCCTGTCTCACTTTGCCTTA  >  NC_000913/2788525‑2788823
                                                                                                                                                     |                                                                                                                                                     
gAAAAATTATCCCCCCTGGTAAAGTATTTACGCTTAGAGGCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATaa                                                                                                                                                       >  1:587261/1‑150 (MQ=255)
 aaaaaTTATCCCCCCTGGGAAAGTATTTCCGCTTAGAGGCGGAACGCTAAATATCAAGATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAAACGCCGAACTATCTCTATTTATAATTTTATGTTATATCTGCCGCGATaaa                                                                                                                                                      >  2:264487/1‑150 (MQ=255)
    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)
          ccccccTGGTAAAGTATTTACGCTTAGAGGCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGTTTTATCGGGGCAg                                                                                                                                             <  2:1324932/150‑1 (MQ=255)
           cccccTGGTAAAGTATTTACGCTTAGAGGCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGTTTTATCGGGGCAGa                                                                                                                                            <  1:376932/150‑1 (MQ=255)
           cccccTGGTAAAGTATTTACGCTTAGAGGCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGa                                                                                                                                            <  1:1382490/150‑1 (MQ=255)
           cccccTGGTAAAGTATTTACGCTTAGAGGCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGa                                                                                                                                            >  2:685486/1‑150 (MQ=255)
           cccccTGGTAAAGTATTTACGCTTAGAGGCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTATTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATTTTATATCTGCCCCGATAAAACGGGGCAGa                                                                                                                                            <  1:1518476/150‑1 (MQ=255)
            ccccTGGTAAAGTATTTACGCTTAGAGGCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGat                                                                                                                                           >  2:141709/1‑150 (MQ=255)
            ccccTGGTAAAGTATTTACGCTTAGAGGCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGat                                                                                                                                           <  2:1225383/150‑1 (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)
                        tATTTACGCTTAGAGGCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATgtttagt                                                                                                                               >  1:782959/1‑150 (MQ=255)
                         aTTTACGCTTAGAGGCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAACACAGGAGAACCGCCGAACTATCTCTATTGATAATTTAATGTTATATCTGCCCCGATAAAACGGGTCAGATAATATgtttagtt                                                                                                                              >  1:133109/1‑150 (MQ=255)
                         aTTTACGCTTAGAGGCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATgtttagtt                                                                                                                              <  1:986698/150‑1 (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:84171/1‑150 (MQ=255)
                                      ggCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCAt                                                                                                                 >  1:1331126/1‑150 (MQ=255)
                                         ggAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCAtttt                                                                                                              >  2:1689475/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                                                                                               <  2:1073967/150‑1 (MQ=255)
                                                          tatTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTAc                                                                                             <  2:1702526/150‑1 (MQ=255)
                                                          tatTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTAc                                                                                             <  1:1629272/150‑1 (MQ=255)
                                                            ttCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTg                                                                                           <  2:10922/150‑1 (MQ=255)
                                                               cTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGttt                                                                                        <  1:573332/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)
                                                                      tAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCg                                                                                 <  2:1645085/150‑1 (MQ=255)
                                                                      actcgcTCTTCTTTTAGGGAAAACAGTAGAACCGCCGACCTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCg                                                                                 <  2:137119/147‑1 (MQ=255)
                                                                          gcTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTTTTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCa                                                                             <  1:1615667/150‑1 (MQ=255)
                                                                          gcTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCa                                                                             >  2:1656612/1‑150 (MQ=255)
                                                                           ctctTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAg                                                                            <  1:379198/150‑1 (MQ=255)
                                                                            tcttctTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGa                                                                           >  1:1624627/1‑150 (MQ=255)
                                                                             cttcttTTAGGTAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGAt                                                                          <  2:1069590/150‑1 (MQ=255)
                                                                                cttTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGTTAAt                                                                       >  2:841294/1‑150 (MQ=255)
                                                                                cttTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAAt                                                                       >  2:1434471/1‑150 (MQ=255)
                                                                                cttTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAAt                                                                       >  1:969796/1‑150 (MQ=255)
                                                                                cttTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAAt                                                                       >  2:308670/1‑150 (MQ=255)
                                                                                 ttttAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAAta                                                                      >  1:416550/1‑149 (MQ=255)
                                                                                  tttAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCt                                                                     >  2:74381/1‑150 (MQ=255)
                                                                                    tAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGt                                                                   <  2:626960/150‑1 (MQ=255)
                                                                                     tgggAAAACAGGAGAACCGCCGAACTATCTATATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGtt                                                                  <  2:1450275/149‑1 (MQ=255)
                                                                                     aGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGtt                                                                  <  2:1729552/150‑1 (MQ=255)
                                                                                      gggAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTtc                                                                 <  2:682286/150‑1 (MQ=255)
                                                                                      gggAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTtc                                                                 <  1:685486/150‑1 (MQ=255)
                                                                                       ggAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTtct                                                                >  2:1371495/1‑150 (MQ=255)
                                                                                       ggAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTtct                                                                >  1:931748/1‑150 (MQ=255)
                                                                                        gAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTtcta                                                               >  2:819908/1‑149 (MQ=255)
                                                                                           aaCAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTaa                                                            <  1:719111/150‑1 (MQ=255)
                                                                                            aCAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAAt                                                           >  2:606891/1‑150 (MQ=255)
                                                                                             cAGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAAtt                                                          >  2:1454892/1‑149 (MQ=255)
                                                                                              aGGAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGt                                                         <  2:719428/150‑1 (MQ=255)
                                                                                               ggAGAACCGCCGAACTATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACCGTTTTTTATCGACCAGATAATCTGTTCTCTAATGtt                                                        <  2:1483795/150‑1 (MQ=255)
                                                                                                              tATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTCCCCCTAACCt                                         >  1:743120/1‑150 (MQ=255)
                                                                                                              tATCTCTATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTCCCACTAACCt                                         >  1:301400/1‑150 (MQ=255)
                                                                                                                 ctctATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTCCCCCTAACCTGtt                                      <  1:1689475/150‑1 (MQ=255)
                                                                                                                  tctATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTAATAAAGTTAATTTTGCAGTGAAGACATTTACTGTTTTTTATCGACCAGATAATCTTTTATCTAATGTTAACTCACCCTAACCTGttt                                     >  1:72576/1‑149 (MQ=255)
                                                                                                                   ctATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTCCCCCTAACCTGTTGc                                    >  2:465663/1‑150 (MQ=255)
                                                                                                                   ctATTTATAATTTAATGTTATATCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTCCCCCTAACCTGTTGc                                    >  2:10877/1‑150 (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)
                                                                                                                          taATTTAATGTTATATCTGCCACGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTCCCCCTAACCTGTTGCTTTAGtt                             <  1:887784/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                 >  2:100827/2‑150 (MQ=255)
                                                                                                                                      ttatCTGCCCCGTTTTATCGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTCCCCCTAACCTGTTGCTTTAGTTATTCATTTCCTg                 >  1:966851/2‑150 (MQ=255)
                                                                                                                                      atatCAGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTCCCCCTAACCTGTTGCTTTAGTTATTCATTTCCTg                 <  1:852125/150‑1 (MQ=255)
                                                                                                                                        atCTGCCCCGTTTTATCGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTCCCCCTAACCTGTTGCTTTAGTTATTCATTTCCTGtc               >  2:247967/1‑150 (MQ=255)
                                                                                                                                        atCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTCCCCCTAACCTGTTGCTTTAGTTATTCATTTCCTGtc               <  2:1331126/150‑1 (MQ=255)
                                                                                                                                        atCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTCCCCCTAACCTGTTGCTTTAGTTATTCATTTCCTGtc               >  2:1455664/1‑150 (MQ=255)
                                                                                                                                        atCTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTAATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATTATCTGTTCTCTAATGTTAACTCCCCCTAACCTGTTGCTTTAGTTATTCATTTCCTGtc               >  2:928121/1‑150 (MQ=255)
                                                                                                                                         tCTGCCCCGTTTTATCGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTACTCTAATGTTAACTCCCCCTAACCTGTTGCTTTAGTTATTCATTTCCTGtct              >  2:1555799/1‑150 (MQ=255)
                                                                                                                                          cTGCCCCGATAAAACGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTCCACATAACCTGTGGCTTTAGTTATTCATTTCCTGtcta             >  2:835631/1‑149 (MQ=255)
                                                                                                                                                     aaaCGGGGCAGATAATATGTTTAGTTTACTAACGGTCATTTTGCAGTGAAGCCATTTACTGTTTTTTATCGACCAGATAATCTGTTCTCTAATGTTAACTCCCCCTAACCTGTTGCTTTAGTTATTCATTTCCTGTCTCACTTTGCCTTa  <  1:1363659/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 < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: