Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,891,504 T→A 17.4% intergenic (+21/‑111) 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,5040TA17.4% 333.1 / 58.7 137intergenic (+21/‑111)tnaB/mdtLtryptophan transporter of low affinity/multidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base T (58/55);  new base A (11/13);  total (69/68)
Fisher's exact test for biased strand distribution p-value = 6.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.66e-01

Only 100 of 138 total aligned reads displayed.

GTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGA  >  NC_000913/3891409‑3891603
                                                                                               |                                                                                                   
gTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGc                                                                                                 >  2:2695694/1‑100 (MQ=255)
  cATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTg                                                                                               <  1:1010832/100‑1 (MQ=255)
  cATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTg                                                                                               <  1:1554568/100‑1 (MQ=255)
   aTTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTgg                                                                                              <  1:2637773/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                                                                                         <  1:450313/100‑1 (MQ=39)
        ctTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCGAATGGCTGGcttct                                                                                         <  2:2356969/100‑1 (MQ=39)
        ctTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCGAATGGCTGGcttct                                                                                         <  2:1943765/100‑1 (MQ=39)
          ttGTTATAACCTTGATTTGCTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTATATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTg                                                                                       <  1:2810271/100‑1 (MQ=255)
          ttGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCGAATGGCTGGCTTCTTg                                                                                       <  1:1415096/100‑1 (MQ=39)
          ttGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTg                                                                                       >  1:2042011/1‑100 (MQ=255)
           tGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCGAATGGCTGGCTTCTTga                                                                                      <  2:1421383/100‑2 (MQ=38)
           tGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCGAATGGCTGGCTTCTTga                                                                                      <  2:1097864/100‑2 (MQ=38)
           tGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGc                                                                                      <  2:1584438/100‑1 (MQ=255)
           tGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGc                                                                                      <  1:2313606/100‑1 (MQ=255)
           tGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGc                                                                                      >  1:616368/1‑100 (MQ=255)
           gatcttaACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCGAATGGCTGGCTTCTTga                                                                                      <  1:222873/95‑2 (MQ=38)
           gatcttaACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCGAATGGCTGGCTTCTTga                                                                                      <  1:2745298/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)
             gTATAACCGTGATTTTGTGCTGGTTAGGCAACGTATTTAACGTGATACCTAAAGTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGCCTTCTTGcc                                                                                     <  1:237788/99‑1 (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)
                taACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCGAATGGCTGGCTTCTTgaagatc                                                                                 >  2:222874/1‑94 (MQ=37)
                taACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCtctc                                                                                 <  2:543240/100‑1 (MQ=255)
                  aCCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAg                                                                               >  2:249935/1‑100 (MQ=255)
                  aCCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAg                                                                               >  1:2137300/1‑100 (MQ=255)
                   ccGTGATTTTGTGCTGGTTCGGCAAGGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAg                                                                               <  2:973238/99‑1 (MQ=255)
                   ccGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAgg                                                                              >  2:2249969/1‑100 (MQ=255)
                   ccGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAgg                                                                              >  2:249025/1‑100 (MQ=255)
                   ccGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAgg                                                                              <  1:1184923/100‑1 (MQ=255)
                    cGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGa                                                                             <  1:2195002/100‑1 (MQ=255)
                    cGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGa                                                                             <  1:1353126/100‑1 (MQ=255)
                     gTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGaa                                                                            >  1:2002231/1‑100 (MQ=255)
                     gTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGaa                                                                            >  2:2797355/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:1489908/1‑100 (MQ=255)
                         ttttGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCa                                                                        <  2:602671/100‑1 (MQ=255)
                          tttGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCAc                                                                       <  2:1404936/100‑1 (MQ=255)
                            tgtgCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACtt                                                                     <  2:459255/100‑1 (MQ=255)
                              tgCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTAt                                                                   >  2:125123/1‑100 (MQ=255)
                               gCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTAt                                                                   >  2:126525/1‑99 (MQ=255)
                               gCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATg                                                                  >  2:1733916/1‑100 (MQ=255)
                               gCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATg                                                                  >  2:126640/1‑100 (MQ=255)
                               gCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCCTCTTGCCTCTCAGGAAATCACTTATg                                                                  <  2:1772520/100‑1 (MQ=255)
                                  ggTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTcc                                                               >  2:980041/1‑100 (MQ=255)
                                   gTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCa                                                              <  2:1348343/100‑1 (MQ=255)
                                    ttCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCaa                                                             >  1:606009/1‑100 (MQ=255)
                                       ggCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATg                                                          >  1:859762/1‑100 (MQ=255)
                                        gCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATgg                                                         <  2:2137313/100‑1 (MQ=255)
                                         cAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGc                                                        >  1:1381310/1‑100 (MQ=255)
                                          aaCGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCa                                                       <  1:777990/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                                                <  2:337843/100‑1 (MQ=255)
                                                 ttAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGcc                                                <  1:963426/100‑1 (MQ=255)
                                                  tAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCt                                               <  2:603299/100‑1 (MQ=255)
                                                       tgtTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGAtcc                                          <  2:2612190/100‑1 (MQ=255)
                                                            ccTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCtt                                     <  1:339127/100‑1 (MQ=255)
                                                            ccTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCtt                                     >  2:1555003/1‑100 (MQ=255)
                                                            ccTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCtt                                     >  1:1821388/1‑100 (MQ=255)
                                                             cTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCtt                                     >  1:1542929/1‑99 (MQ=255)
                                                             cTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTc                                    >  1:2739622/1‑100 (MQ=255)
                                                              tAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTca                                   >  2:1771147/1‑100 (MQ=255)
                                                                aaTTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCGCTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTcacc                                 >  1:1873614/1‑100 (MQ=255)
                                                                aaTTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTcacc                                 <  2:2471526/100‑1 (MQ=255)
                                                                aaTTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTcacc                                 >  2:183905/1‑100 (MQ=255)
                                                                      gCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTAt                           >  1:237433/1‑100 (MQ=255)
                                                                      gCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTAt                           >  1:1597212/1‑100 (MQ=255)
                                                                        tAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGc                         >  2:1636678/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)
                                                                                tCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTc                 <  2:1188765/100‑1 (MQ=255)
                                                                                gcAAGAAGCCAGCGAATGGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTc                 >  1:2290578/2‑100 (MQ=38)
                                                                                gcAAGAAGCCAGCGAATGGATGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTc                 >  2:770762/2‑100 (MQ=38)
                                                                                 cAAGAAGCCAGCGGATGGCTGGCTTCTTGCCTCTCAGGCAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCa                >  1:1795131/1‑100 (MQ=38)
                                                                                 cAAGAAGCCAGCGAATGGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCa                >  1:2251071/1‑100 (MQ=38)
                                                                                 cAAGAAGCCAGCGAATGGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCACCTGATCCTCCTTCACCACGTATGCTTTGCGTCa                >  1:1533820/1‑100 (MQ=38)
                                                                                 cAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACATATGTCCAAATAGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCa                >  2:2170460/1‑100 (MQ=255)
                                                                                   agaagCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCAcc              <  1:887032/100‑1 (MQ=255)
                                                                                    gaagCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCt             >  1:1103161/1‑100 (MQ=255)
                                                                                      agccagcgAATGGCTGGCTTCTGGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACGAGGTATGCTGTGCGTCACCtt            >  1:1695349/1‑99 (MQ=255)
                                                                                         cagccaTTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTa        >  1:942556/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)
                                                                                               |                                                                                                   
GTCATTCTCTTTGGTATAACCGTGATTTTGTGCTGGTTCGGCAACGTCTTTAACGTGTTACCTAAATTTGGCTAAATCCTTCAAGAAGCCAGCCATTCGCTGGCTTCTTGCCTCTCAGGAAATCACTTATGTCCAAATGGCAACTCGCCTGATCCTCCTTCACCACGTATGCTTTGCGTCACCTTACTATCAGGA  >  NC_000913/3891409‑3891603

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: