Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,950,467 G→T 16.0% intergenic (+47/‑40) ilvL → / → ilvX ilvG operon leader peptide/uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,950,4670GT16.0% 311.0 / 34.5 114intergenic (+47/‑40)ilvL/ilvXilvG operon leader peptide/uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (42/53);  new base T (13/6);  total (55/59)
Fisher's exact test for biased strand distribution p-value = 7.75e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 119 total aligned reads displayed.

TGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGAACTAACTA  >  NC_000913/3950368‑3950560
                                                                                                    |                                                                                              
tggtgATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCgggg                                                                                                >  2:805353/1‑100 (MQ=255)
tggtgATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCgggg                                                                                                >  2:2329496/1‑100 (MQ=255)
 ggtgATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCgggg                                                                                                >  1:1695301/1‑99 (MQ=255)
  gtaattattATCCCACAGTGCGGGGCTGCACTTGGCCGAGAAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCgggg                                                                                                <  1:2292972/95‑1 (MQ=255)
   tgATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGttt                                                                                            >  2:2968020/1‑100 (MQ=255)
    gATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTCACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGtttt                                                                                           >  2:1405556/1‑100 (MQ=255)
    gATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGtttt                                                                                           <  2:662167/100‑1 (MQ=255)
      ttattaTCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTACGGGGCtttt                                                                                          >  2:2482748/1‑99 (MQ=255)
       tattatCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTtttttt                                                                                        <  2:595587/100‑1 (MQ=255)
       tattatCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCC‑‑GGGtttttttt                                                                                      <  2:1148548/100‑1 (MQ=255)
         ttatCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTtttttttt                                                                                      <  1:805350/100‑1 (MQ=255)
         ttatCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTg                                                                                     <  2:2460014/100‑1 (MQ=255)
          tatCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTg                                                                                     <  1:2409390/100‑1 (MQ=255)
          tatCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGAGCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTtttttttt                                                                                      >  1:1703706/1‑99 (MQ=255)
           ataccACGGTGCGGGGCTGCACTTGGACGAGGAACGGCTTAGAGATCAAGGCTTAACTAACTAAGCCCCCCGCGCAGAAAGGTCCGCGGTTTTTTCTTGa                                                                                    <  1:1474803/97‑1 (MQ=255)
           atCCCACCGTGCGGGGCTGCACTTAGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGa                                                                                    <  1:300763/100‑1 (MQ=255)
            tCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGAc                                                                                   >  2:3012426/1‑100 (MQ=255)
            tCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGAc                                                                                   >  1:2359614/1‑100 (MQ=255)
            tCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGAc                                                                                   >  2:371941/1‑100 (MQ=255)
             cccACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGAcc                                                                                  <  2:1433289/100‑1 (MQ=255)
              ccACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCt                                                                                 >  2:3005472/1‑100 (MQ=255)
               cACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGTTTTTTTTTTGACCtt                                                                                >  1:1082756/1‑100 (MQ=255)
               cACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCtt                                                                                <  1:1674943/100‑1 (MQ=255)
                aCCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGTTTTTTTTTTGACCTTa                                                                               >  2:1392076/1‑100 (MQ=255)
                    tGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTaaaaa                                                                           <  2:2134888/100‑1 (MQ=255)
                     gCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAAc                                                                          <  2:1626563/100‑1 (MQ=255)
                        gggCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATa                                                                       >  2:491403/1‑100 (MQ=255)
                          gCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATaa                                                                      >  1:1740829/1‑99 (MQ=255)
                          gCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACATAAc                                                                     <  1:1596425/100‑1 (MQ=255)
                           cTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAAcc                                                                    >  2:502529/1‑100 (MQ=255)
                            tGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACATAACCg                                                                   >  1:942250/1‑100 (MQ=255)
                             gCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGTTTTTTTTTTGACCTTAAAAACATAACCga                                                                  >  2:2514956/1‑100 (MQ=255)
                               aCTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACGAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCgagg                                                                <  2:2123023/100‑1 (MQ=255)
                                cTTGGACGAGGAAAGGCGTAGAGATCAAGCCTTAACGAACGAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCgagn                                                                <  2:845411/99‑2 (MQ=255)
                                 ttGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCgaggag                                                              <  2:3082426/100‑1 (MQ=255)
                                 ttGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCgaggag                                                              <  1:143565/100‑1 (MQ=255)
                                  tGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGc                                                             <  2:2709603/100‑1 (MQ=255)
                                  tGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGc                                                             <  2:2049622/100‑1 (MQ=255)
                                   ggACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGATGAGCa                                                            >  1:358113/1‑100 (MQ=255)
                                    gACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGTTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAg                                                           >  2:895938/1‑100 (MQ=255)
                                    gACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAg                                                           >  2:355102/1‑100 (MQ=255)
                                    gACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAg                                                           >  2:1613200/1‑100 (MQ=255)
                                     aggAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGa                                                          <  2:335584/98‑1 (MQ=255)
                                     aCGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGa                                                          <  2:888713/100‑1 (MQ=255)
                                         ggAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGTTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAAt                                                      >  2:2845314/1‑100 (MQ=255)
                                         ggAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAAt                                                      <  2:2911954/100‑1 (MQ=255)
                                          ggaaGGCTTCGAGAGCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATg                                                     <  2:133016/98‑1 (MQ=255)
                                          gAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATg                                                     <  2:1766529/100‑1 (MQ=255)
                                              ggCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGTTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACTATGAAt                                                  >  1:182996/1‑99 (MQ=255)
                                              ggCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGaa                                                   >  1:1699905/1‑98 (MQ=255)
                                              ggCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATa                                                 <  1:1127551/100‑1 (MQ=255)
                                               gCTGAGAGATCATGCCTTAACGAACTAGGACCCCCGCACCGAATGGGCCGGGTTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATa                                                 <  1:2958836/99‑1 (MQ=255)
                                                cgtAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAAc                                               <  2:262345/98‑1 (MQ=255)
                                                cTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAAc                                               <  2:2450674/100‑1 (MQ=255)
                                                 ttAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACa                                              <  2:833518/100‑1 (MQ=255)
                                                 ttAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAAACGAGGAGCAGACAATGAATAAc                                               >  1:800963/1‑99 (MQ=255)
                                                   agagATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAg                                             >  1:1178778/1‑99 (MQ=255)
                                                     agaTCAAGCCTTAACGAACTAAGGCCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGcac                                          <  2:3168544/100‑1 (MQ=255)
                                                      gaTCAAGACTTAACGAACTAAGACCCACGCACCGAAAGGTCCGAAGTTGTCTTTTGACCTTAAAAACATAACCGAGTAGTAGACAATGAATAACAGctta                                         <  1:3171591/100‑4 (MQ=255)
                                                      cttCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGcaca                                         <  2:2620721/98‑1 (MQ=255)
                                                        tCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACaag                                       >  1:2620701/1‑99 (MQ=255)
                                                         cAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACaaaa                                      >  2:1997280/1‑100 (MQ=255)
                                                            gCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACCGCACAAAATTc                                   >  2:710079/1‑100 (MQ=255)
                                                            gCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTc                                   >  2:3246307/1‑100 (MQ=255)
                                                            gCCTTAACGAACCAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTc                                   <  2:21107/100‑1 (MQ=255)
                                                             ccTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTc                                   >  1:1933490/1‑99 (MQ=255)
                                                             ccTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCt                                  >  2:3282265/1‑100 (MQ=255)
                                                              cTTAACGACCGAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCt                                  <  2:1634627/99‑1 (MQ=255)
                                                               ttAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGt                                >  2:1708513/1‑100 (MQ=255)
                                                               ttAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTCTTTTTTTGACCTTAAAAACTTAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGt                                <  1:1841094/100‑1 (MQ=255)
                                                               gaaactaactAAGCCCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGt                                <  2:2644862/98‑1 (MQ=255)
                                                                tAACGATCTAAGACCCCCGAACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGt                                <  1:727755/99‑1 (MQ=255)
                                                                  aCGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTtc                             >  2:2181933/1‑100 (MQ=255)
                                                                  aCGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTtc                             <  1:1118156/100‑1 (MQ=255)
                                                                   tgAATTAAGGCCCCCGCACCGAAAGGTCCGGGTTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTtc                             <  1:764308/98‑1 (MQ=255)
                                                                   cGAACTAAGACCCCCGCACCGAAAGGTCCGGGTTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTtct                            >  2:2184030/1‑100 (MQ=255)
                                                                   cGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGATCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTtct                            >  1:1912429/1‑100 (MQ=255)
                                                                    gAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTtctc                           >  2:316607/1‑100 (MQ=255)
                                                                     aaCTTAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTtc                             <  1:1117059/97‑1 (MQ=255)
                                                                      aCTAAGACCCCCGCACCGAAAGGTCCGGGTTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCaa                         >  2:1550764/1‑100 (MQ=255)
                                                                      aCTAAGACCCCCGCACCGAAAGGTCCGGGTTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCaa                         >  2:1869651/1‑100 (MQ=255)
                                                                      aCTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCaa                         >  1:2898225/1‑100 (MQ=255)
                                                                       cTAAGACCCCCGCACCGAAAGGTCCGGGTTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTaaa                         >  1:1470076/1‑99 (MQ=255)
                                                                       cTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAg                        <  1:39200/100‑1 (MQ=255)
                                                                        tgtgttccccGTACCGACAGGTCCGGGTTTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTtct                            <  1:2040269/90‑1 (MQ=255)
                                                                        tAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGa                       >  2:2595936/1‑100 (MQ=255)
                                                                        gaaGACCCCCCCCCCAAAGGGTCGGGGGTTTTTTTTTGACATCAAAAACATAACCAAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGa                       <  2:3199617/99‑1 (MQ=255)
                                                                             cccccGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAg                  <  1:1110126/100‑1 (MQ=255)
                                                                              ccccGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAgg                 <  2:919527/100‑1 (MQ=255)
                                                                               cccGCACCGAAAGGTCCGGGGTTTTTTTTTTTTTTTTTTTACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAgg                 >  1:2720088/1‑99 (MQ=255)
                                                                               cccGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGa                <  1:2398436/100‑1 (MQ=255)
                                                                               cccGCACCGAAAGGTCCGGGGTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAg                  >  1:338768/1‑97 (MQ=255)
                                                                                     ccGAAAGGTCCGGGTTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGa          >  2:890204/1‑100 (MQ=255)
                                                                                     ccGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGa          <  2:2528656/100‑1 (MQ=255)
                                                                                       gAATGGTCCGGGGTGTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCa                   <  1:1866541/89‑1 (MQ=255)
                                                                                         aaGGTCCGGGGTTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGaact       <  1:2229990/100‑1 (MQ=255)
                                                                                         aaGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGaacta      <  2:1912442/100‑1 (MQ=255)
                                                                                         aaGGGCCGGGTTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGa          <  1:1449860/96‑1 (MQ=255)
                                                                                            gTCCGGGGTTTTTTTTTTTTTTTTATAACATAAACACGGCGCACACCATGGATAACAGCACAAAATTCGTGTTCTCACGCTTCCGGACGGGGaactagct   >  1:501453/1‑100 (MQ=255)
                                                                                             tCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACCGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGAactaacta  >  2:1974507/1‑100 (MQ=255)
                                                                                                    |                                                                                              
TGGTGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGAACTAACTA  >  NC_000913/3950368‑3950560

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