Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 235,082 A→G 71.8% K89K (AAA→AAG yafS → putative S‑adenosyl‑L‑methionine‑dependent methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913235,0820AG71.8% 229.1 / 103.0 151K89K (AAA→AAGyafSputative S‑adenosyl‑L‑methionine‑dependent methyltransferase
Reads supporting (aligned to +/- strand):  ref base A (21/21);  new base G (54/53);  total (76/75)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.12e-01

Only 100 of 152 total aligned reads displayed.

GCGAAGCGTGCGCGGTTTCTCATCAAGTGAATGTTTCTGCGCAAGGAATGCCCGTCCAGGTACAGGCGGACCCACTTCATCTTCCTTTTGCCGATAAATCCGTTGATGTTTGTCTACTGGCACATACATTGCCGTGGTGCACCGATCCGCATCGTTTATTGCGTGAAGCCGATCGGGTATTGATTGATGATGGCT  >  NC_000913/234985‑235179
                                                                                                 |                                                                                                 
gCGAAGCGTGCGCGGTTTCTCATCAAGTGAATGTTTCTGCGCAAGGAATGCCCGTCCAGGTACAGGCGGACCCACTTCATCTTCCTTTTGCCGATAAATc                                                                                                 <  2:1233248/100‑1 (MQ=255)
 cGAAGCGTGCGCGGTTTCTCATCAAGTGAATGTTTCTGCGCAAGGAATGCCCGTCCAGGTACAGGCGGACCCACTTCATCTTCCTTTTGCCGATAAGTcc                                                                                                >  2:1143678/1‑100 (MQ=255)
  gAAGCGTGCGCAGTTTCTCATCAAGTGAATGTTTCTGCGCAAGGAATGCCCGTCCAGGTACAGGCGGACCCACTTCATCTTCCTTTTGCCGATAAATCCg                                                                                               <  1:1545276/100‑1 (MQ=255)
    aGCGTGCGCGGTTTCTCATCAAGTGAATGTTTCTGCGCAAGGAATGCCCGTCCAGGTACAGGCGGACCCACTTCATCTTCCTTTTGCCGATAAGTCCg                                                                                               >  2:1773988/1‑98 (MQ=255)
     gCGTGCGCGGTTTCTCATCAAGTGAATGTTTCTGCGCAAGGAATGCCCGTCCAGGTACAGGCGGACCCACTTCATCTTCCTTTTGCCGATAAGTCCGTTg                                                                                            <  2:4097081/100‑1 (MQ=255)
          cgcgGTTTCTCATCAAGTGAATGTTTCTGCGCAAGGAATGCCCGTCCAGGTACAGGCGGACCCACTTCATCTTCCTTTTGCCGATAAGTCCGTTGATGtt                                                                                       >  1:961270/1‑100 (MQ=255)
          cgcgGTTTCTCATCAAGTGAATGTTTCTGCGCAAGGAATGCCCGTCCAGGTACAGGCGGACCCACTTCATCTTCCTTTTGCCGATAAGTCCGTTGATGtt                                                                                       <  1:858789/100‑1 (MQ=255)
           gcgGTTTCTCATCAAGTGAATGTTTCTGCGCAAGGAATGCCCGTCCAGGTACAGGCGGACCCACTTCATCTTCCTTTTGCCGATAAATCCGTTGATGttt                                                                                      >  1:3690881/1‑100 (MQ=255)
            cgGTTTCTCATCAAGTGTATGTTTCCGCGCAAGGAACGCCCGTCCCGGTACTGGCTGACCCTCTTCGTCTTCCTTTTGCCGATAAGTCCGTTGAGGTTTg                                                                                     <  2:1549873/100‑1 (MQ=255)
            cgGTTTCTCATCAAGTGAATGTTTCTGCGCAAGGAATGCCCGTCCAGGTACAGGCGGACCCACTTCATCTTCCTTTTGCCGATAAGTCCGTTGATGTTTg                                                                                     <  1:1217627/100‑1 (MQ=255)
            cgGTTTCTCATCAAGTGAATGTTTCTGCGCAAGGAATGCCCGTCCAGGTACAGGCGGACCCACTTCATCTTCCTTTTGCCGATAAGTCCGTTGATGTTTg                                                                                     <  2:2102874/100‑1 (MQ=255)
            cgGTTTCTCATCAAGTGAATGTTTCTGCGCAAGGAATGCCCGTCCAGGTACAGGCGGACCCACTTCATCTTCCTTTTGCCGATAAGTCCGTTGATGTTTg                                                                                     <  1:1561141/100‑1 (MQ=255)
            cgGTTTCTCATCAAGTGAATGTTTCTGCGCAAGGAATGCCCGTCCAGGTACAGGCGGACCCACTTCATCTTCCTTTTGCCGATAAATCCGTTGATGTTTg                                                                                     <  1:907241/100‑1 (MQ=255)
             ggTTTCTCATCAAGTGTATGTTTCTGCGCAAGGAATGCCCGTCCAGGTACAGGCGGACCCACTTCATCTTCCTTTTGCCGATAAGTCCGTTGATGTTTGt                                                                                    <  2:2431162/100‑1 (MQ=255)
             ggTTTCTCATCAAGTGAATGTTTCTGCGCAAGGAATGCCCGTCCAGGTACAGGCGGACCCACTTCATCTTCCTTTTGCCGATAAGTCCGTTGATGTTTGt                                                                                    <  1:2328093/100‑1 (MQ=255)
             ggTTTCTCATCAAGTGAATGTTTCTGCGCAAGGAATGCCCGTCCAGGTACAGGCGGACCCACTTCATCTTCCTTTTGCCGATAAGTCCGTTGATGTTTGt                                                                                    >  1:588604/1‑100 (MQ=255)
             ggTTTCTCATCAAGTGAATGTTTCTGCGCAAGGAATGCCCGTCCAGGTACAGGCGGACCCACTTCATCTTCCTTTTGCCGATAAGTCCGTTGATGTTTGt                                                                                    <  2:3550453/100‑1 (MQ=255)
                ttCTCATCAAGTGAATGTTTCTGCGCAAGGAATGCCCGTCCAGGTACAGGCGGACCCACTTCATCTTCCTTTTGCCGATAAATCCGTTGATGTTTGTCTa                                                                                 >  2:3125712/1‑100 (MQ=255)
                 tctcATCAAGTGAATGTTTCTGCGCAAGGAATGCCCGTCCAGGTACAGGCGGACCCACTTCATCTTCCTTTTGCCGATAAGTCCGTTGATGTTTGTCTAc                                                                                >  2:2598175/1‑100 (MQ=255)
                   tcatcaAGTGAATGTTTCTGCGCAAGGATTGCCCGTCCAGGTACAGGCGGACCCACTTCATCTTTCTTTTTCTGATAAGTCTGTTGATGTTTGTCTACTg                                                                              <  2:443604/100‑1 (MQ=255)
                    catcaAGTGAATGTTTCTGCGCAAGGAATGCCCGTCCAGGTACAGGCGGACCCACTTCATCTTCCTTTTGCCGATAAGTCCGTTGATGTTTGTCTACTgg                                                                             >  1:3228216/1‑100 (MQ=255)
                     atcaAGTGAATGTTTCTGCGCAAGGAATGCCCGTCCAGGTACAGGCGGACCCACTTCATCTTCCTTTTGCCGATAAGTCCGTTGATGTTTGTCTACTGGc                                                                            <  2:878568/100‑1 (MQ=255)
                      tcaAGTGAATGTTTCTGCGCAAGGAATGCCCGTCCAGGTACAGGCGGACCCACTTCATCTTCCTTTTGCCGATAAATCCGTTGATGTTTGTCTACTGGca                                                                           <  2:3025469/100‑1 (MQ=255)
                       caAGTGAATGTTTCGTCGCAAGGAATGCCCGTCGAGGTACAGGCGGACTCACTTCATCTTCCTTTTGCCGATAAGTCCGTTGATGTTTGTCTACTGGcac                                                                          <  2:3447412/100‑1 (MQ=255)
                        aaGTGAATGTTTCTGCGCAAGGAATGCCCGTCCAGGTACAGGCGGACCCACTTCATCTTCCTTTTGCCGATAAGTCCGTTGATGTTTGTCTACTGGcaca                                                                         <  2:1139726/100‑1 (MQ=255)
                        aaGTGAATGTTTCTGCGCAAGGAATGCCCGTCCAGGTACAGGCGGACCCACTTCATCTTCCTTTTGCCGATAAGTCCGTTGATGTTTGTCTACTGGcaca                                                                         <  1:3794012/100‑1 (MQ=255)
                        aaGTGAATGTTTCTGCGCAAGGAATGCCCGTCCAGGTACAGGCGGACCCACTTCATCTTCCTTTTGCCGATAAATCCGTTGATGTTTGTCTACTGGcaca                                                                         <  2:1603829/100‑1 (MQ=255)
                          gTGAATGTTTCTGCGCAAGGAATGCCCGTCCAGGTACAGGCGGACCCACTTCATCTTCCTTTTGCCGATAAGTCCGTTGATGTTTGTCTACTGGCacata                                                                       <  2:2833612/100‑1 (MQ=255)
                           tGAATGTTTCTGCGCAAGGAATGCCCGTCCAGGTACAGGCGGACCCACTTCATCTTCCTTTTGCCGATAAGTCCGTTGATGTTTGTCTACTGGCacatac                                                                      <  1:3977782/100‑1 (MQ=255)
                           tGAATGTTTCTGCGCAAGGAATGCCCGTCCAGGTACAGGCGGACCCACTTCATCTTCCTTTTGCCGATAAGTCCGTTGATGTTGGTCTACTGGCacata                                                                       >  2:3625953/1‑99 (MQ=255)
                           tGAATGTTTCTGCGCAAGGAATGCCCGTCCAGGTACAGGCGGACCCACTTCATCTTCCTTTTGCCGATAAATCCGTTGATGTTTGTCTACTGGCacatac                                                                      >  1:2523149/1‑100 (MQ=255)
                           tGAATGTTTCTGCGCAAGGAATGCCCGTCCAGGTACAGGCGGACCCACTTCATCTTCCTTTTGCCGATAAATCCGTTGATGTTTGTCTACTGGCacatac                                                                      >  1:2971618/1‑100 (MQ=255)
                             aaTGTTTCTGCGCAAGGAATGCCCGTCCAGGTACAGGCGGACCCACTTCATCTTCCTTTTGCCGATAAGTCCGTTGATGTTTGTCTACTGGCacatacat                                                                    <  1:3917290/100‑1 (MQ=255)
                             aaTGTTTCTGCGCAAGGAATGCCCGTCCAGGTACAGGCGGACCCACTTCATCTTCCTTTTGCCGATAAATCCGTTGATGTTTGTCTACTGGCacatacat                                                                    >  2:3940295/1‑100 (MQ=255)
                                gTTTCTGCGCAAGGAATGCCCGTCCAGGTACAGGCGGACCCACTTCATCTTCCTTTTGCCGATAAATCCGTTGATGTTTGTCTACTGGCACATACATTGc                                                                 <  1:3581093/100‑1 (MQ=255)
                                  ttCTGCGCAAGGAATGCCCGTCCAGGTACAGGCGGACCCACTTAATCTTCCTTTTGCCGATAAGTCCGTCGATGTTTGTCTACTGGCACATAGATTGCCg                                                               >  2:332841/1‑100 (MQ=255)
                                   tCTGCGCAAGGAATGCCCGTCCAGGTACAGGCGGACCCACTTCATCTTCCTTTTGCCGATAAGTCCGTTGATGTTTGTCTACTGGCACATACATTGCCgt                                                              <  2:588600/100‑1 (MQ=255)
                                    cTGCGCAAGGAATTCCCGTCCAGGTACAGGTGGACCCACTTCATCTTCCTTTTGCCGATAAGTCCGTTGATGTTTGTCTCCTGGCACATACATTGCCgt                                                              <  1:3476906/99‑1 (MQ=255)
                                    cTGCGCAAGGAATGCCCGTCCAGGTACAGGCGGACCCACTTCATCTTCCTTTTGCCGATCAATCCGTTGCTGTTTGTCTACTGGCACATACATTGCGgtg                                                             >  1:1577708/1‑100 (MQ=255)
                                    cTGCGCAAGGAATGCCCGTCCAGGTACAGGCGGACCCACTTCATCTTCCTTTTGCCGATAAGTCCGTTGATGTTTGTCTACTGGCACATACATTGCCgtg                                                             >  1:1148436/1‑100 (MQ=255)
                                    cTGCGCAAGGAATGCCCGTCCAGGTACAGGCGGACCCACTTCATCTTCCTTTTGCCGATAAATCCGTTGATGTTTGTCTACTGGCACATACATTGCCgtg                                                             >  1:1548117/1‑100 (MQ=255)
                                     tGCGCAAGGAATGCCCGTCCAGGTACAGGCGGACCCACTTCATCTTCCTTTTGCCGATAAGTCCGTTGATGTTTGTCTACTGGCACATACATTGCCgtgg                                                            <  2:3652532/100‑1 (MQ=255)
                                     tGCGCAAGGAATGCCCGTCCAGGTACAGGCGGACCCACTTCATCTTCCTTTTGCCGATAAGTCCGTTGATGTTTGTCTACTGGCACATACATTGCCgtgg                                                            <  2:4096755/100‑1 (MQ=255)
                                     tGCGCAAGGAATGCCCGTCCAGGTACAGGCGGACCCACTTCATCTTCCTTTTGCCGATAAATCCGTTGATGTTTGTCTACTGGCACATACATTGCCgtgn                                                            <  2:4074137/100‑2 (MQ=255)
                                      gcgcAAGGAATGCCCGTCCAGGTACAGGCGGACCCACTTCATCTTCCTCTTGCCGATAAGTCCGTTGATGTTTGTCTACTGGCACATACATTGCCgtgg                                                            >  2:3496286/1‑99 (MQ=255)
                                       cgcAAGGAATGCCCGTCCAGGTACAGGCGGACCCACTTCATCTTCCTTTTGCCGATAAGTCCGTTGATGTTTGTCTACTGGCACATACATTGCCgtggtg                                                          <  2:3105934/100‑1 (MQ=255)
                                       cgcAAGGAATGCCCGTCCAGGTACAGGCGGACCCACTTCATCTTCCTTTTGCCGATAAGTCCGTTGATGTTTGTCTACTGGCACATACATTGCCgtggtg                                                          >  2:3511505/1‑100 (MQ=255)
                                       cgcAAGGAATGCCCGTCCAGGTACAGGCGGACCCACTTCATCTTCCTTTTGCCGATAAGTCCGTTGATGTTTGTCTACTGGCACATACATTGCCgt                                                              >  2:1373293/1‑96 (MQ=255)
                                       agcAAGGAATGACCGTCCAGATACAGGCGGACCCACTCCATCTCCGTTTTGCCGATAACTCCGTTGATGTTTGTATACTGGCACATACATTGCCgtggtg                                                          >  2:397739/2‑100 (MQ=255)
                                        gcAAGGAATGCCCGTCCAGGTACAGGCGGACCCACTTCATCTTCCTTTTGCCGATAAGTCCGTTGATGTTTGTCTACTGGCACATACATTGCCGTGGTGc                                                         >  1:20984/1‑100 (MQ=255)
                                        gcAAGGAATGCCCGTCCAGGTACAGGCGGACCCACTTCATCTTCCTTTTGCCGATAAGTCCGTTGATGTTTGTCTACTGGCACATACATTGCCGTGGTGc                                                         >  1:1697748/1‑100 (MQ=255)
                                         cAAGGAATGCCCGTCCAGGTACAGGCGGACCCACTTCATCTTCCTTTTGCCGCAAAGTCCGTTGATGTTTGTCTACTGGCACATACATTGCCGTGGGGCa                                                        >  2:1997108/1‑100 (MQ=255)
                                          aaGGAATGCCCGTCCAGGTACAGGCGGACCCACTTCATCTTCCTTTTGCCGATAAATCCGTTGATGTTTGTCTACTGGCACATACATTGCCGTGGTGc                                                         <  1:2778101/98‑1 (MQ=255)
                                           aGGAATGCCCGTCCAGGTACAGGCGGACCCACTTCATCTTCCTTTTGCCGATAAATCCGTTGATGTTTGTCTACTGGCACATACATTGCCGTGGTGCAcc                                                      <  2:1548108/100‑1 (MQ=255)
                                               aTGCCCGTCCAGGTACAGGCGGACCCACTTCATCTTCCTTTTGCCGATAAGTCCGTTGATGTTTGTCTACTGGCACATACATTGCCGTGGTGCACCGATc                                                  >  1:879333/1‑100 (MQ=255)
                                               aTGCCCGTCCAGGTACAGGCGGACCCACTTCATCTTCCTTTTGCCGATAAGTCCGTTGATGTTTGTCTACTGGCACATACATTGCCGTGGTGCACCGATc                                                  <  1:3793699/100‑1 (MQ=255)
                                                tGCCCGTCCAGGTACAGGCGGACCCACTTCATCTTCCTTTTGCCGATAAGTCCGTTGATGTTTGTCTACTGGCACATACATTGCCGTGGTGCACCGATcc                                                 >  1:2917528/1‑100 (MQ=255)
                                                tGCCCGTCCAGGTACAGGCGGACCCACTTCATCTTCCTTTTGCCGATAAGTCCGTTGATGTTTGTCTACTGGCACATACATTGCCGTGGTGCACCGATcc                                                 <  1:3178538/100‑1 (MQ=255)
                                                tGCCCGTCCAGGTACAGGCGGACCCACTTCATCTTCCTTTTGCCGATAAGTCCGTTGATGTTTGTCTACTGGCACATACATTGCCGTGGTGCACCGATc                                                  >  2:1780673/1‑99 (MQ=255)
                                                  cccGTCCAGGTACAGGCGGACCCACTTCATCTTCCTTTTGCCGATAAGTCCGTTGATGTTTGTCTACTGGCACATACATTCCCGTGGTGCACCGATCCGc                                               >  1:4005541/1‑100 (MQ=255)
                                                  cccGTCCAGGTAAAGGCGGACCCACTTCATCTTCCTTTTGCCGATAAGTCCGTTGATGTTTGTCTACTGGCACATACATTGCCGTGGTGCACCGATCCGc                                               >  1:2735910/1‑100 (MQ=255)
                                                   ccGTCCAGGTCCAGGCGGACCCACTTCATCTTCCTTTTGCCGATAAGTCCGTTTATGTTTGTCTACTGGCACATACATTGCCGTGGTGCACCGATcc                                                 <  2:1423943/97‑1 (MQ=255)
                                                   ccGTCCAGGTACAGGCGGACCCACTTCATCTTCCTTTTGCCGATAAGTCCGTTGATGTTTGTCTACTGGCACATACATTGCCGTGGTGCACCGATCCGcn                                              <  1:39213/100‑2 (MQ=255)
                                                   ccGTCCAGGTACAGGCGGACCCACTTCATCTTCCTTTTGCCGATAAATCCGTTGATGTTTGTCTACTGGCACATACATTGCCGTGGTGCACCGATCCGCa                                              <  1:2627075/100‑1 (MQ=255)
                                                          ggTACAGGCGGACCCACTTCATCTTCCTTTTGCCGATAAATCCGTTGATGTTTGTCTACTGGCACATACATTGCCGTGGTGCACCGATCCGCATCGTTTa                                       >  1:2123112/1‑100 (MQ=255)
                                                            tACAGGCGGACCCACTTCATCTTCCTTTTGCCGATAAGTCCGTTGATGTTTGTCTACTGGCACATACATTGCCGTGGTGCACCGATCCGCATCGTTTAtt                                     >  2:1182425/1‑100 (MQ=255)
                                                             aCAGGCGGACCCACTTCATCTTCCTTTTGCCGATAAGTCCGTTGATGTTTGTCTACTGGCACATACATTGCCGTGGTGCACCGATCCGCATCGTTTATTg                                    <  2:371302/100‑1 (MQ=255)
                                                             aCAGGCGGACCCACTTCATCTTCCTTTTGCCGATAAGTCCGTTGATGTTTGTCTACTGGCACATACATTGCCGTGGTGCACCGATCCGCATCGTTTATTg                                    <  1:436365/100‑1 (MQ=255)
                                                             aCAGGCGGACCCACTTCATCTTCCTTTTGCCGATAAGTCCGTTGATGTTTGTCTACTGGCACATACATTGCCGTGGTGCACCGATCCGCATCGTTTATTg                                    <  2:3203100/100‑1 (MQ=255)
                                                             aCAGGCGGACCCACTTCATCTTCCTTTTGCCGATAAATCCGTTGATGTTTGTCTACTGGCACATACAGTGCCGTGGTGCACCGATCCGCATCGTTTATTg                                    <  2:59059/100‑1 (MQ=255)
                                                               aGGCGGACCCACTTCATCTTCCTTTTGCCGATAAGTCCGTTGATGTTTGTCTACTGGCACATACATTGCCGTGGTGCACCGATCCGCATCGTTTATTGCg                                  <  1:3836350/100‑1 (MQ=255)
                                                                ngCGGACCCACTTCATCTTCCTTTTGCCGATAAATCCGTTGATGTTTGTCTACTGGCACATACATTGCCGTGGTGCACCGATCCGCATCGTTTATTGCGt                                 >  2:3014548/2‑100 (MQ=255)
                                                                ggCGGACCCACTTCATCTTCCTTTTGCCGATAAGTCCGTTGATGTTTGTCTACTGGCACATACATTGCCGTGGTGCACCGATCCGCATCGTTTATTGCGt                                 >  1:1350063/1‑100 (MQ=255)
                                                                ggCGGACCCACTTCATCTTCCTTTTGCCGATAAATCCGTTGATGTTTGTCTACTGGCACATACATTGCCGTGGTGCACCGATCCGCATCGTTTATTGCGt                                 >  1:637232/1‑100 (MQ=255)
                                                                 gCGGACCCACTTCATCTTCCTTTTGCCGATAAGTCCGTTGATGTTTGTCTACTGGCACATACATTGCCGTGGTGCACCGATCCGCATCGTTTATTGCGTg                                <  2:1460493/100‑1 (MQ=255)
                                                                  cGGACCCACTTCATCTTCCTTTTGCCGATAAGTCCGTTGATGTTTGTCTACTGGCACATACATTGCCGTGGTGCACCGATCCGCATCGTTTATTGCGTGa                               >  1:1785672/1‑100 (MQ=255)
                                                                   ggACCCACTTCATCTTCCTTTTGCCGATAAGTCCGTTGATGTTTGTCTACTGGCACATACATTGCCGTGGTGCACCGATCCGCATCGTTTATTGCGTGaa                              <  2:1350055/100‑1 (MQ=255)
                                                                     aCCCACTTCATCTTCCTTTTGCCGATAAGTCCGTTGATGTTTGTCTACTGGCACATACATTGCCGTGGTGCACCGATCCGCATCGTTTATTGCGTGAAGc                            <  1:1535932/100‑1 (MQ=255)
                                                                      cccACTTCATCTTCCTTTTGCCGATAAGTCCGTTGATGTTTGTCTACTGGCACATACATTGCCGTGGTGCACCGATCCGCATCGTTTATTGCGTGAAGcc                           >  2:502539/1‑100 (MQ=255)
                                                                      cccACTTCATCTTCCTTTTGCCGATAAGTCCGTTGATGTTTGTCTACTGGCACATACATTGCCGTGGTGCACCGATCCGCATCGTTTATTGCGTGAAGcc                           >  2:124935/1‑100 (MQ=255)
                                                                      cccACTTCATCTTCCTTTTGCCGATAAATCCGTTTATGTTTGTCTACTGGCACATACATTGCCGTGGTGCACCGATCCGCATCGTTTATTGCGTGAAGcc                           <  2:3252598/100‑1 (MQ=255)
                                                                        cACTTCATCTTCCTTTTGCCGATAAATCCGTTGATGTTTGTCTACTGGCACATACATTGCCGTGGTGCACCGATCCGCATCGTTTATTGCGTGAAGCCGa                         >  2:2073589/1‑100 (MQ=255)
                                                                         aCTTCATCTTCCTTTTGCCGATAAGTCCGTTGATGTTTGTCTACTGGCACATACATTGCCGTGGTGCACCGATCCGCATCGTTTATTGCGTGAAGCCGAt                        >  1:3182810/1‑100 (MQ=255)
                                                                         aCTTCATCTTCCTTTTGCCGATAAGTCCGTTGATGTTTGTCTACTGGCACATACATTGCCGTGGTGCACCGATCCGCATCGTTTATTGCGTGAAGCCGAt                        >  1:1521102/1‑100 (MQ=255)
                                                                             cATCTTCCTTTTGCCGATAAATCCGTTGATGTTTGTCTACTGGCACATACATTGCCGTGGTGCACCGATCCGCATCGTTTATTGCGTGAAGCCGATCggg                    <  1:3298610/100‑1 (MQ=255)
                                                                             cATCTGCCTTTTGCCGATAAATCCGTTGATGTTTGTCTACTGGCACATACATTGCCGTGGTGCACCGATCCGCATCGTTTATTGCGTGAAGCCGATCggg                    <  2:3283083/100‑1 (MQ=255)
                                                                                cTTCCTTTTGCCGATAAGTCCGTTGATGTTTGTCTACTGGCACATACATTGCCGTGGTGCACCGATCCGCATCGTTTATTGCGTGAAGCCGATCGGGTat                 >  2:1294242/1‑100 (MQ=255)
                                                                                 ttCTTTTTACGGATAAGTACGTTGATGTTTGGCTACTGGCACATACATTGCCGGGGTGAACAGATCCGCATCGTTTAATACGTGAACCTGATCGGGTat                 <  2:1590273/99‑1 (MQ=255)
                                                                                    cTTTTGCCGATAAGTCCGTTGATGTTTGTCTACTGGCACATACATTGCCGTGGTGCACCGATCCGCATCGTTTATTGCGTGAAGCCGATCGGGTattgat             <  2:2639663/100‑1 (MQ=255)
                                                                                       ttGCCGATAAGTCCGTTGATGTTTGTCTACTGGCACATACATTGCCGTGGTGCACCGATCCGCATCGTTTATTGCGTGAAGCCGATCGGGTAttgattga          <  1:3827750/100‑1 (MQ=255)
                                                                                       ttGCCGATAAGTCCGTTGATGTTTGTCTACTGGCACATACATTGCCGTGGTGCACCGATCCGCATCGTTTATTGCGTGAAGCCGATCGGGTAttgattga          <  1:4144654/100‑1 (MQ=255)
                                                                                        ttccGATAATTCCGTTGATGTTTGTCTACTGTCACATACATTGCCGTGGTGCACCGATCCGCATCGTTTATTGCGTGAAGCCGATCGGGTATtgattgat         <  1:2311579/98‑1 (MQ=255)
                                                                                        tGCCGATAAGTCCGTTGATGTTTGTCTACTGGCACATACATTGCCGTGGTGCACCGATCCGCATCGTTTATTGCGTGAAGCCGATCGGGTAttgattga          <  2:4084205/99‑1 (MQ=255)
                                                                                        tGCCGATAAATCCGTTGATGTTTGTCTACTGGCACATACATTGCCGTGGTGCACCGATCCGCATCGTTTATTGCGTGAAGCCGATCGGGTATtgattgat         <  1:578475/100‑1 (MQ=255)
                                                                                         gCCGATAAGTCCGTTGATGTTTGTCTACTGGCACATACATTGCCGTGGTGCACCGATCCGCATCGTTTATTGCGTGAAGCCGATCGGGTATTGATtgatg        <  2:2454484/100‑1 (MQ=255)
                                                                                         gCCGATAAGTCCGTTGATGTTTGTCTACTGGCACATACATTGCCGTGGTGCACCGATCCGCATCGTTTATTGCGTGAAGCCGATCGGGTATTGATtgatg        >  1:605670/1‑100 (MQ=255)
                                                                                         gCCGATAAATCCGTTGATGTTTGTCTACTGGCACATACATTGCCGTGGTGCACCGATCCGCATCGTTTATTGCGTGAAGCCGATCGGGTATtgattgat         <  2:1293418/99‑1 (MQ=255)
                                                                                           cGATAAGTCCGTTGATGTTTGTCTACTGGCACATACATTGCCGTGGTGCACCGATCCGCATCGTTTATTGCGTGAAGCCGATCGGGTATTGATTgatgat      >  2:2126952/1‑100 (MQ=255)
                                                                                            gATAAGTCCGTTGATGTTTGTCTACTGGCACATACATTGCCGTGGTGCACCGATCCGCATCGTTTATTGCGTGAAGCCGATCGGGTATTGATTGatgatg     >  2:3613413/1‑100 (MQ=255)
                                                                                                aaTCCGTTGATGTTTGTCTACTGGCACATACATTGCCGTGGTGCACCGATCCGCATCGTTTATTGCGTGAAGCCGATCGGGGATTGATTGATGAtggct  >  1:1262041/1‑99 (MQ=255)
                                                                                                 |                                                                                                 
GCGAAGCGTGCGCGGTTTCTCATCAAGTGAATGTTTCTGCGCAAGGAATGCCCGTCCAGGTACAGGCGGACCCACTTCATCTTCCTTTTGCCGATAAATCCGTTGATGTTTGTCTACTGGCACATACATTGCCGTGGTGCACCGATCCGCATCGTTTATTGCGTGAAGCCGATCGGGTATTGATTGATGATGGCT  >  NC_000913/234985‑235179

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