Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,950,468 G→T 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.8% 283.5 / 0.3 96intergenic (+48/‑39)ilvL/ilvXilvG operon leader peptide/uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (0/3);  new base T (39/51);  total (42/54)
Fisher's exact test for biased strand distribution p-value = 2.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.72e-01

Only 100 of 105 total aligned reads displayed.

TGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑GTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGAACTAACTATG  >  NC_000913/3950371‑3950562
                                                                                                  |                                                                                              
tgATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑ttt                                                                                              <  1:2782383/100‑1 (MQ=255)
tgATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑ttt                                                                                              >  1:2210185/1‑100 (MQ=255)
 gATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑tttt                                                                                             <  1:2515880/100‑1 (MQ=255)
 gATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑tttt                                                                                             >  2:2558485/1‑100 (MQ=255)
  attattATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑Ttttt                                                                                            >  2:3299663/1‑100 (MQ=255)
   ttattaTCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑tttttt                                                                                           <  1:3380949/100‑1 (MQ=255)
   ttattaTCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑tttttt                                                                                           >  2:1225097/1‑100 (MQ=255)
     attatCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑tttttttt                                                                                         <  2:3391236/100‑1 (MQ=255)
     attatCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑tttttttt                                                                                         >  2:3153806/1‑100 (MQ=255)
     attatCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑tttttttt                                                                                         <  2:1154844/100‑1 (MQ=255)
     attatCCCACCGTACGGGGCTACACTTGGACGAGGAAAAGCTTACAGAACAAGCCTTACCGAAGTAAGACCCCCGAACCGAAAGGTCCAGAG‑Ttttt                                                                                            >  1:2150366/1‑97 (MQ=255)
       tatCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTg                                                                                       <  1:1139024/100‑1 (MQ=255)
       tatCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTg                                                                                       <  2:1770533/100‑1 (MQ=255)
       tatCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTg                                                                                       >  2:1072264/1‑100 (MQ=255)
       tatCCCACCGTGCGGGACTGAACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTTAGGCCCCCGCACCGAAAGGTCCGGGG‑TTtttgtttg                                                                                       <  1:774184/100‑1 (MQ=255)
        atCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGa                                                                                      <  1:687047/100‑1 (MQ=255)
        atCCCACAGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGa                                                                                      <  1:1865491/100‑1 (MQ=255)
         tCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGAc                                                                                     <  1:1781577/100‑1 (MQ=255)
          cccACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGAcc                                                                                    >  2:1245520/1‑100 (MQ=255)
           ccACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCt                                                                                   >  1:11782/1‑100 (MQ=255)
           ccACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCt                                                                                   >  2:1597418/1‑100 (MQ=255)
            cACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCtt                                                                                  >  2:133246/1‑100 (MQ=255)
            cACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCt                                                                                   <  1:205246/99‑1 (MQ=255)
              ccGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTaa                                                                                >  1:3445960/1‑100 (MQ=255)
              ccGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTaa                                                                                <  2:1707257/100‑1 (MQ=255)
                gTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTaaaa                                                                              <  2:459569/100‑1 (MQ=255)
                 tGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTaaaaa                                                                             <  2:1116921/100‑1 (MQ=255)
                      ggCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATaa                                                                        >  2:2997977/1‑100 (MQ=255)
                      ggCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATaa                                                                        >  2:1617063/1‑100 (MQ=255)
                       gCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAAc                                                                       <  2:2926495/100‑1 (MQ=255)
                       gCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAAc                                                                       >  2:721118/1‑100 (MQ=255)
                       gCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGAACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAAc                                                                       <  1:834742/100‑1 (MQ=255)
                        cTGCATTTGGACGAGGAAAGGCGTAGAGATCAAGCCAAAACGAAAAAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAAcc                                                                      <  2:1063196/100‑1 (MQ=255)
                          gCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGATCTTAAAAACATaa                                                                        >  1:2635862/1‑96 (MQ=255)
                           cACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCgag                                                                   >  1:1857693/1‑100 (MQ=255)
                           cACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGATACGAGG‑ATTTTTTTTGACATTAAAAACATAACCgag                                                                   >  2:213139/1‑100 (MQ=255)
                            aCTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCgagg                                                                  <  1:721109/100‑1 (MQ=255)
                            aCTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCgagg                                                                  >  2:2105386/1‑100 (MQ=255)
                             cTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCgagga                                                                 >  1:498127/1‑100 (MQ=255)
                               tGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGGCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGc                                                               <  1:2740383/100‑1 (MQ=255)
                                ggACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCa                                                              <  2:478950/100‑1 (MQ=255)
                                 gACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAg                                                             <  1:1926223/100‑1 (MQ=255)
                                  aCGAGGAAAGGCTTAGAGATCGAGCCTTAACGAACGAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGa                                                            <  2:3048379/100‑1 (MQ=255)
                                    gAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACa                                                          >  1:2780160/1‑100 (MQ=255)
                                      ggAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAAt                                                        <  2:1972403/100‑1 (MQ=255)
                                      ggAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAAt                                                        <  1:2558451/100‑1 (MQ=255)
                                      ggAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAAt                                                        <  2:407025/100‑1 (MQ=255)
                                      ggAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCC‑GGG‑GTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATg                                                       <  1:2133333/100‑1 (MQ=255)
                                       gAAAGGCTTAGAGATCAAGCCGTAACGAACTAAGACCCCCGCACAGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATg                                                       <  2:169868/100‑1 (MQ=255)
                                        aaaGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGa                                                      >  2:2266536/1‑100 (MQ=255)
                                        aaaGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGa                                                      >  1:2684719/1‑100 (MQ=255)
                                         aaGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCGCGCACCGAAAGATCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGaa                                                     <  2:2007714/100‑1 (MQ=255)
                                         aaGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGaa                                                     <  2:3446013/100‑1 (MQ=255)
                                          aGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCC‑GGG‑GTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATa                                                   <  1:1956529/100‑1 (MQ=255)
                                           ggCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACTTTAAAAACATAACCAAGGAGTAGACAATGAATa                                                   >  1:3186495/1‑100 (MQ=255)
                                            gCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTTTTTTTAATTACATAAACGAGGGGCAAAAAATg                                                       >  1:2069241/1‑95 (MQ=255)
                                            gCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATaa                                                  >  2:2825709/1‑100 (MQ=255)
                                             cTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCCCCAAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAACACATAACCGAGGACCAGACAATGAATAAc                                                 >  1:2241077/1‑100 (MQ=255)
                                              ttAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACa                                                >  2:1581718/1‑100 (MQ=255)
                                               tAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAg                                               <  2:2460689/100‑1 (MQ=255)
                                               tAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAg                                               <  2:2243844/100‑1 (MQ=255)
                                               tAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGCATCAGACAATGAATAACAg                                               >  1:742541/1‑100 (MQ=255)
                                                agagATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGc                                              >  2:2397922/1‑100 (MQ=255)
                                                agagATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGc                                              <  1:1376201/100‑1 (MQ=255)
                                                 gagaTCAAGCCTTAACGAACTAAGGCCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGca                                             <  2:365891/100‑1 (MQ=255)
                                                     tCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACaaa                                         >  2:2142251/1‑100 (MQ=255)
                                                       aaGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAAt                                       <  1:1650167/100‑1 (MQ=255)
                                                       aaGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAAt                                       <  2:3339336/100‑1 (MQ=255)
                                                       aaGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAAt                                       <  1:138483/100‑1 (MQ=255)
                                                         gCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTc                                     >  2:3143739/1‑100 (MQ=255)
                                                           cTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTg                                   <  2:2225887/100‑1 (MQ=255)
                                                           cTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTg                                   >  2:300411/1‑100 (MQ=255)
                                                            ttAACGAACTAAGACCCCCGCGCCTAATGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCt                                    <  2:1917050/98‑1 (MQ=255)
                                                            ttAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTg                                   <  1:2612107/100‑1 (MQ=255)
                                                            ttAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTg                                   <  2:3364788/99‑1 (MQ=255)
                                                              aaCGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTCAAAACATAACCGAGGAGCAGACACTGAATAACAGCTCAAAATTCTCtt                                 >  1:2415277/1‑99 (MQ=255)
                                                              aaCGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGttt                                >  1:2298300/1‑100 (MQ=255)
                                                              aaCGAACTAAGACCCCCGCACCGAAAGGTCCCGGG‑TTTGTTCTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGtt                                 <  1:2437218/99‑1 (MQ=255)
                                                              aaCGAACTAAGACCCCCGCACCGAAAGGTCC‑GGG‑GTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACTGCACAAAATTCTGTTtc                               <  2:1002799/100‑1 (MQ=255)
                                                                cGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑CTTTTTTTTTTCCTTTAAAAAATAACCGACGCGCAGACAATAAATAACTGCACAACATTCAGTTtct                              >  1:1626792/1‑100 (MQ=255)
                                                                    cTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAg                          <  1:120134/100‑1 (MQ=255)
                                                                       aGACCCCCGCACCAAAAGGCCCGGGC‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCaa                           <  1:2397889/96‑1 (MQ=255)
                                                                        gACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTc                      >  1:3404908/1‑100 (MQ=255)
                                                                        gACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCAAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTc                      <  2:715629/100‑1 (MQ=255)
                                                                          cccccGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCa                     <  1:812012/99‑1 (MQ=255)
                                                                          cccccGCACCAAAAGGTCCGGGG‑TTTTTAAATAACCAAAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTc                      <  1:2825665/98‑1 (MQ=255)
                                                                           ccccGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAgg                   >  1:538956/1‑100 (MQ=255)
                                                                            cccGCACCGAAAGGTCCGGGG‑TTTTTTTTTTTTTTTTAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCa                     >  1:2629691/1‑97 (MQ=255)
                                                                            cccGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGCGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGa                  >  1:538297/1‑100 (MQ=255)
                                                                            cccGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAgg                   <  1:1942405/99‑1 (MQ=255)
                                                                               gCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACgg               <  2:1780223/100‑1 (MQ=255)
                                                                                  ccGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGACTAACAGCACAAAATTCTGTTTCTCAAGATTCAGCAc                 >  1:497206/1‑95 (MQ=255)
                                                                                     aaaGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGaact         <  1:2608331/100‑1 (MQ=255)
                                                                                     aaaGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGaact         <  1:1695645/100‑1 (MQ=255)
                                                                                       aGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGaacta        <  1:3299612/99‑1 (MQ=255)
                                                                                         gTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGGTTCTCAAGATTCAGGACGGGGaactaact     >  2:1835997/1‑100 (MQ=255)
                                                                                           gcGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGaactaact     <  1:1152317/97‑1 (MQ=255)
                                                                                           ccGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGAACTAACTAt   >  2:3306662/1‑100 (MQ=255)
                                                                                           ccGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGAACTAACTAt   <  1:2105355/100‑1 (MQ=255)
                                                                                            cGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGAACTAACTATg  >  2:293354/1‑100 (MQ=255)
                                                                                                  |                                                                                              
TGATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑GTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGAACTAACTATG  >  NC_000913/3950371‑3950562

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: ‑