Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,891,502 C→G 15.1% intergenic (+19/‑113) tnaB → / → mdtL tryptophan transporter of low affinity/multidrug efflux system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,891,5020CG15.1% 300.8 / 45.0 113intergenic (+19/‑113)tnaB/mdtLtryptophan transporter of low affinity/multidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base C (54/42);  new base G (4/13);  total (58/55)
Fisher's exact test for biased strand distribution p-value = 1.72e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.41e-01

Only 100 of 123 total aligned reads displayed.

CGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAG  >  NC_000913/3891404‑3891601
                                                                                                  |                                                                                                   
cGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAgccagcca                                                                                                    <  2:2256233/100‑1 (MQ=255)
     gTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGGAACGTCTTTTACGTGTTACGTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGc                                                                                               >  1:1411120/1‑100 (MQ=255)
       cATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTg                                                                                             >  2:662514/1‑100 (MQ=255)
         ttCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGc                                                                                           <  1:1231399/100‑1 (MQ=255)
          tctcTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGct                                                                                          <  1:459503/100‑1 (MQ=255)
            tctTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCGAATGGCTGGcttc                                                                                        <  1:2531840/100‑1 (MQ=39)
            tctTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGcttc                                                                                        <  2:746500/100‑1 (MQ=255)
             ctTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCGAATGGCTGGcttct                                                                                       <  2:1651192/100‑1 (MQ=39)
             ctTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAAAAGCCAGCCATTCGCTGGcttct                                                                                       >  1:1136155/1‑100 (MQ=255)
              tttGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCGAATGGCTGGcttctt                                                                                      <  1:1299768/100‑1 (MQ=39)
              tttGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCGAATGGCTGGcttctt                                                                                      <  2:1837391/100‑1 (MQ=39)
                tGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCGAATGGCTGGCTTCTTg                                                                                     <  1:632179/99‑1 (MQ=38)
                tGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGc                                                                                    <  2:1377936/100‑1 (MQ=255)
                tGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGc                                                                                    <  2:1482179/100‑1 (MQ=255)
                tGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGc                                                                                    >  2:2301497/1‑100 (MQ=255)
                tGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGc                                                                                    <  2:854268/100‑1 (MQ=255)
                tGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAAGTTGGCTAAATCCTTCAAGAAGCCAGCGAATGGCTGGCTTCTTga                                                                                    <  1:255242/100‑2 (MQ=38)
                 ggTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGcc                                                                                   <  1:266493/100‑1 (MQ=255)
                 ggTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGcc                                                                                   <  2:2038289/100‑1 (MQ=255)
                  gTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCGAATGGCTGGCTTCTTgaag                                                                                  <  2:2337609/100‑4 (MQ=38)
                  gTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCct                                                                                  >  1:145008/1‑100 (MQ=255)
                   tataACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTGCAAGAAGCCAGCGAATGGCTGGCTTCTTgaagg                                                                                 <  2:1761294/100‑5 (MQ=37)
                      aaCCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCa                                                                              <  2:2128898/100‑1 (MQ=255)
                       aCCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCGTCTCAg                                                                             >  1:612271/1‑100 (MQ=255)
                       aCCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAg                                                                             <  1:2041129/100‑1 (MQ=255)
                       aCCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAg                                                                             >  2:779746/1‑100 (MQ=255)
                        ccGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAgg                                                                            >  1:1983474/1‑100 (MQ=255)
                         cGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGACTCTCAGGa                                                                           >  1:989310/1‑100 (MQ=255)
                         cGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACCTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGa                                                                           <  1:1363700/100‑1 (MQ=255)
                             aTTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATc                                                                       >  2:49197/1‑100 (MQ=255)
                             aTTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATc                                                                       >  2:299742/1‑100 (MQ=255)
                              ttttGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCa                                                                      >  2:15679/1‑100 (MQ=255)
                                ttGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGC‑‑CTCAGGAAATCACTTa                                                                  <  1:20136/100‑1 (MQ=255)
                                      tGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTc                                                              <  1:1092065/100‑1 (MQ=255)
                                      tGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTc                                                              >  1:2178020/1‑100 (MQ=255)
                                       ggatCGGCAACGTCTTTAACGTGATACCTAAATTTGGCTATATCCTTCAAGAAGACCGCCAATCGCGCTCTTCTTGCCTCCCAGGAAATCACt                                                                    <  1:311370/90‑1 (MQ=255)
                                        gTTCGGCAACGTCGTTAACGTGTTACCTAAATTTGGCTACAAACGTCACGAAGACACGCATACGGTGGTCTATTGACTCTCAGGACATGAAATATGt                                                               >  1:2230268/1‑97 (MQ=255)
                                         ttCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCaa                                                           >  2:208353/1‑100 (MQ=255)
                                          tCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCaaa                                                          >  1:507957/1‑100 (MQ=255)
                                          tCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCaaa                                                          <  2:1320306/100‑1 (MQ=255)
                                             gCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATgg                                                       >  2:1488393/1‑100 (MQ=255)
                                             gCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATgg                                                       <  1:2468885/100‑1 (MQ=255)
                                             gCAACGTCTGGAACGTGTTATCTAAATTTTGCTAAATCCGTCAATAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCTGGAAATCACTTATGTCCAAAtt                                                        <  1:204832/99‑2 (MQ=255)
                                              cAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGc                                                      <  2:2443828/100‑1 (MQ=255)
                                                 cGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAAc                                                   >  1:2465772/1‑100 (MQ=255)
                                                   tCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTc                                                 >  1:1615787/1‑100 (MQ=255)
                                                         aCGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGa                                           <  1:299739/100‑1 (MQ=255)
                                                          cGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGAt                                          <  2:1121703/100‑1 (MQ=255)
                                                          cGTGTGACCTAGGTTTGGCTAAATCCTTCAAGAATCCAGCCATTCGCTGGCTGCTTGCCTCTCAGGAAACCACTTATGTCCAAATGGCAAGTCGCCTGa                                           <  1:933404/99‑1 (MQ=255)
                                                            tgtTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGAtcc                                        >  1:2045235/1‑100 (MQ=255)
                                                            tgtTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGAtcc                                        <  2:2311691/100‑1 (MQ=255)
                                                             gtTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATTTCCAAATGGCAACTCGCCTGAtcct                                       >  2:765230/1‑100 (MQ=255)
                                                             gtTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGAtcct                                       >  2:1082932/1‑100 (MQ=255)
                                                               tACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGAtcctcc                                     >  1:1689065/1‑100 (MQ=255)
                                                                 ccTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCtt                                   >  1:1883099/1‑100 (MQ=255)
                                                                 ccTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCtt                                   >  1:322350/1‑100 (MQ=255)
                                                                 ccTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCtt                                   >  1:512696/1‑100 (MQ=255)
                                                                 ccTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCtt                                   >  2:2121173/1‑100 (MQ=255)
                                                                  cTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTc                                  >  2:758587/1‑100 (MQ=255)
                                                                  cTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTc                                  >  1:79519/1‑100 (MQ=255)
                                                                  cTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTc                                  >  2:1599009/1‑100 (MQ=255)
                                                                   tAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTca                                 <  1:806237/100‑1 (MQ=255)
                                                                    aaaTTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTcac                                <  2:2045335/100‑1 (MQ=255)
                                                                        ttGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACg                            <  1:341947/100‑1 (MQ=255)
                                                                          ggCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTa                          <  2:512698/100‑1 (MQ=255)
                                                                          ggCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTa                          <  1:2254341/100‑1 (MQ=255)
                                                                           gCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGTCAACTCGCCTGATCCTCCTTCACCACGTAt                         >  1:53111/1‑100 (MQ=255)
                                                                            cTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATg                        >  2:372697/1‑100 (MQ=255)
                                                                              aaaTCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCt                      >  1:679413/1‑100 (MQ=255)
                                                                                 tCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTg                   <  2:2550607/100‑1 (MQ=255)
                                                                                 tCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTg                   <  2:164996/100‑1 (MQ=255)
                                                                                   cTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCTCCACGTATGCTTTGCg                 <  2:495177/100‑1 (MQ=255)
                                                                                   cTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGc                       <  1:746995/94‑1 (MQ=255)
                                                                                   cTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCg                 >  2:1693538/1‑100 (MQ=255)
                                                                                   cTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCg                 >  2:900094/1‑100 (MQ=255)
                                                                                   cTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCg                 >  1:724091/1‑100 (MQ=255)
                                                                                   cTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCg                 >  2:1216260/1‑100 (MQ=255)
                                                                                   acgcAAGAAGCCAGCGAATGGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCg                 <  2:1120634/97‑1 (MQ=38)
                                                                                     tCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTc               >  2:2401229/1‑100 (MQ=255)
                                                                                     tCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTc               <  1:1691344/100‑1 (MQ=255)
                                                                                     gcAAGAAGCCAACGAATGGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTc               >  1:6459/2‑100 (MQ=38)
                                                                                      cAAGAAGCCAGCGAATGGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCa              >  2:2114114/1‑100 (MQ=38)
                                                                                      cAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCa              <  2:1509957/100‑1 (MQ=255)
                                                                                       aagaagCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCAc             >  2:687513/1‑100 (MQ=255)
                                                                                        agaagCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCTAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCAcc            <  2:1983485/100‑1 (MQ=255)
                                                                                        agaagCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCAcc            >  2:1048519/1‑100 (MQ=255)
                                                                                        agaagCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCAcc            >  1:337110/1‑100 (MQ=255)
                                                                                        agaagCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCAcc            >  1:1830054/1‑100 (MQ=255)
                                                                                         gaagCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCt           <  1:758584/100‑1 (MQ=255)
                                                                                          aagCCAGCGAATGGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCtt          >  2:1446848/1‑100 (MQ=255)
                                                                                          aagACAGCCATTCGCTCGCTTCTTGCCTCTCAGGAAATCACTTATGTCCCAATGGCAACTCGCCTGATCCTCCTTCACCACTTATGCTTTGCGTCACCt           <  1:1578762/99‑1 (MQ=255)
                                                                                            gccagccaTTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTAc        >  1:1165590/1‑100 (MQ=255)
                                                                                              cagccaTTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTa      >  2:972532/1‑100 (MQ=255)
                                                                                              cagccaTTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTa      <  2:1217061/100‑1 (MQ=255)
                                                                                              cagccaTTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTa      <  1:879125/100‑1 (MQ=255)
                                                                                              cagacaTTAGCTAGATTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCGTCACCACGTATGCTTTGCGTCACCTTAc        <  1:467615/98‑1 (MQ=255)
                                                                                                gccaTTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATc    >  1:732896/1‑100 (MQ=255)
                                                                                                gccaTTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATc    <  1:15679/100‑1 (MQ=255)
                                                                                                gccaTTCGCTGGCTTCTTGCCTCCCCGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTCTACCACGTATGCTTTCCGTCACCTTACTATc    >  1:1795654/1‑100 (MQ=255)
                                                                                                  caTTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAg  >  1:1139528/1‑100 (MQ=255)
                                                                                                  |                                                                                                   
CGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAG  >  NC_000913/3891404‑3891601

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