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.

TCTTCTGCCGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCA‑GTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCTGATTGTG  >  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)
                      cACCTTACCGATTGATAATAAAAGACATTCCCCGTCGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGGTCCCCTCACCCTAACCCTCTCCCCAACGGGGCGAGGGGACCGTTCA‑g                                                                                                                                 <  2:718748/150‑1 (MQ=38)
                      cACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCA‑g                                                                                                                                 <  2:105278/150‑1 (MQ=255)
                        ccTTAACTATTGATAATAAAAGACCTTAACCGTAGGCCGGAGAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCA‑GTa                                                                                                                               <  2:315975/150‑1 (MQ=255)
                         cTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCA‑GTAc                                                                                                                              <  2:1896036/150‑1 (MQ=255)
                         cTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCA‑GTAc                                                                                                                              <  2:4279330/150‑1 (MQ=255)
                         cTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCACCACACCCCAAACCACAACCCAAAGGGGCGGGGGGCCCCGTAA‑GTAc                                                                                                                              >  2:2108781/1‑150 (MQ=37)
                          ttAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCAAACCCTCACCCCAAAGGGGCGAGGGGCCCGGTCA‑GTACg                                                                                                                             >  2:1729717/1‑150 (MQ=255)
                          ttAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTACCCACACCCAAAACCACAACCCAAAGGGGCGAGGGGGCCCGTCA‑GTAcc                                                                                                                             >  2:1478322/1‑149 (MQ=37)
                          gtcaCTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGAGGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCA‑GTACg                                                                                                                             <  1:2912922/147‑1 (MQ=255)
                           tAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCA‑GTACgg                                                                                                                            <  1:4333458/150‑1 (MQ=255)
                             aCTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCA‑GTACGGaa                                                                                                                          <  2:2046001/150‑1 (MQ=255)
                                  tGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCACACCCTAAACCTCACCCCAAAGGGGCGAGGGGACCGGTCA‑GGACGGACTTTAg                                                                                                                     >  1:4703009/1‑150 (MQ=39)
                                   gATAATAAAAGACATTACCCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCA‑GTACGGAATTTAGc                                                                                                                    <  1:861357/150‑1 (MQ=255)
                                     taataaAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTACCCACACCCCAAACCACACCCCCAAGGGGCGAGGGGGCCGGTCA‑GTAACGAATGTAGCtt                                                                                                                  >  1:2348958/1‑150 (MQ=37)
                                        taaAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCA‑GTACGGAATTTAGCTTGAc                                                                                                               <  1:4593310/150‑1 (MQ=255)
                                        taaAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGGACGTTCCCCACCCCCAAAACCACAACCCAAAAGGGCGAGGGGACCCGTCA‑GTACCGAAAATAGCTTGAc                                                                                                               >  1:2886797/1‑150 (MQ=25)
                                          aaaGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTGCCCCCACCCAAAACCGAAACCCAAAGGGGCGGGGGGACCCTTCA‑GTTCCGAAGTTATCTTgtagg                                                                                                             >  1:1991632/1‑146 (MQ=25)
                                             gACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCA‑GTACGGAATTTAGCTTGACGGTGc                                                                                                          <  1:2823855/150‑1 (MQ=255)
                                                 ttAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCA‑GTACGGAATTTAGCTTGACGGTGCTATc                                                                                                      >  2:2900803/1‑150 (MQ=255)
                                                 ttAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCA‑GTACGGAATTTAGCTTGACGGTGCTATc                                                                                                      <  2:3297730/150‑1 (MQ=255)
                                                 ttAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTACCCTCACCCTAAACCTCTACCCAAAGGGGCGTGGGGACCGTTGA‑GTACGGAATTTTGCTTGTAGGTGCTAtg                                                                                                      >  2:3967097/1‑149 (MQ=38)
                                                   aaCCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCA‑GTACGGAATTTAGCTTGACGGTGCTATCCa                                                                                                    <  2:4288062/150‑1 (MQ=255)
                                                    aCCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCA‑GTACGGAATTTAGCTTGACGGTGCTATCCAg                                                                                                   <  1:3376944/150‑1 (MQ=255)
                                                     ccGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCCCCCGCACCCCCTCCCCCCCGGGGCGAGGGGACCGTTCA‑GTACGGAATTTAGCTTGACGGTGCTATCCAGc                                                                                                  <  2:1293862/150‑1 (MQ=38)
                                                      cGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTACCCACACCCAAACCCTCTCCCCAAAGGGGCGAGGGGACCGGTCA‑GTACGGAAATTAGCTTGACGGTGCTATCCAGCa                                                                                                 >  2:2723436/1‑150 (MQ=39)
                                                      cGTAGGCCGGATAAGATGCGCCAGCATCACATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCACACCCAAAACCACACCCCCAAGGGGCGCGGGGCCCCGTAA‑GTAAGGAATTTATCTTGGAGGTGCGATACaccc                                                                                                 >  1:3819054/1‑147 (MQ=25)
                                                           gCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCA‑GTACGGAATTTAGCTTGACGGTGCTATCCAGCATTaaa                                                                                            <  1:1263688/150‑1 (MQ=255)
                                                               gATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTGCCCTCCCCCAAACCCGCTACCCCAAGGGGCGGGGGGACCGGGGA‑GAACGGCATTTATCGTGAAGGGGCGATCCAGCATTAAATaat                                                                                        >  2:1391296/1‑148 (MQ=255)
                                                               gATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGGTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCA‑GTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACt                                                                                        <  1:3283169/150‑1 (MQ=255)
                                                                aTAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCA‑GTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTg                                                                                       <  2:3994681/150‑1 (MQ=255)
                                                                 tAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCA‑GTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTgg                                                                                      <  1:3501560/150‑1 (MQ=255)
                                                                      tGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTGCCCGCACCCAAAACCACAACCCAAAGCGGCCGGGGGGCCGGTGAGGAACGGAATTT‑TATTGTCGGGGCTATCCACCAGTAAAAACAGGccaaa                                                                                 >  1:3044486/1‑146 (MQ=255)
                                                                       gcgcCAGCATCGCCGGATGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGGTCCCCCCACCCAAACCCCCACCCCCAAGGGGCGTGGGGCACGTGCA‑GTACGGAAATTAGCTTGACGGGGCTATCCAGCCATAAAAAAAGGCCATAc                                                                                >  2:1829427/1‑150 (MQ=255)
                                                                        cgcCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCACAGGGGCGAGGGGACCGTTCA‑GTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATAct                                                                               <  2:3103734/150‑1 (MQ=255)
                                                                        cgcCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCA‑GTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATAct                                                                               <  1:1060671/150‑1 (MQ=255)
                                                                        cgcCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCA‑GTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATAct                                                                               <  1:3196486/150‑1 (MQ=255)
                                                                         gcCAGCATCGCATCCCGCGATGCTGGCGCGAAGAAATTAACATCCCGTACCTTCCCCTCACCCAAACCCTCACCCCAAACGGGCGAGGGGACCGTTCA‑GTACGGAATTTAGCATGACGGTGCTATCCAGCATTAAATACTGGCAATActc                                                                              <  2:824877/150‑1 (MQ=255)
                                                                            aGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCA‑GTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCa                                                                           <  1:4458167/150‑1 (MQ=255)
                                                                             gcatcgcatcCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCCAAGGGGCGCGGGGACCCTTCA‑GTACGGAATTTAGCTTGCCGGTGCTATCCAGCATTAAATACTGGCCATtctcgccc                                                                          >  2:2529204/1‑146 (MQ=39)
                                                                             gcatcgcatcCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCA‑GTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCAc                                                                          <  1:1411696/150‑1 (MQ=255)
                                                                                tcgcatcCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTGCCCCCACCCAAAACCACACCCCCAAGGGGCGGGGGGGCCGGTCA‑GTACCGAAATTATCGTGGAGGTGCTAACCAGCAATAAAAAAAGGCAATACAATCACCCa                                                                       >  2:4551242/1‑150 (MQ=255)
                                                                                     tcCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCGAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCA‑GTACCGAATTTAGCTTGACGGTGCTATCCAGCAATAAATAATTGCAAAAATCACACGCAAAAta                                                                  >  2:1548508/1‑149 (MQ=38)
                                                                                     tcCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCAAAACCTATACCCAAAGGGGCGAGGGGACCCTTCA‑GTAAGGAATTTATCTTTACCGTGCTATACAGCATTAAAAACTGGCAAAACACTCAAGCaaaaaa                                                                  >  2:3785047/1‑148 (MQ=255)
                                                                                      ccGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCA‑GTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATcc                                                                 <  1:4321280/150‑1 (MQ=255)
                                                                                        ggCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCA‑GTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCgc                                                               <  2:3556441/150‑1 (MQ=255)
                                                                                        ggCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCA‑GTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCgc                                                               <  1:3967097/150‑1 (MQ=255)
                                                                                          cGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAATGGGGCGAGGGGACCGTTCA‑GTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCgcag                                                             <  1:4461402/150‑1 (MQ=255)
                                                                                          cGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCA‑GTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCgcag                                                             <  2:3450220/150‑1 (MQ=255)
                                                                                                tGGCGCGTTGAATTTTACATCCCGTACGTTCCCCGCACCCAAACCCACACCCCCAACGGGCGCGGGGCCCGGTCA‑GTAAGGAAATTAGCATGAAGGGGCTATACAGCAATAAAAACAGGCAAAACTCACACCCAAAATACGCaccaccca                                                       >  2:584178/1‑148 (MQ=255)
                                                                                                   cgcgTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCA‑GTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGtctc                                                    <  1:389737/150‑1 (MQ=255)
                                                                                                   cgcgTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCA‑GTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGtctc                                                    <  1:180272/150‑1 (MQ=255)
                                                                                                      gTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCA‑GTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTc                                                 <  1:1363303/150‑1 (MQ=255)
                                                                                                        tGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCA‑GTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCAc                                               <  1:1566316/150‑1 (MQ=255)
                                                                                                            ttttACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCA‑GTACGGAATTTAGCTTGGCGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCt                                           <  2:87610/150‑1 (MQ=255)
                                                                                                             tttACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCA‑GTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTg                                          <  2:1159920/150‑1 (MQ=255)
                                                                                                             tttACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCA‑GTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTg                                          <  1:1584579/150‑1 (MQ=255)
                                                                                                                          cGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCA‑GTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTaa                             <  1:1478322/150‑1 (MQ=255)
                                                                                                                                ccTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCA‑GTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGaa                       >  1:4769075/1‑150 (MQ=255)
                                                                                                                                ccTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCA‑GTACGGAATTTAGCTTGACGGTGCTATCCAGCATAAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGaa                       >  1:2740658/1‑150 (MQ=255)
                                                                                                                                 cTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCA‑GTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAAt                      <  1:1953602/150‑1 (MQ=255)
                                                                                                                                    aCCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCA‑GTACGGAATTTAGCTTGACGGTGCAATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGc                   >  1:1728466/1‑150 (MQ=255)
                                                                                                                                     cccTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCA‑GTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGca                  <  1:4650107/150‑1 (MQ=255)
                                                                                                                                      ccTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCA‑GTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCACAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGcac                 >  2:514977/1‑150 (MQ=255)
                                                                                                                                      ccTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCA‑GTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGcac                 >  2:3811042/1‑150 (MQ=255)
                                                                                                                                      ccTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCA‑GTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGcac                 >  1:4766073/1‑150 (MQ=255)
                                                                                                                                       cTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCA‑GTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGcaca                >  2:3866835/1‑150 (MQ=255)
                                                                                                                                       cTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCA‑GTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGcaca                >  2:250981/1‑150 (MQ=255)
                                                                                                                                        tAACCCTCTCCCCAAAGGGGCGAGGGGACCGTGCA‑GTACGGAATTTAGCTTGACGGTGCTATCCAGCAGTAAATACTGGCAATACTCTCACGCAAAATCCGCCGCACCGTCTCTGCACCGCGGCGATGCGGGGTCATGCGACGCGcaccc               >  2:2700830/1‑148 (MQ=255)
                                                                                                                                           cccTCTCCCCAAAGGGGCGAGGGGACCGTTCA‑GTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGt            >  2:1526513/1‑150 (MQ=255)
                                                                                                                                           cccTCGCCCCTTTGGGGAGAGGGGACCGTTCA‑GTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGt            >  2:2402585/1‑150 (MQ=255)
                                                                                                                                             ctctCCCCAAAGGGGCGAGGGGACCGTTCA‑GTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTc          >  1:1535716/1‑150 (MQ=255)
                                                                                                                                                tcCCCACAGGGGCGAGGGGACCGTTCA‑GTACGGACTTTAGCTTGACGGTGCTATCCAGCAGTCAATACTGGCAATACTCGCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGAGTCGGGTTAATGCGAATCGCACACTGGTCTGa       >  2:4636547/1‑150 (MQ=255)
                                                                                                                                                 ccccAAAGGGGCGAGGGGACCGTTCA‑GTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCTGAt      >  1:4057748/1‑150 (MQ=255)
                                                                                                                                                 ccccAAAGGGGCGAGGGGACCGTTCA‑GTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCTGAt      >  1:3551282/1‑150 (MQ=255)
                                                                                                                                                  cccAAAGGGGCGAGGGGACCGTTCA‑GTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCTGAtt     >  1:1660620/1‑150 (MQ=255)
                                                                                                                                                  cccAAAGGGGCGAGGGGACCGTTCA‑GTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCTGAtt     >  2:2592369/1‑150 (MQ=255)
                                                                                                                                                   ccAAAGGGGCGAGGGGACCGTTCA‑GTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACACTGTTCTGAGtg    >  2:3386422/1‑150 (MQ=255)
                                                                                                                                                   ccAAAGGGGCGAGGGGACCGTTCA‑GTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCTGATtg    >  1:3131392/1‑150 (MQ=255)
                                                                                                                                                   ccAAAGGGGCGAGGGGACCGTTCA‑GTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCTGATtg    >  1:2231357/1‑150 (MQ=255)
                                                                                                                                                   ccAAAGGGGCGAGGGGACCGTTCA‑GTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCTGATtg    >  2:319512/1‑150 (MQ=255)
                                                                                                                                                   ccAAAGGGGCGAGGGGACCGTTCA‑GTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCTGATtg    >  2:1251720/1‑150 (MQ=255)
                                                                                                                                                    cAAAGGGGCGAGGGGACCGTTCA‑GTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCTGATtgt   >  2:2332970/1‑150 (MQ=255)
                                                                                                                                                    cAAAGGGGCGAGGGGACCGTTCA‑GTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCTGATtgt   >  2:2919357/1‑150 (MQ=255)
                                                                                                                                                    cAAAGGGGCGAGGGGACCGTTCA‑GTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCTGATtgt   <  2:1615882/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
TCTTCTGCCGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACCGTTCA‑GTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCTGATTGTG  >  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: ‑