Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,950,468 G→T 100% intergenic (+48/‑39) 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,4680GT96.6% 465.3 / 4.9 156intergenic (+48/‑39)ilvL/ilvXilvG operon leader peptide/uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (2/4);  new base T (57/91);  total (61/95)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.51e-01

Only 100 of 173 total aligned reads displayed.

TGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGAACTAACTATGAA  >  NC_000913/3950371‑3950564
                                                                                                 |                                                                                                
tgATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGttt                                                                                                >  1:3853898/1‑100 (MQ=255)
tgATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGttt                                                                                                >  2:1928438/1‑100 (MQ=255)
tgATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGttt                                                                                                >  1:3581273/1‑100 (MQ=255)
 gATTATTATCCCACCGTGCGGGGCTGCGCTTGGACGAGGAACGGCTTAGAGATCAAGCCCTAACGAACTAAGACCCCCGCTCCGAAAGGTCCGGGGtttt                                                                                               >  2:3738280/1‑100 (MQ=255)
 gATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGtttt                                                                                               >  1:185680/1‑100 (MQ=255)
 gATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGtttt                                                                                               >  1:3313389/1‑100 (MQ=255)
   ttattaTCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGtttttt                                                                                             >  2:1398072/1‑100 (MQ=255)
   ttattaTCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTtttt                                                                                              >  2:3602069/1‑99 (MQ=255)
    tattatCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTtttttt                                                                                            >  1:3421309/1‑100 (MQ=255)
     attatCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGtttttttt                                                                                           <  2:3157918/100‑1 (MQ=255)
     attatCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGtttttttt                                                                                           >  1:2270269/1‑100 (MQ=255)
     attatCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGtttttttt                                                                                           <  1:1831994/100‑1 (MQ=255)
      ttatCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTtttttttt                                                                                          >  2:3729620/1‑100 (MQ=255)
      ttatCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTtttttttt                                                                                          >  1:2756887/1‑100 (MQ=255)
       tatCCCACCGTGCGGGGCTGCACTTAGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTg                                                                                         <  2:1365679/100‑1 (MQ=255)
        atCCCACCGTGCGGGGCTCCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGa                                                                                        <  2:1026804/100‑1 (MQ=255)
        atCCAACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGa                                                                                        <  2:1023003/100‑1 (MQ=255)
         tCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGAc                                                                                       <  1:4131589/100‑1 (MQ=255)
           ccACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGAcc                                                                                      <  2:2324144/99‑1 (MQ=255)
           ccACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCt                                                                                     >  1:3565777/1‑100 (MQ=255)
           ccACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCATTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCt                                                                                     >  1:3512540/1‑100 (MQ=255)
             aCCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTa                                                                                   <  2:1949651/100‑1 (MQ=255)
              ccGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTATCGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTa                                                                                   <  1:468225/99‑1 (MQ=255)
                  gCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAAc                                                                              <  2:3334586/100‑1 (MQ=255)
                   cGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACa                                                                             <  1:3729647/100‑1 (MQ=255)
                   cGGGGCGGCACTTGGACGAGGAAAGGCTTAGAGATCAAGGCGTAACGAACTAAGACCCCCGCACCGAAAGGGCCGGGGTTTTTTTTTGACCTTAAAAACa                                                                             <  2:2064943/100‑1 (MQ=255)
                       gCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAAc                                                                         <  1:2650093/100‑1 (MQ=255)
                       gCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAAc                                                                         >  1:4175518/1‑100 (MQ=255)
                       gCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAAc                                                                         >  1:1498316/1‑100 (MQ=255)
                         tGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCg                                                                       <  2:1129263/100‑1 (MQ=255)
                            aCTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGATTTTTTTTGACCTTAAAAACATAACCgcg                                                                     >  2:970996/1‑97 (MQ=255)
                            aCTTGGAACAGGAAAGGCTTAGAGAACAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCgagg                                                                    <  1:1589018/100‑1 (MQ=255)
                             cTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCgagga                                                                   >  2:2769617/1‑100 (MQ=255)
                             cTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCgagga                                                                   <  2:54353/100‑1 (MQ=255)
                              ttGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAATCATAACCgaggag                                                                  <  2:2582587/100‑1 (MQ=255)
                                 gACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAg                                                               <  2:3853869/100‑1 (MQ=255)
                                    gAGGAAAGGCTTAGAGATCAAGCCTTAACTAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTGTTTGACCTTAAAAACATAACCGAGGAGCAGACa                                                            <  2:2510596/100‑1 (MQ=255)
                                     aGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGAACTTAAAAACATAACCGAGGAGCAGACaa                                                           >  2:3272719/1‑100 (MQ=255)
                                      ggAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAAt                                                          >  1:1727140/1‑100 (MQ=255)
                                       gAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATg                                                         <  2:2533924/100‑1 (MQ=255)
                                       gAAAGGCTTAGAGATCAAGCCTTAACGAACGAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATg                                                         <  1:2446224/100‑1 (MQ=255)
                                         aaGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGCGCAGACACTGaa                                                       >  2:2891803/1‑100 (MQ=255)
                                         aaGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGaa                                                       >  1:902224/1‑100 (MQ=255)
                                         aaGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAAt                                                      <  1:3455182/100‑1 (MQ=255)
                                          aGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAAt                                                      >  1:615613/1‑100 (MQ=255)
                                           ggCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGTTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATa                                                     >  1:2125096/1‑100 (MQ=255)
                                           ggCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATa                                                     <  2:1484524/100‑1 (MQ=255)
                                           ggCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATa                                                     <  1:2336378/100‑1 (MQ=255)
                                             cTTAGAGATCAAGCCTTAACGAACTGAGCCCCCCGCACCGAAAGGTCCGGGGTCTTTGTTTTACCTTACAAAACTAACGGCGGCGCAGAAAATGAATaa                                                    >  2:2625375/1‑99 (MQ=255)
                                              ttAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACa                                                  <  2:162049/100‑1 (MQ=255)
                                               tAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAg                                                 >  1:2191850/1‑100 (MQ=255)
                                               tAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAg                                                 <  2:3957635/100‑1 (MQ=255)
                                                 gagaTCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGca                                               >  1:2951648/1‑100 (MQ=255)
                                                    aTCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGcacan                                            <  1:2348478/100‑2 (MQ=255)
                                                       aaGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAAt                                         <  2:2613263/100‑1 (MQ=255)
                                                         gCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTc                                       <  1:1519872/100‑1 (MQ=255)
                                                         gCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTc                                       >  2:4095412/1‑100 (MQ=255)
                                                          ccTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCt                                      >  1:354809/1‑100 (MQ=255)
                                                          ccTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGCTTTTTTTTGCCCTTAACACCATAACCGACGAGCCGACAATGAATAACACcaca                                             >  2:3247990/1‑93 (MQ=255)
                                                           cTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTg                                     <  2:216772/100‑1 (MQ=255)
                                                           cTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTg                                     <  1:941479/100‑1 (MQ=255)
                                                           cTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTg                                     <  1:1041370/100‑1 (MQ=255)
                                                            ttAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTg                                     >  1:3649442/1‑99 (MQ=255)
                                                            ttAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGt                                    <  2:2943064/100‑1 (MQ=255)
                                                            ttAACGAACTAAGACCCCCGCACAGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGt                                    <  2:4038041/100‑1 (MQ=255)
                                                            ttAACGAACTAACACCCCCGGACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGt                                    <  1:660576/100‑1 (MQ=255)
                                                             tAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTTTTCTTTTAAACATAACCGAGGAGCAAACAATGAATAACAGCACAACATTCTGtt                                   >  2:654768/1‑100 (MQ=255)
                                                              aaCGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTg                                     <  1:2413937/97‑1 (MQ=255)
                                                               aCGAACTAAGACCCCCGCACCGGAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTtc                                 <  1:3133545/100‑1 (MQ=255)
                                                               aCGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTtc                                 <  2:1897323/100‑1 (MQ=255)
                                                               aCGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTtc                                 <  1:3977717/100‑1 (MQ=255)
                                                                 gAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTtctc                               <  1:918112/100‑1 (MQ=255)
                                                                 gAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTtc                                 <  1:501272/98‑1 (MQ=255)
                                                                 gAACTAAGACCACCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAACCATAAACGGCGAGAAGCCAATGACTAACCGCACAACAtac                                       >  2:3442769/1‑90 (MQ=255)
                                                                   aCTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCaa                             <  1:3638079/100‑1 (MQ=255)
                                                                     tAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAGCCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAg                            <  2:864581/99‑1 (MQ=255)
                                                                          cccccGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCa                       <  2:2710530/99‑1 (MQ=255)
                                                                          cccccGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAg                      >  1:307222/1‑100 (MQ=255)
                                                                          cccccGCACCGAAAGGTCCGGGGTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAgg                     <  2:3797681/100‑1 (MQ=255)
                                                                           ccccGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAgg                     <  1:3161281/100‑1 (MQ=255)
                                                                           ccccGAGCCGAAAGGCCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAgg                     <  1:885701/100‑1 (MQ=255)
                                                                            cccGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGa                    >  1:2140136/1‑100 (MQ=255)
                                                                             ccGCATCGAAAGGGCCGGGGTTTTTTTTTGACCTTAAAAACATACCCGGGTAGTAGACAATGAAGAACAACACAAAATGCTGTTTCTCAAGTTTCAGGAc                   <  2:901153/100‑1 (MQ=255)
                                                                                cACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGCGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACggg                >  1:3842509/1‑100 (MQ=255)
                                                                                cACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACggg                <  1:1050027/100‑1 (MQ=255)
                                                                                 agcGAAGGGTCCGGGGTTTTTTTTTGACCTTAGAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAGATTCTGTTTCTCAAGATTCAGGACgggg               <  2:1165103/98‑1 (MQ=255)
                                                                                 aCCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACgggg               <  2:2962989/100‑1 (MQ=255)
                                                                                  ccGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACggggn              <  1:1685162/100‑2 (MQ=255)
                                                                                  ccGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACgggg               >  1:1256744/1‑99 (MQ=255)
                                                                                  ccGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGa              <  1:1380132/100‑1 (MQ=255)
                                                                                  ccGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGa              >  1:1252349/1‑100 (MQ=255)
                                                                                    gAAAGGTCCGGGGTTTTTTTTTTTTTTTTTATTCTTAACCGAGGAGCAGTCAATGAATAATTGCACAAGATTCTGTTTCTCAAGATTCCGGACGGGGaa             >  2:3542231/1‑99 (MQ=255)
                                                                                      aaGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACCTTCAATAACAGCTCAAAATTGTGTTTCTCAAGCTTCAGCACgcg                >  2:681050/1‑92 (MQ=255)
                                                                                      aaGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGaact           <  2:1125001/99‑1 (MQ=255)
                                                                                       aGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGaactaa         >  1:1058688/1‑100 (MQ=255)
                                                                                           ccGGGGTTTTTTTTTTTTTTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGAACTAACTAt     >  2:242390/1‑100 (MQ=255)
                                                                                            ccggggTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGAACTAACTATg    >  1:1521468/3‑100 (MQ=255)
                                                                                              gggtTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGAACTAACTATGaa  <  2:3260148/100‑1 (MQ=255)
                                                                                              ggggTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGAACTAACTATGa   <  2:2682128/99‑1 (MQ=255)
                                                                                               gggTTTTTTTTTTTTTTTAAAACATAACCGAGGAGCAGACAATGCATAACAGCACAACAGTCTGGTTCTCTAGTTTCAGGACGGGGAACTTACTAt     >  2:1132432/1‑96 (MQ=255)
                                                                                                 |                                                                                                
TGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGAACTAACTATGAA  >  NC_000913/3950371‑3950564

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: