Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,827,811 G→C 6.1% G223A (GGG→GCG)  mtgA → monofunctional biosynthetic peptidoglycan transglycosylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,827,8110GC 6.1% 184.4 / 2.4 60G223A (GGG→GCG) mtgAmonofunctional biosynthetic peptidoglycan transglycosylase
Reads supporting (aligned to +/- strand):  ref base G (26/29);  new base C (2/2);  total (29/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCGTCGACGCCAGCCGGCTCAGCCGCCAGCAGGCAGCGCAGCTGGCAGCGGTGCTGCCAAGCCCAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAGGCCCGGCCCTTTCGCGGGTAAACCCGCTCCTACGGGCAACCGCATCCCCGGCCTCCTCGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCC  >  NC_002947/5827663‑5827953
                                                                                                                                                    |                                                                                                                                              
ggCGTCGACGCCAGCCGGCTCAGCCGCCAGCAGGCAGCGCAGCTGGCAGCGGTGCTGCCAAGCCCAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTgggcggg                                                                                                                                               <  2:350899/150‑1 (MQ=255)
ggCGTCGACGCCAGCCGGCTCAGCCGCCAGCAGGCAGCGCAGCTGGCAGCGGTGCTGCCAAGCCCAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTgggcggg                                                                                                                                               <  1:642506/150‑1 (MQ=255)
 gCGTCGACGCCAGCCGGCTCAGCCGCCAGCAGGCAGCGCAGCTGGCAGCGGTGCTGCCAAGCCCAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTgggcgggc                                                                                                                                              <  2:1253751/150‑1 (MQ=255)
         gCCAGCCGGCTCAGCCGCCAGCAGGCAGCGCAGCTGGCAGCGGTGCTGCCAAGCCCAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCCCGCCCAGCtgg                                                                                                                                      <  2:988727/150‑3 (MQ=255)
          ccAGCCGGCTCAGCCGCCAGCAGGCAGCGCAGCTGGCAGCGGTGCTGCCAAGCCCAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCCGCTAcc                                                                                                                                     <  2:216497/150‑1 (MQ=255)
              ccGGCTCAGCCGCCAGCAGGCAGCGCAGCTGGCAGCGGTGCTGCCAAGCCCAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTAGCCGGCAGATCATCCATCGCGGCGTACCCAGCCACCTGAt                                                                                                                                 >  2:462483/1‑150 (MQ=255)
                   tCAGCCGCCAGCAGGCAGCGCAGCTGGCAGCGGTGCTGCCAAGCCCAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTAc                                                                                                                                      <  2:559284/140‑1 (MQ=255)
                   tCAGCCGCCAGCAGGCAGCGCAGCTGGCAGCGGTGCTGCCAAGCCCAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTAc                                                                                                                                      >  1:559284/1‑140 (MQ=255)
                   tCAGCCGCCAGCAGGCAGCGCAGCTGGCAGCGGTGCTGCCAAGCCCAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGc                                                                                                                            <  1:961617/150‑1 (MQ=255)
                    cAGCCGCCAGCAGGCAGCGCAGCTGGCAGCGGTGCTGCCAAGCCCAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCt                                                                                                                           <  2:626913/150‑1 (MQ=255)
                     aGCCGCCAGCAGGCAGCGCAGCTGGCAGCGGTGCTGCCAAGCCCAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTg                                                                                                                          <  1:1219924/150‑1 (MQ=255)
                       ccgccAGCAGGCAGCGCAGCTGGCAGCGGTGCTGCCAAGCCCAATCAAGTGGAGTGCCAGCCGACCGAGAGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGa                                                                                                                        <  1:961164/150‑1 (MQ=255)
                             gcaggcagCGCAGCTGGCAGCGGTGCTGCCAAGCCCAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGActcctc                                                                                                                  <  2:224497/150‑1 (MQ=255)
                                 gcagcgcagcTGGCAGCGGTGCTGCCAAGCCCAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGa                                                                                                              >  1:408479/1‑150 (MQ=255)
                                 gcagcgcagcTGGCAGCGGTGCTGCCAAGCCCAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGCGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGa                                                                                                              >  2:1066112/1‑150 (MQ=255)
                                      gcagcTGGCAGCGGTGCTGCCAAGCCCAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCt                                                                                                         >  2:880654/1‑150 (MQ=255)
                                                 cGGTGCTGCCAAGCCCAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAgg                                                                                              >  1:781035/1‑150 (MQ=255)
                                                        gcCAAGCCCAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTAcc                                                                                                                                     >  2:843138/1‑104 (MQ=255)
                                                        gcCAAGCCCAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTAcc                                                                                                                                     <  1:843138/104‑1 (MQ=255)
                                                         ccAAGCCCAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTgggcgggc                                                                                                                                              <  2:890090/94‑1 (MQ=255)
                                                         ccAAGCCCAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTgggcgggc                                                                                                                                              >  1:890090/1‑94 (MQ=255)
                                                         ccAAGCCCAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCAAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGCCTCCTCCAGAAAGCTTTGAATGCAggcccggccc                                                                                      >  2:411351/1‑150 (MQ=255)
                                                           aaGCCCAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCttt                                                                                                       >  2:1141314/1‑131 (MQ=255)
                                                           aaGCCCAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCttt                                                                                                       <  1:1141314/131‑1 (MQ=255)
                                                              cccAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAGGCCCGGCCCTTTcg                                                                                 >  2:709131/1‑150 (MQ=255)
                                                               ccAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAggcccggcctcttcgc                                                                                <  2:408479/150‑8 (MQ=255)
                                                                cAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCa                                                                                                                >  2:129183/1‑117 (MQ=255)
                                                                cAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCa                                                                                                                <  1:129183/117‑1 (MQ=255)
                                                                cAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAGGCCCGGCCCTTTcgcg                                                                               <  1:462483/150‑1 (MQ=255)
                                                                cAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTAGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAGGCCCGGCCCTTTcgcg                                                                               <  2:901270/150‑1 (MQ=255)
                                                                 aaTCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAGGCCCGGCCCTTTCGCgg                                                                              <  1:950209/150‑1 (MQ=255)
                                                                   tCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCt                                                                                                                           <  1:1052820/103‑1 (MQ=255)
                                                                   tCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCt                                                                                                                           >  2:1052820/1‑103 (MQ=255)
                                                                      aGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAGGCCCGGCCCTTTCGCgg                                                                              <  2:725513/145‑1 (MQ=255)
                                                                      aGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAGGCCCGGCCCTTTCGCgg                                                                              >  1:725513/1‑145 (MQ=255)
                                                                      aGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAGGCCCGGCCCTTTCGCGGGTaaa                                                                         <  1:880654/150‑1 (MQ=255)
                                                                       gTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAGGCCCGGCCCTTTCGCGGGTAAAc                                                                        <  1:709131/150‑1 (MQ=255)
                                                                          gAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAGGCCCGGCCCTTTCGCgg                                                                              <  1:925149/141‑1 (MQ=255)
                                                                          gAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAGGCCCGGCCCTTTCGCgg                                                                              >  2:925149/1‑141 (MQ=255)
                                                                           aGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGTCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAGGCCCGGCCCTTTCGCGGGTAAACCCGc                                                                    <  1:913299/150‑1 (MQ=255)
                                                                                  gCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAGGCCCGGCCCTTTCGCGGGTAAACCCGCTCCTACg                                                             >  1:555649/1‑150 (MQ=255)
                                                                                  gCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAGGCCCGGCCCTTTCGCGGGTAAACCCCCTCCTCCg                                                             >  2:417299/1‑150 (MQ=255)
                                                                                   ccgaccgaGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAGGCCCGGCCCTTTCGCGGGTAAACCCGCTCCTACgg                                                            <  2:907388/150‑1 (MQ=255)
                                                                                   ccgaccgaGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAGGCCCGGCCCTTTCGCGGGTAAACCCGCTCCTACgg                                                            <  1:316279/150‑1 (MQ=255)
                                                                                     gaccgaGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAGGCCCGGCCCTTTCGCGGGTAAACCCGCTCCTACGGGc                                                          >  2:956623/1‑150 (MQ=255)
                                                                                      accgaGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAGGCCCGGCCCTTTCGCGGGTAAACCCGCTCCTACGGGCa                                                         >  1:16896/1‑150 (MQ=255)
                                                                                         gaGTGCCTACTTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAggc                                                                                             <  1:200661/111‑1 (MQ=255)
                                                                                         gaGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAggc                                                                                             >  2:200661/1‑111 (MQ=255)
                                                                                               cTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAGGCCCGGCCCTTTCGCGGGTAAACCCGCTCCTACGGGCAACCGCATcc                                                <  1:1073014/150‑1 (MQ=255)
                                                                                                  cGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAGGCCCGGCCCTTTCGCGGGTAAACCCGCTCCTACGGGCAACCGCATCCCCg                                             >  2:620387/1‑150 (MQ=255)
                                                                                                     atccAGCCGGCGCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAGGCCCGGCCCTTTCGCGGGTAAACCCGCTCCTACGGGCAACCGCATCCCCGGcc                                          >  2:282890/3‑150 (MQ=255)
                                                                                                     atccAGCCGGCCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAGGCCCGGCCCTTTCGCGGGTAAACCCGCTCCTACGGGCAACCGCATCCCCGGcc                                          >  2:582847/3‑150 (MQ=255)
                                                                                                      tccAGCCGGCCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAGGCCCGGCCCTTTCGCGGGTAAACCCGCTCCTACGGGCAACCGCATCCCCGGcct                                         >  2:507523/2‑150 (MQ=255)
                                                                                                      tccAGCCGGCCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAGGCCCGGCCCTTTCGCGGGTAAACCCGCTCCTACGGGCAACCGCATCCCCGGcct                                         >  1:995969/2‑150 (MQ=255)
                                                                                                      gCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAGGCCCGGCCCTTTCGCGGGTAAACCCGCTCCTACGGGCAACCGCATCCCCGGcct                                         >  1:990909/1‑150 (MQ=255)
                                                                                                                            gCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAGGCCCGGCCCTTTCGCGGGTAAACCCGCTCCTACGGGCAACCGCATCCCCGGCCTCCTCGCGGGCAAACCCGCGCCt                   >  1:264391/1‑150 (MQ=255)
                                                                                                                                         ccAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAGGCCCGGCCCTTTCGCgg                                                                              <  2:357815/78‑1 (MQ=255)
                                                                                                                                         ccAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAGGCCCGGCCCTTTCGCGGGTAAACCCGCTCCTACGGGCAACCGCATCCCCGGCCTCCTCGCGGGCAAACCcgc                       >  1:1053410/1‑133 (MQ=255)
                                                                                                                                         ccAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAGGCCCGGCCCTTTCGCGGGTAAACCCGCTCCTACGGGCAACCGCATCCCCGGCCTCCTCGCGGGCAAACCcgc                       <  2:1053410/133‑1 (MQ=255)
                                                                                                                                         ccAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGATGTGAATGCAGGCCCGGCCCTTTCGCgg                                                                              >  1:357815/1‑78 (MQ=255)
                                                                                                                                          cAGCTGGGCGGGCCCAGCGACCTAAAGCAGCTCGACTCCTCCAGAAAGCTTTGAATGCAGGCCCGGCCCTTTCGCGCTTAAACACGCTCCTACGGGCAACCGCATCCCCGGCCACCTCGCGGGCAAACCCGCGCCTACAGGCAACCGCAt     <  1:411351/150‑1 (MQ=255)
                                                                                                                                           aGCTGGGCCCGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAGGCCCGGCCCTTTCGCGGGTAAACCCGCTCCTACGGGCAACCGCATCCCCGGCCTCCTCGCGGGCAAACCCGCGCCTAc                 >  2:768611/1‑137 (MQ=255)
                                                                                                                                           aGCTGGGCCCGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAGGCCCGGCCCTTTCGCGGGTAAACCCGCTCCTACGGGCAACCGCATCCCCGGCCTCCTCGCGGGCAAACCCGCGCCTAc                 <  1:768611/137‑1 (MQ=255)
                                                                                                                                             cTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAGGCCCGGCCCTTTCGCGGGTAAACCCGCTCCTACGGGCAACCGCATCCCCGGCCTCCTCGCGGGCAAACCCGCGCCTACAGGCAACCGCATccc  <  2:825579/150‑1 (MQ=255)
                                                                                                                                             cTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAGGCCCGGCCCTTTCGCGGGTAAACCCGCTCCTACGGGCAACCGCATCCCCGGCCTCCTCGCGGGCAAACCCGCGCCTACAGGCAACCGCATccc  <  2:995969/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                              
GGCGTCGACGCCAGCCGGCTCAGCCGCCAGCAGGCAGCGCAGCTGGCAGCGGTGCTGCCAAGCCCAATCAAGTGGAGTGCCAGCCGACCGAGTGCCTACGTGGCCAGCCGGGCCGGCTGGATTCGCCGGCAGATGAGCCAGCTGGGCGGGCCCAGCTACCTGATGCAGCTGGACTCCTCCAGAAAGCTTTGAATGCAGGCCCGGCCCTTTCGCGGGTAAACCCGCTCCTACGGGCAACCGCATCCCCGGCCTCCTCGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCC  >  NC_002947/5827663‑5827953

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