Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,504,880 A→C 6.3% intergenic (+97/+55) frlR → / ← yhfS putative DNA‑binding transcriptional regulator/FNR‑regulated pyridoxal phosphate‑dependent aminotransferase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,504,8800AC 6.3% 563.3 / 4.8 198intergenic (+97/+55)frlR/yhfSputative DNA‑binding transcriptional regulator/FNR‑regulated pyridoxal phosphate‑dependent aminotransferase family protein
Reads supporting (aligned to +/- strand):  ref base A (84/94);  new base C (12/5);  total (97/101)
Fisher's exact test for biased strand distribution p-value = 7.82e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 206 total aligned reads displayed.

TCTTCTGCCGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCTGATTGTG  >  NC_000913/3504731‑3505029
                                                                                                                                                     |                                                                                                                                                      
tcttctGCCGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCGAACCCACACCCCa                                                                                                                                                        >  1:2578664/1‑150 (MQ=255)
 cttctGCCGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCaa                                                                                                                                                       <  1:4242460/150‑1 (MQ=255)
   tctGCCGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAg                                                                                                                                                     <  1:1297915/150‑1 (MQ=255)
   tctGCCGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAg                                                                                                                                                     <  1:3721504/150‑1 (MQ=255)
   tctGCCGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAg                                                                                                                                                     <  1:1880084/150‑1 (MQ=255)
        cgcgCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGTCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAACTTTACATCCCGTACGTTCCCCTAACCCAAACCCTCTCACCAAAgggggg                                                                                                                                                >  1:4010807/1‑148 (MQ=39)
         gcgcCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTGCCCACACCCAAAACCACACCCCAAAagggggg                                                                                                                                               >  1:1008492/1‑147 (MQ=38)
         gcgcCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCACACCCAAACCCACACCCCAAAAGGGCGa                                                                                                                                               >  1:1769010/1‑150 (MQ=39)
         gcgcCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGGAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCgg                                                                                                                                               >  2:2488043/1‑149 (MQ=255)
           gcCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCAAACCCACTCCCCCAAGGGGCGAgg                                                                                                                                             >  1:806082/1‑150 (MQ=255)
            ccAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCTTTGGGGAGAggg                                                                                                                                            <  2:4443162/150‑1 (MQ=255)
                    aTCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCAACCAAAACCACTCCCCAAACGGGCGACGGGGCCGGTc                                                                                                                                    >  2:1123259/1‑150 (MQ=38)
                    aTCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTACCCTCACCCTAAACCACACCCCAAAGGGGCGAGGGGACCGTTc                                                                                                                                    >  1:2784600/1‑150 (MQ=255)
                     tCACCTTAACTATTGGTAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCa                                                                                                                                   <  1:3289713/150‑1 (MQ=255)
                     tCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCa                                                                                                                                   <  2:3236909/150‑1 (MQ=255)
                     tCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCACACCCAAAACCACACCCCAAACGGGCGCGGGGACCGGTCa                                                                                                                                   >  1:2760468/1‑150 (MQ=38)
                      cACCTTACCGATTGATAATAAAAGACATTCCCCGTCGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGGTCCCCTCACCCTAACCCTCTCCCCAACGGGGCGAGGGGACCGTTCAg                                                                                                                                 <  2:718748/150‑1 (MQ=38)
                      cACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAg                                                                                                                                 <  2:105278/150‑1 (MQ=255)
                        ccTTAACTATTGATAATAAAAGACCTTAACCGTAGGCCGGAGAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTa                                                                                                                               <  2:315975/150‑1 (MQ=255)
                         cTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTAc                                                                                                                              <  2:1896036/150‑1 (MQ=255)
                         cTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTAc                                                                                                                              <  2:4279330/150‑1 (MQ=255)
                         cTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCACCACACCCCAAACCACAACCCAAAGGGGCGGGGGGCCCCGTAAGTAc                                                                                                                              >  2:2108781/1‑150 (MQ=37)
                          ttAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCAAACCCTCACCCCAAAGGGGCGAGGGGCCCGGTCAGTACg                                                                                                                             >  2:1729717/1‑150 (MQ=255)
                          ttAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTACCCACACCCAAAACCACAACCCAAAGGGGCGAGGGGGCCCGTCAGTAcc                                                                                                                             >  2:1478322/1‑149 (MQ=37)
                          gtcaCTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGAGGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACg                                                                                                                             <  1:2912922/147‑1 (MQ=255)
                           tAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACgg                                                                                                                            <  1:4333458/150‑1 (MQ=255)
                             aCTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGaa                                                                                                                          <  2:2046001/150‑1 (MQ=255)
                                  tGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCACACCCTAAACCTCACCCCAAAGGGGCGAGGGGACCGGTCAGGACGGACTTTAg                                                                                                                     >  1:4703009/1‑150 (MQ=39)
                                   gATAATAAAAGACATTACCCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGc                                                                                                                    <  1:861357/150‑1 (MQ=255)
                                     taataaAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTACCCACACCCCAAACCACACCCCCAAGGGGCGAGGGGGCCGGTCAGTAACGAATGTAGCtt                                                                                                                  >  1:2348958/1‑150 (MQ=37)
                                        taaAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGAc                                                                                                               <  1:4593310/150‑1 (MQ=255)
                                        taaAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGGACGTTCCCCACCCCCAAAACCACAACCCAAAAGGGCGAGGGGACCCGTCAGTACCGAAAATAGCTTGAc                                                                                                               >  1:2886797/1‑150 (MQ=25)
                                          aaaGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTGCCCCCACCCAAAACCGAAACCCAAAGGGGCGGGGGGACCCTTCAGTTCCGAAGTTATCTTgtagg                                                                                                             >  1:1991632/1‑146 (MQ=25)
                                             gACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGc                                                                                                          <  1:2823855/150‑1 (MQ=255)
                                                 ttAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATc                                                                                                      >  2:2900803/1‑150 (MQ=255)
                                                 ttAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATc                                                                                                      <  2:3297730/150‑1 (MQ=255)
                                                 ttAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTACCCTCACCCTAAACCTCTACCCAAAGGGGCGTGGGGACCGTTGAGTACGGAATTTTGCTTGTAGGTGCTAtg                                                                                                      >  2:3967097/1‑149 (MQ=38)
                                                   aaCCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCa                                                                                                    <  2:4288062/150‑1 (MQ=255)
                                                    aCCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAg                                                                                                   <  1:3376944/150‑1 (MQ=255)
                                                     ccGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCCCCCGCACCCCCTCCCCCCCGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGc                                                                                                  <  2:1293862/150‑1 (MQ=38)
                                                      cGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTACCCACACCCAAACCCTCTCCCCAAAGGGGCGAGGGGACCGGTCAGTACGGAAATTAGCTTGACGGTGCTATCCAGCa                                                                                                 >  2:2723436/1‑150 (MQ=39)
                                                      cGTAGGCCGGATAAGATGCGCCAGCATCACATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCACACCCAAAACCACACCCCCAAGGGGCGCGGGGCCCCGTAAGTAAGGAATTTATCTTGGAGGTGCGATACaccc                                                                                                 >  1:3819054/1‑147 (MQ=25)
                                                           gCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTaaa                                                                                            <  1:1263688/150‑1 (MQ=255)
                                                               gATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTGCCCTCCCCCAAACCCGCTACCCCAAGGGGCGGGGGGACCGGGGAGAACGGCATTTATCGTGAAGGGGCGATCCAGCATTAAATaat                                                                                        >  2:1391296/1‑148 (MQ=255)
                                                               gATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGGTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACt                                                                                        <  1:3283169/150‑1 (MQ=255)
                                                                aTAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTg                                                                                       <  2:3994681/150‑1 (MQ=255)
                                                                 tAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTgg                                                                                      <  1:3501560/150‑1 (MQ=255)
                                                                      tGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTGCCCGCACCCAAAACCACAACCCAAAGCGGCCGGGGGGCCGGTGAGGAACGGAATTTTATTGTCGGGGCTATCCACCAGTAAAAACAGGccaaa                                                                                 >  1:3044486/1‑146 (MQ=255)
                                                                       gcgcCAGCATCGCCGGATGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGGTCCCCCCACCCAAACCCCCACCCCCAAGGGGCGTGGGGCACGTGCAGTACGGAAATTAGCTTGACGGGGCTATCCAGCCATAAAAAAAGGCCATAc                                                                                >  2:1829427/1‑150 (MQ=255)
                                                                        cgcCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCACAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATAct                                                                               <  2:3103734/150‑1 (MQ=255)
                                                                        cgcCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATAct                                                                               <  1:1060671/150‑1 (MQ=255)
                                                                        cgcCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATAct                                                                               <  1:3196486/150‑1 (MQ=255)
                                                                         gcCAGCATCGCATCCCGCGATGCTGGCGCGAAGAAATTAACATCCCGTACCTTCCCCTCACCCAAACCCTCACCCCAAACGGGCGAGGGGACCGTTCAGTACGGAATTTAGCATGACGGTGCTATCCAGCATTAAATACTGGCAATActc                                                                              <  2:824877/150‑1 (MQ=255)
                                                                            aGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCa                                                                           <  1:4458167/150‑1 (MQ=255)
                                                                             gcatcgcatcCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCCAAGGGGCGCGGGGACCCTTCAGTACGGAATTTAGCTTGCCGGTGCTATCCAGCATTAAATACTGGCCATtctcgccc                                                                          >  2:2529204/1‑146 (MQ=39)
                                                                             gcatcgcatcCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCAc                                                                          <  1:1411696/150‑1 (MQ=255)
                                                                                tcgcatcCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTGCCCCCACCCAAAACCACACCCCCAAGGGGCGGGGGGGCCGGTCAGTACCGAAATTATCGTGGAGGTGCTAACCAGCAATAAAAAAAGGCAATACAATCACCCa                                                                       >  2:4551242/1‑150 (MQ=255)
                                                                                     tcCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCGAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACCGAATTTAGCTTGACGGTGCTATCCAGCAATAAATAATTGCAAAAATCACACGCAAAAta                                                                  >  2:1548508/1‑149 (MQ=38)
                                                                                     tcCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCAAAACCTATACCCAAAGGGGCGAGGGGACCCTTCAGTAAGGAATTTATCTTTACCGTGCTATACAGCATTAAAAACTGGCAAAACACTCAAGCaaaaaa                                                                  >  2:3785047/1‑148 (MQ=255)
                                                                                      ccGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATcc                                                                 <  1:4321280/150‑1 (MQ=255)
                                                                                        ggCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCgc                                                               <  2:3556441/150‑1 (MQ=255)
                                                                                        ggCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCgc                                                               <  1:3967097/150‑1 (MQ=255)
                                                                                          cGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAATGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCgcag                                                             <  1:4461402/150‑1 (MQ=255)
                                                                                          cGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCgcag                                                             <  2:3450220/150‑1 (MQ=255)
                                                                                                tGGCGCGTTGAATTTTACATCCCGTACGTTCCCCGCACCCAAACCCACACCCCCAACGGGCGCGGGGCCCGGTCAGTAAGGAAATTAGCATGAAGGGGCTATACAGCAATAAAAACAGGCAAAACTCACACCCAAAATACGCaccaccca                                                       >  2:584178/1‑148 (MQ=255)
                                                                                                   cgcgTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGtctc                                                    <  1:389737/150‑1 (MQ=255)
                                                                                                   cgcgTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGtctc                                                    <  1:180272/150‑1 (MQ=255)
                                                                                                      gTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTc                                                 <  1:1363303/150‑1 (MQ=255)
                                                                                                        tGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCAc                                               <  1:1566316/150‑1 (MQ=255)
                                                                                                            ttttACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGGCGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCt                                           <  2:87610/150‑1 (MQ=255)
                                                                                                             tttACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTg                                          <  2:1159920/150‑1 (MQ=255)
                                                                                                             tttACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTg                                          <  1:1584579/150‑1 (MQ=255)
                                                                                                                          cGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTaa                             <  1:1478322/150‑1 (MQ=255)
                                                                                                                                ccTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGaa                       >  1:4769075/1‑150 (MQ=255)
                                                                                                                                ccTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATAAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGaa                       >  1:2740658/1‑150 (MQ=255)
                                                                                                                                 cTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAAt                      <  1:1953602/150‑1 (MQ=255)
                                                                                                                                    aCCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCAATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGc                   >  1:1728466/1‑150 (MQ=255)
                                                                                                                                     cccTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGca                  <  1:4650107/150‑1 (MQ=255)
                                                                                                                                      ccTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCACAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGcac                 >  2:514977/1‑150 (MQ=255)
                                                                                                                                      ccTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGcac                 >  2:3811042/1‑150 (MQ=255)
                                                                                                                                      ccTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGcac                 >  1:4766073/1‑150 (MQ=255)
                                                                                                                                       cTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGcaca                >  2:3866835/1‑150 (MQ=255)
                                                                                                                                       cTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGcaca                >  2:250981/1‑150 (MQ=255)
                                                                                                                                        tAACCCTCTCCCCAAAGGGGCGAGGGGACCGTGCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCAGTAAATACTGGCAATACTCTCACGCAAAATCCGCCGCACCGTCTCTGCACCGCGGCGATGCGGGGTCATGCGACGCGcaccc               >  2:2700830/1‑148 (MQ=255)
                                                                                                                                           cccTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGt            >  2:1526513/1‑150 (MQ=255)
                                                                                                                                           cccTCGCCCCTTTGGGGAGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGt            >  2:2402585/1‑150 (MQ=255)
                                                                                                                                             ctctCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTc          >  1:1535716/1‑150 (MQ=255)
                                                                                                                                                tcCCCACAGGGGCGAGGGGACCGTTCAGTACGGACTTTAGCTTGACGGTGCTATCCAGCAGTCAATACTGGCAATACTCGCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGAGTCGGGTTAATGCGAATCGCACACTGGTCTGa       >  2:4636547/1‑150 (MQ=255)
                                                                                                                                                 ccccAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCTGAt      >  1:4057748/1‑150 (MQ=255)
                                                                                                                                                 ccccAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCTGAt      >  1:3551282/1‑150 (MQ=255)
                                                                                                                                                  cccAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCTGAtt     >  1:1660620/1‑150 (MQ=255)
                                                                                                                                                  cccAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCTGAtt     >  2:2592369/1‑150 (MQ=255)
                                                                                                                                                   ccAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACACTGTTCTGAGtg    >  2:3386422/1‑150 (MQ=255)
                                                                                                                                                   ccAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCTGATtg    >  1:3131392/1‑150 (MQ=255)
                                                                                                                                                   ccAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCTGATtg    >  1:2231357/1‑150 (MQ=255)
                                                                                                                                                   ccAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCTGATtg    >  2:319512/1‑150 (MQ=255)
                                                                                                                                                   ccAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCTGATtg    >  2:1251720/1‑150 (MQ=255)
                                                                                                                                                    cAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCTGATtgt   >  2:2332970/1‑150 (MQ=255)
                                                                                                                                                    cAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCTGATtgt   >  2:2919357/1‑150 (MQ=255)
                                                                                                                                                    cAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCTGATtgt   <  2:1615882/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
TCTTCTGCCGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCTGATTGTG  >  NC_000913/3504731‑3505029

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