Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 838,599 G→T 6.3% D33E (GAC→GAA rplY ← 50S ribosomal protein L25/general stress protein Ctc

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947838,5990GT 6.3% 230.7 / 3.5 73D33E (GAC→GAArplY50S ribosomal protein L25/general stress protein Ctc
Reads supporting (aligned to +/- strand):  ref base G (31/37);  new base T (2/3);  total (33/40)
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.43e-01

AACGTGGCTGAAGGCAGCCTCGTTTTCGAACAGTTTGGTGATTTCCTTGGCCAGGATGGTCAGGGACTGGGCGTCCTGGTTACCACCGTAGACAACGGCTGGGATGTTCGCGGAGTGACGCAGGCGGCGGCTCGCACCTTTCCCCAGGTCAGTACGCGCTTGGGCGTTCAGAGTGAATTCAGTCATTTTGTTTCTCCAAAGTAGCCCCCCGAGGGCGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATCAACA  >  NC_002947/838452‑838745
                                                                                                                                                   |                                                                                                                                                  
aaCGTGGCTGAAGGCAGCCTCGTTTTCGAACAGTTTGGTGATTTCCTTGGCCAGGATGGTCAGGGACTGGGCGTCCTGGTTACCACCGTAGACAACGGCTGGTATGTTCGCGGAGTGACGCAGGCGGCGGCTCGCACCTTTCCCCAGGt                                                                                                                                                   <  1:81041/149‑1 (MQ=255)
   gTGGCTGAAGGCAGCCTCGTTTTCGAACAGTTTGGTGATTTCCTTGGCCAGGATGGTCAGGGACTGGGCGTCCTGGTTACCACCGTAGACAACGGCTGGGATGTTCGCGGAGTGACGCAGGCGGCGGCTCGCACCTTTCCCCAGGTCAg                                                                                                                                                >  6:320198/1‑149 (MQ=255)
      gCTGAAGGCAGCCTCGTTTTCGAACAGTTTGGTGATTTCCTTGGCCAGGATGGTCAGGGACTGGGCGTCCTGGTTACCACCGTAGACAACGGCTGGGATGTTCGCGGAGTGACGCAGGCGGCGGCTCGCACCTTTCCCCAGGTCAGTAc                                                                                                                                             <  2:106049/149‑1 (MQ=255)
        tGAAGGCAGCCTCGTTTTCGAACAGTTTGGTGATTTCCTTGGCCAGGATGGTCAGGGACTGGGCGTCCTGGTTACCACCGTAGACAACGGCTGGGATGTTCGCGGAGTGACGCAGGCGGCGGCTCGCACCTTTCCCCAGGTCAGTAcgc                                                                                                                                           >  8:38717/1‑149 (MQ=255)
            ggCAGCCTCGTTTTCGAACAGTTTGGTGATTTCCTTGGCCAGGATGGTCAGGGACTGGGCGTCCTGGTTACCACCGTAGACAACGGCTGGGATGTTCGCGGAGTGACGCAGGCGGCGGCTCGCACCTTTCCCCAGGTCAGTACGCGCtt                                                                                                                                       <  4:349007/149‑1 (MQ=255)
               aGCCTCGTTTTCGAACAGTTTGGTGATTTCCTTGGCCAGGATGGTCAGGGACTGGGCGTCCTGGTTACCACCGTAGACAACGGCTGGGATGTTCGCGGAGTGACGCAGGCGGCGGCTCGCACCTTTCCCCAGGTCAGTACGCGCTTggg                                                                                                                                    <  8:251411/149‑1 (MQ=255)
               aGCCTCGTTTTCGAACAGTTTGGTGATTTCCTTGGCCAGGATGGTCAGGGACTGGGCGTCCTGGTTACCACCGTAGACAACGGCTGGGATGTTCGCGGAGTGACGCAGGCGGCGGCTCGCACCTTACCCCAGGTCAGTACGCGCTTggg                                                                                                                                    >  8:88286/1‑149 (MQ=255)
                gCCTCGTTTTCGAACAGTTTGGTGATTTCCTTGGCCAGGATGGTCAGGGACTGGGCGTCCTGGTTACCACCGTAGACAACGGCTGGGATGTTCGCGGAGTGACGCAGGCGGCGGCTCGCACCTTTCCCCATGTCAGTACGCGCTTGGGc                                                                                                                                   <  7:206176/149‑1 (MQ=255)
                 ccTCGTTTTCGAACAGTTTGGTGATTTCCTTGGCCAGGATGGTCAGGGACTGGGCGTCCTGGTTACCACCGTAGACAACGGCTGGGATGTTCGCGGAGTGACGCAGGCGGCGGCTCGCACCTTTCCCCAGGTCAGTACGCGCTTGGGCg                                                                                                                                  >  1:223349/1‑149 (MQ=255)
                    cGTTTTCGAACAGTTTGGTGATTTCCTTGGCCAGGATGGTCAGGGACTGGGCGTCCTGGTTACCACCGTAGACAACGGCTGGGATGTTCGCGGAGTGACGCAGGCGGCGGCTCGCACCTTTCCCCAGGTCAGTACGCGCTTgg                                                                                                                                     <  4:261408/143‑1 (MQ=255)
                    cGTTTTCGAACAGTTTGGTGATTTCCTTGGCCAGGATGGTCAGGGACTGGGCGTCCTGGTTACCACCGTAGACAACGGCTGGGATGTTCGCGGAGTGACGCAGGCGGCGGCTCGCACCTTTCCCCAGGTCAGTACGCGCTTgg                                                                                                                                     >  3:261408/1‑143 (MQ=255)
                     gTTTTCGAACAGTTTGGTGATTTCCTTGGCCAGGATGGTAAGGGACTGGGCGTCCTGGTTACCACCGTAGACAACGGCTGGGATGTTCGCGGAGTGACGCAGGAGGAGGCTCGCACCTTTCCCCAGGTCAGTACGCGCTTGGGCGTTc                                                                                                                               >  6:230102/1‑148 (MQ=255)
                       tttCGAACAGTTTGGTGATTTCCTTGGCCAGGATGGTCAGGGACTGGGCGTCCTGGTTACCACCGTAGACAACGGCTGGGATGTTCGCGGAGTGACGCAGGCGGCGGCTCGCACCTTTCCCCAGGTCAGTACGCGCTTGGGCGTTCa                                                                                                                              >  7:87370/1‑147 (MQ=255)
                       ctacGAACAGTTTGGTGATTTCCTTGGCCAGGATGGTCAGGGACTGGGCGTCCTGGTTACCACCGTAGACAACGGCTGGGATGTTCGCGGAGTGACGCAGGCGGCGGCTCGCACCTTTCCCCAGGTCAGTACGCGCTTGGGCGTTCag                                                                                                                             <  1:132056/145‑1 (MQ=255)
                            aaCAGTTTGGTGATTTCCTTGGCCAGGATGGTCAGGGACTGGGCGTCCTGGTTACCACCGTAGACAACGGCTGGGATGTTCGCGGAGTGACGCAGGCGGCGGCTCGCACCTTTCCCCAGGTCAGTACGCGCTTGGGCGTTCAGAGTGaa                                                                                                                       >  3:303207/1‑149 (MQ=255)
                                 tttGGTGATTTCCTTGGCCAGGATGGTCAGGGACTGGGCGTCCTGGTTACCACCGTAGACAACGGCTGGGATGTACGCGGAGTGACGCAGGCGGCGGCTCGCCCCGTTCCCCAGGACAGTCCGCGCGTGGGCGTTCAGAGGGAATACAg                                                                                                                  >  2:77998/1‑149 (MQ=255)
                                     gTGATTTCCTTGGCCAGGATGGTCAGGGACTGGGCGTCCTGGTTACCACCGTAGACAACGGCTGGGATGTTCGCGGAGTGACGCAGGCGGCGGCTCGCACCTTTCCCCAGGTCAGTACGCGCTTGGGCGTTCAGAGTGAATTCAGTCa                                                                                                               <  4:303207/148‑1 (MQ=255)
                                       gATTTCCTTGGCCAGGATGGTCAGGGACTGGGCGTCCTGGTTACCACCGTAGACAACGGCTGGGATGTTCGCGGAGTGACGCAGGCGGCGGCTCGCACCTTTCCCCAGGTCAGTACGCGCTTGGGCGTTCAGAGTGAATTCAGTCAttt                                                                                                            >  3:3780/1‑149 (MQ=255)
                                       gATTTCCTTGGCCAGGATGGTCAGGGACTGGGCGTCCTGGTTACCACCGTAGACAACGGCTGGGATGTTCGCGGAGTGACGCAGGCGGCGGCTCGCACCTTTCCCCAGGTCAGTACGCGCTTGGGCGTTCAGAGTGAATTCAGTCAttt                                                                                                            >  3:281466/1‑149 (MQ=255)
                                          ttCCTTGGCCAGGATGGTCAGGGACTGGGCGTCCTGGTTACCACCGTAGACAACGGCTGGGATGTTCGCGGAGTGACGCAGGCGGCGGCTCGCACCTTTCCCCAGGTCAGTACGCGCTTGGGCGTTCAGAGTGAATTCAGTCATTTTGt                                                                                                         >  5:361953/1‑149 (MQ=255)
                                             cTTGGCCAGGATGGTCAGGGACTGGGCGTCCTGGTTACCACCGTAGACAACGGCTGGGATGTTCGCGGAGTGACGCAGGCGGCGGCTCGCACCTTTCCCCAGGTCAGTACGCGCTTGGGCGTTCAGAGTGAATTCAGTCATTTTGTTtc                                                                                                      <  6:203432/149‑1 (MQ=255)
                                                 gCCAGGATGGTCAGGGACTGGGCGTCCTGGTTACCACCGTAGACAACGGCTGGGATGTTCGCGGAGTGACGCAGGCGGCGGCTCGCACCTTTCACCAGGTCAGTACGCGCTTGTGCGTTCAGAGTGAATTACGTCATTTTGTTTCTcc                                                                                                   >  7:351976/1‑148 (MQ=255)
                                                   cAGGATGGTCAGGGACTGGGCGTCCTGGTTAGCACCGTAGACAACGGCTGGGATGTTCGCGGAGTGACGCAGTCGGCGGCTCGGACCTTTCCCCAGGTCAGTACGCGCTTGGGCGTTCAGAGTGAATTCAGTCATTTTGTTTCTCCaaa                                                                                                <  5:230102/149‑1 (MQ=255)
                                                   cAGGATGGTCAGGGACTGGGCGTCCTGGTTACCACAGTAGACAACGGCTGGGATGTTCGCGGAGTGACGCAGGCGGCGGCTCGCACCTTTCCCCAGGTCAGTACGCGCTTGGGCGTTCAGAGTGAATTCAGTCATTTTGTTTCTCCa                                                                                                  >  2:82748/1‑147 (MQ=255)
                                                   cAGGATGGGCAGGGACTGGGCGTCCTGGTTACCAGCGTAGACAACGGCTGGGAGGTTCGCGGAGTGACGCAGGCGGCGGCTCGCACCTTTCCCCAGGTCAGTACGCGCTTGGGCGTTCAGAGTGAATTCAGTCATTTTGTTTCTCCaaa                                                                                                <  1:77998/149‑1 (MQ=255)
                                                      gATGGTCAGGGACTGGGCGTCCTGGTTACCACCGTAGACAACGGCTGGGATGTTCGCGGAGTGACGCAGGCGGCGGCTCGCACCTTTCCCCAGGTCAGTACGCGCTTGGGCGTTCAGAGTGAATTCAGTCATTTTGTTTCTCCAAAGTa                                                                                             <  4:281466/149‑1 (MQ=255)
                                                      gATGGTCAGGGACTGGGCGTCCTGGTTACCACCGTAGACAACGGCTGGGATGTTCGCGGAGTGACGCAGGCGGCGGCGCGCACCTTTCCCCAGGTCAGTACGCGCTTGGGCGTTCAGAGTGAATTCAGTCATTTTGTTTCTCCAAAg                                                                                               >  5:376180/1‑147 (MQ=255)
                                                        tGGTCAGGGACTGGGCGTCCTGGTTACCACCGTAGACAACGGCTGGGATGTTCGCGGAGTGACGCAGGCGGCGGCTCGCACCTTTCCCCAGGTCAGTACGCGCTTGGGCGTTCAGAGTGAATTCAGTCATTTTGTTTCTCCAAAGTAGc                                                                                           <  4:359423/149‑1 (MQ=255)
                                                          gTCAGGGACTGGGCGTCCTGGTTACCACCGTAGACAACGGCTGGGATGTTCGCGGAGTGACGCAGGCGGCGGCTCGCACCTTTCCCCAGGTCAGTACGCGCTTGGGCGTTCAGAGTGAATTCAGTCATTTTGTTTCTCCAAAGTAGccc                                                                                         <  3:90415/149‑1 (MQ=255)
                                                          gTCAGGGACTGGGCGTCCTGGTTACCACCGTAGACAACGGCTGGGATGTTCGCGGAGTGACGCAGGCGGCGGCTCGCACCTTTCCCCAGGTCAGTACGCGCTTGGGCGTTCAGAGTGAATTCAGTCATTTTGTTTCTCCAAAGTAGcc                                                                                          <  8:87370/148‑1 (MQ=255)
                                                            cAGGGACTGGGCGTCCTGGTTACCACCGTAGACAACGGCTGGGATGTTCGCGGAGTGACGCAGGCGGCGGCTCGCACCTTTCCCCAGGTCAGTACGCGCTTGGGCGTTCAGAGTGAATTCAGTCATTTTGTTTCTCCAAAGTAGccccc                                                                                       >  2:90646/1‑149 (MQ=255)
                                                                 aCTGGGCGTCCTGGTTACCACCGTAGACAACGGCTGGGATGTTCGCGGAGTGACGCAGGCGGCGGCTCGCACCTTTCCCTAGGTCAGTACGCGCTTGGGCGTTCAGAGTGAATTCAGTCATTTTGTTTCTCCAAAg                                                                                               <  8:352659/136‑1 (MQ=255)
                                                                 aCTGGGCGTCCTGGTTACCACCGTAGACAACGGCTGGGATGTTCGCGGAGTGACGCAGGCGGCGGCTCGCACCTTTCCCGAGGTCAGTACGCGCTTGGGCGTTCAGAGTGAATTCAGTCATTTTGTTTCTCCAAAg                                                                                               >  7:352659/1‑136 (MQ=255)
                                                                 aCTGGGCGTCCTGGTTACCACCGTAGACAACGGCTGGGATGTTCGCGGAGTGACGCAGGCGGCGGCTCGCACCTTTCCCCAGGTCAGTACGCGCTTGGGCGTTCAGAGTGAATTCAGTCATTTTGTTTCTCCAAAGTAGCCCCCCGAgg                                                                                  <  7:38717/149‑1 (MQ=255)
                                                                 aCTGGGCGTCCTGGTTACCACCGTAGACAACGGCTGGGATGTTCGCGGAGTGACGCAGGCGGCGGCTCGCACCTTTCCCCAGGTCAGTACGCGCTTGGGCGTTCAGAGTGAATTCAGTCATTTTGTTTCTCCAAAGTAGCCCCCCGAgg                                                                                  <  4:3780/149‑1 (MQ=255)
                                                                 aCTGGGCGTCCTGGTTACCACCGTAGACAACGGCTGGGATGTTCGCGGAGTGACGCAGGCGGCGGCTCGCACCTTTCCCCAGGTCAGTACGCGCTTGGGCGTTCAGAGTGAATTCAGTCATTTTGTTTCTCCAAAGTAGCCCCCCGAgg                                                                                  <  1:249774/149‑1 (MQ=255)
                                                                         tCCTGGTTACCACCGTAGACAACGGCTGGGATGTTCGCGGAGTGACGCAGGCGGCGGCTCGCACCTTTCCCCAGGTCAGTACGCGCTTGGGCGTTCAGAGTGAATTCAGTCATTTTGTTTCTCCAAAGTAGCCCCCCGAGGGCGTTTGc                                                                          <  2:210342/149‑1 (MQ=255)
                                                                                 accaccGTAGACAACGGCTGGGATGTTCGCGGAGTGACGCAGGCGGCGGCTCGCACCTTTCCCCAGGTCAGTACGCGCTTGGGCGTTCAGAGTGAATTCAGTCATTTTGTTTCTCCAAAGTAGCCCCCCGAGGGCGTTTGCGACCAgc                                                                   <  3:66234/148‑1 (MQ=255)
                                                                                   caccGTAGACAACGGCTGGGATGTTCGCGGAGTGACGCAGGCGGCGGCTCGCACCTTTCCCCAGGTCAGTACGCGCTTGGGCGTTCAGAGTGAATTCAGTCATTTTGTTTCTCCAAAGTAGCCCCCCGAGGGCGTTTGCGACCAgcgc                                                                 >  6:145475/1‑148 (MQ=255)
                                                                                        tAGACAACGGCTGGGATGTTCGCGGAGTGACGCAGGCGGCGGAGCGCACCTTTCCCCAGGTCAGTACGCGCTTGGGCGTTCAGAGTGAATTCAGTCATTTTGTTTCTCCAAAGTAGCCCCCCGAGGGCGTTTTCGACCAGCGCCAACgg                                                           <  6:376180/149‑1 (MQ=255)
                                                                                          gACAACGGCTGTGATGTTCGCGGAGTGACTCAGGCGGCGGCTAGCACCTTTACCCAGTTCATTACGCGCTTGGGCGTTCAGAGTGAATTCAGTAATTTTGTTTATCAAAAGTAGCCCCCCGAGGGCGTTTGCGACAAGCGCCAACggg                                                          >  2:106734/1‑148 (MQ=255)
                                                                                          gACAACGGCTGGGATGTTCGCGGAGTGACGCAGGCGGCGGCTCGCACCTTTCCCCAGTTCAGTACGCGCTTGGGCGTTCAGAGTGAATTCAGTCATTTTGTTTCTCCAAAGTAGCCCCCCGAGGGCGTTTGCGACCAGCGCCAACggg                                                          >  2:111366/1‑148 (MQ=255)
                                                                                           aCAACGGCTGGGATGTTCGCGGAGTGACGCAGGCGGCGGCTCGCACCTTTCCCCAGGTCAGTACGCGCTTGGGCGTTCAGAGTGAATTCAGTCATTTTGTTTCTCCAAAGTAGCCCCCCGAGGGCGTTTGCGCCCAGCGCCAACGGGGa                                                        >  6:87228/1‑149 (MQ=255)
                                                                                            ccaCGGCTGGGATGTTCGCGGAGTGTCGCAGGGGGCGGCTCGCACCTTTCCCCAGTTCAGTACGCGCTTGGGCGTTCAGAGTGAATTCAGTCATTTTGTTTCTCCATAGTAGCCCCCCGAGGGCGTTTGCGACCAGCGCCATCGGGGAt                                                       <  4:32964/147‑1 (MQ=255)
                                                                                            cAACGGCTGGGATGTTCGCGGAGTGACGCAGGCGGCGGCTCGCACCTTTCCCCAGGTCAGTACGCGCTTGGGCGTTCAGAGTGAATTCAGTCATTTTGTTTCTCCAAAGTAGCCCCCCGAGGGCGTTTGCGACCAGCGCCAACGGGGAt                                                       <  5:320198/149‑1 (MQ=255)
                                                                                                     ggATGTTCGCGGAGTGACGCAGGCGGCGGCTCGCACCTTTCCCCAGGTCAGTACGCGCTTGGGCGTTCAGAGTGAATTCAGTCATTTTGTTTCTCCAAAGTAGCCCCCCGAGGGCGTTTGCGACCAGCGCCAACGGGGATGGCaaaaaa                                              <  1:90646/149‑1 (MQ=255)
                                                                                                             gcgGAGTGACGCAGGCGGCGGCTCGCACCTTTCCCCAGGTCAGTACGCGCTTGGGCGTTCAGAGTGAATTCAGTCATTTTGTTTCTCCAAAGTAGCCCCCCGAGGGTGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAgccccgcc                                      >  5:23810/1‑149 (MQ=255)
                                                                                                               ggAGTGACGCATGCGGCGGCTCGCACATTTCCCCATGTCAGTACGAGCTTTGGAGTTCATAGTGAATTCAGTCATTTTGTTTCTCCAAAGTAGCCCCCCGAGGGCGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAgccccgc                                       <  8:351976/146‑1 (MQ=255)
                                                                                                               ggAGTGACGCAGGCGGCGGCTCGCACCTTTCCCCAGGTCAGTACGCGCTTGGGCGTTCAGAGTGAATTCAGTCATTTTGTTTCTCCAAAGTAGCCTCCCGAGGGCGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAgccccgccc                                     >  2:149789/1‑148 (MQ=255)
                                                                                                                     aCGCAGGCGGCGGCTCGCACCTTTCCCCAGGTCAGTACGCGCTTGGGCGTTCAGAGTGAATTCAGTCATTTTGTTTCTCCAAAGTAGCCCCCCGAGGGCGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAaca                                >  5:121285/1‑147 (MQ=255)
                                                                                                                      cGCAGGCGGCGGCTCGCACCTTTCCCCAGGTCAGTACGCGCTTGGGCGTTCAGAGTGAATTCAGTCATTTTGTTTCTCCAAAGTAGCCCCCCGAGGGCGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAAcaca                              <  2:223349/148‑1 (MQ=255)
                                                                                                                        cAGGCTGCGGCTCGCACCTTTCCCCAGTTCAGTACGCTCTTGTGCGTTCAGAGTGAATTCAGTCATTTTGTTTCTCCAAAGTAGCTTCCCGCGGGCGTTTGCGAACAGCGCCAACGGGGATGGCAAAAAAGCCCAGCCCCAACACAGGt                           <  1:106734/149‑1 (MQ=255)
                                                                                                                        cAGGCGGCGGCTCGCACCTTTCCCCAGTTCAGTACGCGCTTGGGCGTTCAGAGTGAATTCAGTCATTTTGTTTCTCCAAAGTAGCTTCCCGAGGGCGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACACAGGt                           <  1:111366/149‑1 (MQ=255)
                                                                                                                         aGGCGGCGGCTCGCACCTTTCCCCAGGTCAGTACGCGCTTGGGCGTTCAGAGTGAATTCAGTCATTTTGTTTCTCCAAAGTAGCCCCCCGAGGGCGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACACAGGtt                          >  1:313175/1‑149 (MQ=255)
                                                                                                                           gcggcggcTCGCACCTTTCCCCAGGTCAGTACGCGCTTGGGCGTTCAGAGTGAATTCAGTCATTTTGTTTCTCCAAAGTAGCCCCCCGAGGGCGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACACAGGTTgg                        >  2:199805/1‑149 (MQ=255)
                                                                                                                             ggcggcTCGCACCTTTCCCCAGGTCAGTACGCGCTTGGGCGTTCAGAGTGAATTCAGTCATTTTGTTTCTCCAAAGTAGCCCCCCGAGGGCGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACACAGGTTgggg                      >  7:190552/1‑149 (MQ=255)
                                                                                                                                ggcTCGCACCTTTCCCCAGGTCAGTACGCGCTTGGGCGTTCAGAGTGAATTCAGTCATTTTGTTTCTCCAAAGTAGCCTCCCGAGGGCGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACACAGGTTggggcgg                   <  1:149789/149‑1 (MQ=255)
                                                                                                                                 gcTCGCACCTTTCCCCCTGTCAGTACGCGCTTGGGCGTTCATATTGAATTCATTCTTTTTGTTTCTCCAAAGTATCCCCCCTATGGTGTTTGCGACCAGCGCCAACGGGGTTGTCAAAAAACCCCCTCCCCAACCTGTGTTggggcggg                  <  7:235373/149‑1 (MQ=255)
                                                                                                                                 gcTCGCACCTTTCCCCAGGTCAGTACGCGCTTGGGCGTTCAGAGTGAATTCAGTCATTTTGTTTCTCCAAAGTAGCCCCCCGAGGGTGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACACAGGTTggggcggg                  <  6:23810/149‑1 (MQ=255)
                                                                                                                                   tCGCACCTTTCCCCAGGTCAGTACGCGCTTGGGCGTTCAGAGTGAATTCAGTCATTTTGTTTCTCAAAAGTAGCCCCCCGAGGGCGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACACAGGTTggggcggggc                >  4:120212/1‑149 (MQ=255)
                                                                                                                                       aCCTTTCCCCAGGTCAGTACGCGCTTGGGCGTTCAGAGTGAATTCAGTCATTTTGTTTCTCCAAAGTAGCCCCCCGAGGGCGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCtt            <  5:66473/149‑1 (MQ=255)
                                                                                                                                       aCCTTTCCCCAGGTCAGTACGCGCTTGGGCGTTCAGAGTGAATTCAGTCATTTTGTTTCTCCAAAGTAGCCCCCCGAGGGCGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCtt            >  8:189625/1‑149 (MQ=255)
                                                                                                                                         cTTTCCCCAGGTCAGTACGCGCTTGGGCGTTCAGAGTGAATTCAGTCATTTTGTTTCTCCAAAGTAGCCCCCCGAGGGCGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTca          >  6:319771/1‑149 (MQ=255)
                                                                                                                                           ttCCCCAGGTCAGTACGCGCTTGGGCGTTCAGAGTGAATTCAGTCATTTTGTTTCTCCAAAGTAGCCCCCCGAGGGCGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTcaca        >  4:298033/1‑149 (MQ=255)
                                                                                                                                            tCCCCAGGTCAGTACGCGCTTGGGCGTTCAGAGTGAATTCAGTCATTTTGTTTCTCCAAAGTAGCCCCCCGAGGGCGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTcac         <  6:121285/147‑1 (MQ=255)
                                                                                                                                              cccAGGTCAGTACGCGCTTGGGCGTTCAGAGTGAATTCAGTCATTTTGTTTCTCCAAAGTAGCCCCCCGAGGGCGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATc      <  5:145475/148‑1 (MQ=255)
                                                                                                                                               ccAGGTCAGTACGCGCTTGGGCGTTCAGAGTGAATTCAGTCATTTTGTTTCTCCAAAGTAGGCCCCCGAGGGCGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACACAGGtt                          <  2:128076/127‑1 (MQ=255)
                                                                                                                                                cAGGTCAGTACGCGCTTGGGCGTTCAGAGTGAATTCAGTCATTTTGTTTCTCCAAAGTAGCCCCCCGAGGGCGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATCAac   <  1:199805/149‑1 (MQ=255)
                                                                                                                                                 aGGTCAGTACGCGCTTGGGCGTTCAGAGTGAATTCAGTCATTTTGTTtctc                                                                                                    >  7:310253/1‑51 (MQ=255)
                                                                                                                                                 aGGTCAGTACGCGCTTGGGCGTTCAGAGTGAATTCAGTCATTTTGTTtctc                                                                                                    <  8:310253/51‑1 (MQ=255)
                                                                                                                                                 aGGTCAGTACGCGCTTGGGCGTTCAGAGTGAATTCAGTCATTTTGTTTCTCCAAAGTAGCCCCCCGAGGGCGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCaa                                  >  7:320222/1‑117 (MQ=255)
                                                                                                                                                 aGGTCAGTACGCGCTTGGGCGTTCAGAGTGAATTCAGTCATTTTGTTTCTCCAAAGTAGCCCCCCGAGGGCGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCaa                                  <  8:320222/117‑1 (MQ=255)
                                                                                                                                                 aGGTCAGTACGCGCTTGGGCGTTCAGAGTGAATTCAGTCATTTTGTTTCTCCAAAGTAGCCCCCCGAGGGCGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATCAaca  <  5:87228/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
AACGTGGCTGAAGGCAGCCTCGTTTTCGAACAGTTTGGTGATTTCCTTGGCCAGGATGGTCAGGGACTGGGCGTCCTGGTTACCACCGTAGACAACGGCTGGGATGTTCGCGGAGTGACGCAGGCGGCGGCTCGCACCTTTCCCCAGGTCAGTACGCGCTTGGGCGTTCAGAGTGAATTCAGTCATTTTGTTTCTCCAAAGTAGCCCCCCGAGGGCGTTTGCGACCAGCGCCAACGGGGATGGCAAAAAAGCCCCGCCCCAACACAGGTTGGGGCGGGGCGCTTCACATCAACA  >  NC_002947/838452‑838745

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