Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,827,774 G→C 5.8% A211P (GCC→CCC)  mtgA → monofunctional biosynthetic peptidoglycan transglycosylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,827,7740GC 5.8% 188.6 / 2.4 62A211P (GCC→CCC) mtgAmonofunctional biosynthetic peptidoglycan transglycosylase
Reads supporting (aligned to +/- strand):  ref base G (30/27);  new base C (2/2);  total (32/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

GCCCAGGCAGCGGCGCGTTATCACTTTGGCGTCGACGCCAGCCGGCTCAGCCGCCAGCAGGCAGCGCAGCTGGCAGCGGTGCTGCCAAGCCCAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAGGCCCGGCCCTTTCGCGGGTAAACCCGCTCCTACGGGCAACCGCATCCCCGGCCTCCTCGCGG  >  NC_002947/5827636‑5827922
                                                                                                                                           |                                                                                                                                                    
gCCCAGGCAGCGGCGCGTTATCACTTTGGCGTCGACGCCAGCCGGCTCAGCCGCCAGCAGGCAGCGCAGCTGGCAGCGGTGCTGCCAAGCCCAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCGGCTGGAtga                                                                                                                                          <  2:821513/150‑3 (MQ=255)
   cAGGCAGCGGCGCGTTATCACTTTGGCGTCGACGCCAGCCGGCTCAGCCGCCAGCAGGCAGCGCAGCTGGCAGCGGTGCTGCCAAGCCCAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGACGGGCCGGCTGGATTCGc                                                                                                                                        >  1:754667/1‑150 (MQ=255)
       cAGCGGCGCGTTATCACTTTGGCGTCGACGCCAGCCGGCTCAGCCGCCAGCAGGCAGCGCAGCTGGCAGCGGTGCTGCCAAGCCCAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGc                                                                                                                                    >  2:1282778/1‑150 (MQ=255)
           ggcgCGTTATCACTTTGGCGTCGACGCCAGCCGGCTCAGCCGCCAGCAGGCAGCGCAGCTGGCAGCGGTGCTGCCAAGCCCAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGAt                                                                                                                                <  1:532749/150‑1 (MQ=255)
               cgTTATCACTTTGGCGTCGACGCCAGCAGGATCAGCGGCTAGCAGGAGATCCGCTGGGCAGCGGTGCTGCCAAGCACAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGGGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGc                                                                                                                            <  2:457221/150‑1 (MQ=255)
                        tttGGCGTCGACGCCAGCCGGCTCAGCCGCCAGCAGGCAGCGCAGCTGGCAGCGGTGCTGCCAAGCCCAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTgggc                                                                                                                   <  2:953358/150‑1 (MQ=255)
                        tttGGCGTCGACGCCAGCCGGCTCAGCCGCCAGCAGGCAGCGCAGCTGGCAGCGGTGCTGCCAAGCCCAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCAGCAGGTGAGGCAGCTgggc                                                                                                                   >  2:819517/1‑150 (MQ=255)
                          tGGCGTCGACGCCAGCCGGCTCAGCCGCCAGCAGGCAGCGCAGCTGGCAGCGGTGCTGCCAAGCCCAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTgggcgg                                                                                                                 <  2:883454/150‑1 (MQ=255)
                             cGTCGACGCCAGCCGGCTCAGCCGCCAGCAGGCAGCGCAGCTGGCAGCGGTGCTGCCAAGCCCAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGcc                                                                                                              >  2:241541/1‑150 (MQ=255)
                              gTCGACGCCAGCCGGCTCAGCCGCCAGCAGGCAGCGCAGCTGGCAGCGGTGCTGCCAAGCCCAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGccc                                                                                                             >  2:638790/1‑150 (MQ=255)
                               tCGACGCCAGCCGGCTCAGCCGCCAGCAGGCAGCGCAGCTGGCAGCGGTGCTGCCAAGCCCAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCa                                                                                                            <  2:214478/150‑1 (MQ=255)
                                  aCGCCAGCCGGCTCAGCCGCCAGCAGGCAGCGCAGCTGGCAGCGGTGCTGCCAAGCCCAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCt                                                                                                         <  2:656887/150‑1 (MQ=255)
                                  aCGCCAGCCGGCTCAGCCGCCAGCAGGCAGCGCAGCTGGCAGCGGTGCTGCCAAGCCCAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCt                                                                                                         <  1:241541/150‑1 (MQ=255)
                                    gCCAGCCGGCTCAGCCGCCAGCAGGCAGCGCAGCTGGCAGCGGTGCTGCCAAGCCCAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGa                                                                                                                                             <  1:343920/112‑1 (MQ=255)
                                    gCCAGCCGGCTCAGCCGCCAGCAGGCAGCGCAGCTGGCAGCGGTGCTGCCAAGCCCAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGa                                                                                                                                             >  2:343920/1‑112 (MQ=255)
                                    gCCAGCCGGCTCAGCCGCCAGCAGGCAGCGCAGCTGGCAGCGGTGCTGCCAAGCCCAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGTTAAGACAGCGGGGGGGGCGCAGCTAc                                                                                                       >  2:785455/1‑150 (MQ=255)
                                       aGCCGGCTCAGCCGCCAGCAGGCAGCGCAGCTGGCAGCGGTGCTGCCAAGCCCAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTg                                                                                                    >  1:161442/1‑150 (MQ=255)
                                        gCCGGCTCAGCCGCCAGCAGGCAGCGCAGCTGGCAGCGGTGCTGCCAAGCCCAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGa                                                                                                   <  1:638790/150‑1 (MQ=255)
                                            gCTCAGCCGCCAGCAGGCAGCGCAGCTGGCAGCGGTGCTGCCAAGCCCAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCa                                                                                               >  1:2592/1‑150 (MQ=255)
                                             cTCAGCCGCCAGCAGGCAGCGCAGCTGGCAGCGGTGCTGCCAAGCCCAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAg                                                                                              <  1:233969/150‑1 (MQ=255)
                                               cAGCCGCCAGCAGGCAGCGCAGCTGGCAGCGGTGCTGCCAAGCCCAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATNCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCt                                                                                            <  1:819517/150‑1 (MQ=255)
                                                aGCCGCCAGCAGGCAGCGCAGCTGGCAGCGGTGCTGCCAAGCCCAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTg                                                                                           <  1:135390/150‑1 (MQ=255)
                                                          aggcagCGCAGCTGGCAGCGGTGCTGCCAAGCCCAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCa                                                                                 <  2:1602900/150‑1 (MQ=255)
                                                          aggcagCGCAGCTGGCAGCGGTGCTGCCAAGCCCAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCCGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCa                                                                                 >  1:1593308/1‑150 (MQ=255)
                                                            gcagcgcagcTGGCAGCGGTGCTGCCAAGCCCAATCAAGTGGAGTGCCAGCCGACCGAGTGCCCACGTGGCCAGCCGGGCCGGCTGGAGTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGa                                                                               >  2:938559/1‑150 (MQ=255)
                                                            gcagcgatgcTGGCAGCGGTGCTGCCAAGCCCAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGa                                                                               >  2:1257757/1‑150 (MQ=255)
                                                              agcgcagcTGGCAGCGGTGCTGCCAAGCCCAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAgccgggc                                                                                                                                                     <  2:940113/78‑1 (MQ=255)
                                                              agcgcagcTGGCAGCGGTGCTGCCAAGCCCAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAgccgggc                                                                                                                                                     >  1:940113/1‑78 (MQ=255)
                                                              agcgcagcTGGCAGCGGTGCTGCCAAGCCCAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGaaa                                                                             >  2:1426500/1‑150 (MQ=255)
                                                              agcgcagcTGGCAGCGGTGCTGCCAAGCCCAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGCCCGGCTggccac                                                                                                                                          <  2:55052/89‑5 (MQ=255)
                                                              agcgcagcTGGCAGCGGTGCTGCCAAGCCCAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGCCCGGCTggccac                                                                                                                                          >  1:55052/1‑85 (MQ=255)
                                                                     ttGGCAGCACCGCTGCCAAGCCCAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGa                                                                      <  2:1155952/149‑1 (MQ=255)
                                                                                   tctcAAGCCCAATCAAGTGGAGTGCCAGCCGACCGCGTTCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAggcccggc                                                         <  2:794084/147‑1 (MQ=255)
                                                                                      cAAGCCCAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAGGCCCGGCCCt                                                      <  1:1282778/150‑1 (MQ=255)
                                                                                       aaGCCCAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCAGTCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAGGCCCGGCCCtt                                                     <  1:918723/150‑1 (MQ=255)
                                                                                       aaGCCCAATCAAGTGGAGTGCCAGCCGACCGAGTAGCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAggcccggccaga                                                     >  1:794084/1‑147 (MQ=255)
                                                                                         gCCCAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGAt                                                                                                  >  1:1632417/1‑103 (MQ=255)
                                                                                         gCCCAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGAt                                                                                                  <  2:1632417/103‑1 (MQ=255)
                                                                                            cAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAGGCCCGGCCCTTTcgcg                                                >  2:398358/1‑150 (MQ=255)
                                                                                             aaTCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGc                                                                                                          >  2:989915/1‑91 (MQ=255)
                                                                                             aaTCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGc                                                                                                          <  1:989915/91‑1 (MQ=255)
                                                                                              aTCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAggcccgg                                                          <  1:282030/138‑1 (MQ=255)
                                                                                              aTCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAggcccgg                                                          >  2:282030/1‑138 (MQ=255)
                                                                                              aTCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAGGCCCGGCCCTTTCGCggg                                              >  1:1420667/1‑150 (MQ=255)
                                                                                                  aGTGGAGTTCCAGCCGACCGATTCCCTACGTACCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATCCAGGCCCGGCCCTTTCGCGGGTaaa                                          <  1:785455/150‑1 (MQ=255)
                                                                                                        gTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAGGCCCGGCCCTTTCGCGGGTAAACCCGCt                                    >  2:575741/1‑150 (MQ=255)
                                                                                                        gTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAGGCCCGGCCCTTTCGCGGGTAAACCCGCt                                    >  2:575825/1‑150 (MQ=255)
                                                                                                         tGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAGGCCCGGCCCTTTCGCGGGTAAACCCGCTc                                   >  1:647086/1‑150 (MQ=255)
                                                                                                           ccAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGGGGGCGCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAGGCCCGGCCCTTTTGCGGGTAAACTCgtttct                                 >  2:781334/1‑145 (MQ=255)
                                                                                                            cAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACtcctcc                                                                                  >  2:539434/1‑100 (MQ=255)
                                                                                                            cAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACtcctcc                                                                                  <  1:539434/100‑1 (MQ=255)
                                                                                                            cAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAGGCCCGGCCCtt                                                     <  2:263245/129‑1 (MQ=255)
                                                                                                            cAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAGGCCCGGCCCtt                                                     >  1:263245/1‑129 (MQ=255)
                                                                                                            cAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAGGCCCGGCCCTTTCGCGGGTTTACCCGcgaaag                                <  1:398358/150‑6 (MQ=255)
                                                                                                             aGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAGGCCCGGCCCTTTCGCGGGTAAACCCGCTCCTAc                               <  1:1426500/150‑1 (MQ=255)
                                                                                                                          ccTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAGGCCCGGCCCTTTCGCGGGTAAACCCGCTCCTACggg                            <  1:931076/140‑1 (MQ=255)
                                                                                                                          ccTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAGGCCCGGCCCTTTCGCGGGTAAACCCGCTCCTACggg                            >  2:931076/1‑140 (MQ=255)
                                                                                                                           cTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAGGCCCGGCCCTTTCGCGGGTAAACCCGCTCCTACGGGCAACCGCATcc                 >  2:1052316/1‑150 (MQ=255)
                                                                                                                           cTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAGGCCCGGCCCTTTCGCGGGTAAACCCGCTCCTACGGGCAACCGCATcc                 >  2:419074/1‑150 (MQ=255)
                                                                                                                             aCGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAGGCCCGGCCCTTTCGCGGGTAAACCCGCTCCTACGGGCAACCGCATcccc               >  1:992066/1‑150 (MQ=255)
                                                                                                                               gaatccAGCCGGCCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAGGCCCGGCCCTTTCGCGGGTAAACCCGCTCCTACgg                             >  2:1481264/5‑134 (MQ=255)
                                                                                                                               gaatccAGCCGGCCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAGGCCCGGCCCTTTCGCGGGTAAACCCGCTCCTACgg                             <  1:1481264/130‑1 (MQ=255)
                                                                                                                                tGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAGGCCCGGCCCTTTCGCGGGTAAACCCGCTCCTACGGGCAACCGCATCCCCGGc            <  1:1257757/150‑1 (MQ=255)
                                                                                                                                 ggCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAGGCCCGGCCCTTTCGCGGGTAAACCCGCTCCTACGGGCAACCGCATCCCCGGcc           >  2:1393946/1‑150 (MQ=255)
                                                                                                                                        cgggccggCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAGGCCCGGCCGTTTCGCGGGTAAACCCGCTCCTACGGGCAACCGCATCCCCGGCCTCCTcgc    <  1:938559/150‑1 (MQ=255)
                                                                                                                                          ggccggCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAGGCCCGGCCCTTTCGCGGGTAAACCCGCTCCTACGGGCAACCGCATCCCCGGCCTCCTCGCgg  <  2:161442/150‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                                    
GCCCAGGCAGCGGCGCGTTATCACTTTGGCGTCGACGCCAGCCGGCTCAGCCGCCAGCAGGCAGCGCAGCTGGCAGCGGTGCTGCCAAGCCCAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAGGCCCGGCCCTTTCGCGGGTAAACCCGCTCCTACGGGCAACCGCATCCCCGGCCTCCTCGCGG  >  NC_002947/5827636‑5827922

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: