Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,891,502 C→G 15.5% 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.5% 533.0 / 85.2 195intergenic (+19/‑113)tnaB/mdtLtryptophan transporter of low affinity/multidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base C (65/98);  new base G (13/17);  total (80/115)
Fisher's exact test for biased strand distribution p-value = 8.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.01e-01

Only 100 of 210 total aligned reads displayed.

CGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAG  >  NC_000913/3891404‑3891601
                                                                                                  |                                                                                                   
cGATTGTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAgccagcca                                                                                                    <  2:1961392/100‑1 (MQ=255)
   ttGTCATTCTCTTTGGTATAACCGTGAATTTGTGTTGGTTCGGCAACGTCTTTAACGTGTTACCTAACTTTGGCTCTATATTTCTATAAGGCAGCAATTc                                                                                                 >  2:2725792/1‑100 (MQ=255)
       cATTCTTTTTGGTATAACCGTGATTTTGTGATGGTTGGGCAACGTCTTTAACGTGTTACCTATATTTGGCTAAATCCTTCGAGAACCCAGCCATTCGCTg                                                                                             <  2:2071028/100‑1 (MQ=255)
       cATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTg                                                                                             <  1:594078/100‑1 (MQ=255)
       cATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTg                                                                                             <  1:3438706/100‑1 (MQ=255)
       cATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTg                                                                                             >  2:2680741/1‑100 (MQ=255)
       cATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAAACTTTCAAGAAGCGAGCCGTTc                                                                                                 >  2:4158885/1‑96 (MQ=255)
         ttCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGc                                                                                           <  1:1412054/100‑1 (MQ=255)
          tctcTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGct                                                                                          <  1:2242693/100‑1 (MQ=255)
            tctTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCGAATGGCTGGcttc                                                                                        <  1:3359585/100‑1 (MQ=39)
            tctTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCGAATGGCTGGcttc                                                                                        <  1:2764706/100‑1 (MQ=39)
            tctTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGcttc                                                                                        >  1:308485/1‑100 (MQ=255)
             ctTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCGAATGGCTGGcttct                                                                                       <  2:891431/100‑1 (MQ=39)
             ctTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCGAATGGCTGGcttct                                                                                       <  1:828671/100‑1 (MQ=39)
             ctTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGcttct                                                                                       <  2:252999/100‑1 (MQ=255)
             ctTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGcttct                                                                                       <  2:715321/100‑1 (MQ=255)
              tttGGTATATCCGTGATTTTGTGCTGGTTCGGCAACGTTTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCGAATGGCTGGcttct                                                                                       <  2:3562500/99‑1 (MQ=38)
              tttGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCGAATGGCTGGcttct                                                                                       <  2:1933856/99‑1 (MQ=39)
               ttGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCGAATGGCTGGcttctt                                                                                      <  1:2092345/99‑1 (MQ=38)
               ttGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTg                                                                                     <  1:3201933/100‑1 (MQ=255)
                tGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGc                                                                                    >  2:1765372/1‑100 (MQ=255)
                tGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGACTTCTTGc                                                                                    >  1:3265127/1‑100 (MQ=255)
                 ggTATACCCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCGAATGGCTGGCTTCTTg                                                                                     <  2:3882559/98‑1 (MQ=38)
                 ggTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCGAATGGCTGGCTTCTTgaa                                                                                   <  1:4059480/100‑3 (MQ=38)
                 ggTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGcc                                                                                   <  2:1287781/100‑1 (MQ=255)
                  gTATAACCGTGATTTTGTGCTGGTTCGGCAACGTATTTAACGTGTTACCTAAATTTGGCTTAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTg                                                                                     <  2:1514209/97‑1 (MQ=255)
                   tataACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCGAAATTTGGCTAAATCCTTCAAGAAGCCAGCGAATGGCTGGCTTCTTgaag                                                                                  <  2:1725647/99‑4 (MQ=37)
                     taACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCtctc                                                                               <  1:1576665/100‑1 (MQ=255)
                       aCCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAg                                                                             <  2:585144/100‑1 (MQ=255)
                       aCCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTACTTGCCTCTCAg                                                                             <  1:2739608/100‑1 (MQ=255)
                        ccGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAgg                                                                            >  2:2403091/1‑100 (MQ=255)
                         cGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGa                                                                           <  2:3931409/100‑1 (MQ=255)
                          gTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGaa                                                                          >  1:3316909/1‑100 (MQ=255)
                           tGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGACAGCCATTCGCTGGCTTCTTGCCTCTCAGGaaa                                                                         <  2:1669252/100‑1 (MQ=255)
                            gATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAAt                                                                        <  2:2444706/100‑1 (MQ=255)
                                ttGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACt                                                                    <  2:1359488/100‑1 (MQ=255)
                                 tgtgCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACtt                                                                   <  1:4091246/100‑1 (MQ=255)
                                   tgCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTAt                                                                 <  1:1215221/100‑1 (MQ=255)
                                   tgCTGGTTCGGCAACGTCTTTAACGTGGTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTAt                                                                 >  2:3056671/1‑100 (MQ=255)
                                    gCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATg                                                                >  1:187863/1‑100 (MQ=255)
                                       ggTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTcc                                                             >  2:3855194/1‑100 (MQ=255)
                                        gTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCa                                                            <  2:22871/100‑1 (MQ=255)
                                         ttCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTTTCTTCTTGCCTCTCAGGAAATCACTTATGTCCa                                                            <  2:4081532/99‑1 (MQ=255)
                                           cGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAAt                                                         >  2:3304870/1‑100 (MQ=255)
                                            ggCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATg                                                        >  1:2588511/1‑100 (MQ=255)
                                            ggCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATg                                                        <  1:3533686/100‑1 (MQ=255)
                                            ggCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATg                                                        >  1:1800280/1‑100 (MQ=255)
                                             gCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATgg                                                       <  2:1178838/100‑1 (MQ=255)
                                              caatgTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGc                                                      <  1:1821296/96‑1 (MQ=255)
                                              cAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGc                                                      <  2:811877/100‑1 (MQ=255)
                                                aCGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCaa                                                    <  1:3644847/100‑1 (MQ=255)
                                                    cTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCg                                                <  1:3616254/100‑1 (MQ=255)
                                                    cTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCg                                                <  2:1565167/100‑1 (MQ=255)
                                                      ttAACGTGTTACCTAAATTTGGCTAAATCCTTTAAGAAGGCATCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGcc                                              <  2:2924367/100‑1 (MQ=255)
                                                      ttAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGcc                                              <  1:1952393/100‑1 (MQ=255)
                                                      ttAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGcc                                              <  2:3361042/100‑1 (MQ=255)
                                                       tAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCt                                             <  1:2031170/100‑1 (MQ=255)
                                                       tAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCt                                             >  1:2512448/1‑100 (MQ=255)
                                                           gtgtTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGAtc                                         <  1:3843930/100‑1 (MQ=255)
                                                            tgtTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGAtcc                                        <  1:1709600/100‑1 (MQ=255)
                                                            tgtTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGAtcc                                        <  2:3328350/100‑1 (MQ=255)
                                                               tACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGAtcctcc                                     <  1:985698/100‑1 (MQ=255)
                                                               tACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGAtcctcc                                     >  2:1427963/1‑100 (MQ=255)
                                                                 ccTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCTGGAAATCACTTATGTGCAAATGGCAACTCGCCTGATcctcct                                    <  1:2531063/99‑1 (MQ=255)
                                                                 ccTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCtt                                   >  2:1404549/1‑100 (MQ=255)
                                                                 ccTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCtt                                   <  1:3566362/100‑1 (MQ=255)
                                                                  cTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTc                                  <  2:3789923/100‑1 (MQ=255)
                                                                  cTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTc                                  <  1:1549463/100‑1 (MQ=255)
                                                                  cTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTc                                  >  2:1118580/1‑100 (MQ=255)
                                                                   tAAATTTGGCTAACTCCTTCAAGAAGCCAGCCATCCGCTCGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTca                                 <  2:2877541/100‑1 (MQ=255)
                                                                   tAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTca                                 >  1:192382/1‑100 (MQ=255)
                                                                    aaaTTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTcac                                <  2:2315252/100‑1 (MQ=255)
                                                                      aTTTGGCTAAATCCTTCAAGTAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTCATGTCCAAATGGCAGCTCGCCTGATCCTCCTTcacca                              <  2:1839237/100‑1 (MQ=255)
                                                                       tttGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTcaccac                             <  2:44880/100‑1 (MQ=255)
                                                                        ttGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTGCACCACg                            <  2:3097796/100‑1 (MQ=255)
                                                                           gCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACGCGCCTGATCCTCCTTCACCACGTAt                         >  1:3508972/1‑100 (MQ=255)
                                                                             tAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGc                       >  1:507316/1‑100 (MQ=255)
                                                                              aaaTGCTGCAAGAAGCCAGGCATTCGCTGGCGTCTTGCCTCTCAGGAAATCACTTAGGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGc                       <  2:325680/99‑1 (MQ=255)
                                                                                 tCCTTCAAGAAGCCAGCCATTCGCTGGGTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTg                   <  2:783301/100‑1 (MQ=255)
                                                                                 tCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTg                   >  1:2896891/1‑100 (MQ=255)
                                                                                  gagccAAGAAGCCAGCTATTCGCTGGCTTCTTGCCTCTCAGGACATCCCTTATGTCCAAATATCAACTCCCTTGCTCCCCCTTCACTCCGTATGCCTTGc                  >  2:3947287/5‑100 (MQ=255)
                                                                                   cTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGcn                 <  2:727181/100‑2 (MQ=255)
                                                                                   cTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCg                 >  1:4136755/1‑100 (MQ=255)
                                                                                   cTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCg                 >  1:866366/1‑100 (MQ=255)
                                                                                   aggcAAGAAGCCAGCGAATGGCTGGCTTCTTGCCTCTCAGGAAATCACTTCTGTCCAACTGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGc                  >  2:2156718/4‑99 (MQ=37)
                                                                                    ttCAAGAAGCCAGCCATTCGCTGGCTTCTGGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTGTGCg                 <  2:2521374/99‑1 (MQ=255)
                                                                                     tCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATACTTTGCGt                >  1:2385749/1‑99 (MQ=255)
                                                                                      cAAGAAGCCAGCGAATGGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCa              >  1:3805966/1‑100 (MQ=38)
                                                                                      cAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCa              <  1:421403/100‑1 (MQ=255)
                                                                                      cAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCa              <  1:2008144/100‑1 (MQ=255)
                                                                                       aagaagCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCAc             >  2:186609/1‑100 (MQ=255)
                                                                                        agaagCCAGCGAATGGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCAcc            >  1:1080778/1‑100 (MQ=39)
                                                                                        agaagCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCAcc            <  1:3023229/100‑1 (MQ=255)
                                                                                         gaagCCAGCGAATGGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCt           >  1:3884331/1‑100 (MQ=39)
                                                                                         gaagCCAGCGAATGGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCt           >  1:3663278/1‑100 (MQ=39)
                                                                                          aagCCATCCATTCGCTGGCTTCTTGCCTTTCAGGTAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCGCCACGTATGCTTTGCGTCACCtt          <  2:3224849/100‑1 (MQ=255)
                                                                                            gccagccaTTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTAc        >  1:1192415/1‑100 (MQ=255)
                                                                                              cagccaTTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTa      >  2:662084/1‑100 (MQ=255)
                                                                                                gccaTTCACTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATc    >  1:2186664/1‑100 (MQ=255)
                                                                                                  caTTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTAACCTTACTATCAg  >  2:4002461/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: