Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,891,506 C→G 17.3% intergenic (+23/‑109) 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,5060CG17.3% 368.8 / 65.9 138intergenic (+23/‑109)tnaB/mdtLtryptophan transporter of low affinity/multidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base C (58/56);  new base G (11/13);  total (69/69)
Fisher's exact test for biased strand distribution p-value = 8.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.68e-01

Only 100 of 139 total aligned reads displayed.

GTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACG  >  NC_000913/3891409‑3891605
                                                                                                 |                                                                                                   
gTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGc                                                                                                   >  2:2695694/1‑100 (MQ=255)
  cATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTg                                                                                                 <  1:1010832/100‑1 (MQ=255)
   aTTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTgg                                                                                                >  2:339568/1‑100 (MQ=255)
    ttCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGc                                                                                               <  2:741004/100‑1 (MQ=255)
      ctctTTGGTATAACCGTGATTTTGTGCTGGTTCGGCATCGTCTTTAACGTGTTACCTAACTTTTGCTAAATCCTTCAAACAGCCAGCCATTCATTGGct                                                                                              >  1:1550371/1‑99 (MQ=255)
      ctctTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGctt                                                                                             >  2:1882245/1‑100 (MQ=255)
       tctTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCGAATGGCTGGctt                                                                                             <  1:2673881/99‑1 (MQ=39)
        ctTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCGAATGGCTGGcttct                                                                                           <  2:2356969/100‑1 (MQ=39)
        ctTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCGAATGGCTGGcttct                                                                                           <  2:1943765/100‑1 (MQ=39)
        ctTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCGAATGGCTGGcttct                                                                                           <  1:450313/100‑1 (MQ=39)
          ttGTTATAACCTTGATTTGCTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTATATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTg                                                                                         <  1:2810271/100‑1 (MQ=255)
          ttGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCGAATGGCTGGCTTCTTg                                                                                         <  1:1415096/100‑1 (MQ=39)
           tGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCGAATGGCTGGCTTCTTga                                                                                        <  2:1421383/100‑2 (MQ=38)
           tGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGc                                                                                        <  2:1584438/100‑1 (MQ=255)
           tGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGc                                                                                        <  1:2313606/100‑1 (MQ=255)
           gatcttaACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCGAATGGCTGGCTTCTTga                                                                                        <  1:2745298/95‑2 (MQ=38)
           gatcttaACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCGAATGGCTGGCTTCTTga                                                                                        <  1:222873/95‑2 (MQ=38)
            ggTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCGAATGGCTGGCTTCTTgaa                                                                                       <  1:824852/100‑3 (MQ=38)
            ggTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGcc                                                                                       >  1:2581000/1‑100 (MQ=255)
             gTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCGAGAAGCCAGCCATTCGCTGGCTTCTTGCct                                                                                      >  1:797078/1‑100 (MQ=255)
             gTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCGAATGGCTGGCTTCTTgaag                                                                                      <  1:1237723/100‑4 (MQ=38)
             gTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTgaag                                                                                      >  1:1096434/1‑97 (MQ=255)
              tataACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGCAGCCAGCCATTCGCTCGCTTCTTCCctc                                                                                     >  1:2054611/1‑100 (MQ=255)
              tataACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCGAATGGCTGGCTTCTTgaagg                                                                                     <  2:2615520/100‑5 (MQ=38)
              tataACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCctc                                                                                     <  2:388858/100‑1 (MQ=255)
               ataACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAGATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTgaagga                                                                                    <  2:1096439/100‑6 (MQ=255)
                taACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCtctc                                                                                   <  2:543240/100‑1 (MQ=255)
                  aCCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAg                                                                                 >  1:2137300/1‑100 (MQ=255)
                  aCCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAg                                                                                 >  2:249935/1‑100 (MQ=255)
                   ccGTGATTTTGTGCTGGTTCGGCAAGGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAg                                                                                 <  2:973238/99‑1 (MQ=255)
                   ccGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAgg                                                                                >  2:249025/1‑100 (MQ=255)
                   ccGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAgg                                                                                <  1:1184923/100‑1 (MQ=255)
                   ccGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAgg                                                                                >  2:2249969/1‑100 (MQ=255)
                    cGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGa                                                                               <  1:1353126/100‑1 (MQ=255)
                    cGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGa                                                                               <  1:2195002/100‑1 (MQ=255)
                     gTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGaa                                                                              >  2:2797355/1‑100 (MQ=255)
                     gTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGaa                                                                              >  1:2002231/1‑100 (MQ=255)
                        aTTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATc                                                                           >  2:2816548/1‑100 (MQ=255)
                        aTTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCAAAACATTCGCTGGCTTCTTGCCTCTCAGGAAATc                                                                           >  2:2233227/1‑100 (MQ=255)
                         ttttGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATc                                                                           <  1:1172200/99‑1 (MQ=255)
                         ttttGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCa                                                                          <  2:602671/100‑1 (MQ=255)
                          tttGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCAc                                                                         <  2:1404936/100‑1 (MQ=255)
                            tgtgCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACtt                                                                       <  2:1765275/100‑1 (MQ=255)
                            tgtgCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACtt                                                                       <  2:459255/100‑1 (MQ=255)
                              tgCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTAt                                                                     >  2:125123/1‑100 (MQ=255)
                               gCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATg                                                                    >  2:1733916/1‑100 (MQ=255)
                               gCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATg                                                                    <  2:1605175/100‑1 (MQ=255)
                               gCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATg                                                                    >  2:126640/1‑100 (MQ=255)
                               gCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCCTCTTGCCTCTCAGGAAATCACTTATg                                                                    <  2:1772520/100‑1 (MQ=255)
                                cTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGt                                                                   >  1:822396/1‑100 (MQ=255)
                                  ggTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTcc                                                                 >  2:980041/1‑100 (MQ=255)
                                    ttCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCaa                                                               >  2:2226793/1‑100 (MQ=255)
                                    ttCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCaa                                                               >  1:606009/1‑100 (MQ=255)
                                      cGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAAt                                                             >  1:743681/1‑100 (MQ=255)
                                        gCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATgg                                                           <  2:2137313/100‑1 (MQ=255)
                                        gCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATgg                                                           <  2:1804271/100‑1 (MQ=255)
                                         cAATGTCTGTAACGGGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGc                                                          <  1:2372624/100‑1 (MQ=255)
                                         cAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGc                                                          >  1:1381310/1‑100 (MQ=255)
                                            cGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAAc                                                       <  2:131199/100‑1 (MQ=255)
                                            cGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAAc                                                       <  1:273536/100‑1 (MQ=255)
                                              tCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTc                                                     >  2:178773/1‑100 (MQ=255)
                                                tttAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGc                                                   <  2:542185/100‑1 (MQ=255)
                                                 ttAACGTGTTACCTAAATTTGGCTAAATCCTTCATGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCTACTCGcc                                                  <  1:1529441/100‑1 (MQ=255)
                                                 ttAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGcc                                                  <  1:963426/100‑1 (MQ=255)
                                                  tAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCt                                                 <  2:603299/100‑1 (MQ=255)
                                                     cGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGAt                                              <  2:2367176/100‑1 (MQ=255)
                                                          tACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGAtcctcc                                         <  2:1166248/100‑1 (MQ=255)
                                                           aCCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATcctcct                                        <  2:1038632/100‑1 (MQ=255)
                                                            ccTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCtt                                       >  2:1555003/1‑100 (MQ=255)
                                                             cTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCtt                                       >  1:1542929/1‑99 (MQ=255)
                                                                aaTTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCGCTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTcacc                                   >  1:1873614/1‑100 (MQ=255)
                                                                    tGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGt                               >  1:2285272/1‑100 (MQ=255)
                                                                      gCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTAt                             >  1:1597212/1‑100 (MQ=255)
                                                                      gCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTAt                             >  1:237433/1‑100 (MQ=255)
                                                                       cTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTAt                             <  1:2797337/99‑1 (MQ=255)
                                                                        tAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGc                           >  1:1583521/1‑100 (MQ=255)
                                                                         aaaTCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCt                          >  2:781238/1‑100 (MQ=255)
                                                                         aaaTCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCt                          <  2:1125307/100‑1 (MQ=255)
                                                                              cTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGCAATCACTTATGTCCAAATGGCAACTCGCCTGATGCTCCTTCACCAGGTATGAtttcag                     >  1:1773222/1‑97 (MQ=255)
                                                                              cTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCg                     <  2:2327284/100‑1 (MQ=255)
                                                                               ggcAAGAAGCCAGCGAATGGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCCCCACGTATGCTTTGCGt                    >  1:2572750/3‑100 (MQ=38)
                                                                                gcAAGAAGCCAGCGAATGGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTc                   >  1:2290578/2‑100 (MQ=38)
                                                                                gcAAGAAGCCAGCGAATGGATGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTc                   >  2:770762/2‑100 (MQ=38)
                                                                                 cAAGAAGCCAGCGGATGGCTGGCTTCTTGCCTCTCAGGCAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCa                  >  1:1795131/1‑100 (MQ=38)
                                                                                 cAAGAAGCCAGCGAATGGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCACCTGATCCTCCTTCACCACGTATGCTTTGCGTCa                  >  1:1533820/1‑100 (MQ=38)
                                                                                 cAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCa                  >  2:1256113/1‑100 (MQ=255)
                                                                                 cAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACATATGTCCAAATAGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCa                  >  2:2170460/1‑100 (MQ=255)
                                                                                   agaagCCAGCGAATGGCTGGCTTCTTGCCTCTCAGGAACTCACTTCTGTCCACATGCCAACTCGCCTGCTCCTCCTTCACCACGCCTGCTTTGCGTccc                 >  1:1835290/1‑97 (MQ=25)
                                                                                   agaagCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCAcc                <  1:887032/100‑1 (MQ=255)
                                                                                    gaagCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCt               >  1:1103161/1‑100 (MQ=255)
                                                                                    gaagCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCt               >  2:1423469/1‑100 (MQ=255)
                                                                                      agccagcgAATGGCTGGCTTCTGGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACGAGGTATGCTGTGCGTCACCtt              >  1:1695349/1‑99 (MQ=255)
                                                                                        ccagccaTTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCGTGATCCACCTTCACCACGTATGCTTTGCGTCACCTTACt           >  1:2049137/1‑100 (MQ=255)
                                                                                        ccagccaTTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACt           <  2:2536641/100‑1 (MQ=255)
                                                                                         cagccaTTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCATCACCTTACTa          >  1:934415/1‑100 (MQ=255)
                                                                                           gccaTTCGCTGGCTTCTTGCCCCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATc        >  1:732714/1‑100 (MQ=255)
                                                                                           gccaTGCGCTGGCTTCTTGGGTCTCAGGAAATCGCTTATGTCCAAATGGTAAGTCGCCTGATCCTCCTTCACCAGGTATGCTTTGCGTCGCCTTACTATc        >  1:2434806/1‑100 (MQ=255)
                                                                                            ccaTTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCTAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCa       <  2:1190271/100‑1 (MQ=255)
                                                                                             caTTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCa       >  2:318108/1‑99 (MQ=255)
                                                                                              aTTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAgg     <  1:497369/100‑1 (MQ=255)
                                                                                                 |                                                                                                   
GTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGACG  >  NC_000913/3891409‑3891605

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