Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,650,018 T→A 8.4% intergenic (+1495/+196) PP_3214 → / ← PP_3216 hypothetical protein/cysteine dioxygenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,650,0180TA 8.4% 231.1 / 8.0 80intergenic (+1495/+196)PP_3214/PP_3216hypothetical protein/cysteine dioxygenase
Reads supporting (aligned to +/- strand):  ref base T (39/34);  new base A (2/5);  total (41/39)
Fisher's exact test for biased strand distribution p-value = 2.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.72e-01

TTTGCTTGGCACGTTGACATCACGAAACAACTCGCCGCCCGCAGCGCGCCAGCCGGCAAGGCCGCCATCAAGCAAGGCAACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGAGCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTATGTTGGCCGCGAACAGCGGATGGGCCTGGGCGAAAGGGTCTTCTTCGCGTACGTCGATCAGGGCCAGCTCCTCCTTGGC  >  NC_002947/3649873‑3650166
                                                                                                                                                 |                                                                                                                                                    
tttGCTTGGCACGTTGACATCACGAAACAACTCGCCGCCCGCAGCGCGCCAGCCGGCAAGGCCGCCATCAAGCAAGGCAACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGAGCCCCTCACCGTCATCGTAGa                                                                                                                                                   >  6:18308/1‑149 (MQ=255)
  tGCTTGGCACGTTGACATCACGAAACAACTCGCCGCCCGCAGCGCGCCAGCCGGCAAGGCCGCCATCAAGCAAGGCAACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGAGCCCCTCACCGTCATCGTAGAc                                                                                                                                                  <  2:233745/148‑1 (MQ=255)
     ttGGCACGTTGACATCACGAAACAACTCGCCGCCCGCAGCGCGCCAGCCGGCAAGGCCGCCATCAAGCAAGGCAACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGAGCCCCTCACCGTCATCGTAGACGGt                                                                                                                                               >  3:21387/1‑148 (MQ=255)
       ggCACGTTGACATCACGAAACAACTCGCCGCCCGCAGCGCGCCAGCCGGCAAGGCCGCCATCAAGCAAGGCAACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGAGCCCCTCACCGTCATCGTAGACGGTGAg                                                                                                                                            >  4:275605/1‑149 (MQ=255)
        gCACGTTGACATCACGAAACAATTCGCCGCCCGCAGCGCGCCAGCCGGCAAGGCCGCCATCAAGCAAGGCAACATTGACATAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGAGCCCCTCACCGTCTACGATGACGGTGAg                                                                                                                                            <  6:203889/148‑1 (MQ=255)
        gCACGTTGACATCACGAAACAACTCGCCGCCCGCAGCGCGCCAGCCGGCAAGGCCGCCATCAAGCAAGGCAACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGAGCCCCTCACCGTCTACGATGACGGTGAg                                                                                                                                            <  2:226312/148‑1 (MQ=255)
        gCACGTTGACATCACGAAACAACTCGCCGCCCGCAGCGCGCCAGCCGGCAAGGCCGCCATCAAGCAAGGCAACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGAGCCCCTCACCGTCTACGATGACGGTGAg                                                                                                                                            <  2:232131/148‑1 (MQ=255)
        gCACGTTGACATCACGAAACAACTCGCCGCCCGCAGCGCGCCAGCCGGCAAGGCCGCCATCAAGCAAGGCAACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGAGCCCCTCACCGTCTACGATGACGGTGAg                                                                                                                                            <  6:235110/148‑1 (MQ=255)
        gCACGTTGACATCACGAAACAACTCGCCGCCCGCAGCGCGCCAGCCGGCAAGGCCGCCATCAAGCAAGGCAACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGAGCCCCTCACCGTCTACGATGACGGTGAg                                                                                                                                            <  5:99527/148‑1 (MQ=255)
            gTTGACATCACGAAACAACTCGCCGCCCGCAGCGCGCCAGCCGGCAAGGCCGCCATCAAGCAAGGCAACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGAGCCCCTCACCGTCATCGAACACGGTGAGTGCgg                                                                                                                                       >  7:248278/1‑149 (MQ=255)
                     aCGAAACAACTCGCCGCCCGCAGCGCGCCAGCCGGCAAGGCCGCCATCAAGCAAGGCAACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGAGCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCccgcc                                                                                                                              >  7:216448/1‑149 (MQ=255)
                       gAAACAACTCGCCGCCCGCAGCGCGCCAGCCGGCAAGGCCGCCATCAAGCAAGGCAACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGAGCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCgg                                                                                                                            >  4:76975/1‑149 (MQ=255)
                       gAAACAACTCGCCGCCCGCAGCGCGCCAGCCGGCAAGGCCGCCATCAAGCAAGGCAACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGAGCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCgg                                                                                                                            >  2:38183/1‑149 (MQ=255)
                       gAAACAACTCGCCGCCCGCAGCGCGCCAGCCGGCAAGGCCGCCATCAAGCAAGGCAACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGAGCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCgg                                                                                                                            <  4:226450/149‑1 (MQ=255)
                        aaaCAACTCGCCGCCCGCAGCGCGCCAGCCGGCAAGGCCGCCATCAAGCAAGGCAACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGAGCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCggg                                                                                                                           >  2:325/1‑149 (MQ=255)
                         aacaacTCGCCGCCCGCAGCGCGCCAGCCGGCAAGGCCGCCATCAAGCAAGGCAACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGAGCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCggg                                                                                                                           <  4:205814/148‑1 (MQ=255)
                            aacTCGCCGCCCGCAGCGCGCCAGCCGGCAAGGCCGCCATCAAGCAAGGCAACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGAGCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCccgcccgcgttc                                                                                                                       >  5:218188/1‑142 (MQ=255)
                            aacTCGCCGCCCGCAGCGCGCCAGCCGGCAAGGCCGCCATCAAGCAAGGCAACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGAGCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAAc                                                                                                                       >  7:269237/1‑149 (MQ=255)
                               tCGCCGCCCGCAGCGCGCCAGCCGGCAAGGCCGCCATCAAGCAAGGCAACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGAGCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAAcgcg                                                                                                                    >  7:49618/1‑149 (MQ=255)
                                 gccgccCGCAGCGCGCCAGCCGGCAAGGCCGCCATCAAGCAAGGCAACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGAGCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCggg                                                                                                                  >  4:114946/1‑149 (MQ=255)
                                 gccgccCGCAGCGCGCCAGCCGGCAAGGCCGCCATCAAGCAAGGCAACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGAGCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCgg                                                                                                                   <  4:65939/148‑1 (MQ=255)
                                 gccgccCGCAGCGCGCCAGCCGGCAAGGCCGCCATCAAGCAAGGCAACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCAGCGAGCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCgg                                                                                                                   <  6:190294/148‑1 (MQ=255)
                                     cccgcAGCGCGCCAGCCGGCAAGGCCGCCATCAAGCAAGGCAACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGAGCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGaa                                                                                                              <  3:114946/149‑1 (MQ=255)
                                     cccgcAGCGCGCCAGCCGGCAAGGCCGCCATCAAGCAAGGCAACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGAGCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGaa                                                                                                              <  2:337609/149‑1 (MQ=255)
                                       cgcAGCGCGCCAGCCGGCAAGGCCGCCATCAAGCAAGGCAACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGAGCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCg                                                                                                                >  5:85288/1‑145 (MQ=255)
                                          aGCGCGCCAGCCGGCAAGGCCGCCATCAAGCAAGGCAACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGAGCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGAttt                                                                                                         >  8:296395/1‑149 (MQ=255)
                                          aGCGCGCCAGCCGGCAAGGCCGCCATCAAGCAAGGCAACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGAGCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGAttt                                                                                                         <  1:325/149‑1 (MQ=255)
                                           gcgcgcCAGCCGGCAAGGCCGCCATCAAGCAAGGCAACATTGTCGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGAGCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGa                                                                                                               >  4:88755/1‑142 (MQ=255)
                                           gcgcgcCAGCCGGCAAGGCCGCCATCAAGCAAGGCAACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGAGCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGa                                                                                                               <  3:88755/142‑1 (MQ=255)
                                           gcgcgcCAGCCGGCAAGGCCGCCATCAAGCAAGGCAACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGAGCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTc                                                                                                        <  5:18308/149‑1 (MQ=255)
                                               gcCAGCCGGCAAGGCCGCCATCAAGCAAGGCAACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGAGCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGc                                                                                                    <  8:216448/149‑1 (MQ=255)
                                                       gCAAGGCCGCCATCAAGCAAGGCAACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGAGCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGA‑TTCCAGCTCCAGTTTc                                                                                           <  3:275605/149‑1 (MQ=255)
                                                               gccATCAAGCAAGGCAACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGAGCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAggg                                                                                    <  1:180801/149‑1 (MQ=255)
                                                                  aTCAAGCAAGGCAACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGAGCCCCTCACCGTCATCGTAGATGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTa                                                                                 >  2:320110/1‑149 (MQ=255)
                                                                     aagcaagGCAACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGAGCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTATGt                                                                              >  3:294360/1‑149 (MQ=255)
                                                                      agcaagGCAACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGAGCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTATGtt                                                                             >  2:243047/1‑149 (MQ=255)
                                                                       gcaagTCAACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGTCGAGCCCCTCACCGACATCGTAGAAGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTATGTTg                                                                            >  5:268303/1‑149 (MQ=255)
                                                                       gcaagGCAACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGAGCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTATGTTg                                                                            <  2:355393/149‑1 (MQ=255)
                                                                         aagGCAACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGAGCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTATGTTGGc                                                                          <  3:76975/149‑1 (MQ=255)
                                                                           ggCAACATTAACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGAGCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGa                                                                                         <  4:110970/132‑1 (MQ=255)
                                                                           ggCAACATTAACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGAGCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGa                                                                                         >  3:110970/1‑132 (MQ=255)
                                                                               aCATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGAGCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTATGTTGGCCGCGaa                                                                    <  6:385096/149‑1 (MQ=255)
                                                                                 aTTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGAGCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTATGTTGGCCGCGaa                                                                    >  7:318281/1‑147 (MQ=255)
                                                                                 aTTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGAGCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTATGTTGGCCGCGAACa                                                                  <  8:49618/149‑1 (MQ=255)
                                                                                       gTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGAGCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTATGTTGGCCGCGAACAGCGGAt                                                            >  8:47326/1‑149 (MQ=255)
                                                                                          gCCCATGGCCAGTAAACGCTGCTCGGCGACTGCGGCGAGCCCCTCAGCGTCATCGTAGACGGTGAGTGCGGTGTCCCGGCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTATGTTGGCCGCGAACAGCGGATgg                                                          <  8:248278/148‑1 (MQ=255)
                                                                                          gCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGAGCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTATGTTGGCCGCGAACAGCGGATggg                                                         >  7:19625/1‑149 (MQ=255)
                                                                                           cccAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGAGCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTATGTTGGCCGCGAACAGCGGATGGGc                                                        >  7:278574/1‑149 (MQ=255)
                                                                                            ccAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGAGCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTATGTTGGCCGCGAACAGCGGATGGGcc                                                       >  4:319845/1‑149 (MQ=255)
                                                                                             cAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGAGCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTATGTTGGCCGCGAACAGCGGATGGGCCt                                                      >  5:151783/1‑149 (MQ=255)
                                                                                                 gCCAGTAAACGCTGCGCGGCGACTGCGGCGAGCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTATGTTGGCCGCGAACAGCGGATGGGCCTGGGc                                                  >  1:61543/1‑149 (MQ=255)
                                                                                                   cAGTAAACGCTGCGCGGCGACTGCGGCGAGCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGACCGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTATGTTGGCCGCGAACCGCGGATGGTCCTGGGCg                                                 >  8:311537/1‑148 (MQ=255)
                                                                                                     gTAAACGCTGCGCGGCGACTGCGGCGAGCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTATGTTGGCCGCGAACAGCGGATGGGCCTGGGCGaaa                                              <  7:296395/149‑1 (MQ=255)
                                                                                                     gTAAACGCTGCGCGGCGACTGCGGCGAGCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTATGTTGGCCGCGAACAGCGGATGGGCCTGGGCGaaa                                              <  8:269237/149‑1 (MQ=255)
                                                                                                      tAAACGCTGCGCGGCGACTGCGGCGAGCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTATGTTGGCCGCGAACAGCGGATGGGCCTGGGCGaaa                                              <  4:27298/148‑1 (MQ=255)
                                                                                                      tAAACGCTGCGCGGCGACTGCGGCGAGCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTATGTTGGCCGCGAACAGCGGATGGGCCTGGGCGa                                                >  3:27298/1‑146 (MQ=255)
                                                                                                      tAAACGCTGCGCGGCGACTGCGGCGAGCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTATGTTGGCCGCGAACAGCGGATGGGCCTGGGCGAAAg                                             >  1:113179/1‑149 (MQ=255)
                                                                                                        atcGCTTCGCGGAGTCTGCGTCGTGCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTATGTTGGCCGCGAACAGCGGATGGGCCTGGGCGAAAgg                                            <  6:268303/146‑1 (MQ=255)
                                                                                                        aaCGCTGCGCGGCGACTGCGGCGAGCCCCTGACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTATGTTGGCCGCGAACAGCGGATGGGCCTGGGCGAAAgg                                            <  4:334961/148‑1 (MQ=255)
                                                                                                        aaCGCTGCGCGGCGACTGCGGCGAGCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTATGTTGGCCGCGAACAGCGGATGGGCCTGGGCGAAAggg                                           >  7:116020/1‑149 (MQ=255)
                                                                                                           gCTGCGCGGCGACTGCGGCGAGCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTATGTTGGCCGCGAACAGCGGATGGGCCTGGGCGAAAGGGtct                                        <  7:47326/149‑1 (MQ=255)
                                                                                                                gcggcgACTGCGGCGAGCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGt                                                                                  >  8:218136/1‑102 (MQ=255)
                                                                                                                gcggcgACTGCGGCGAGCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGt                                                                                  <  7:218136/102‑1 (MQ=255)
                                                                                                                gcggcgACTGCGGCGAGCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTATGTTGGCCGCGAACAGCGGATGGGCCTGGGCGAAAGGGTCttcttc                                   <  3:319845/149‑1 (MQ=255)
                                                                                                                     gACTGCGGCGAGCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTATGTTGGCCGCGAACAGCGGATGGGCCTGGGCGAAAGGGTCTTCTTCGCGTa                              >  1:359852/1‑149 (MQ=255)
                                                                                                                          cggcgAGCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTATGTTGGCCGCGAACAGCGGATGGGCCTGGGCGAAAGGGTCTTCTTCGCGTACGTc                          <  7:119606/148‑1 (MQ=255)
                                                                                                                           ggcgAGCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTATGTTGGCCGCGAACAGCGGATGGGCCTGGGCGAAAGGGTCTTCTTCGCGTACGTCGa                        >  4:290617/1‑149 (MQ=255)
                                                                                                                            gcgAGCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTATGTTGGCCGCGAACAGCGGATGGGCCTGGGCGAAAGGGTCTTCTTCGCGTACGTCGAt                       <  3:290617/149‑1 (MQ=255)
                                                                                                                             cgAGCCTCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTATGTTGGCCGCGAACAGCGGATGGGCCTGGGCGAAAGGGTCTTCTTCGCGTACGt                           <  1:134373/144‑1 (MQ=255)
                                                                                                                             cgAGCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTATGTTGGCCGCGAACAGCGGATGGGCCTGGGCGAAAGGGTCTTCTTCGCGTACGt                           >  2:134373/1‑144 (MQ=255)
                                                                                                                             cgAGCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTATGTTGGCCGCGAACAGCGGATGGGCCTGGGCGAAAGGGTCTTCTTCGCGTACGTCGATc                      <  7:337902/149‑1 (MQ=255)
                                                                                                                               ccgcacTCACCGTCTACGATGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTATGTTGGCCGCGAACAGCGGATGGGCCTGGGCGAAAGGGTCTTCTTCGCGTACGTCGATCAg                    >  2:129240/6‑149 (MQ=255)
                                                                                                                                gCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTATGTTGGCCGCGAACAGCGGATGGGCCTGGGCGAAAGGGTCTTCTTCGCGTACGTCGATCAgg                   <  2:113179/149‑1 (MQ=255)
                                                                                                                                  cccTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTATGTTGGCCGCGAACAGCGGATGGGCCTGGGCGAAAGGGTCTTCTTCGCGTACGTCGATCAGGGc                 <  8:19625/149‑1 (MQ=255)
                                                                                                                                   ccTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTATGTTGGCCGCGAACAGCGGATGGGCCTGGGCGAAAGGGTCTTCTTCGCGTACGTCGATCAGGGc                 >  1:223289/1‑148 (MQ=255)
                                                                                                                                         cGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGGATGTTGGCCGCGAACAGCGGATGGGCCTGGGCGAAAGGGGCTTCTTCGCGTACGTCGATCAGGGCCAGctcc          >  5:121691/1‑149 (MQ=255)
                                                                                                                                         cGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCACCTCCAGTTTCGACAGGGGTATGTTGGCCGCGAACAGCGGATGGGCCTGGGCGAAAGGGTCTTCTTCGCGTACGTCGATCAGGGCCAGctcc          >  3:219251/1‑149 (MQ=255)
                                                                                                                                          gTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTATGTTGGCCGCGAACAGCGGATGGGCCTGGGCGAAAGGGTCTTCTTCGCGTACGTCGATCAGGGCCAGctcct         >  3:314726/1‑149 (MQ=255)
                                                                                                                                               cGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTATGTTGGCCGCGAACAGCGGATGGGCCTGGGCGAAAGGGTCTTCTTCGCGTACGTCGATCTGGGCCAGGTGCTCCTTg    >  8:341735/1‑149 (MQ=255)
                                                                                                                                                gTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTATGTTGGCCGCGAACAGCGGATGGGCCTGGGCGAAAGGGTCTTCTTCGCGTACGTCGATCAGGGCCAGCTCCTCCTTgg   <  5:283658/149‑1 (MQ=255)
                                                                                                                                                 tAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCATTTTCGACAGGGGTATGTTGGCCGCGAACAGCGGATGGGCCTGGGCGAAAGGGTCTTCTTCGCGTACGTCGATCAGGGCCAGCTCCTCCTTGGc  <  8:278574/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
TTTGCTTGGCACGTTGACATCACGAAACAACTCGCCGCCCGCAGCGCGCCAGCCGGCAAGGCCGCCATCAAGCAAGGCAACATTGACGTAGCCCAAGGCCAGTAAACGCTGCGCGGCGACTGCGGCGAGCCCCTCACCGTCATCGTAGACGGTGAGTGCGGTGTCCCGCCGGGGAACGCGGGCGAAGATTTCCAGCTCCAGTTTCGACAGGGGTATGTTGGCCGCGAACAGCGGATGGGCCTGGGCGAAAGGGTCTTCTTCGCGTACGTCGATCAGGGCCAGCTCCTCCTTGGC  >  NC_002947/3649873‑3650166

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