Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,351,779 A→T 6.1% intergenic (‑19/+49) ycjD ← / ← sapF DUF559 domain‑containing protein YcjD/putrescine ABC exporter ATP binding protein SapF

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,351,7790AT 6.1% 380.8 / 8.1 122intergenic (‑19/+49)ycjD/sapFDUF559 domain‑containing protein YcjD/putrescine ABC exporter ATP binding protein SapF
Reads supporting (aligned to +/- strand):  ref base A (52/62);  new base T (5/3);  total (57/65)
Fisher's exact test for biased strand distribution p-value = 4.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 127 total aligned reads displayed.

CCCACTGGATGTTGACGGCGAAATTTGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCGTTTGGTGAGTTCATGCAGCGGCGATGCCAGCACATC  >  NC_000913/1351630‑1351926
                                                                                                                                                       |                                                                                                                                                   
cccACTGGATGTTGACGGCGAAATTTGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTa                                                                                                                                                     >  2:1388807/1‑150 (MQ=255)
cccACTGGATGTTGACGGCGAAATTTGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTa                                                                                                                                                     >  2:285515/1‑150 (MQ=255)
cccACTGGATGTTGACGGCGAAATTTGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATCGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTa                                                                                                                                                     <  1:1219916/150‑1 (MQ=255)
 ccACTGGATGTTGACGGCGAAATTTGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTaa                                                                                                                                                    >  2:953885/1‑150 (MQ=255)
        aTGTTGACGGCGAAATTTGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTg                                                                                                                                             >  1:578373/1‑150 (MQ=255)
          gTTGACGGCGAAATTTGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCATTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGta                                                                                                                                           >  2:342340/1‑149 (MQ=255)
            tGACGGCGAAATTTGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAg                                                                                                                                         >  2:902569/1‑150 (MQ=255)
            tGACGGCGAAATTTGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAg                                                                                                                                         >  1:1292159/1‑150 (MQ=255)
            tGACGGCGAAATTTGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAg                                                                                                                                         >  1:44113/1‑150 (MQ=255)
            tGACGGCGAAATTTGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAg                                                                                                                                         >  1:1120367/1‑150 (MQ=255)
               cGGCGAAATTTGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAAt                                                                                                                                      <  2:1542677/150‑1 (MQ=255)
               cGGCGAAATTTGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAAt                                                                                                                                      <  1:902569/150‑1 (MQ=255)
                ggCGAAATTTGAAATCACTAAAACGTCGGCTGCGAATATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAAAGATTAAATCTTGTCAGAATg                                                                                                                                     <  2:1538427/150‑1 (MQ=38)
                 gCGAAATTTGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGa                                                                                                                                    <  1:177181/150‑1 (MQ=255)
                  cGAAATTTGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAg                                                                                                                                   <  2:578373/150‑1 (MQ=255)
                  cGAAATTTGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAg                                                                                                                                   >  1:645936/1‑150 (MQ=255)
                  cGAAATTTCAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTACTGCAAAGTGAGATTGCGTCTTAAATCACGCCCATTTCATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAg                                                                                                                                   >  2:1602759/1‑150 (MQ=255)
                    aaaTTTGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGt                                                                                                                                 >  1:1366742/1‑150 (MQ=255)
                     aaTTTGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTg                                                                                                                                <  2:587406/150‑1 (MQ=255)
                     aaTTTGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTg                                                                                                                                <  1:1388807/150‑1 (MQ=255)
                      aTTTGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGa                                                                                                                               >  1:1242935/1‑150 (MQ=255)
                       tttGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGaa                                                                                                                              <  2:332944/150‑1 (MQ=255)
                       tttGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTAAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAAAGATTAAATCTTGTCAGAATGAGGTGat                                                                                                                              <  2:810463/150‑2 (MQ=37)
                         tGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGaat                                                                                                                             >  2:1002765/1‑149 (MQ=255)
                         tGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGaat                                                                                                                             <  1:1002765/149‑1 (MQ=255)
                                aCTAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGa                                                                                                                     <  1:526008/150‑1 (MQ=255)
                                aCTAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGa                                                                                                                     <  1:531451/150‑1 (MQ=255)
                                aCTAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCAATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGa                                                                                                                     <  2:1542672/150‑1 (MQ=255)
                                 cTAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGaa                                                                                                                    <  1:1388071/150‑1 (MQ=255)
                                 cTAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGaa                                                                                                                    >  2:645966/1‑150 (MQ=255)
                                  tAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAag                                                                                                                   <  1:285515/150‑1 (MQ=255)
                                     aaCGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGAATTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAaggag                                                                                                                <  1:540009/150‑1 (MQ=255)
                                       cGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAggagga                                                                                                              <  1:1294167/150‑1 (MQ=255)
                                         tCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGAtt                                                                                                            >  2:64638/1‑150 (MQ=255)
                                           ggCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTTCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCa                                                                                                          <  2:404040/150‑1 (MQ=255)
                                           ggCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCa                                                                                                          <  1:1191324/150‑1 (MQ=255)
                                           ggCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCa                                                                                                          <  1:455655/150‑1 (MQ=255)
                                              tGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGaa                                                                                                       <  1:953885/150‑1 (MQ=255)
                                                     atatCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGc                                                                                                <  1:64638/150‑1 (MQ=255)
                                                      tatCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCg                                                                                               >  2:53787/1‑150 (MQ=255)
                                                          gCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCt                                                                                           >  1:896459/1‑150 (MQ=255)
                                                          gCCAGAGCTTTCGTTCCTGCAAACTGACATTGCGTCTTAAATCACGCGAATATGATTTAATTTTATCAATCACCTCATTCTGAC‑‑AAGATTTTATCTTTTCTCAGAATGATGTGAATAATCTGGAAGGAGGATTCAGAAACTTAGCGATCt                                                                                           <  2:1003553/150‑1 (MQ=255)
                                                                 cTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGcc                                                                                    >  2:506721/1‑150 (MQ=255)
                                                                 cTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGcc                                                                                    <  1:829011/150‑1 (MQ=255)
                                                                 cTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGcc                                                                                    >  1:868526/1‑150 (MQ=255)
                                                                 cTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGcc                                                                                    >  2:539160/1‑150 (MQ=255)
                                                                    tCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACg                                                                                 <  1:23352/150‑1 (MQ=255)
                                                                      gTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCa                                                                               >  1:109910/1‑150 (MQ=255)
                                                                      gTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCa                                                                               <  1:910115/150‑1 (MQ=255)
                                                                         ccTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATcc                                                                            >  2:564837/1‑150 (MQ=255)
                                                                           tGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAaggagg                                                                                                               <  2:804392/113‑1 (MQ=255)
                                                                           tGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAaggagg                                                                                                               >  1:804392/1‑113 (MQ=255)
                                                                             cAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATccgccg                                                                        >  2:1216518/1‑150 (MQ=255)
                                                                               aaGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTc                                                                      >  1:367463/1‑150 (MQ=255)
                                                                                aGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCa                                                                     <  1:1216518/150‑1 (MQ=255)
                                                                                    agaTTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGc                                                                 <  2:790825/150‑1 (MQ=255)
                                                                                      aTTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCt                                                               >  1:1277388/1‑150 (MQ=255)
                                                                                             cTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCaaa                                                        >  2:1418877/1‑150 (MQ=255)
                                                                                              ttAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCaaaa                                                       >  2:405565/1‑150 (MQ=255)
                                                                                               tAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAAt                                                      <  1:365393/150‑1 (MQ=255)
                                                                                               tAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAAt                                                      >  1:1108931/1‑150 (MQ=255)
                                                                                               tAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAAt                                                      <  1:166695/150‑1 (MQ=255)
                                                                                                   tCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGcc                                                           <  1:581506/141‑1 (MQ=255)
                                                                                                   tCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGcc                                                           >  2:581506/1‑141 (MQ=255)
                                                                                                   tCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGAcc                                                  <  1:101155/150‑1 (MQ=255)
                                                                                                      cgcgCATTTGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGc                                               <  2:877737/150‑1 (MQ=255)
                                                                                                      cgcgCATTTGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGc                                               <  2:1120367/150‑1 (MQ=255)
                                                                                                      cgcgCATTTGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGc                                               <  1:1124805/150‑1 (MQ=255)
                                                                                                           aTTTGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCa                                          >  2:42539/1‑150 (MQ=255)
                                                                                                             ttGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGc                                        <  2:1284888/150‑1 (MQ=255)
                                                                                                             ttGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGc                                        >  2:1094906/1‑150 (MQ=255)
                                                                                                               gATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCAACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCg                                      >  2:553631/1‑150 (MQ=255)
                                                                                                                aTTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCGt                                     >  2:1215732/1‑150 (MQ=255)
                                                                                                                   tAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCGTTTg                                  >  1:40156/1‑150 (MQ=255)
                                                                                                                    aaTTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCGTTTgg                                 <  2:645936/150‑1 (MQ=255)
                                                                                                                      ttttATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCGTTTGGTg                               <  1:1094906/150‑1 (MQ=255)
                                                                                                                      ttttATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCGTTTGGTg                               <  1:145707/150‑1 (MQ=255)
                                                                                                                      ttttATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCGTTTGGTg                               <  2:109910/150‑1 (MQ=255)
                                                                                                                      ttttATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCGTTTGGTg                               <  1:53787/150‑1 (MQ=255)
                                                                                                                       tttATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCGTTTGGTGa                              >  2:1088352/1‑150 (MQ=255)
                                                                                                                        ttATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCGTTTGGTGAg                             <  2:425067/150‑1 (MQ=255)
                                                                                                                              ttcACCTCATTCTGAC‑‑AAAAGATTAAATCTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCGTTTGGTGAGTTCATg                       >  1:450902/2‑150 (MQ=38)
                                                                                                                              ttcACCTCATTCTGAC‑‑AAAAGATTAAATCTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCGTTTGGTGAGTTCATg                       >  1:215381/2‑150 (MQ=38)
                                                                                                                               tcACCTCATTCTGAC‑‑AAAAGATTAAATCTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCGTTTGGTGAGTTCATGc                      >  2:384355/1‑150 (MQ=38)
                                                                                                                                cACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCGTTTGGTGAGTTCATGCa                     <  2:983473/150‑1 (MQ=255)
                                                                                                                                 aCCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCGTTTGGTGAGTTCATGCAg                    >  2:1277955/1‑150 (MQ=255)
                                                                                                                                  ccTCATTCTGAC‑‑AAGATTTAATCTTTTCTCAGAATGAGCTCAATAATCTGGAAGCAGGATTCAGAAAATTACCGATCTTTACGCCACGCATCCGCCGTCAATGCAACGCAAAAATGAACGGCAATAACCCGAATCGAGACTTCATGCacc                   >  2:705790/1‑148 (MQ=255)
                                                                                                                                  ccTCATTCTGAC‑‑AAAAGATTAAATCTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGAGTCAGAAAATTAGCGATCTTTAGGCCACGCATCCGCAGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCGTTTGGTGAGTTCCTGCAgc                   >  2:334518/1‑150 (MQ=37)
                                                                                                                                   cTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCGTTTGGTGAGTTCATGCAgcg                  <  2:1083009/150‑1 (MQ=255)
                                                                                                                                     cATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCGTTTGGTGAGTTCATGCAgcggc                >  2:105193/1‑150 (MQ=255)
                                                                                                                                     cATTCTGAC‑‑AAAAGATTAAATCTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCGTTTGGTGAGTTCATGCAgcggc                >  1:431213/1‑150 (MQ=38)
                                                                                                                                      aTTCTGACAAAAGATTAAATC‑‑TTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCGTTTGGTGAGTTCATGCAgcggcg               >  1:1444491/1‑150 (MQ=38)
                                                                                                                                         cTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCGTTTGGTGAGTTCATGCAGCGGCGATg            <  1:1088352/150‑1 (MQ=255)
                                                                                                                                           gAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCGTTTGGTGAGTTCATGCAGCGGCGATGcc          <  1:420836/150‑1 (MQ=255)
                                                                                                                                           gAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCGTTTGGTGAGTTCATGCAGCGGCGATGcc          <  2:1266832/150‑1 (MQ=255)
                                                                                                                                             c‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCGTTTGGTGAGTTCATGCAGCGGCGATGCCAg        >  2:1410469/1‑150 (MQ=255)
                                                                                                                                             c‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCGTTTGGTGAGTTCATGCAGCGGCGATGCCAg        >  2:1228683/1‑150 (MQ=255)
                                                                                                                                                aaGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCGTTTGGTGAGTTCATGCAGCGGCGATGCCAGc       <  2:1366742/150‑1 (MQ=255)
                                                                                                                                                 aGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCGTTTGGTGAGTTCATGCAGCGGCGATGCCAGca      >  1:1387156/1‑150 (MQ=255)
                                                                                                                                                     ttAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCGTTTGGTGAGTTCATGCAGCGGCGATGCCAGCACATc  >  2:422469/1‑150 (MQ=255)
                                                                                                                                                       |                                                                                                                                                   
CCCACTGGATGTTGACGGCGAAATTTGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCGTTTGGTGAGTTCATGCAGCGGCGATGCCAGCACATC  >  NC_000913/1351630‑1351926

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