Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,504,880 A→C 5.5% 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 5.5% 1403.0 / 8.6 470intergenic (+97/+55)frlR/yhfSputative DNA‑binding transcriptional regulator/FNR‑regulated pyridoxal phosphate‑dependent aminotransferase family protein
Reads supporting (aligned to +/- strand):  ref base A (180/243);  new base C (31/6);  total (218/252)
Fisher's exact test for biased strand distribution p-value = 1.18e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 486 total aligned reads displayed.

TTCTGCCGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTC‑TCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCTGATTGTG  >  NC_000913/3504733‑3505029
                                                                                                                                                    |                                                                                                                                                     
ttctGCCGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGGTCCCCTCACCCAAACCCGC‑ACCCCaaa                                                                                                                                                     >  1:515921/1‑150 (MQ=255)
  ctGCCGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTC‑TCCCCAAAgg                                                                                                                                                   <  2:6051635/150‑1 (MQ=255)
    gCCGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTGCCCTCACCCAAACCCAC‑ACACCAAAcggg                                                                                                                                                 >  2:2655159/1‑150 (MQ=39)
     ccGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTC‑TCCCCAAAGGGGc                                                                                                                                                <  1:5316632/150‑1 (MQ=255)
        cgcCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTATGTTCCCCTCACCCAAACCCAC‑TACCCAAAGGGGCGAg                                                                                                                                             >  1:6184548/1‑150 (MQ=255)
         gcCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTC‑TCCCCAAAGGGGCGAgg                                                                                                                                            <  2:4450497/150‑1 (MQ=255)
         gcCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTC‑TCCCCAAAGGGGCGAgg                                                                                                                                            <  2:2390820/150‑1 (MQ=255)
          ccAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTC‑TCCCCAAAGGGGCGAggg                                                                                                                                           <  1:1012199/150‑1 (MQ=255)
          ccAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGGTCCCCTCACCCTAACCCTC‑TCCCCAAAGGGGCGAggg                                                                                                                                           <  1:220092/150‑1 (MQ=255)
           cAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTC‑TCCCCAAAGGGGCGAgggg                                                                                                                                          <  1:7411051/150‑1 (MQ=255)
           cAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCGCCCCCGAACCCACAACCCCAAAgggcgccggg                                                                                                                                           >  1:4308791/1‑143 (MQ=38)
            aaTCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTC‑TCCCCAAAGGGGCGAGGGGa                                                                                                                                         <  1:1967878/150‑1 (MQ=255)
             aTCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTC‑TCCCCAAAGGGGCGAGGGGAc                                                                                                                                        <  1:3829341/150‑1 (MQ=255)
             aTCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGGAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTACCCTCACCCAAAACCAC‑TCCCCAAAGGGCCCAGGGGAc                                                                                                                                        >  2:2135893/1‑150 (MQ=39)
              tCGCATCACCTTAACTATTGATAATAAAAGACCTTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTC‑TCCCCAAAGGGGCGAGGGGAcc                                                                                                                                       <  2:1434742/150‑1 (MQ=255)
              tCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTC‑TCCCCAAAGGGGCGAGGGGAcc                                                                                                                                       <  1:800177/150‑1 (MQ=255)
              tCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTC‑TCCCCAAAGGGGCGAGGGGAcc                                                                                                                                       <  1:8132845/150‑1 (MQ=255)
                gCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTC‑TCCCCAAAGGGGCGAGGGGACCGt                                                                                                                                     <  1:3977195/150‑1 (MQ=255)
                 cATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTC‑TACCCACAGGGGCGAGGGGCCCGtt                                                                                                                                    >  2:5854113/1‑150 (MQ=255)
                   tCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGTCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTC‑TGCCCAAATGGGCGAGGGGACCGTTCa                                                                                                                                  <  2:4852688/150‑1 (MQ=255)
                      ccTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTC‑TCCCCAAAGGGGCGAGGGGACCGTTCAGTa                                                                                                                               <  1:6074392/150‑1 (MQ=255)
                       cTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTC‑TCCCCAAAGGGGCGAGGGGACCGTTCAGTAc                                                                                                                              <  1:4968215/150‑1 (MQ=255)
                       cTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTC‑TCCCCAAAGGGGCGAGGGGACCGTTCAGTAc                                                                                                                              <  1:5597737/150‑1 (MQ=255)
                         tAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCACACCCAAACCCAC‑ACCCCCCACGGGCGCGGGGCCCCGTCAGTaacg                                                                                                                            >  1:3167036/1‑147 (MQ=37)
                               ttGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTC‑TCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTa                                                                                                                      <  1:2699929/150‑1 (MQ=255)
                                tGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTC‑TCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAg                                                                                                                     <  1:7821200/150‑1 (MQ=255)
                                 gagaataaCAGACAGGAACCGTAGGCCGGATAAGAGGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACAGCCCGTACGTTCCCCTCACCCGAACCCTC‑TCCCCAACGGGGCGAGGGGACCGTTCAGTACGGAATTTAGc                                                                                                                    <  1:4487828/147‑1 (MQ=255)
                                 gATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTC‑TCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGTATTTAGc                                                                                                                    <  1:7633205/150‑1 (MQ=255)
                                 gATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTC‑TCCCCAAAGGGGCGAGGGGACCGGTCAGTACGGAATTTAGc                                                                                                                    >  2:6143526/1‑150 (MQ=255)
                                  ataataAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTC‑TCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCt                                                                                                                   >  2:6465602/1‑150 (MQ=255)
                                   taataaAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTC‑TCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCtt                                                                                                                  <  2:3030814/150‑1 (MQ=255)
                                      taaAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTC‑TCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGAc                                                                                                               <  2:8176469/150‑1 (MQ=255)
                                        aaaGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACCTCCCGTACGTTCCCCTCACCCTAACCCTC‑TCCCCAAAGGGGCGGGGGGCCCGTTCAGTACGGAATTTAGCTTGACgg                                                                                                             >  2:9649817/1‑150 (MQ=255)
                                        aaaGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTACCCGCCCCCAAACCCAC‑ACCCCCAAGGGGCGAGGGGGCCGGTCAGTACGGGCTTTATCTTGGCgg                                                                                                             >  1:1367480/1‑150 (MQ=37)
                                        aaaGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATACCGTACGTTCCCCTCACCCAAACCCTC‑TCCCCAAAGGGGCGAGGGGACCGTTAAGTACGGAATTTATCTTGAAgg                                                                                                             >  1:7057164/1‑150 (MQ=39)
                                               ttAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCAACACCGAACCCGC‑TACCCCAAGTGGCCCGGGGCCCCTGCACTACCGAATTTAGCTTGACGGTGCtctc                                                                                                      >  2:8530177/1‑150 (MQ=25)
                                                    cGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTC‑TCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCa                                                                                                 <  2:1881101/150‑1 (MQ=255)
                                                    cGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTC‑TCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCa                                                                                                 <  1:1523048/150‑1 (MQ=255)
                                                         gCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTC‑TCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTaaa                                                                                            <  2:11643/150‑1 (MQ=255)
                                                         gCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTC‑TCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTaaa                                                                                            <  2:8932747/150‑1 (MQ=255)
                                                             gATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTC‑TCCCCAAAGGGGCGCGGGGACCGGTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACt                                                                                        >  2:5593037/1‑150 (MQ=255)
                                                              aTAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTC‑TCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTg                                                                                       <  1:5419788/150‑1 (MQ=255)
                                                              aTAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTC‑TCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTg                                                                                       <  1:9164520/150‑1 (MQ=255)
                                                               tAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTC‑TCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTgg                                                                                      <  1:5824390/150‑1 (MQ=255)
                                                               tAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTC‑TCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTgg                                                                                      <  1:8971263/150‑1 (MQ=255)
                                                                  gATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTGCCCGCACCCAAACCCGC‑ACCCCCAAGGGGCGGGGGGACCGGTCACTACCGCATTTAGCGTGACGGGGCTATCCCGCCGTAAAAACTGGCaa                                                                                   >  2:9168730/1‑150 (MQ=255)
                                                                      cgcCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTC‑TCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATAct                                                                               <  1:652987/150‑1 (MQ=255)
                                                                      cgcCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTC‑TCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATAct                                                                               <  2:4340682/150‑1 (MQ=255)
                                                                       gcCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTGCCCACACCCAAAACCAC‑AACCCAAAGGGGCGAGGGGTCCGTTCAGTAAGGAAATTATCTTGACGGTGCTATACAGCAATAAATAATGGCAAAAcac                                                                              >  2:2870095/1‑148 (MQ=255)
                                                                       gcCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTGTACATCCCGTACGGTCCCCTCACCCTAACCCTC‑TCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATActc                                                                              <  2:5164288/150‑1 (MQ=255)
                                                                        ccAGCATCGCATCCGGCGATGCTGGCGCGATGAATTTTACATCCCCTACGTACCCCTCACCCTAACCCAC‑ACCCCAAACCGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGAGCTATCCAGCATTAAATACTGGCAATActct                                                                             <  1:6684095/150‑1 (MQ=255)
                                                                          aGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGGTCCCCTCACCCTAACCCTC‑TCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCa                                                                           <  1:3460431/150‑1 (MQ=255)
                                                                                 catcCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCGCACCCAAACCCAC‑ACCCCAAAGGGGCGAGGGGGCCGGTCAGTACGGCAATTAGCGTGACGGTGCGATCCACCAATAAAAACAGGCAAAACACACACGCaaaa                                                                    >  2:7467808/1‑150 (MQ=255)
                                                                                     cGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTC‑TCCCCAACGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCg                                                                <  1:6262113/150‑1 (MQ=255)
                                                                                     cGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTC‑TCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCg                                                                <  1:1281944/150‑1 (MQ=255)
                                                                                      ggCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTC‑TCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCgc                                                               <  2:1051951/150‑1 (MQ=255)
                                                                                      ggCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCACACCCGAACCCTC‑ACCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCCGCCTTAAATACTGGCAATACTCTCACGCAAAATCCgc                                                               >  2:4678426/1‑150 (MQ=39)
                                                                                      ggCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCACACCCAAAACCGC‑TCCCCCAAGGGGCGAGGGGCCCCTTCAGTACCGAAAATATCTTGGCGGGGCTCTCCACCATAAAATACTGGCCATAATATCACCCAAAATCCgc                                                               >  2:8419821/1‑150 (MQ=255)
                                                                                         gATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCAAAACCTC‑ACCCCAAAGGGGCGGGGGGACCCTTCACTACGGAATTTAGCTTGACGGTGCTAACCAGCATTAAATACTGGCAAAACTCTCACGCAAAAACCgcccc                                                            >  2:7609724/1‑147 (MQ=255)
                                                                                          aTGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCGAAACCAC‑TACCCAAAGGGGCGAGGGGACCCTTCAGTACCGCAATTAGCATGACGGTGCTATCCAGCCTTAAATAATGGCAATAAACTCACGCAAAATCCgcagcc                                                           >  2:3869901/1‑149 (MQ=255)
                                                                                            gCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTC‑TCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCAcc                                                         <  1:488869/150‑1 (MQ=255)
                                                                                            gCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCGAACCCTC‑TCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCAcc                                                         <  1:8220179/150‑1 (MQ=255)
                                                                                              tGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTC‑TCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGt                                                       <  2:6965829/150‑1 (MQ=255)
                                                                                              tGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTC‑TCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGt                                                       <  1:2852487/150‑1 (MQ=255)
                                                                                               ggCGCGATGAATTTTACATCCCGTACGTTCCCCACACCCAAACCCTC‑TCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGCTGCTATCCAGCATTAAATACTGGCAATACTCTCACCCAAAATCCGCAGCACCGtc                                                      <  1:2514998/150‑1 (MQ=39)
                                                                                                     ttGAATTTTACATCCCGTACGTTCCCCACACCCAAACCCTC‑TACCCCAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCAATCCACCATTAAAAAATGGCAAAAATCACACGCAAAAAACGCAGCCCCGTCTCttaa                                                >  2:969837/1‑148 (MQ=255)
                                                                                                     ttGAATTTTACATCCCGTACGTTCCCCACACCCAAACCCAC‑ACCCCAAAGGGGCGGGGGGGCCCGTCAGTACGGAATTTATCTTTACGGTGCGATCCACCAATAAATAAAGGCCAAACACGAACCCCAAAACACCAGCACCCACTCTTCa                                                >  2:4541016/1‑150 (MQ=255)
                                                                                                     ttGAATTTTACATCCCGTACGGTCCCCTCACCCTAACCCTC‑TCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCa                                                <  1:113618/150‑1 (MQ=255)
                                                                                                      tGAAATTTACATCCCGTACGTTCCCCACACCCTAACCCTC‑TCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCAc                                               <  1:3313175/150‑1 (MQ=255)
                                                                                                          ttttACATCCCGTACGTTCCCCTCACCCTAAACCTC‑TCCCCAAAGGGGCGAGGGGACCGGTAAGTAACGAAATTAGCATGAAGGTGCTAACCAGCATTAAAAACAGGCAATAAACACACGCAAAATCCGCAGCAACCTCGCTTCACCGCt                                           >  2:9424971/1‑150 (MQ=25)
                                                                                                          ttttACATCCCGTACGTTCCCCGCACCCTAACCCGC‑TCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAAATTAGCTTGACGGTGCGATCCAGCAATAAAAACGGGCCATACTCTCAAGCAAAATCCGCCGCACCGTCTCATCACCGCt                                           >  2:669131/1‑150 (MQ=37)
                                                                                                              aCATCCCGTACGTTCCCCTCACCCTAACCCTC‑TCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGa                                       <  1:458184/150‑1 (MQ=255)
                                                                                                                   ccGTACGTTCCCCTCACCCTAAACCTC‑ACCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACCCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCgg                                  >  2:5505738/1‑150 (MQ=255)
                                                                                                                     gTACGTGCCCCTCCCCCTACCCCTC‑GCCCCCACGGGGCGAGGGGACCGTTCAGTACGGACTTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGt                                <  2:3545466/150‑1 (MQ=255)
                                                                                                                         gTTCCCCTCACCCTAACCCTC‑TCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAAt                            >  2:8876932/1‑150 (MQ=255)
                                                                                                                          ttCCCCTCACCCTAACCCTC‑TCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATg                           <  2:324428/150‑1 (MQ=255)
                                                                                                                               cTCACCCTAACCCTC‑TCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAAt                      >  1:6853433/1‑150 (MQ=255)
                                                                                                                                tCACCCTAACCCTC‑TCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATc                     >  1:9624183/1‑150 (MQ=255)
                                                                                                                                  aCCCTAACCCTC‑TCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGc                   >  2:3970042/1‑150 (MQ=255)
                                                                                                                                   cccTAACCCTC‑TCCCCCCCGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGca                  <  2:2291570/150‑1 (MQ=255)
                                                                                                                                    cggtccccTC‑GCCCCTTTGGGGAGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGcac                 >  2:1255263/6‑150 (MQ=255)
                                                                                                                                     cTAACCCTC‑TCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCAGTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTGTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGcaca                >  1:8341825/1‑150 (MQ=255)
                                                                                                                                       tccccTC‑GCCCCTTTGGGGAGAGGGTACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAAt              >  1:727445/3‑150 (MQ=255)
                                                                                                                                       tccccTC‑GCCCCTTTGGGGAGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAAt              >  1:8707491/3‑150 (MQ=255)
                                                                                                                                        aCCCTC‑TCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATg             >  2:8746031/1‑150 (MQ=255)
                                                                                                                                         cccTC‑TCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGt            >  2:7745564/1‑150 (MQ=255)
                                                                                                                                         cccTC‑TCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGt            >  2:1709116/1‑150 (MQ=255)
                                                                                                                                          ccTC‑TCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGtt           >  2:3397134/1‑150 (MQ=255)
                                                                                                                                          ccTC‑TCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGtt           >  2:5688359/1‑150 (MQ=255)
                                                                                                                                           ctc‑tCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTc          <  1:9336367/150‑1 (MQ=255)
                                                                                                                                           ctc‑tCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTc          >  1:3657444/1‑150 (MQ=255)
                                                                                                                                            tc‑tcCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCt         >  2:2567262/1‑150 (MQ=255)
                                                                                                                                            tc‑tcCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCt         >  2:7888800/1‑150 (MQ=255)
                                                                                                                                                ccccAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCTGAt      <  2:6184747/150‑1 (MQ=255)
                                                                                                                                                 cccAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCTGAtt     >  1:6448198/1‑150 (MQ=255)
                                                                                                                                                  ccAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACCATGTTCTGATtg    >  2:5006838/1‑150 (MQ=255)
                                                                                                                                                  ccAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCTGATtg    >  1:4025990/1‑150 (MQ=255)
                                                                                                                                                  ccAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCTGATtg    <  2:6850485/150‑1 (MQ=255)
                                                                                                                                                    aaaGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCTGATtgtg  >  1:5093860/1‑150 (MQ=255)
                                                                                                                                                    aaaGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCTGATtgtg  >  1:8890513/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TTCTGCCGCGCCAATCGCATCACCTTAACTATTGATAATAAAAGACATTAACCGTAGGCCGGATAAGATGCGCCAGCATCGCATCCGGCGATGCTGGCGCGTTGAATTTTACATCCCGTACGTTCCCCTCACCCTAACCCTC‑TCCCCAAAGGGGCGAGGGGACCGTTCAGTACGGAATTTAGCTTGACGGTGCTATCCAGCATTAAATACTGGCAATACTCTCACGCAAAATCCGCAGCACCGTCTCTTCACCGCTGCGATTCGGGTTAATGCGAATCGCACAATGTTCTGATTGTG  >  NC_000913/3504733‑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: ‑