Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 401,802 T→G 5.9% intergenic (+329/+100) PP_0334 → / ← PP_0335 transposase/glutathione S‑transferase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947401,8020TG 5.9% 210.8 / 3.6 68intergenic (+329/+100)PP_0334/PP_0335transposase/glutathione S‑transferase family protein
Reads supporting (aligned to +/- strand):  ref base T (33/31);  new base G (2/2);  total (35/33)
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.55e-01

CCTGATACATAGATGCAACCGGGTTCCTATGTTCAAAAAAGCAGCTAGACAGTGTGGGCGAGTCCATACATAGACAGTGTGGGCGAGTCCATACATAGGAGCGGGCTTGCCCGCGAACACACATAATGGGGTGTCAGGCA‑‑CCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGA  >  NC_002947/401657‑401942
                                                                                                                                                   |                                                                                                                                             
ccTGATACATAGATGCAACCGGGTTCCTATGTTCAAAAAAGCAGCTAGACAGTGTGGGCGAGTCCATACATAGACAGTGTGGGCGAGTCCATACATAGGAGCGGGCTTGCCCGCGAACACACATAATGGGGTGTCAGGCA‑‑CCGGCTTCg                                                                                                                                            <  8:107843/149‑1 (MQ=255)
    aTACATAGATGCAACCGGGTTCCTATGTTCAAAAAAGCAGCTAGACAGTGTGGGCGAGTCCATACATAGACAGTGTGGGCGAGTCCATACATAGGAGCGGGCTTGCCCGCGAACACACATAATGGGGTGTCAGGCA‑‑CCGGCTTCGCCg                                                                                                                                         <  3:165763/148‑1 (MQ=255)
        aTAGATGCAACCGGGTTCCTATGTTCAAAAAAGCAGCTAGACAGTGTGGGCGAGTCCATACATAGACAGTGTGGGCGAGTCCATACATAGGAGCGGGCTTGCCCGCGAAAACACATAATGGGGTGTCAGGCA‑‑CCGGCTTCGCCGgtgt                                                                                                                                     >  6:246385/1‑148 (MQ=255)
         tAGATGCAACCGGGTTCCTATGTTCAAAAAAGCAGCTAGACAGTGTGGGGGAGTCCATACATAGACAGTGTGGGCGAGTCCATACATAGGAGCGGGCTTGCCCGCGAACACACATAATGGGGTGTCAGGCA‑‑CCGGCTTCGCCGGTGTTc                                                                                                                                   >  1:72854/1‑149 (MQ=255)
                aaCCGGGTTCCTATGTTCAAAAAAGCAGCTAGACAGTGTGGGCGAGTCCATACATAGACAGTGTGGGCGAGTCCATACATAGGAGCGGGCTTGCCCGCGAACACACATAATGGGGTGTCAGGCA‑‑CCGGCTTCGCCGGTGTTCGCGGGCa                                                                                                                            >  5:230704/1‑149 (MQ=255)
                aaCCGGGTTCCTATGTTCAAAAAAGCAGCTAGACAGTGTGGGCGAGTCCATACATAGACAGTGTGGGCGAGTCCATACATAGGAGCGGGCTTGCCCGCGAACACACATAATGGGGTGTCAGACA‑‑CCGGCGAAGCCGgtgcct                                                                                                                                   >  3:285946/1‑139 (MQ=38)
                aaCCGGGTTCCTATGTTCAAAAAAGCAGCTAGACAGTGTGGGCGAGTCCATACATAGACAGTGTGGGCGAGTCCATACATAGGAGCGGGCTTGCCCGCGAACACACATAATGGGGTGTCAGACA‑‑CCGGCGAAGCCGgtgcct                                                                                                                                   <  4:285946/142‑4 (MQ=38)
                 aCCGGGTTCCTATGTTCAAAAAAGCAGCTAGACAGTGTGGGCGAGTCCATACATAGACAGTGTGGGCGAGTCCATACATAGGAGCGGGCTTGCCCGCGAACACACATAATGGGGTGTCAGGCA‑‑CCGGCGAAGCCGgt                                                                                                                                       <  6:148112/137‑1 (MQ=255)
                 aCCGGGTTCCTATGTTCAAAAAAGCAGCTAGACAGTGTGGGCGAGTCCATACATAGACAGTGTGGGCGAGTCCATACATAGGAGCGGGCTTGCCCGCGAACACACATAATGGGGTGTCAGGCA‑‑CCGGCGAAGCCGgt                                                                                                                                       >  5:148112/1‑137 (MQ=255)
                   cGGGTTCCTATGTTCAAAAAAGCAGCTAGACAGTGTGGGCGAGTCCATACATAGACAGTGTGGGCGAGTCCATACATAGGAGCGGGCTTGCCCGCGAACACACATAATGGGGTGTCAGGCA‑‑CCGGCTTCGCCGGTGTTCGCGGGCACGc                                                                                                                         <  6:178741/149‑1 (MQ=255)
                          cTATGTTCAAAAAAGCAGCTAGACAGTGTGGGCGAGTCCATACATAGACAGTGTGGGCGAGTCCATACATAGGAGCGGGCTTGCCCGCGAACACACATAATGGGGTGTCAGGCA‑‑CCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTcc                                                                                                                  >  1:24756/1‑149 (MQ=255)
                           tATGTTCAAAAAAGCAGCTAGACAGTGTGGGCGAGTCCATACATAGACAGTGTGGGCGAGTCCATACATAGGAGCGGGCTTGCCCGCGAACACACATAATGGGGTGTCAGGCA‑‑CCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTcc                                                                                                                  >  2:12991/1‑148 (MQ=255)
                               ttCAAAAAAGCAGCTAGACAGTGTGGGCGAGTCCATACATAGACAGTGTGGGCGAGTCCATACATAGGAGCGGGCTTGCCCGCGAACACACATAATGGGGTGTCAGGCA‑‑CCGGCTTCGCCGGTGTTCGCGTGCACGCCCGCTACCaaa                                                                                                              >  8:209597/1‑148 (MQ=255)
                                  aaaaaaGCAGCTAGACAGTGTGGGCGAGTCCATACATAGACAGTGTGGGCGAGTCCATACATAGGAGCGGGCTTGCCCGCGAACACACATAATGGGGTGTCAGGCA‑‑CCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGcc                                                                                                          >  4:18049/1‑149 (MQ=255)
                                  aaaaaaGCAGCTAGACAGTGTGGGCGAGTCCATACATAGACAGTGTGGGCGAGTCCATACATAGGAGCGGGCTTGCCCGCGAACACACATAATGGGGTGTCAGGCA‑‑CCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGcc                                                                                                          <  1:51652/149‑1 (MQ=255)
                                    aaaaGCAGCTAGACAGTGTGGGCGAGTCAATACATAGACAGTGTGGGCGAGTCCATACATAGGAGCGGGCTTGCCCGCGAACACACATAATGGGGTGTCAGGCA‑‑CCGGCTTCGCCGGTGTTAGCGGGCACGCCCGCTCCCAAAAGCcca                                                                                                        >  1:154633/1‑149 (MQ=255)
                                     aaaGCAGCTAGACAGTGTGGGCGAGTCCATACATAGACAGTGTGGGCGAGTCCATACATAGGAGCGGGCTTGCCCGCGAACACACATAATGGGGTGTCAGGCA‑‑CCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCccac                                                                                                       <  2:24756/149‑1 (MQ=255)
                                             tAGACAGTGTGGGCGAGTCCATACATAGACAGTGTGGGCGAGTCCATACATAGGAGCGGGCTTGCCCGCGAACACACATAATGGGGTGTCAGGCA‑‑CCGGCTTCGCCGgtg                                                                                                                                      >  7:181704/1‑110 (MQ=255)
                                             tAGACAGTGTGGGCGAGTCCATACATAGACAGTGTGGGCGAGTCCATACATAGGAGCGGGCTTGCCCGCGAACACACATAATGGGGTGTCAGGCA‑‑CCGGCTTCGCCGgtg                                                                                                                                      <  8:181704/110‑1 (MQ=255)
                                             tAGACAGTGTGGGCGAGTCCATACATAGACAGTGTGGGCGAGTCCATACATAGGAGCGGGCTTGCCCGCGAACACACATAATGGGGTGTCAGGCA‑‑CCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGa                                                                                                 >  1:250465/1‑147 (MQ=255)
                                                aCAGTGTGGGCGAGTCCATACATAGACAGTGTGGGCGAGTCCATACATAGGAGCGGGCTTGCCCGCGAACACACATAATGGGGTGTCAGGCA‑‑CCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACcaca                                                                                                   <  8:228143/142‑1 (MQ=255)
                                                aCAGTGTGGGCGAGTCCATACATAGACAGTGTGGGCGAGTCCATACATAGGAGCGGGCTTGCCCGCGAACACACATAATGGGGTGTCAGGCA‑‑CCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGagatag                                                                                            >  7:228143/1‑144 (MQ=255)
                                                        ggCGAGTCCATACATAGACAGTGTGGGCGAGTCCATACATAGGAGCGGGCTTGCCCGCGAACACACATAATGGGGTGTCAGGCA‑‑CCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACtgt                                                                                    <  6:101418/149‑1 (MQ=255)
                                                             gTCCATACATAGACAGTGTGGGCGAGTCCATACATAGGAGCGGGCTTGCCCGCGAACACACATAATGGGGTGTCAGGCA‑‑CCGGCTTCGCCg                                                                                                                                         <  5:31196/91‑1 (MQ=255)
                                                             gTCCATACATAGACAGTGTGGGCGAGTCCATACATAGGAGCGGGCTTGCCCGCGAACACACATAATGGGGTGTCAGGCA‑‑CCGGCTTCGCCg                                                                                                                                         >  6:31196/1‑91 (MQ=255)
                                                                    cttaGACAGTGTGGGCGTGTCGATACATAGGTGCGGGCTTGCCCTCGAACAGAGATAATGGGGTGTCAGGCA‑‑GCGGCTTCGCCGGTGTTCGCTGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCTTACACTGTGGAAAACCTTg                                                                         <  2:233908/146‑1 (MQ=38)
                                                                    cataGACAGTGTGGGCGAGTCCATACATAGGAGCGGGCTTGCCCGCGAACACACATAATGGGGTGTCAGGCA‑‑CCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTg                                                                         <  2:250465/148‑1 (MQ=255)
                                                                    cataGACAGTGTGGGCGAGTCCATACATAGGAGCGGGCTTGCCCGCGAACACACATAATGGGGTGTCAGGCA‑‑CCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGt                                                                        <  6:230704/149‑1 (MQ=255)
                                                                     ataGACAGTGTGGGCGAGTCCAAACATAGGACCGCGCTTCCCCCCGAACACACATAATGGCGTGTCAGGCA‑‑CCGGCTTCGCCGGTGTTCGCGCGCACGCCCGCTCCCAAAAGCCCACCGCAGAAAGCGAACACTGTGGAAAACCtag                                                                         >  7:179632/1‑145 (MQ=38)
                                                                      taGACAGTGTGGGCGAGTCCATACATAGGAGCGGGCTTGCCCGCGAACACACATAATGGGGTGTCAGGCA‑‑CCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAg                                                                      <  2:72854/149‑1 (MQ=255)
                                                                        gACAGTGTGGGCGAGTCCATACATAGGAGCGGGCTTGCCCGCGAACACACATAATGGGGTGTCAGGCA‑‑CCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAg                                                                      >  8:33285/1‑147 (MQ=255)
                                                                         aCAGTGTGGGCGAGTCCATACATAGGAGCGGGCTTGCCCGCGAACACACATAATGGGGTGTCAGGCA‑‑CCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAg                                                                   >  4:43091/1‑149 (MQ=255)
                                                                             tgtgGGCGAGTCCATACATAGGAGCGGGCTTGCCCGCGAACACACATAATGGGGTGTCAGGCAGGGCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCg                                                                 >  1:246386/1‑149 (MQ=255)
                                                                                  gCGAGTACATACATAGGAGGGGGCTTGCGCGGGAACACACATAATGGGGTGTCAGGCA‑‑CCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACAGTGTGGAAAACCTTGTAGGAGCGGGTTTAc                                                         <  8:179632/149‑1 (MQ=39)
                                                                                     aGTCCATACATAGGAGCGGGCTTGCCCGCGAACACAAATAATGTGATGTCAGGCA‑‑CCGGCTTCGCCGGTGTTCGCGGTCACGCCCGCTCCCAAAAGCCCACCACAGATAGGGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCcg                                                      <  5:246385/149‑1 (MQ=39)
                                                                                          atacataGGAGCGGGCTTGCCCGCGAACACACATAATGGGGTGTCAGGCA‑‑CCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGt                                                                        >  5:207906/1‑127 (MQ=255)
                                                                                          atacataGGAGCGGGCTTGCCCGCGAACACACATAATGGGGTGTCAGGCA‑‑CCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGt                                                                        <  6:207906/127‑1 (MQ=255)
                                                                                          atacataGGAGCGGGCTTGCCCGCGAACACACATAATGGGGTGTCAGGCA‑‑CCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAg                                                 <  1:12991/149‑1 (MQ=255)
                                                                                            acataGGAGCGGGCTTGCCCGCGAACACACATAATGGGGTGTCAGGCA‑‑CCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAggg                                               >  3:92052/1‑149 (MQ=255)
                                                                                             cataGGAGCGGGCTTGCCCGCGAACACACATAATGGGGTGTCAGGCA‑‑CCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTAccc                                                       >  5:4102/1‑140 (MQ=255)
                                                                                             cataGGAGCGGGCTTGCCCGCGAACACACATAATGGGGTGTCAGGCA‑‑CCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTAccc                                                       <  6:4102/140‑1 (MQ=255)
                                                                                                  gAGCGGGCTTGCCCGCGAACACACATAATGGGGTGTCAGGCA‑‑CCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCt                                          >  6:301850/1‑148 (MQ=255)
                                                                                                      gggCTTGCCCGCGAACACACATAATGGGGTGTCAGGCA‑‑CCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGc                                     <  3:18049/149‑1 (MQ=255)
                                                                                                        gCTTGCCCGCGAACACACATAATGGGGTGTCAGGCA‑‑CCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCa                                   >  2:227003/1‑149 (MQ=255)
                                                                                                        gCTTGCCCGCGAACACACATAATGGGGTGTCAGGCA‑‑CCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCa                                   >  1:47114/1‑149 (MQ=255)
                                                                                                        gCTTGCCCGCGAACACACATAATGGGGTGTCAGGCA‑‑CCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCa                                   <  2:47114/149‑1 (MQ=255)
                                                                                                          ttGCCCGCGAACACACATAATGGGGTGTCAGGCA‑‑CCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACt                                 >  1:244851/1‑149 (MQ=255)
                                                                                                            gCCCGCGAACACACATAATGGGGTGTCAGGCA‑‑CCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTg                                <  2:246386/148‑1 (MQ=255)
                                                                                                            gCCCGCGAACACACATAATGGGGTGTCAGGCA‑‑CCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTg                                <  5:301850/148‑1 (MQ=255)
                                                                                                            ccccGCGAACACACATATTGGGGTGTCAGGCA‑‑CGGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGGGGGTTTACCCGCGAAGGGGCCCTCAGGCCACt                                 <  2:154633/146‑1 (MQ=255)
                                                                                                                        acaTAATGGGGTGTCAGGCA‑‑CCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACtgtg                                                                                   <  6:221083/86‑1 (MQ=255)
                                                                                                                        acaTAATGGGGTGTCAGGCA‑‑CCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACtgtg                                                                                   >  5:221083/1‑86 (MQ=255)
                                                                                                                        acaTAATGGGGTGTCAGGCA‑‑CCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGAc                   <  2:244851/149‑1 (MQ=255)
                                                                                                                          aTAATGGGGTGTCAGGCA‑‑CCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAggg                                               <  7:49664/119‑1 (MQ=255)
                                                                                                                          aTAATGGGGTGTCAGGCA‑‑CCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAggg                                               >  8:49664/1‑119 (MQ=255)
                                                                                                                            aaTGGGGTGTCAGGCA‑‑CCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATc               >  8:274701/1‑149 (MQ=255)
                                                                                                                               ggggTGTCAGGCA‑‑CCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGAt            >  5:163572/1‑149 (MQ=255)
                                                                                                                               ggggTGTCAGGCA‑‑CCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGAt            >  3:78942/1‑149 (MQ=255)
                                                                                                                                gggTGTCTGGCA‑‑CCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTCTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGAt            <  1:227003/149‑1 (MQ=255)
                                                                                                                                gggTGTCAGGCA‑‑CCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATg           <  3:43091/149‑1 (MQ=255)
                                                                                                                                     tCAGGCA‑‑CCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGc      <  7:33285/149‑1 (MQ=255)
                                                                                                                                       aGGCA‑‑CCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATCGCCGCtt    >  1:191748/1‑149 (MQ=255)
                                                                                                                                       aGGCA‑‑CCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCAAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCtt    >  4:69226/1‑149 (MQ=255)
                                                                                                                                         gCA‑‑CCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGa  >  5:153068/1‑149 (MQ=255)
                                                                                                                                         gCA‑‑CCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGAGAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCAGCGATGGCCGCTTgc  >  5:116425/1‑148 (MQ=255)
                                                                                                                                          cA‑‑CCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGa                    <  4:45197/130‑1 (MQ=255)
                                                                                                                                          cA‑‑CCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGa                    >  3:45197/1‑130 (MQ=255)
                                                                                                                                              ccGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAg                                       <  6:159775/109‑1 (MQ=255)
                                                                                                                                                   ttCGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGAc                   <  6:56452/124‑1 (MQ=255)
                                                                                                                                                   ttCGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGAc                   >  5:56452/1‑124 (MQ=255)
                                                                                                                                                   |                                                                                                                                             
CCTGATACATAGATGCAACCGGGTTCCTATGTTCAAAAAAGCAGCTAGACAGTGTGGGCGAGTCCATACATAGACAGTGTGGGCGAGTCCATACATAGGAGCGGGCTTGCCCGCGAACACACATAATGGGGTGTCAGGCA‑‑CCGGCTTCGCCGGTGTTCGCGGGCACGCCCGCTCCCAAAAGCCCACCACAGATAGCGAACACTGTGGAAAACCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCTCAGGCCACTGGCGCTGTAGCGACCATCGATGGCCGCTTGA  >  NC_002947/401657‑401942

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: