Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,351,784 T→C 5.7% intergenic (‑24/+44) 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,7840TC 5.7% 379.7 / 9.1 123intergenic (‑24/+44)ycjD/sapFDUF559 domain‑containing protein YcjD/putrescine ABC exporter ATP binding protein SapF
Reads supporting (aligned to +/- strand):  ref base T (52/63);  new base C (5/2);  total (58/65)
Fisher's exact test for biased strand distribution p-value = 2.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.65e-01

Only 100 of 124 total aligned reads displayed.

ATGTTGACGGCGAAATTTGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCGTTTGGTGAGTTCATGCAGCGGCGATGCCAGCACATCCGCGGT  >  NC_000913/1351638‑1351932
                                                                                                                                                    |                                                                                                                                                    
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:1120367/1‑150 (MQ=255)
    tGACGGCGAAATTTGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAg                                                                                                                                               >  1:44113/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)
        ggCGAAATTTGAAATCACTAAAACGTCGGCTGCGAATATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAAAGATTAAATCTTGTCAGAATg                                                                                                                                           <  2:1538427/150‑1 (MQ=38)
        ggCGAAATTTGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATg                                                                                                                                           <  1:180227/150‑1 (MQ=255)
         gCGAAATTTGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGa                                                                                                                                          <  1:177181/150‑1 (MQ=255)
          cGAAATTTGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAg                                                                                                                                         <  2:578373/150‑1 (MQ=255)
          cGAAATTTGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAg                                                                                                                                         >  1:645936/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)
              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)
                      tCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTg                                                                                                                             <  1:1144242/150‑1 (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)
                            aaaCGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAagga                                                                                                                       <  1:871971/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:1321759/150‑1 (MQ=255)
                                   ggCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCa                                                                                                                <  1:455655/150‑1 (MQ=255)
                                    gCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAg                                                                                                               <  2:1055027/150‑1 (MQ=255)
                                      tGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGaa                                                                                                             <  1:953885/150‑1 (MQ=255)
                                           aGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTa                                                                                                        <  2:1103469/150‑1 (MQ=255)
                                              tatCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCg                                                                                                     >  2:53787/1‑150 (MQ=255)
                                              tatCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCg                                                                                                     <  1:645966/150‑1 (MQ=255)
                                                  gCCAGAGCTTTCGTTCCTGCAAACTGACATTGCGTCTTAAATCACGCGAATATGATTTAATTTTATCAATCACCTCATTCTGAC‑‑AAGATTTTATCTTTTCTCAGAATGATGTGAATAATCTGGAAGGAGGATTCAGAAACTTAGCGATCt                                                                                                 <  2:1003553/150‑1 (MQ=255)
                                                   ccAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATGTGGAAGGATGATTCAGAAAATTTGCGATCtt                                                                                                <  2:893502/150‑1 (MQ=255)
                                                       agCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACg                                                                                            >  1:1096189/1‑150 (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)
                                                         cTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGcc                                                                                          >  1:868526/1‑150 (MQ=255)
                                                            tCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACg                                                                                       <  1:23352/150‑1 (MQ=255)
                                                              gTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCa                                                                                     <  1:910115/150‑1 (MQ=255)
                                                                 ccTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATcc                                                                                  >  2:564837/1‑150 (MQ=255)
                                                                   tGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAaggagg                                                                                                                     >  1:804392/1‑113 (MQ=255)
                                                                   tGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAaggagg                                                                                                                     <  2:804392/113‑1 (MQ=255)
                                                                   tGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATccgc                                                                                >  2:1251067/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)
                                                                                      ttAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCaaaa                                                             >  2:405565/1‑150 (MQ=255)
                                                                                       tAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAAt                                                            <  2:1292159/150‑1 (MQ=255)
                                                                                       tAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAAt                                                            >  1:1108931/1‑150 (MQ=255)
                                                                                       tAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAAt                                                            >  1:983473/1‑150 (MQ=255)
                                                                                       tAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAAt                                                            <  1:365393/150‑1 (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                                                     <  1:1124805/150‑1 (MQ=255)
                                                                                              cgcgCATTTGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGc                                                     <  2:1120367/150‑1 (MQ=255)
                                                                                                     ttGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGc                                              <  2:1284888/150‑1 (MQ=255)
                                                                                                      tGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGcc                                             <  1:1414701/150‑1 (MQ=255)
                                                                                                      tGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGcc                                             >  2:734813/1‑150 (MQ=255)
                                                                                                        aTTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCGt                                           >  2:1215732/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)
                                                                                                            aaTTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCGTTTgg                                       <  2:645936/150‑1 (MQ=255)
                                                                                                              ttttATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCGTTTGGTg                                     <  1:145707/150‑1 (MQ=255)
                                                                                                              ttttATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCGTTTGGTg                                     <  2:109910/150‑1 (MQ=255)
                                                                                                              ttttATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCGTTTGGTg                                     <  1:1094906/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:215381/2‑150 (MQ=38)
                                                                                                                      ttcACCTCATTCTGAC‑‑AAAAGATTAAATCTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCGTTTGGTGAGTTCATg                             >  1:450902/2‑150 (MQ=38)
                                                                                                                      atcACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCGTTTGGTGAGTTCATg                             <  2:1506785/150‑1 (MQ=255)
                                                                                                                        cACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCGTTTGGTGAGTTCATGCa                           <  2:983473/150‑1 (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)
                                                                                                                                  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)
                                                                                                                                             ttAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCGTTTGGTGAGTTCATGCAGCGGCGATGCCAGCACATc        >  2:422469/1‑150 (MQ=255)
                                                                                                                                                  cTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCGTTTGGTGAGTTCATGCAGCGGCGATGCCAGCACATCCGCgg   >  2:982694/1‑150 (MQ=255)
                                                                                                                                                   ttttGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCGTTTGGTGAGTTCATGCAGCGGCGATGCCAGCACATCCGCGGt  <  2:1242935/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
ATGTTGACGGCGAAATTTGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCGTTTGGTGAGTTCATGCAGCGGCGATGCCAGCACATCCGCGGT  >  NC_000913/1351638‑1351932

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: