Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,351,776 T→G 5.6% intergenic (‑16/+52) 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,7760TG 5.6% 403.4 / 12.0 126intergenic (‑16/+52)ycjD/sapFDUF559 domain‑containing protein YcjD/putrescine ABC exporter ATP binding protein SapF
Reads supporting (aligned to +/- strand):  ref base T (55/64);  new base G (5/2);  total (60/66)
Fisher's exact test for biased strand distribution p-value = 2.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01

Only 100 of 129 total aligned reads displayed.

TCCCCACTGGATGTTGACGGCGAAATTTGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCGTTTGGTGAGTTCATGCAGCGGCGATGCCAGCA  >  NC_000913/1351628‑1351922
                                                                                                                                                      |                                                                                                                                                  
tCCCCACTGGATGTTGACGGCGAAATTTGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGAtt                                                                                                                                                   >  1:1103469/1‑150 (MQ=255)
tCCCCACTGGATGTTGACGGCGAAATTTGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGAtt                                                                                                                                                   >  2:1124805/1‑150 (MQ=255)
 ccccACTGGATGTTGACGGCGAAATTTGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGAttt                                                                                                                                                  <  1:1617023/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                                                                                                                                     >  1:44113/1‑150 (MQ=255)
              tGACGGCGAAATTTGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAg                                                                                                                                     >  2:902569/1‑150 (MQ=255)
               gACGGCGAAATTTGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGa                                                                                                                                    >  2:23352/1‑150 (MQ=255)
                 cGGCGAAATTTGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAAt                                                                                                                                  <  2:1542677/150‑1 (MQ=255)
                 cGGCGAAATTTGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAAt                                                                                                                                  <  1:902569/150‑1 (MQ=255)
                 cGGCGAAATTTGAAACCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAAt                                                                                                                                  <  2:512792/150‑1 (MQ=255)
                  ggCGAAATTTGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATg                                                                                                                                 <  1:180227/150‑1 (MQ=255)
                   gCGAAATTTGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGa                                                                                                                                <  1:177181/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                                                                                                                            <  1:1265129/150‑1 (MQ=255)
                       aaTTTGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTg                                                                                                                            <  1:1388807/150‑1 (MQ=255)
                       aaTTTGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTg                                                                                                                            <  2:587406/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                                                                                                                         <  1:1002765/149‑1 (MQ=255)
                           tGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGaat                                                                                                                         >  2:1002765/1‑149 (MQ=255)
                                tCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTg                                                                                                                   <  1:1144242/150‑1 (MQ=255)
                                 cACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTgg                                                                                                                  >  1:567088/1‑150 (MQ=255)
                                  acaaaacCGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGAGTGCGTCTTAAATCACGCCCATTTGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGa                                                                                                                 <  2:570887/147‑1 (MQ=255)
                                  aCTAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGa                                                                                                                 <  1:531451/150‑1 (MQ=255)
                                  aCTAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGa                                                                                                                 <  1:526008/150‑1 (MQ=255)
                                  aCTAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCAATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGa                                                                                                                 <  2:1542672/150‑1 (MQ=255)
                                   cTAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGaa                                                                                                                >  2:645966/1‑150 (MQ=255)
                                    tAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAag                                                                                                               <  1:285515/150‑1 (MQ=255)
                                      aaaCGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAagga                                                                                                             <  1:871971/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)
                                             ggCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCa                                                                                                      <  1:1321759/150‑1 (MQ=255)
                                              gCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAg                                                                                                     <  2:1055027/150‑1 (MQ=255)
                                                     aGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTa                                                                                              <  2:1103469/150‑1 (MQ=255)
                                                       atatCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGc                                                                                            <  1:64638/150‑1 (MQ=255)
                                                           cGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGAta                                                                                        >  1:425067/1‑149 (MQ=255)
                                                            gCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCt                                                                                       >  1:896459/1‑150 (MQ=255)
                                                            gCCAGAGCTTTCGTTCCTGCAAACTGACATTGCGTCTTAAATCACGCGAATATGATTTAATTTTATCAATCACCTCATTCTGAC‑‑AAGATTTTATCTTTTCTCAGAATGATGTGAATAATCTGGAAGGAGGATTCAGAAACTTAGCGATCt                                                                                       <  2:1003553/150‑1 (MQ=255)
                                                                   cTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGcc                                                                                <  1:829011/150‑1 (MQ=255)
                                                                   cTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGcc                                                                                >  2:506721/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)
                                                                             tGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATccgc                                                                      >  2:1251067/1‑150 (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)
                                                                                        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:983473/1‑150 (MQ=255)
                                                                                                 tAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAAt                                                  <  1:166695/150‑1 (MQ=255)
                                                                                                 tAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAAt                                                  <  2:1292159/150‑1 (MQ=255)
                                                                                                 tAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAAt                                                  >  1:1108931/1‑150 (MQ=255)
                                                                                                     tCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGcc                                                       >  2:581506/1‑141 (MQ=255)
                                                                                                     tCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGcc                                                       <  1:581506/141‑1 (MQ=255)
                                                                                                     tCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGAcc                                              <  1:101155/150‑1 (MQ=255)
                                                                                                        cgcgCATTTGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGc                                           <  2:877737/150‑1 (MQ=255)
                                                                                                        cgcgCATTTGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGc                                           <  1:1124805/150‑1 (MQ=255)
                                                                                                        cgcgCATTTGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGc                                           <  2:1120367/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)
                                                                                                                tGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGcc                                   >  2:734813/1‑150 (MQ=255)
                                                                                                                tGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGcc                                   <  1:1414701/150‑1 (MQ=255)
                                                                                                                 gATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCAACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCg                                  >  2:553631/1‑150 (MQ=255)
                                                                                                                    ttAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATATGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCGttt                               >  2:608887/1‑150 (MQ=255)
                                                                                                                     tAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCGTTTg                              >  1:40156/1‑150 (MQ=255)
                                                                                                                     tAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCGTTTg                              <  1:342340/150‑1 (MQ=255)
                                                                                                                        ttttATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCGTTTGGTg                           <  1:53787/150‑1 (MQ=255)
                                                                                                                        ttttATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCGTTTGGTg                           <  1:1094906/150‑1 (MQ=255)
                                                                                                                         tttATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCGTTTGGTGa                          >  2:1088352/1‑150 (MQ=255)
                                                                                                                                ttcACCTCATTCTGAC‑‑AAAAGATTAAATCTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCGTTTGGTGAGTTCATg                   >  1:450902/2‑150 (MQ=38)
                                                                                                                                atcACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCGTTTGGTGAGTTCATg                   <  2:1506785/150‑1 (MQ=255)
                                                                                                                                 tcACCTCATTCTGAC‑‑AAAAGATTAAATCTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCGTTTGGTGAGTTCATGc                  >  2:384355/1‑150 (MQ=38)
                                                                                                                                  cACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCGTTTGGTGAGTTCATGCa                 <  2:983473/150‑1 (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)
                                                                                                                                        aTTCTGACAAAAGATTAAATC‑‑TTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCGTTTGGTGAGTTCATGCAgcggcg           >  1:1444491/1‑150 (MQ=38)
                                                                                                                                           cTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCGTTTGGTGAGTTCATGCAGCGGCGATg        <  1:1088352/150‑1 (MQ=255)
                                                                                                                                            tGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCGTTTGGTGAGTTCATGCAGCGGCGATGc       >  2:412957/1‑150 (MQ=255)
                                                                                                                                             gAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCGTTTGGTGAGTTCATGCAGCGGCGATGcc      <  1:420836/150‑1 (MQ=255)
                                                                                                                                             gAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCGTTTGGTGAGTTCATGCAGCGGCGATGcc      <  2:1266832/150‑1 (MQ=255)
                                                                                                                                               c‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCGTTTGGTGAGTTCATGCAGCGGCGATGCCAg    >  2:1228683/1‑150 (MQ=255)
                                                                                                                                               c‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCGTTTGGTGAGTTCATGCAGCGGCGATGCCAg    >  2:1410469/1‑150 (MQ=255)
                                                                                                                                                  aaGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCGTTTGGTGAGTTCATGCAGCGGCGATGCCAGc   <  2:1366742/150‑1 (MQ=255)
                                                                                                                                                   aGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCGTTTGGTGAGTTCATGCAGCGGCGATGCCAGca  >  1:1387156/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                  
TCCCCACTGGATGTTGACGGCGAAATTTGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCGTTTGGTGAGTTCATGCAGCGGCGATGCCAGCA  >  NC_000913/1351628‑1351922

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