Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 95,418 T→G S257A (TCG→GCG)  murF → UDP‑N‑acetylmuramoyl‑tripeptide:D‑alanyl‑D‑ alanine ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091395,4180TG99.3% 482.1 / ‑4.7 138S257A (TCG→GCG) murFUDP‑N‑acetylmuramoyl‑tripeptide:D‑alanyl‑D‑ alanine ligase
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base G (70/67);  total (71/67)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.97e-01

Only 100 of 143 total aligned reads displayed.

GGCTGAACTGGCAGAGCGTAATTGGCTCACGCAAAGTGTGGCGTTTCTCACCCAATGCCGCCAACAGCGATTTCACCGCCACCAATATCCATGTGACCTCGCACGGTACGGAATTTACCCTACAAACCCCAACCGGTAGCGTCGATGTTCTGCTGCCGTTGCCGGGGCGTCACAATATTGCGAATGCGCTGGCAG  >  NC_000913/95320‑95514
                                                                                                  |                                                                                                
ggCTGAACTGGCAGAGCGTAATTGGCTCACGCAAAGTGTGGCGTTTCTCACCCAATGCCGCCAACAGCGATTTCACCGCCACCAATATCCATGTGACCTc                                                                                                 >  2:3428732/1‑100 (MQ=255)
 gCTGAACTGGCAGAGCGTAATTGGCTCACGCAAAGTGTGGCGTTTCTCACCCAATGCCGCCAACAGCGATTTCACCGCCACCAATATCCATGTGACcgcg                                                                                                <  2:1450302/100‑1 (MQ=255)
  cTGAACTGGCAGAGCGTAATTGGCTCACGCAAAGTGTGGCGTTTCTCACCCAATGCCGCCAACAGCGATTTCACCGCCACCAATATCCATGTGACCgcgc                                                                                               <  2:3174536/100‑1 (MQ=255)
   tGAACTGGCAGAGCGTAATTGGCTCACGCAAAGTGTGGCGTTTCTCACCCAATGCCGCCAACAGCGATTTCACCGCCACCAATATCCATGTGACCGCGCa                                                                                              <  1:3409415/100‑1 (MQ=255)
   tGAACTGGCAGAGCGTAATTGGCTCACGCAAAGTGTGGCGTTTCTCACCCAATGCCGCCAACAGCGATTTCACCGCCACCAATATCCATGTGACCGCGCa                                                                                              <  1:3316385/100‑1 (MQ=255)
     aaCTGGCAGAGCGTAATTGGCTCACGCAAAGTGTGGCGTTTCTCACCCAATGCCGCCAACAGCGATTTCACCGCCACCAATATCCATGTGACCGCGCACg                                                                                            <  1:3432356/100‑1 (MQ=255)
      aCTGGCAGAGCGTAATTGGCTCACGCAAAGTGTGGCGTTTCTCACCCAATGCCGCCAACAGCGATTTCACCGCCACCAATATCCATGTGACCGCGCACgg                                                                                           >  1:1324096/1‑100 (MQ=255)
      aCTGGCAGAGCGTAATTGGCTCACGCAAAGTGTGGCGTTTCTCACCCAATCCCGCCAACAGCGATTTCACCGCCACCAATATCCATGTGACCGCGCACgg                                                                                           >  2:2436591/1‑100 (MQ=255)
       cTGGCAGAGCGTAATTGGCTCACGCAAAGTGTGGCGTTTCTCACCCAATGCCGCCAACAGCGATTTCACCGCCACCAATATCCATGTGACCGCGCACGGt                                                                                          >  1:1087366/1‑100 (MQ=255)
       cTGGCAGAGCGTAATTGGCTCACGCAAAGTGTGGCGTTTCTCACCCAATGCCGCCAACAGCGATTTCACCGCCACCAATATCCATGTGACCGCGCACGGt                                                                                          <  1:1362175/100‑1 (MQ=255)
         ggCAGAGCGTAATTGGCTCACGCAAAGTGTGGCGTTTCTCACCCAATGCCGCCAACAGCGATTTCACCGCCACCAATATCCATGTGACCGCGCACGGTAc                                                                                        >  2:2713130/1‑100 (MQ=255)
         ggCAGAGCGTAATTGGCTCACGCAAAGTGTGGCGTTTCTCACCCAATGCCGCCAACAGCGATTTCACCGCCACCAATATCCATGTGACCGCGCACGGTAc                                                                                        >  2:1061557/1‑100 (MQ=255)
         ggCAGAGCGTAATTGGCTCACGCAAAGTGTGACGTTTCTCGCCCAAGGCCGCCAACAGCGATTTCACCGCCACCAATATCCATGTGACCGCGCACGGTAc                                                                                        >  2:2565725/1‑100 (MQ=255)
          gCAGAGCGTAATTGGCTCACGCAAAGTGTGGCGTTTCTCACCCAATGCCGCCAACAGCGATTTCACCGCCACCAATATCCATGTGACCGCGCACGGTACg                                                                                       >  2:3520701/1‑100 (MQ=255)
            agagCGTAATTGGCTCACGCAAAGTGTGGCGTTTCTCACCCAATGCCGCCAACAGCGATTTCACCGCCACCAATATCCATGTGACCGCGCACGGTACGGa                                                                                     >  2:3757404/1‑100 (MQ=255)
            agagCGTAATTGGCTCACGCAAAGTGTGGCGTTTCTCACCCAATGCCGCCAACAGCGATTTCACCGCCACCAATATCCATGTGACCGCGCACGGTACGGa                                                                                     >  1:3995466/1‑100 (MQ=255)
              agCGTAATTGGCTCACGCAAAGTGTGGCGTTTCTCACCCAATGCCGCCAACAGCGATTTCACCGCCACCAATATCCATGTGACCGCGCACGGTACGGAAt                                                                                   <  1:624629/100‑1 (MQ=255)
                 gTAATTGGCTCACGCAAAGTGTGGCGTTTCTCACCCAATGCCGCCAACAGCGATTTCACCGCCACCAATATCCATGTGACCGCGCACGGTACGGAATTTa                                                                                <  2:2609480/100‑1 (MQ=255)
                   aaTTGGCTCACGCAAAGTGTGGCGTTTCTCACCCAATGCCGCCAACAGCGATTTCACCGCCACCAATATCCATGTGACCGCGCACGGTACGGAATTTAcc                                                                              <  1:2104708/100‑1 (MQ=255)
                   aaTTGGCTCACGCAAAGTGTGGCGTTTCTCACCCAATGCCGCCAACAGCGATTTCACCGCCACCAATATCCATGTGACCGCGCACGGTACGGAATTTAcc                                                                              <  1:3507899/100‑1 (MQ=255)
                    aTTGGCTCACGCAAAGTGTGGCGTTTCTCACCCAATGCCGCCAACAGCGATTTCACCGCCACCAATATCCATGTGACCGCGCACGGTACGGAATTTAccc                                                                             <  1:1027627/100‑1 (MQ=255)
                     ttGGCTCACGCAAAGTGTGGCGTTTCTCACCCAATGCCGCCAACAGCGATTTCACCGCCACCAATATCCATGTGACCGCGCACGGTACGGAATTTACCCt                                                                            <  2:2569846/100‑1 (MQ=255)
                     ttGGCTCACGCAAAGTGTGGCGTTTCTCACCCAATGCCGCCAACAGCGATTTCACCGCCACCAATATCCATGTGACCGCGCACGGTACGGAATTTACCCt                                                                            <  2:698543/100‑1 (MQ=255)
                            aCGCAAAGTGTGGCGTTTCTCACCCAATGCCGCCAACAGCGATTTCACCGCCACCAATATCCATGTGACCGCGCACGGTACGGAATTTACCCTACAAAcc                                                                     <  1:931830/100‑1 (MQ=255)
                            aCGCAAAGTGTGGCGTTTCTCACCCAATGCCGCCAACAGCGATTTCACCGCCACCAATATCCATGTGACCGCGCACGGTACGGAATTTACCCTACAAAcc                                                                     <  1:1186374/100‑1 (MQ=255)
                            aCGCAAAGTGTGGCGTTTCTCACCAAATGGCGCCAACAGCGATTTCACCGCCACCAATATCCATGTGACCGCGCACGGTACGGAATTTACCCTACAAAcc                                                                     <  2:397657/100‑1 (MQ=255)
                              gCAAAGTGTGCCGTTTCTCACCCAATGCCGCCAACAGCGATTTCACCGCCACCAATATCCATGTGACCGCGCACGGTACGGAATTTACCCTACAAAcccc                                                                   <  2:3960564/100‑1 (MQ=255)
                               cAAAGTGTGGCGTTTCTCACCCAATGCCGCCAACAGCGATTTCACCGCCACCAATATCCATGTGACCGCGCACGGTACGGAATTTACCCTACAAACCCCa                                                                  <  2:1625186/100‑1 (MQ=255)
                               cAAAGTGTGGCGTTTCTCACCCAATGCCGCCAACAGCGATTTCACCGCCACCAATATCCATGTGACCGCGCACGGTACGGAATTTACCCTACAAACCCCa                                                                  >  2:299534/1‑100 (MQ=255)
                                  aGTGTGGCGTTTCTCACCCAATGCCGCCAACAGCGATTTCACCGCCACCAATATCCATGTGACCGCGCACGGTACGGAATTTACCCTACAAACCCCAAcc                                                               <  1:3994912/100‑1 (MQ=255)
                                   gtgtGGCGTTTCTCACCCAATGCCGCCAACAGCGATTTCACCGCCACCAATATCCATGTGACCGCGCACGGTACGGAATTTACCCTACAAACCCCAACCg                                                              <  2:2263049/100‑1 (MQ=255)
                                     gtgGCGTTTCTCACCCAATGCCGCCAACAGCGATTTCACCGCCACCAATATCCATGTGACCGCGCACGGTACGGAATTTACCCTACAAACCCCAACCGGt                                                            >  2:749419/1‑100 (MQ=255)
                                      tgGCGTTTCTCACCCAATGCCGCCAACAGCGATTTCACCGCCACCAATATCCATGTGACCGCGCACGGTACGGAATTTACCCTACAAACCCCAACCGGTa                                                           <  2:3928048/100‑1 (MQ=255)
                                       ggCGTTTCTCACCCAATGCCGCCAACAGCGATTTCACCGCCACCAATATCCATGTGACCGCGCACGGTACGGAATTTACCCTACAAACCCCAACCGGTAg                                                          <  1:3103463/100‑1 (MQ=255)
                                        gCGTTTCTCACCCAATGCCGCCAACAGCGATTTCACCGCCACCAATATCCATGTGACCGCGCACGGTACGGAATTTACCCTACAAACCCCAACCGGTAGc                                                         >  1:1648796/1‑100 (MQ=255)
                                         cGTTTCTCACCCAATGCCGCCAACAGCGATTTCACCGCCACCAATATCCATGTGACCGCGCACGGTACGGAATTTACCCTACAAACCCCAACCGGTAGCg                                                        >  2:2706008/1‑100 (MQ=255)
                                              ctcACCCAATGCCGCCAACAGCGATTTCACCGCCACCAATATCCATGTGACCGCGCACGGTACGGAATTTACCCTACAAACCCCAACCGGTAGCGTCGAt                                                   <  2:64135/100‑1 (MQ=255)
                                              ctcAACCAATGCCGCCAACAGCGATTTCACCGCCACCAATATCCATGTGACCGCGCACGGTACGGAATTTACCCTACAAACCCCCACGGGTAGCGTCGAt                                                   <  2:1033778/100‑1 (MQ=255)
                                               tcACCCAATGCCGCCAACAGCGATTTCACCGCCACCAATATCCATGTGACCGCGCACGGTACGGAATTTACCCTACAAACCCCAACCGGTAGCGTCGATg                                                  >  1:841604/1‑100 (MQ=255)
                                               tcACCCAATGCCGCCAACAGCGATTTCACCGCCACCAATATCCATGTGACCGCGCACGGTACGGAATTTACCCTACAAACCCCAACCGGTAGCGTCGATg                                                  <  1:3369265/100‑1 (MQ=255)
                                                cACCCAATGCCGCCAACAGCGATTTCACCGCCACCAATATCCATGTGACCGCGCACGGTACGGAATTTACCCTACAAACCCCAGCCGGTAGCGTCGATGt                                                 <  1:407517/100‑1 (MQ=255)
                                                 aCCCAATGCCGCCAACAGCGATTTCACCGCCACCAATATCCCTGTCACCGCGCACGGTACGGACTTTACCCTACAAACCCCAACCGGTAGCGACCATGtt                                                >  2:531398/1‑100 (MQ=255)
                                                 aCCCAATGCCGCCAACAGCGATTTCACCGCCACCAATATCCATGTGACCGCGCACGGTACGGAATTTACCCTACAAACCCCAACCGGTAGCGTCGATGtt                                                >  1:3489627/1‑100 (MQ=255)
                                                  cccAATGCCGCCAACAGCGATTTCACCGCCACCAATATCCATGTGACCGCGCACGGTACGGAATTTACCCTACAAACCCCAACCGGTAGCGTCGATGTTc                                               <  1:963124/100‑1 (MQ=255)
                                                  cccAATGCCGCCAACAGCGATTTCACCGCCACCAATATCCATGTGACCGCGCACGGTACGGAATTTACCCTACAAACCCCAACCGGTAGCGTCGATGTTc                                               >  2:1925504/1‑100 (MQ=255)
                                                   ccACTGCCGCCAACAGCGATTTCACCGCCACCAATATCCATGTGACCGCGCACGGTACGGAATTTACCCTACAAACCCCAACCGGTAGCGTCGATGTTct                                              <  1:2299185/100‑1 (MQ=255)
                                                    cAATGCCGCCAACAGCGATTTCACCGCCACCAATATCCATGTGACCGCGCACGGTACGGAATTTACCCTACAAACCCCAACCGGTAGCGTCGATGTTctg                                             >  1:3793788/1‑100 (MQ=255)
                                                     aaTGCCGCCAACAGCGATTTCACCGCCACCAATATCCATGTGACCGCGCACGGTACGGAATTTACCCTACAAACCCCAACCGGTAGCGTCGATGTTctgc                                            <  1:171810/100‑1 (MQ=255)
                                                     aaTGCCGCCAACAGCGATTTCACCGCCACCAATATCCATGTGACCGCGCACGGTACGGAATTTACCCTACAAACCCCAACCGGTAGCGTCGATGTTctgc                                            <  1:1995099/100‑1 (MQ=255)
                                                        gccgccAACAGCGATTTCACCGCCACCAATATCCATGTGACCGCGCACGGTACGGAATTTACCCTACAAACCCCAACCGGTAGCGTCGATGTTCtgctgc                                         >  1:759337/1‑100 (MQ=255)
                                                         ccgccAACAGCGATTTCACCGCCACCAATATCCATGTGACCGCGCACGGTACGGAATTTACCCTACAAACCCCAACCGGTAGCGTCGATGTTCTGCTGcc                                        >  1:1827838/1‑100 (MQ=255)
                                                          cgccAACAGCGATTTCACCGCCACCAATATCCATGTGACCGCGCACGGTACGGAATTTACCCTACAAACCCCAACCGGTAGCGTCGATGTTCTGCTGCCg                                       >  1:563203/1‑100 (MQ=255)
                                                          cgccAACAGCGATTTCACCGCCACCAATATCCATGTGACCGCGCACGGTACGGAATTTACCCTACAAACCCCAACCGGTAGCGTCGATGTTCTGCTGCCg                                       >  1:218937/1‑100 (MQ=255)
                                                           gccAACAGCGATTTCACCGCCACCAATATCCATGTGACCGCGCACGGTACGGAATTTACCCTACAAACCCCAACCGGTAGCGTCGATGTTCTGCTGCCGt                                      <  1:1952916/100‑1 (MQ=255)
                                                               aCAGCGATTTCACCGCCACCAATATCCATGTGACCGCGCACGGTACGGAATTTACCCTACAAACCCCAACCGGTAGCGTCGATGTTCTGCTGCCGTTGcc                                  <  1:2206140/100‑1 (MQ=255)
                                                                cAGCGATTTCACCGCCACCAATATCCATGTGACCGCGCACGGTACGGAATTTACCCTACAAACCCCAACCGGTAGCGTCGATGTTCTGCTGCCGTTGCCg                                 >  2:1592190/1‑100 (MQ=255)
                                                                cAGCGATTTCACCGCCACCAATATCCATGTGACCGCGCACGGTACGGAATTTACCCTACAAACCCCAACCGGTAGCGTCGATGTTCTGCTGCCGTTGCCg                                 >  2:3017196/1‑100 (MQ=255)
                                                                cAGCGATTTCACCGCCACCAATATCCATGTGACCGCGCACGGTACGGAATTTACCCTACAAACCCCAACCGGTAGCGTCGATGTTCTGCTGCCGTTGCCg                                 >  2:254534/1‑100 (MQ=255)
                                                                cAGCGATTTCACCGCCACCAATATCCATGTGACCGCGCACGGTACGGAATTTACCCTACAAACCCCAACCGGTAGCGTCGATGTTCTGCTGCCGTTGCCg                                 >  2:2266836/1‑100 (MQ=255)
                                                                  gCGATTTCACCGCCACCAATATCCATGTGACCGCGCACGGTACGGAATTTACCCTACAAACCCCAACCGGTAGCGTCGATGTTCTGCTGCCGTTGCCggg                               >  1:3188781/1‑100 (MQ=255)
                                                                  gCGATTTCACCGCCACCAATATCCATGTGACCGCGCACGGTACGGAATTTACCCTACAAACCCCAACCGGTAGCGTCGATGTTCTGCTGCCGTTGCCggg                               <  2:562728/100‑1 (MQ=255)
                                                                  gCGATTTCACCGCCACCAATATCCATGTGACCGCGCACGGTACGGAATTTACCCTACAAACCCCAACCGGTAGCGTCGATGTTCTGCTGCCGTTGCCggg                               >  2:2608253/1‑100 (MQ=255)
                                                                  gCGATTTCACCGCCACCAATATCCATGTGACCGCGCACGGTACGGAATTTACCCTACAAACCCCAACCGGTAGCGTCGATGTTCTGCTGCCGTTGCCggg                               >  1:3642786/1‑100 (MQ=255)
                                                                   cGATTTCACCGCCACCAATATCCATGTGACCGCGCACGGTACGGAATTTACCCTACAAGCCCCAGCCGGTAGCGTCGATGTTCTGCTGCCGTTGCCgggg                              <  2:3140171/100‑1 (MQ=255)
                                                                   cGATTTCACCGCCACCAATATCCATGTGACCGCGCACGGTACGGAATTTACCCTACAAACCCCAACCGGTAGCGTCGATGTTCTGCTGCCGTTGCCgggg                              >  1:2212003/1‑100 (MQ=255)
                                                                    gATTTCACCGCCACCAATATCCATGTGACCGCGCACGGTACGGAATTTACCCTACAAACCCCAACCGGTAGCGTCGATGTTCTGCTGCCGTTGCCGGGGc                             >  1:3781228/1‑100 (MQ=255)
                                                                    gATTTCACCGCCACCAATATCCATGTGACCGCGCACGGTACGGAATTTACCCTACAAACCCCAACCGGTAGCGTCGATGTTCTGCTGCCGTTGCCGGGGc                             >  2:2469936/1‑100 (MQ=255)
                                                                    gATTTCACCGCCACCAATATCCATGTGACCGCGCACGGTACGGAATTTACCCTACAAACCCCAACCGGTAGCGTCGATGTTCTGCTGCCGTTGCCGGGGc                             >  2:1517986/1‑100 (MQ=255)
                                                                     aTTTCACCGCCACCAATATCCATGTGAGCGCGCACCTTACGCAATAAACCCTACAAACCCCAACCGGCAGCGTCGATGTTCTGCGGCCCTTGCCGGGGCg                            >  2:2184132/1‑100 (MQ=255)
                                                                     aTTTCACCGCCACCAATATCCATGTGACCGCGCACGGTACGGAATTTACCCTACAAACCCCAGCCGGTAGCGTCGATGTTCTGCTGCCGTTGCCGGGGCg                            <  2:1050425/100‑1 (MQ=255)
                                                                     aTTTCACCGCCACCAATATCCATGTGACCGCGCACGGTACGGAATTTACCCTACAAACCCCACCTGGTAGCGTCGATGTTCTGCTGCCGTTGCCGGGGCg                            <  2:3905336/100‑1 (MQ=255)
                                                                     aTTTCACCGCCACCAATATCCATGTGACCGCGCACGGTACGGAATTTACCCTACAAACCCCAACCGGTAGCGTCGATGTTCTGCTGCCGTTGCCGGGGCg                            <  2:1284267/100‑1 (MQ=255)
                                                                     aTTTCACCGCCACCAATATCCATGTGACCGCGCACGGTACGGAATTTACCCTACAAACCCCAACCGGTAGCGTCGATGTTCTGCTGCCGTTGCCGGGGCg                            <  2:3174624/100‑1 (MQ=255)
                                                                     aTTTCACCGCCACCAATATCCATGTGACCGCGCACGGTACGGAATTTACCCTACAAACCCCAACCGGTAGCGTCGATGTTCTGCTGCCGTTGCCGGGGCg                            <  2:3042175/100‑1 (MQ=255)
                                                                      tttCACCGCCACCAATATCCATGTGACCGCGCACGGTACGGAATTTACCCTACAAACCCCAACCGGTAGCGTCGATGTTCTGCTGCCGTTGCCGGGGCGt                           >  1:2420977/1‑100 (MQ=255)
                                                                      tttCACCGCCACCAATAGCCATGTGACCGCGCACGGTACGGAATTTACCCTAGAAACCCCAACCGGTAGCGTCGATGTTCTGCTGCCGTTGCCGGGGCGt                           <  2:2234194/100‑1 (MQ=255)
                                                                       ttCACCGCCACCAATATCCATGTGACCGCGCACGGTACGGAATTTACCCTACAAACCCCAACCGGTAGCGTCGATGTTCTGCTGCCGTTGCCGGGGCGTc                          >  2:2018695/1‑100 (MQ=255)
                                                                       ttCACCGCCACCAATATCCATGTGACCGCGCACGGTACGGAATTTACCCTACAAACCCCAACCGGTAGCGTCGATGTTCTGCTGCCGTTGCCGGGGCGTc                          <  2:2889074/100‑1 (MQ=255)
                                                                       ttCACCGCCACCAATATCCATGTGACCGCGCACGGTACGGAATTTACCCTACAAACCCCAACCGGTAGCGTCGATGTTCTGCTGCCGTTGCCGGGGCGTc                          >  2:2019612/1‑100 (MQ=255)
                                                                        tCACCGCCACCAATATCCATGTGACCGCGCACGGTACGGAATTTACCCTACAAACCCCAACCGGTAGCGTCGATGTTCTGCTGCCGTTGCCGGGGCGTca                         >  1:3547181/1‑100 (MQ=255)
                                                                        tCACCGCCACCAATATCCATGTGACCGCGCACGGTACGGAATTTACCCTACAAACCCCAACCGGTAGCGTCGATGTTCTGCTGCCGTTGCCGGGGCGTca                         >  2:2165007/1‑100 (MQ=255)
                                                                         cACGGCCACCAATATCCATGTGACCGCGCACGGTACGGCATTTACCCTACAAGTTCCAACCGGCAGAGTCGATGTTCTGATGCCGTTGCCGGGGCGCcac                        >  2:441100/1‑100 (MQ=255)
                                                                         cACCGCCACCAATATCCATGTGACCGCGCACGGTACGGAATTTACCCTACAAACCCCAACCGGTAGCGTCGATGTTCTGCCGCCGTTGCCGGGGCGTcac                        >  1:3364606/1‑100 (MQ=255)
                                                                          aCCGCCACCAATATCCATGTGACCGCGCACGGTACGGAATTTACCCTACAAACCCCAACCGGTAGCGTCGATGTTCTGCTGCCGTTGCCGGGGCGTcaca                       >  2:2081533/1‑100 (MQ=255)
                                                                          aCCGCCACCAATATCCATGTGACCGCGCACGGTACGGAATTTACCCTACAAACCCCAACCGGTAGCGTCGATGTTCTGCTGCCGTTGCCGGGGCGTcaca                       <  2:1621252/100‑1 (MQ=255)
                                                                          aCCGCCACCAATATCCATGTGACCGCGCACGGTACGGAATTTACCCTACAAACCCCAACCGGTAGCGTCGATGTTCTGCTGCCGTTGCCGCGGCGCcaca                       >  2:2585265/1‑100 (MQ=255)
                                                                          aCCGCCACCAATATCCATGTCACCGCGCACGGTACGGAATTTACCCTACAAACCCCAACCGGTAGCGTCGATGTTCTCCTGCCGTTGCCGGGGCGTcaca                       >  2:2181716/1‑100 (MQ=255)
                                                                             agcaccaATATCCATGTGACCGCGCACGGTACGGAATTTACCCTACAAACCCCAACCGGTAGCGTCGATGTTCTGCTGCCGTTGCCGGGGCGTCACAata                    <  2:3602975/98‑1 (MQ=255)
                                                                               caccaATATCCATGTGACCGCGCACGGTACGGAATTTACCCTACAAACCCCAACCGGTAGCGTCGATGTTCTGCTGCCGTTGCCGGGGCGTCACAATAtt                  >  2:3773098/1‑100 (MQ=255)
                                                                                 ccaATATCCATGTGACCGCGCACGGTACGGAATTTACCCTACAAACCCCAACCGGTAGCGTCGATGTTCTGCTGCCGTTGCCGGGGCGTCACAATATTGc                >  2:1806878/1‑100 (MQ=255)
                                                                                 ccaATATCCATGTGACCGCGCACGGTACGGAATTTACCCTACAAACCCCAACCGGTAGCGTCGATGTTCTGCTGCCGTTGCCGGGGCGTCACAATATTGc                <  2:1798349/100‑1 (MQ=255)
                                                                                  caATATCCATGTGACCGCGCACGGTACGGAATTTACCCTACAAACCCCAACCGGTAGCGTCGATGTTCTGCTGCCGTTGCCGGGGCGTCACAATATTGCg               <  2:1339968/100‑1 (MQ=255)
                                                                                   natatCCATGTGACCGCGCACGGTACGGAATTTACCCTACAAACCCCAACCGGTAGCGTCGATGTTCTGCTGCCGTTGCCGGGGCGTCACAATATTGCGa              >  1:3231957/2‑100 (MQ=255)
                                                                                      atCCATGTGACCGCGCACGGTACGGAATTTACCCTACAAACCCCAACCGGTAGCGTCGATGTTCTGCTGCCGTTGCCGGGGCGTCACAATATTGCGAATg           <  1:2330547/100‑1 (MQ=255)
                                                                                       tCCATGTGACCGCGCACGGTACGGAATTTACCCTACAAACCCCAACCGGTAGCGTCGATGTTCTGCTGCCGTTGCCGGGGCGTCACAATATTGCGAATgc          >  1:2353692/1‑100 (MQ=255)
                                                                                        ccATGTGACCGCGCACGGTACGGAATTTACCCTACAAACTCCAACCGGTAGCGTCCATGTTCTGCTGCCGTTGCCGGGGCGTCACAATATTGAGAATgcg         >  1:3748345/1‑100 (MQ=255)
                                                                                             tgACCGCGCACGGTACGGAATTTACCCTACAAACCCCAACCGGTAGCGTCGATGTTCTGCTGCCGTTGCCGGGGCGTCACAATATTGCGAATGCGCTGGc    >  1:2694291/1‑100 (MQ=255)
                                                                                              gACCGCGCACGGTACGGAATTTACCCTACAAACCCCAACCGGTAGCGTCGATGTTCTGCTGCCGTTGCCGGGGCGTCACAATATTGCGAATGCGCTGGCa   <  2:637764/100‑1 (MQ=255)
                                                                                              gACCGCGCACGGTACGGAATTTACCCTACAAACCCCAACCGGTAGCGTCGATGTTCTGCTGCCGTTGCCGGGGCGTCACAATATTGCGAATGCGCTGGCa   >  1:2395916/1‑100 (MQ=255)
                                                                                               aCCGCGCACGGTACGGAATTTACCCTACAAACCCCAACCGGTAGCGTCGATGTTCTGCTGCCGTTGCCGGGGCGTCACAATATTGCGAGTGCGCTCGCAg  >  2:2532191/1‑100 (MQ=255)
                                                                                                  |                                                                                                
GGCTGAACTGGCAGAGCGTAATTGGCTCACGCAAAGTGTGGCGTTTCTCACCCAATGCCGCCAACAGCGATTTCACCGCCACCAATATCCATGTGACCTCGCACGGTACGGAATTTACCCTACAAACCCCAACCGGTAGCGTCGATGTTCTGCTGCCGTTGCCGGGGCGTCACAATATTGCGAATGCGCTGGCAG  >  NC_000913/95320‑95514

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: