Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,023,216 T→C 9.8% intergenic (+75/+69) PP_5273 → / ← PP_5274 oxidoreductase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,023,2160TC 9.8% 163.3 / 7.5 61intergenic (+75/+69)PP_5273/PP_5274oxidoreductase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (25/30);  new base C (2/4);  total (27/34)
Fisher's exact test for biased strand distribution p-value = 6.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.48e-01

CCTGGCGGCAAGCACCTGCGCTCGCCACTGCTGGCGCTGGGAATGCTGTGGCATCGCCTCAAAGAGCTGACTTGAAAGTCCCCATCCTTGTGGCAGCAACTGTCTTGCTTAAAATCTAAAAGCTGACGCGCCCCCTGTGGGAGCGGGCGTGCCCGCTCCCACAGGGTCGGTGCAAGTTTCAAGTCTGGGGAACAACAGGCCCTCCAGCGACCATCGTGTCAGTCCTTCCAGAAGGGTTTCCGCCCTTCGATCCGGGCCTGCTCCCAGCTCAGCCCGACATCGCGCAACTCGTCATCGG  >  NC_002947/6023067‑6023364
                                                                                                                                                     |                                                                                                                                                    
ccTGGCGGCAAGCACCTGCGCTCGCCACTGCTGGCGCTGGGAATGCTGTGGCATCGCCTCAAAGAGCTGACTTGAAAGTCCCCATCCTTGTGGCAGCAACTGTCTTGCTTAAAATCTAAAAGCTGACGCGCCCCCTGTGGGAGCGGGCGt                                                                                                                                                      <  1:1037676/150‑1 (MQ=255)
 cTGGCGGCAAGCACCTGCGCTCGCCACTGCTGGCGCTGGGAATGCTGTGGCATCGCCTCAAAGAGCTGACTTGAAAGTCCCCATCCTTGTGGCAGCAACTGTCTTGCTTAAAATCTAAAAGCTGACGCGCCCCCTGTGGGAGCGGGCGTg                                                                                                                                                     >  2:736431/1‑150 (MQ=255)
    gcggcAAGCACCTGCGCTCGCCACTGCTGGCGCTGGGAATGCTGTGGCATCGCCTCAAAGAGCTGACTTGAAAGTCCCCATCCTTGTGGCAGCAACTGTCTTGCTTAAAATCTAAAAGCTGACGCGCCCCCTGTGGGAGCGGGCGTGccc                                                                                                                                                  >  2:63135/1‑150 (MQ=255)
       gcAAGCACCTGCGCTCGCCACTGCTGGCGCTGGGAATGCTGTGGCATCGCCTCAAAGAGCTGACTTGAAAGTCCCCATCCTTGTGGCAGCAACTGTCTTGCTTAAAATCTAAAAGCTGACGCGCCCCCTGTGGGAGCGGGCACGCCCGCt                                                                                                                                               <  2:1194281/150‑1 (MQ=255)
        cAAGCACCTGCGCTCGCCACTGCTGGCGCTGGGAATGCTGTGGCATCGCCTCAAAGAGCTGACTTGAAAGTCCCCATCCTTGTGGCAGCAACTGTCTTGCTTAAAATCTAAAAGCTGACGCGCCCCCTGTGGGAGCGGGCACGCCCGCTc                                                                                                                                              <  2:310668/150‑1 (MQ=255)
         aaGCACCTGCGCTCGCCACTGCTGGCGCTGGGAATGCTGTGGCATCGCCTCAAAGAGCTGACTTGAAAGTCCCCATCCTTGTGGCAGCAACTGTCTTGCTTAAAATCTAAAAGCTGACGCGCCCCCTGTGGGAGCGGGCACGCCCGCTcc                                                                                                                                             <  2:8054/150‑1 (MQ=255)
            cACCTGCGCTCGCCACTGCTGGCGCTGGGAATGCTGTGGCATCGCCTCAAAGAGCTGACTTGAAAGTCCCCATCCTTGTGGCAGCAACTGTCTTGCTTAAAATCTAAAAGCTGACGCGCCCCCTGTGGGAGCGGGCGTGCCCGCTCCcac                                                                                                                                          <  1:736431/150‑1 (MQ=255)
             aCCTGCGCTCGCCACTGCTGGCGCTGGGAATGCTGTGGCATCGCCTCAAAGAGCTGACTTGAAAGTCCCCATCCTTGTGGCAGCAACTGTCTTGCTTAAAATCTAAAAGCTGACGCGCCCCCTGTGGGAGCGGGCGTGCCCGCTCCcaca                                                                                                                                         <  2:470382/150‑1 (MQ=255)
             aCCTGCGCTCGCCACTGCTGGCGCTGGGAATGCTGTGGCATCGCCTCAAAGAGCTGAAATGAAAGTCCGCATCCTTGTGGCAGCAACTGTCTTGCTTAAAATCTAAAAGCTGACGCGCCCCCTGTGGGAGCGGGGGTGGCAGCTCCcaca                                                                                                                                         >  1:392796/1‑150 (MQ=255)
                 gcgcTCGCCACTGCTGGCGCTGGGAATGCTGTGGCATCGCCTCAAAGAGCTGACTTGAAAGTCCCCATCCTTGTGGCAGCAACTGTCTTGCTTAAAATCTAAAAGCTGACGCGCCCCCTGTGGGAGCGGGCGTGCCCGCTCCCACAGGGt                                                                                                                                     <  1:280395/150‑1 (MQ=255)
                     tCGCCACTGCTGGCGCTGGGAATGCTGTGGCATCGCCTCAAAGAGCTGACTTGAAAGTCCCCATCCTTGTGGCAGCACCTGTCTTGCTTACAATCTAAAAGCTGATGCGCCCCCTGTGGGAGCGGGCGTGCCCGCTCCCACAGGGTCGGt                                                                                                                                 <  1:885732/150‑1 (MQ=255)
                     tCGCCACTGCTGGCGCTGGGAATGCTGTGGCATCGCCTCAAAGAGCTGACTTGAAAGTCCCCATCCTTGTGGCAGCAACTGTCTTGCTTAAAATCTAAAAGCTGACGCGCCCCCTGTGGGAGCGGGCGTGCCCGCTCCCACAGGGTCGGt                                                                                                                                 >  2:563001/1‑150 (MQ=255)
                        ccACTGCTGGCGCTGGGAATGCTGTGGCATCGCCTCAAAGAGCTGACTTGAAAGTCCCCATCCTTGTGGCAGCAACTGTCTTGCTTAAAATCTAAAAGCTGACGCGCCCCCTGTGGGAGCGGGCACGCCCGCTCCcaca                                                                                                                                         <  1:315518/139‑1 (MQ=255)
                        ccACTGCTGGCGCTGGGAATGCTGTGGCATCGCCTCAAAGAGCTGACTTGAAAGTCCCCATCCTTGTGGCAGCAACTGTCTTGCTTAAAATCTAAAAGCTGACGCGCCCCCTGTGGGAGCGGGCACGCCCGCTCCcaca                                                                                                                                         >  2:315518/1‑139 (MQ=255)
                          aCTGCTGGCGCTGGGAATGCTGTGGCATCGCCTCAAAGAGCTGACTTGAAAGTCCCCATCCTTGTGGCAGCAACTGTCTTGCTTAAAATCTAAAAGCTGACGCGCCCCCTGTGGGAGCGGGCGTGCCCGCTCCCACAGGGTCGGTGCAAg                                                                                                                            >  1:809733/1‑150 (MQ=255)
                             gctgGCGCTGGGAATGCTGTGGCATCGCCTCAAAGAGCTGACTTGAAAGTCCCCATCCTTGTGGCAGCAACTGTCTTGCTTAAAATCTAAAAGCTGACGCGCCCCCTGTGGGAGCGGGCGTGCCCGCTCCCACAGGGTCGGTGCAAGttt                                                                                                                         <  2:310884/150‑1 (MQ=255)
                                ggCGCTGGGAATGCTGTGGCATCGCCTCAAAGAGCTGACTTGAAAGTCCCCATCCTTGTGGCAGCAACTGTCTTGCTTAAAATCTAAAAGCTGACGCGCCCCCTGTGGGAGCGGGCGTGCCCGCTCCCACAGGGTCGGTGCAAGTTTCaa                                                                                                                      <  1:413396/150‑1 (MQ=255)
                                ggCGCTGGGAATGCTGTGGCATCGCCTCAAAGAGCTGACTTGAAAGTCCCCATCCTTGTGGCAGCAACTGTCTTGCTTAAAATCTAAAAGCTGACGCGCCCCCTGTGGGAGCGGGCGTGCCCGCTCCCACAGGGTCGGTGCAAGTTTCaa                                                                                                                      <  1:751159/150‑1 (MQ=255)
                                ggCGCTGGGAATGCTGTGGCATCGCCTCAAAGAGCTGACTTGAAAGTCCCCATCCTTGTGGCAGCAACTGTCTTGCTTAAAATCTAAAAGCTGACGCGCCCCCTGTGGGAGCGGGCGTGCCCGCTCCCACAGGGTCGGTGCAAGTTTCaa                                                                                                                      <  2:578721/150‑1 (MQ=255)
                                ggCGCTGGGAATGCTGTGGCATCGCCTCAAAGAGCTGACTTGAAAGTCCCCATCCTTGTGGCAGCAACTGTCTTGCTTAAAATCTAAAAGCTGACGCGCCCCCTGTGGGAGCGGGCGTGCCCGCTCCCACAGGGTCGGTGCAAGTTTCaa                                                                                                                      <  2:1037737/150‑1 (MQ=255)
                                ggCGCTGGGAATGCTGTGGCATCGCCTCAAAGAGCTGACTTGAAAGTCCCCATCCTTGTGGCAGCAACTGTCTTGCTTAAAATCTAAAAGCTGACGCGCCCCCTGTGGGAGCGGGCGTGCCCGCTCCCACAGGGTCGGTGCAAGTTTCaa                                                                                                                      <  2:809733/150‑1 (MQ=255)
                                     tGGGAATGCTGTGGCATCGCCTCAAAGAGCTGACTTGAAAGTCCCCATCCTTGTGGCAGCAACTGTCTTGCTTAAAATCTAAAAGCTGACGCGCCCCCTGTGGGAGCGGGCGTGCCCGCTCCCACAGGGTCGGTGCAAGTTTCAAGTCTg                                                                                                                 <  1:886591/150‑1 (MQ=255)
                                     tGGGAATGCTGTGGCATCGCCTCAAAGAGCTGACTTGAAAGTCCCCATCCTTGTGGCAGCAACTGTCTTGCTTAAAATCTAAAAGCTGACGCGCCCCCTGTGGGAGCGGGCGTGCCCGCTCCCACAGGGTCGGTGCAAGTTTCAAGTCTg                                                                                                                 >  2:640627/1‑150 (MQ=255)
                                     tGGGAATGCTGTGGCATCGCCTCAAAGAGCTGACTTGAAAGTCCCCATCCTTGTGGCAGCAACTGTCTTGCTTAAAATCTAAAAGCTGACGCGCCCCCTGTGGGAGCGGGCGTGCCCGCTCCCACAGGGGCGGTGCAAGTTTCAAGTCTg                                                                                                                 >  1:637484/1‑150 (MQ=255)
                                      gggAATGCTGTGGCATCGCCTCAAAGAGCTGACTTGAAAGTCCCCATCCTTGTGGCAGCAACTGTCTTGCTTAAAATCTAAAAGCTGACGCGCCCCCTGTGGGAGCGGGCGTGCCCGCTCCCACAGGGTCGGTGCAAGTTTCAAGTCTgg                                                                                                                <  2:987454/150‑1 (MQ=255)
                                       ggAATGCTGTGGCATCGCCTCAAAGAGCTGACTTGAAAGTCCCCATCCTTGTGGCAGCAACTGTCTTGCTTAAAATCTAAAAGCTGACGCGCCCCCTGTGGGAGCGGGCGTGCCCGCTCCCACAGGGTCGGTGCAAGTTTCAAGTCTggg                                                                                                               <  1:598204/150‑1 (MQ=255)
                                       ggAATGCTGTGGCATCGCCTCAAAGAGCTGACTTGAAAGTCCCCATCCTTGTGGCAGCAACTGTCTTGCTTAAAATCTAAAAGCTGACGCGCCCCCTGTGGGAGCGGGCGTGCCCGCTCCCACAGGGTCGGTGCAAGTTTCAAGTCTggg                                                                                                               <  1:497932/150‑1 (MQ=255)
                                        gAATGCTGTGGCATCGCCTCAAAGAGCTGACTTGAAAGTCCCCATCCTTGTGGCAGCAACTGTCTTGCTTAAAATCTAAAAGCTGACGCGCCCCCTGTGGGAGCGGGCGTGCCCGCTCCCACAGGGTCGGTGCAAGTTTCAAGTCGgggg                                                                                                              >  1:632020/1‑150 (MQ=255)
                                              tgtgGCATCGCCTCAAAGAGCTGACTTGAAAGTCCCCATCCTTGTGGCAGCAACTGTCTTGCTTAAAATCTAAAAGCTGACGCGCCCCCTGTGGGAGCGGGCGTGCCCGCTCCCACAGGGTCGGTGCAAGGTTCAAGTCTGGGGaacaac                                                                                                        >  2:1409910/1‑150 (MQ=255)
                                                tgGCATCGCCTCAAAGAGCTGACTTGAAAGTCCCCATCCTTGTGGCAGCAACTGTCTTGCTTAAAATCTAAAAGCTGACGCGCCCCCTGTGGGAGCGGGCGTGCCCGCTCCCACAGGGTCGGTGCAAGTTTCAAGTCTGGGGAACAACAg                                                                                                      >  1:713275/1‑150 (MQ=255)
                                                tgGCATCGCCTCAAAGAGCTGACTTGAAAGTCCCCATCCTTGTGGCAGCAACTGTCTTGCTTAAAATCTAAAAGCTGACGCGCCCCCTGTGGGAGCGGGCGTGCCCGCTCCCACAGGGTCGGTGCAAGTTTCAAGTCTGGGGAACAACAg                                                                                                      <  2:1204727/150‑1 (MQ=255)
                                                  ccagtgCCTCAAAGTGCTTACGTGAAAGTCCCCAACCTTGTGGCAGCAACTGTCTTGCTTAAAATCTAACACCTGACGCGCCCCCTCTGGGAGCAGGCGTTCACGCTCCCACAGGGTCGGTGCAAGTTTCAAGTCTGGGGAACAACAGGc                                                                                                    <  2:392796/145‑1 (MQ=255)
                                                     tCGCCTCAAAGAGCTGACTTGAAAGTCCCCATCCTTGTGGCAGCAACTGTCTTGCTTAAAATCTAAAAGCTGACGCGCCCCCTGTGGGAGCGGGCGTGCCCGCTCCCACAGGGTCGGTGCAAGTTTCAAGTCTGGGGAACAACAGGCCCt                                                                                                 <  2:713275/150‑1 (MQ=255)
                                                         cTCAAAGAGCTGACTTGAAAGTCCCCATCCTTGTGGCAGCAACTGTCTTGCTTAAAATCTAAAAGCTGACGCGCCCCCTGTGGGAGCGGGCGTGCCCGCTCCcaca                                                                                                                                         >  1:1166932/1‑106 (MQ=255)
                                                         cTCAAAGAGCTGACTTGAAAGTCCCCATCCTTGTGGCAGCAACTGTCTTGCTTAAAATCTAAAAGCTGACGCGCCCCCTGTGGGAGCGGGCGTGCCCGCTCCcaca                                                                                                                                         <  2:1166932/106‑1 (MQ=255)
                                                                       ttGAAAGTCCCCATCCTTGTGGCAGCAACTGTCTTGCTTAAAATCTAAAAGCTGACGCGCCCCCTGTGGGAGCGGGCGTGCCCGCTCCCACAGGGTCGGTGCAAGTTTCAAGTCTGGGGAACAACAGGCCCTCCAGCGACCATc                                                                                     <  1:108729/144‑1 (MQ=255)
                                                                       ttGAAAGTCCCCATCCTTGTGGCAGCAACTGTCTTGCTTAAAATCTAAAAGCTGACGCGCCCCCTGTGGGAGCGGGCGTGCCCGCTCCCACAGGGTCGGTGCAAGTTTCAAGTCTGGGGAACAACAGGCCCTCCAGCGACCATc                                                                                     >  2:108729/1‑144 (MQ=255)
                                                                       ttGAAAGTCCCCATCCTTGTGGCAGCAACTGTCTTGCTTAAAATCTAAAAGCTGACGCGCCCCCTGTGGGAGCGGGCGTGCCCGCTCCCACAGGGTCGGTGCAAGTTTCAAGTCTGGGGAACAACAGGCCCTCCAGCGACCATCgt                                                                                   >  1:1057407/1‑146 (MQ=255)
                                                                       ttGAAAGTCCCCATCCTTGTGGCAGCAACTGTCTTGCTTAAAATCTAAAAGCTGACGCGCCCCCTGTGGGAGCGGGCGTGCCCGCTCCCACAGGGTCGGTGCAAGTTTCAAGTCTGGGGAACAACAGGCCCTCCAGCGACCATCgt                                                                                   <  2:1057407/146‑1 (MQ=255)
                                                                       ttGAAAGTCCCCATCCTTGTGGCAGCAACTGTCTTGCTTAAAATCTAAAAGCTGACGCGCCCCCTGTGGGAGCGGGCGTGCCCGCTCCCACAGGGTCGGTGCAAGTTTCAAGTCTGGGGAACAACAGGCCCTCCAGCGACCATCGTGTCa                                                                               <  1:63135/150‑1 (MQ=255)
                                                                                     ccTTGTGGCAGCAAATGTCTTGCTTAAAATCTAAAAGCTGACGCGCCCCCTGTGGGAGCGGGCGTGCCCGCTCCCACAGGGTCGGTGCAAGTTTAAAGTCTGGGGAACAACAGGCCCTCCCGCGACCATCGTGTCAGTCCTTCCAGAAgg                                                                 <  2:1173802/150‑1 (MQ=255)
                                                                                       ttGTGGCAGCAACTGTCTTGCTTAAAATCTAAAAGCTGACGCGCCCCCTGTCGGAGCGGGCGTGCCCGCTCCCAAAGGGGCGGTGCACGTTTAATGACTGGGGAAGAACAGGCCCTCCAGCGACCGTCGTGTCAGTCCTTCCAGAAGGGt                                                               >  2:338428/1‑150 (MQ=255)
                                                                                         gtgGCAGCAACTGTCTTGCTTAAAATCTAAAAGCTGACGCGCCCCCTGTGGGAGCGGGCGTGCCCGCTCCCACAGGGTCGGTGCAAGTTTCAAGTCTGGGGAACAACAGGCCCTCCAGCGACCATCGTGTCAGTCCTTCCAGAAGGGttt                                                             >  1:498625/1‑150 (MQ=255)
                                                                                            gcagcaACTGTCTTGCTTAAAATCTAAAAGCTGACGCGCCCCCTGTGGGAGCGGGCGTGCCCGCTCCCACAGGGTCGGTGCAAGTTTCAAGTCTggg                                                                                                               <  2:817530/97‑1 (MQ=255)
                                                                                            gcagcaACTGTCTTGCTTAAAATCTAAAAGCTGACGCGCCCCCTGTGGGAGCGGGCGTGCCCGCTCCCACAGGGTCGGTGCAAGTTTCAAGTCTggg                                                                                                               >  1:817530/1‑97 (MQ=255)
                                                                                            gcagcaACTGTCTTGCTTAAAATCTAAAAGCTGACGCGCCCCCTGTGGGAGCGGGCGTGCCCGCTCCCACAGGGTCGGTGCAAGTTTCAAGTCTGGGGaacaa                                                                                                         >  1:355847/1‑103 (MQ=255)
                                                                                            gcagcaACTGTCTTGCTTAAAATCTAAAAGCTGACGCGCCCCCTGTGGGAGCGGGCGTGCCCGCTCCCACAGGGTCGGTGCAAGTTTCAAGTCTGGGGaacaa                                                                                                         <  2:355847/103‑1 (MQ=255)
                                                                                               gcaACTGTCTTGCTTAAAATCTAAAAGCTGACGCGCCCCCTGTGGGAGCGGGCGTGCCCGCTCCCACAGGGTCGGTGCAAGTTTCAAGTCTGGGGAACAACAGGCCCTCCAGCGACCATCGTGTCAGTCCTTCCAGAAGGGTTTCCGccc                                                       <  1:640627/150‑1 (MQ=255)
                                                                                                  aCTGTCTTGCTTAAAATCTAAAAGCTGACGCGCCCCCTGTGGGAGCGGGCGTGCCCGCTCCCACAGGGTCGGTGCAAGTTTCAAGTCTGGGGAACAACAGGCCCTCCAGCGACCATCGTGTCAGTCCTTCCCGTAGGGTTTCCGCCCTTc                                                    <  1:409719/150‑1 (MQ=255)
                                                                                                    tGTCTTGCTTAAAATCTAAAAGCTGACGCGCCCCCTGTGGGCGCGGGCGTGCCCGCTCCCACAGGGTCGGTGCAAGTTTCAAGTCTGGGGAACAACAGGCCCTCCAGCGACCATCGTGTCAGTCCTTCCAGAAGGGTTTCCGCCCTTCGa                                                  >  1:997652/1‑150 (MQ=255)
                                                                                                        ttGCTTAAAATCTAAAAGCTGACGCGCCCCCTGTGGGAGCGGGCGTGCCCGCTCCCACAGGGTCGGTGCAAGTTTCAAGTCTGGGGAACAACAGGCCCTCCAGCGACCATCGTGTCAGTCCTTCCAGAAGGGTTTCCGCCCTTCGATCCg                                              >  2:768358/1‑150 (MQ=255)
                                                                                                            ttAAAATCTAAAAGCTGACGCGCCCCCTGTGGGAGCGGGCGTTCCCGCTCCCCCAGGGTCGGTGCAAGTTTCAAGTCTGGCGAACAACAGGCCCTCCAGCGACCATCGTGTCAGTCCTTCCAGAAGGGTTTCCGCCCTTCGATCCGGGcc                                          >  1:1111989/1‑150 (MQ=255)
                                                                                                            ttAAAATCTAAAAGCTGACGCGCCCCCTGTGGGAGCGGGCGTGCCCGCTCCCACAGGGTCGGTGCAAGTTTCAAGTCTGGGGAACAACAGGCCCTCCAGCGACCATCGTGTCAGTCCTTCCAGAAGGGTTTCCGCCCTTCGATCCGGGcc                                          >  1:1251439/1‑150 (MQ=255)
                                                                                                                aaTCTAAAAGCTGACGCGCCCCCTGTGGGAGCGGGCGTGCCCGCTCCCACAGGGTCGGTGCAAGTTTCAAGTCTGGGGaacaac                                                                                                        <  2:1265851/84‑1 (MQ=255)
                                                                                                                aaTCTAAAAGCTGACGCGCCCCCTGTGGGAGCGGGCGTGCCCGCTCCCACAGGGTCGGTGCAAGTTTCAAGTCTGGGGaacaac                                                                                                        >  1:1265851/1‑84 (MQ=255)
                                                                                                                 aTCTAAAAGCTGACGCGCCCCCTGTGGGAGCGGGCGTGCCCGCTCCCACAGGGTCGGTGCAAGTTTCAAGTCTGGGGAACAACAGGCCCTCCAGCGACCATCGTGTCAGTCCTTCCAGAAGGGTTTCCGCCCTTCGATCCGGGCCTGCTc                                     <  2:997652/150‑1 (MQ=255)
                                                                                                                   cTAAAAGCTGACGCGCCCCCTGTGGGAGCGGGCGTGCCCGCTCCCACAGGGTCGGTGCAAGTTTCAAGTCTGGGGAACAACAGGCCCTCCAGCGACCATCGTGTCAGTCCTTCCAGAAGGGTTTCCGCCCTTCGATCCGGGCCTGCTccc                                   <  1:1266445/150‑1 (MQ=255)
                                                                                                                                    cccTGTGGGAGCGGGCGTGCCCGCTCCCACAGGGTCGGTGCAAGTTTCAAGTCTGGGGAACAACAGGCCCTCCAGCGACCATCGTGTCAGACGTTCCGCAGGTGCAGGCGTGGTGTGATCCGGGCCTGCTCCCAGCTCAGCCCGACATcg                  >  1:1385654/1‑150 (MQ=255)
                                                                                                                                     ccTGTGGGAGCGGGCGTGCCCGCTCCCACAGGGTCGGTGCAAGTTTCAAGTCTGGGGAACAACAGGCCCTCCAGCGACCATCGTGTCAGTCCTTCCAGAAGGGTTTCCGCCCTTCGATCCGGGCCTGCTCCCAGCTCAGCCCGACATcgc                 <  1:563001/150‑1 (MQ=255)
                                                                                                                                      cTGTGGGAGCGGGCACGCCCGCTCCCACAGGGTCGGTGCAAGTTTCAAGTCTGGGGAACAACAGGCCCTCCAGCGACCATCGTGTCAGTCCTTCCAGAAGGGTTTCCGCCCTTCGATCCGGGCCTGCTCCCAGCTCAGCCCGACATcgcg                >  2:1094223/1‑150 (MQ=255)
                                                                                                                                                 ggCGTGCCCGCTCCCACAGGGTCGGTGCAAGTTTCAAGTCTGGGGAACAACAGGCCCTCCAGCGACCATCGTGTCAGTCCTTCCAGAAGGGTTTCCGCCCTTCGATCCGGGCCTGCTCCCAGCTCAGCCCGACATCGCGCAACTCGTCAt     >  2:1112447/1‑150 (MQ=255)
                                                                                                                                                    gTGCCCGCTCCCACAGGGTCGGTGCAAGTTTCAAGTCTGGGGAACAACAGGCCCTCCAGCGACCATCGTGTCAGTCCTTCCAGAAGGGTTTCCGCCCTTCGATCCGGGCCTGCTCCCAGCTCAGCCCGACATCGCGCAACTCGTCATCgg  >  2:925720/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CCTGGCGGCAAGCACCTGCGCTCGCCACTGCTGGCGCTGGGAATGCTGTGGCATCGCCTCAAAGAGCTGACTTGAAAGTCCCCATCCTTGTGGCAGCAACTGTCTTGCTTAAAATCTAAAAGCTGACGCGCCCCCTGTGGGAGCGGGCGTGCCCGCTCCCACAGGGTCGGTGCAAGTTTCAAGTCTGGGGAACAACAGGCCCTCCAGCGACCATCGTGTCAGTCCTTCCAGAAGGGTTTCCGCCCTTCGATCCGGGCCTGCTCCCAGCTCAGCCCGACATCGCGCAACTCGTCATCGG  >  NC_002947/6023067‑6023364

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