Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,351,781 T→A 5.1% intergenic (‑21/+47) 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,7810TA 5.1% 397.9 / 5.3 122intergenic (‑21/+47)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 A (5/2);  total (57/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 = 1.00e+00

Only 100 of 122 total aligned reads displayed.

ATGTTGACGGCGAAATTTGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCGTTTGGTGAGTTCATGCAGCGGCGATGCCAGCACATC  >  NC_000913/1351638‑1351926
                                                                                                                                                 |                                                                                                                                                 
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:1120367/1‑150 (MQ=255)
     gACGGCGAAATTTGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGa                                                                                                                                        >  2:23352/1‑150 (MQ=255)
       cGGCGAAATTTGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAAt                                                                                                                                      <  1:902569/150‑1 (MQ=255)
       cGGCGAAATTTGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAAt                                                                                                                                      <  2:1542677/150‑1 (MQ=255)
         gCGAAATTTGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGa                                                                                                                                    <  1:177181/150‑1 (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                                                                                                                                <  2:587406/150‑1 (MQ=255)
             aaTTTGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTg                                                                                                                                <  1:1388807/150‑1 (MQ=255)
              aTTTGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGa                                                                                                                               >  1:1242935/1‑150 (MQ=255)
                 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                                                                                                                    <  1:1388071/150‑1 (MQ=255)
                          tAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAag                                                                                                                   <  1:285515/150‑1 (MQ=255)
                             aaCGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGAATTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAaggag                                                                                                                <  1:540009/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:1191324/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)
                                             atatCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGc                                                                                                <  1:64638/150‑1 (MQ=255)
                                              tatCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCg                                                                                               >  2:53787/1‑150 (MQ=255)
                                              tatCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCg                                                                                               <  1:645966/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)
                                                   ccAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATGTGGAAGGATGATTCAGAAAATTTGCGATCtt                                                                                          <  2:893502/150‑1 (MQ=255)
                                                       agCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACg                                                                                      >  1:1096189/1‑150 (MQ=255)
                                                         cTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGcc                                                                                    >  2:506721/1‑150 (MQ=255)
                                                         cTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGcc                                                                                    >  1:868526/1‑150 (MQ=255)
                                                         cTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGcc                                                                                    <  1:829011/150‑1 (MQ=255)
                                                         cTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGcc                                                                                    >  2:539160/1‑150 (MQ=255)
                                                            tCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACg                                                                                 <  1:23352/150‑1 (MQ=255)
                                                              gTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCa                                                                               >  1:457515/1‑150 (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)
                                                                            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:166695/150‑1 (MQ=255)
                                                                                       tAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAAt                                                      >  1:1108931/1‑150 (MQ=255)
                                                                                       tAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAAt                                                      <  2:1292159/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                                               <  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:1094906/1‑150 (MQ=255)
                                                                                                      tGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGcc                                       >  2:734813/1‑150 (MQ=255)
                                                                                                       gATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCAACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCg                                      >  2:553631/1‑150 (MQ=255)
                                                                                                        aTTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCGt                                     >  2:1215732/1‑150 (MQ=255)
                                                                                                          ttAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATATGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCGttt                                   >  2:608887/1‑150 (MQ=255)
                                                                                                           tAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCGTTTg                                  <  1:342340/150‑1 (MQ=255)
                                                                                                           tAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCGTTTg                                  >  1:40156/1‑150 (MQ=255)
                                                                                                            aaTTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCGTTTgg                                 <  2:645936/150‑1 (MQ=255)
                                                                                                              ttttATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCGTTTGGTg                               <  2:109910/150‑1 (MQ=255)
                                                                                                              ttttATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCGTTTGGTg                               <  1:145707/150‑1 (MQ=255)
                                                                                                              ttttATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCGTTTGGTg                               <  1:1094906/150‑1 (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)
                                                                                                                       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: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)
                                                                                                                                                 |                                                                                                                                                 
ATGTTGACGGCGAAATTTGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGAC‑‑AAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCGTTTGGTGAGTTCATGCAGCGGCGATGCCAGCACATC  >  NC_000913/1351638‑1351926

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: