Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 772,766 T→C 5.4% A142A (GCT→GCC PP_0664 → homoserine dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947772,7660TC 5.4% 340.7 / 3.9 109A142A (GCT→GCCPP_0664homoserine dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (47/56);  new base C (3/3);  total (50/59)
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.90e-01

Only 100 of 112 total aligned reads displayed.

GACATCATCGCCGAGGCCACATTCACCAACCCCAAGGACGGTGAGCCGGCGGCAACCTTCTGCCGCTGGGCGCTGGAGGCGGGCAAGCACGTGGTCACCACCAACAAGGGGCCGATTGCCTTGCATGCCCAGGCACTCAAGGCCTTGGCTAAGGCCAACAATGTGTCGTTCGAATACGAAGGCGCCGTCATGAGCGGCACACCGGTGCTGCGTATGGCGCGCCAGACC‑‑CTGGCGGGCGCGGGCTTGGTAGGTTTCGAGGGCATTCTCAACGGTACCTCTAACTATGTGCTCACGCGCAT  >  NC_002947/772617‑772915
                                                                                                                                                     |                                                                                                                                                        
  catcatCGCCGAGGCCACATTCACCAACCCCAAGGACGGTGAGCCGGCGGCAACCTTCTGCCGCTGGGCGCTGGAGGCGGGCAAGCACGTGGTCACCACCAACAAGGGGCCGATTGCCTTGCATGCCCAGGCACTCAAGGCCTTGGCTaa                                                                                                                                                        <  2:412013/150‑1 (MQ=255)
   atcatcGCCGAGGCCACATTCACCAACCCCAAGGACGGTGAGCCGGCGGCAACCTTCTGCCGCTGGGCGCTGGAGGCGGGCAAGCACGTGGTCACCACCAACAAGGGGCCGATTGCCTTGCATGCCCAGGCACTCAAGGCCTTGGCTAAg                                                                                                                                                       >  2:797172/1‑150 (MQ=255)
      atagCCGAGGCCACATTCACCAACCCCAAGGACGGTGAACCGGCGGCAACCTTCTGCCGATGGGCGCTGGAGGCGGGCAAGCACGTGGTCACCACCAACAAGGGGCCGATTGCCTTGCATGCCCAGGCCCTCAAGGCCTTGGCTAAGGcc                                                                                                                                                    >  2:133295/4‑150 (MQ=255)
            gAGGCCACATTCACCAAGCCCAAGGACGGTGAGCCGGCGGCAACCTTCTGCCGCTGGGCGCTGGAGGCGGGCAAGCACGTGGTCACCACCAACAAGGGGCCGATTGCCTTGCATGCCCAGGCACTCAAGGCCTTGGCTAAGGCCAACAAt                                                                                                                                              <  1:425105/150‑1 (MQ=255)
            gAGGCCACATTCACCAACCCCAAGGACGGTGAGCCGGCGGCAACCTTCTGCCGCTGGGCGCTGGAGGCGGGCAAGCACGTGGTCACCACCAACAAGGGGCCGATTGCCTTGCATGCCCAGGCACTCAAGGCCTTAGCCAAGGcctt                                                                                                                                                  <  2:912049/146‑3 (MQ=255)
            gAGGCCACATTCACCAACCCCAAGGACGGTGAGCCGGCGGCAACCTTCTGCCGCTGGGCGCTGGAGGCGGGCAAGCACGTGGTCACCACCAACAAGGGGCCGATTGCCTTGCATGCCCAGGCACTCAAGGCCTTAGCCAAGGcctt                                                                                                                                                  >  1:912049/1‑144 (MQ=255)
              ggCCACATTCACCAACCCCAAGGACGGTGAGCCGGCGGCAACCTTCTGCCGCTGGGCGCTGGAGGCGGGCAAGCACGTGGTCACCACCAACAAGGGGCCGATTGCCTTGCATGCCCAGGCACTCAAGGCCTTGGCTAAGGCCAACAAtgt                                                                                                                                            <  1:797172/150‑1 (MQ=255)
               gCCACATTCACCAACCCCAAGGACGGTGAGCCGGCGGCAACCTTCTGCCGCTGGGCGCTGGAGGCGGGCAAGCACGTGGTCACCACCAACAAGGGGCCGATTGCCTTGCATGCCCAGGCACTCAAGGCCTTGGCTAAGGCCAACAAtgtg                                                                                                                                           <  1:332093/150‑1 (MQ=255)
                 cacaTTCACCAACCCCAAGGACGGTGAGCCGGCGGCAACCTTCTGCCGCTGGGCGCTGGAGGCGGGCAAGCACGTGGTCACCACCAACAAGGGGCCGATTGCCTTGCATGCCCAGGCACTCAAGGCCTTGGCTAAGGCCAACAATGTGTc                                                                                                                                         >  2:1804587/1‑150 (MQ=255)
                     ttCACCAACCCCAAGGACGGTGAGCCGGCGGCAACCTTCTGCCGCTGGGCGCTGGAGGCGGGCAAGCACGTGGTCACCACCAACAAGGGGCCGATTGCCTTGCATGCCCAGGCACTCAAGGCCTTGGCTAAGGCCAACAATGTGTCGTTc                                                                                                                                     >  2:885971/1‑150 (MQ=255)
                        aCCAACCCCAAGGACGGTGAGCCGGCGGCAACCTTCTGCCGCTGGGCGCTGGAGGCGGGCAAGCACGTGGTCACCACCAACAAGGGGCCGATTGCCTTGCATGCCCAGGCACTCAAGGCCTTGGCTAAGGCCAACAATGTGTCGTTCGaa                                                                                                                                  <  2:2137881/150‑1 (MQ=255)
                        aCCAACCCCAAGGACGGTGAGCCGGCGGCAACCTTCTGCCGCTGGGCGCTGGAGGCGGGCAAGCACGTGGTCACCACCAACAAGGGGCCGATTGCCTTGCATGCCCAGGCACTCAAGGCCTTGGCTAAGGCCAACAATGTGTCGTTCGaa                                                                                                                                  <  1:1729460/150‑1 (MQ=255)
                        aCCAACCCCAAGGACGGTGAGCCGGCGGCAACCTTCTGCCGCTGGGCGCTGGAGGCGGGCAAGCACGTGGTCACCACCAACAAGGGGCCGATTGCCTTGCATGCCCAGGCACTCAAGGCCTTGGCTAAGGCCAACAATGTGTCGTTCGaa                                                                                                                                  >  1:2227042/1‑150 (MQ=255)
                         ccAACCCCAAGGACGGTGAGCCGGCGGCAACCTTCTGCCGCTGGGCGCTGGAGGCGGGCAAGCACGTGGTCACCACCAACAAGGGGCCGATTGCCTTGCATGCCCAGGCACTCAAGGCCTTTGCTAAGGCCAACAATGTGTCGTTCGAAt                                                                                                                                 <  2:7670/150‑1 (MQ=255)
                          cAACCCCAAGGACGGTGAGCCGGCGGCAACCTTCTGCCGCTGGGCGCTGGAGGCGGGCAAGCACGTGGTCACCACCAACAAGGGGCCGATTGCCTTGCATGCCCAGGCACTCAAGGCCTTGGCTAAGGCCAACAATGTGTCGTTCGaa                                                                                                                                  >  1:1678310/1‑148 (MQ=255)
                          cAACCCCAAGGACGGTGAGCCGGCGGCAACCTTCTGCCGCTGGGCGCTGGAGGCGGGCAAGCACGTGGTCACCACCAACAAGGGGCCGATTGCCTTGCATGCCCAGGCACTCAAGGCCTTGGCTAAGGCCAACAATGTGTCGTTCGaa                                                                                                                                  <  2:1678310/148‑1 (MQ=255)
                          cAACCCCAAGGACGGTGAGCCGGCGGCAACCTTCTGCCGCTGGGCGCTGGAGGCGGGCAAGCACGTGGTCACCACCAACAAGGGGCCGATTGCCTTGCATGCCCAGGCACTCAAGGCCTTGGCTAAGGCCAACAATGTGTCGTTCGAATa                                                                                                                                >  1:1063493/1‑150 (MQ=255)
                              cccAAGGACGGTGAGCCGGCGGCAACCTTCTGCCGCTGGGCGCTGGAGGCGGGCAAGCACGTGGTCACCACCAACAAGGGGCAGCTCGGTTTGCTGGACGAGGCGCTCAAGGCCTTGGCTCAGGCCAACAATGTGTCGTTCGAATACGaa                                                                                                                            >  1:2090122/1‑150 (MQ=255)
                               ccAAGGACGGTGAGCCGGCGGCAACCTTCTGCCGCTGGGCGCTGGAGGCGGGCAAGCACGTGGTCACCACCAACAAGGGGCCGATTGCCTTGCATGCCCAGGCACTCAAGGCCTTGGCTAAGGCCAACAATGTGTCGTTCGAATACGAAg                                                                                                                           <  1:109993/150‑1 (MQ=255)
                               ccAAGGACGGTGAGCCGGCGGCAACCTTCTGCCGCTGGGCGCTGGAGGCGGGCAAGCACGTGGTCACCACCAACAAGGGGCCGATTGCCCTGCATGCCGAGGCACTCAAGGCCTTGGCTAAGGCCAACAATGTGTCGTTCGAATACGAAg                                                                                                                           <  2:214323/150‑1 (MQ=255)
                               ccAAGGACGGTGAGCCGGCGGCAACCTTCTGCCGCTGGGCGCCGGAGGCGGGCAAGTACGGGGTCACCACCAACAAGGGGCCGATTGCCTTGCATGCCCAGGCACTCAAGGCCTTGGCTAAGGCCAACAATGTGTCGTTCGAATACGAAg                                                                                                                           >  1:1253335/1‑150 (MQ=255)
                                cAAGGACGGTGAGCCGGCGGCAACCTTCTGCCGCTGGGCGCTGGAGGCGGGCAAGCACGTGGTCACCACCAACAAGGGGCCGATTGCCTTGCATGCCCAGGCACTCAAGGCCTTGGCTAAGGCCAACAATGTGTCGTTCGAATACGAAgg                                                                                                                          >  2:1439012/1‑150 (MQ=255)
                                cAAGGACGGTGAGCCGGCGGCAACCTTCTGCCGCTGGGCGCTGGAGGCGGGCAAGCACGTGGTCACCACCAACAAGGGGCCGATTGCCTTGCATGCCCAGGCACTCAAGGCCTTGGCTAAGGCCAACAATGTGTCGTTCGAATACGAAgg                                                                                                                          >  1:1965988/1‑150 (MQ=255)
                                 aaGGACGGTGAGCCGGCGGCAACCTTCTGCCGCTGGGCGCTGGAGGCGGGCAAGCACGTGGTCACCACCAACAAGGGGCCGATTGCCTTGCATGCCCAGGCACTCAAGGCCTTGGCTAAGGCCAACAATGTGTCGTTCGAATACGAAGgc                                                                                                                         <  1:1728094/150‑1 (MQ=255)
                                       ggTGAGCCGGCGGCAACCTTCTGCCGCTGGGCGCTGGAGGCGGGCAAGCACGTGGTCACCACCAACAAGGGGCCGATTGCCTTGCATGCCCAGGCACTCAAGGCCTTGGCTAAGGCCAACAATGTGTCGTTCGAATACGAAGGCGCCGTc                                                                                                                   >  2:2047640/1‑150 (MQ=255)
                                       ggTGAGCCGGCGGCAACCTTCTGCCGCTGGGCGCTGGAGGCGGGCAAGCACGTGGTCACCACCAACAAGGGGCCGATTGCCTTGCATGCCCAGGCACTCAAGGCCTTGGCTAAGGCCAACAATGTGTCGTTCGAATACGAAGGCGCCGTc                                                                                                                   <  1:1366886/150‑1 (MQ=255)
                                            gCCGGCGGCAACCTTCTGCCGCTGGGCGCTGGAGGCGGGCAAGCACGTGGTCACCACCAACAAGGGGCCGATTGCCTTGCATGCCCAGGCACTCAAGGCCTTGGCTAAGGCCAACAATGTGTCGTTCGAATACGAAGGCGCCGTCATGAg                                                                                                              <  1:721661/150‑1 (MQ=255)
                                            gCCGGCGGCAACCTTCTGCCGCTGGGCGCGGGAGGCGGGCAAGCACGTGGTCACCACCAACAAGGGGCCGATTGCCTTGCATGCCCAGGCACTCAAGGCCTTGGCTAAGGCCAACAATGTGTCGTTCGAATACGAAGGCGCCGTCATGAg                                                                                                              >  1:1529736/1‑150 (MQ=255)
                                              ggccggcAACCTTCTGCCGCTGGGCGCTGGAGGCGGGCAAGCACGTGGTCACCACCAACAAGGGGCCGATTGCCTTGCATGCCCAGGCACTCAAGGACTTGGCTAAGGCCAACAATGTGTCGTTCGAATACGAAGGCGCCGTCATGAGCg                                                                                                            >  1:2269215/4‑150 (MQ=255)
                                                      aCCTTCTGCCGCTGGGCGCTGGAGGCGGGCAAGCACGTGGTCACCACCAACAAGGGGCCGATTGCCTTGCATGCCCAGGCACTCAAGGCCTTGGCTAAGGCCAACAATGTGTCGTTCGAATACGAAGGCGCCGTCATGAGCGGCACACCg                                                                                                    <  2:1755640/150‑1 (MQ=255)
                                                              ccGCTGGGCGCTGGAGGCGGGCAAGCACGTGGTCACCACCAACAAGGGGCCGATTGCCTTGCATGCCCAGGCACTCAAGGCCTTGACTAAGGCCAACAATGTGTCGTTCGAATACGAAGGCGCCGTCATGAGCGGCACACCGGTGCTGCg                                                                                            >  1:1804307/1‑150 (MQ=255)
                                                                gCTGGGCGCTGGAGGCGGGCAAGCACGTGGTCACCACCAACAAGGGGCCGATTGCCTTGCATGCCCAGGCACTCAAGGCCTTGGCTAAGGCCAACAATGTGTCGTTCGAATACGAAGGCGCCGTCATGAGCGGCACACCGGTGCTGCGTa                                                                                          <  2:1965988/150‑1 (MQ=255)
                                                                    ggCGCTGGAGGCGGGCAAGCACGTGGTCACCACCAACAAGGGGCCGATTGCCTTGCATGCCCAGGCACTCAAGGCCTTAGCCAAGGccttga                                                                                                                                                <  1:1739337/92‑5 (MQ=255)
                                                                    ggCGCTGGAGGCGGGCAAGCACGTGGTCACCACCAACAAGGGGCCGATTGCCTTGCATGCCCAGGCACTCAAGGCCTTAGCCAAGGccttga                                                                                                                                                >  2:1739337/1‑88 (MQ=255)
                                                                      cgcTGGAGGCGGGCAAGCACGTGGTCACCACCAACAAGGGGCCGATTGCCTTGCATGCCCAGGCACTCAAGGCCTTGGCTAAGGCCAACAATGTGTCGTTCGAATACGAAGGCGCCGTCATGAGCGGCACACCGGTGCTGCGTATGgcgc                                                                                    <  2:1458214/150‑1 (MQ=255)
                                                                        cTGGAGGCGGGCAAGCACGTGGTAACCACCAAAAAGTGGCCGATTGTCTTGCATGCCCAGGCACTCAAGGCCTTGGCTAAGGCCAACAATGTGTCGTTCGAATACGAAGGCGCCGTCATGAGCGGCACACCGGTGCTGCGTATGgcgcgc                                                                                  <  1:2338717/150‑1 (MQ=255)
                                                                                gggCAAGCACGTGGTCACCACCAACAAGGGGCCGATTGCCTTGCATGCCCAGGCACTCAAGGCCTTGGCTAAGGCCAACAATGTGTCGTTCGAATACGAAGGCGCCGTCATGAGCGGCACACCGGTGCTGCGTATGGCGCGCCAGACC‑‑Ct                                                                        >  2:1252898/1‑150 (MQ=255)
                                                                                  gCAAGCACGTGGTCACCACCAACAAGGGGCCGATTGCCTTGCATGCCCAGGCACTCAAGGCCTTGGCTAAGGCCAACAATGTGTCGTTCGAATACGAAGGCGCCGTCATGAGCGGCACACCGGTGCTGCGTATGGCGCGCCAGACC‑‑CTgg                                                                      >  2:2247624/1‑150 (MQ=255)
                                                                                   cAAGCACGTGGTCACCACCAACAAGGGGCCGATTGCCTTGCATGCCCAGGCACTCAAGGCCTTGGCTAAGGCCAACAATGTGTCGTTCGAATACGAAGGCGCCGTCATGAGCGGCACACCGGTGCTGCGTATGGCGCGCCAGACC‑‑CTggc                                                                     <  1:2116641/150‑1 (MQ=255)
                                                                                    aaGCACGTGGTCACCACCAACAAGGGGCCGATTGCCTTGCATGCCCAGGCACTCAAGGCCTTGGCTAAGGCCAACAATGTGTCGTTCGAATACGAAGGCGCCGTCATGAGCGGCACACCGGTGCTGCGTATGGCGCGCCAGAc                                                                             >  2:1413600/1‑143 (MQ=255)
                                                                                    aaGCACGTGGTCACCACCAACAAGGGGCCGATTGCCTTGCATGCCCAGGCACTCAAGGCCTTGGCTAAGGCCAACAATGTGTCGTTCGAATACGAAGGCGCCGTCATGAGCGGCACACCGGTGCTGCGTATGGCGCGCCAGAc                                                                             <  1:1413600/143‑1 (MQ=255)
                                                                                     aGCACGTGGTCACCACCAACAAGGGGCCGATTGCCTTGCATGCCCAGGCACTCAAGGCTTTGGCTAAGGCCAACAATGTGTCGTTCGAATACGAAGGCGCCGTCATGAGCGGCACACCGGTGCTGCGTATGGCGCGCCAGACC‑‑CTggcgg                                                                   <  1:826554/150‑1 (MQ=255)
                                                                                     aGCACGTGGTCACCACCAACAAGGGGCCGATTGCCTTGCATGCCCAGGCACTCAAGGCCTTGGCTAAGGCCAACAATGTGTCGTTCGAATACGAAGGCGCCGTCATGAGCGGCACACCGGTGCTGCGTATGGCGCGCCAGACC‑‑CTggcgg                                                                   >  1:1013817/1‑150 (MQ=255)
                                                                                      gCACGTGGTCACCACCAACAAGGGGCCGATTGCCTTGCATGCCCAGGCACTCAAGGCCTTGGCTAAGGCCAACAATGTGTCGTTCGAATACGAAGGCGCCGTCATGAGCGGCACACCGGTGCTGCGTATGGCGCGCCAGACC‑‑CTggcggg                                                                  >  1:764240/1‑150 (MQ=255)
                                                                                      gCACGTGGTCACCACCAACAAGGGGCCGATTGCCTTGCATGCCCAGGCACTCAAGGCCTTGGCTAAGGCCAACAATGTGTCGTTCGAATACGAAGGCGCCGTCATGAGCGGCACACCGGTGCTGCGTATGGCGCGCCAGACC‑‑CTggcggg                                                                  >  1:765306/1‑150 (MQ=255)
                                                                                       cACGTGGTCACCACCAACAAGGGGCCGATTGCCTTGCATGCCCAGGCACTCAAGGCCTTGGCTAAGGCCAACAATGTGTCGTTCGAATACGAAGGCGCCGGCATGAGCGGCACACCGGTGCTGCGTATGGCGCGCCAGACC‑‑CTggcgggc                                                                 <  1:62229/150‑1 (MQ=255)
                                                                                          gTGGTCACCACCAACAAGGGGCCGATTGCCTTGCATGCCCAGGCACTCAAGGCCTTGGCTAAGGCCAACAATGTGTCGTTCGAATACGAAGGCGCCGTCATGAGCGGCACACCgg                                                                                                   >  2:1909607/1‑115 (MQ=255)
                                                                                          gTGGTCACCACCAACAAGGGGCCGATTGCCTTGCATGCCCAGGCACTCAAGGCCTTGGCTAAGGCCAACAATGTGTCGTTCGAATACGAAGGCGCCGTCATGAGCGGCACACCgg                                                                                                   <  1:1909607/115‑1 (MQ=255)
                                                                                          gTGGTCACCACCAACAAGGGGCCGATTGCCTTGCATGCCCAGGCACTCAAGGCCTTGGCTAAGGCCAACAATGTGTCGTTCGAATACGAAGGCGCCGTCATGAGCGGCACACCGGTGCTGCGTATGGCGCGCCAGACC‑‑CTGgcgggcgcg                                                              >  1:1216520/1‑150 (MQ=255)
                                                                                            ggTCACCACCAACAAGGGGCCGATTGCCTTGCATGCCCAGGCACTCAAGGCCTTGGCTAAGGCCAACAATGTGTCGTTCGAATACGAAGGCGCCGTCATGAGCGGCACACCGGTGCTGCGTATGGCGCGCCAGACC‑‑Ct                                                                        <  2:1781275/138‑1 (MQ=255)
                                                                                            ggTCACCACCAACAAGGGGCCGATTGCCTTGCATGCCCAGGCACTCAAGGCCTTGGCTAAGGCCAACAATGTGTCGTTCGAATACGAAGGCGCCGTCATGAGCGGCACACCGGTGCTGCGTATGGCGCGCCAGACC‑‑Ct                                                                        >  1:1781275/1‑138 (MQ=255)
                                                                                             gTCACCACCAACAAGGGGCCGATTGCCTTGCATGCCCAGGCACTCAAGGCCTTGGCTAAGGCCAACAATGTGTCGTTCGAATACGAAGGCGCCGTCATGAGCGGCACACCGGTGCTGCGTATGGCGCGCCAGACC‑‑CTGgcgggcgcgggc                                                           >  2:16105/1‑150 (MQ=255)
                                                                                              tCACCACCAACAAGGGGCCGATTGCCTTGCATGCCCAGGCACTCAAGGCGTTGGCTAAGGACAACAATGTGTCGTTCGCATACGAAGGCACCGTAATGAGCGGCACACCGGTGCTGCGTATGGCGCGCCAGACC‑‑CTGGCGGGCGCGGGCt                                                          >  2:494514/1‑150 (MQ=255)
                                                                                               caccacCAACAAGGGGCCGATTGCCTTGCATGCCCAGGCACTCAAGGCCTTGGCTAAGGCCAACAATGTGTCGTTCGAATACGAAGGCGCCGTCATGAGCGGCACACCGGTGCTGCGTATGGCGCGCCAGACC‑‑CTGGCGGGCGCGGGCtt                                                         <  2:475114/150‑1 (MQ=255)
                                                                                                 ccaccaACAAGGGGCCGATTGCCTTGCATGCCCAGGCACTCAAGGCCTTGGCTAAGGCCAACAATGTGTCGTTCGAATACGAAGGCGCCGTCATGAGCGGCACACCGGTGCTGCGTATGGCGCGCCAGACC‑‑CTGGCGGGCGCGGGCTTgg                                                       >  1:1361880/1‑150 (MQ=255)
                                                                                                    ccaACAAGGGGCCGATTGCCTTGCATGCCCAGGCACTCAAGGCCTTGGCTAAGGCCAACAATGTGTCGTTCGAATACGAAGGCGCCGTCATGAGCGGCACACCGGTGCTGCGTATGGCGCGCCAGACC‑‑CTGGCGGGCGCGGGCTTGGTAg                                                    <  1:2047640/150‑1 (MQ=255)
                                                                                                      aacaaGGGGCCGATTGCCTTGCATGCCCAGGCACTCAAGGCCTTGGCTAAGGCCAACAATGTGTCGTTCGAATACGAAGGCGCCGTCATGAGCGGCACACCGGTGCTGCGTATGGCGCGCCAGACC‑‑CTGGCGGGCGCGGGCTTGGTAGGt                                                  <  2:2227042/150‑1 (MQ=255)
                                                                                                         aaGGGGCCGATTGCCTTGCATGCCCAGGCACTCAAGGCCTTGGCTAAGGCCAACAATGTGTCGTTCGAATACGAAGGCGCCGTCATGAGCGGCACACCGGTGCTGCGTATGGCGCGCCAGACC‑‑CTGGCGGGCGCGGGCTTGGTAGGTTTc                                               <  1:16105/150‑1 (MQ=255)
                                                                                                            ggtccGATTGCCTTGCATGCCCAGGCACTCAAGGCCTTGGCTAAGGCCAACAATGTGTCGTTCGAATACGAAGGCGCCGTCATGAGCGGCACACCGGTGCTGCGTATGGCGCGCCAGACC‑‑CTGGCGGGCGCGGGCTTGGTAGGTTTCGAg                                            <  2:1672634/147‑1 (MQ=255)
                                                                                                            gggCCGATTGCCTTGCATGCCCAGGCACTCAAGGCCTTGGCTAAGGCCAACAATGTGTCGTTCGAATACGAAGGCGCCGTCATGAGCgg                                                                                                           <  2:1936227/89‑1 (MQ=255)
                                                                                                            gggCCGATTGCCTTGCATGCCCAGGCACTCAAGGCCTTGGCTAAGGCCAACAATGTGTCGTTCGAATACGAAGGCGCCGTCATGAGCgg                                                                                                           >  1:1936227/1‑89 (MQ=255)
                                                                                                            gggCCGATTGCCTTGCATGCCCAGGCACTCAAGGCCTTGGCTAAGGCCAACAATGTGTCGTTCGAATACGAAGGCGCCGTCATGAGCGGCACACCGGTGCTGCGTATGGCGCGCCAGACC‑‑CTGGCGGGCGCGGGCTTGGTAGGTTTCGAg                                            >  2:1461059/1‑150 (MQ=255)
                                                                                                            cgagatcttGCCTTTTATTCCCAGTCACTCAAGGCCTTGGCTAAGGCCAACAATGTGTCGTTCGAATACGAAGGCGCCGTCATGAGCGGCACACCGGTGCTGCGTATGGCGCGCCAGACC‑‑CTGGCGGGCGCGGGCTTGGTAGGTTTCGAg                                            <  2:1118033/143‑1 (MQ=255)
                                                                                                             ggCCGATTGCCTTGCATGCCCAGGCACTCAAGGCCTTGGCTAAGGCCCACACTGTTTCTCTCGAATACGAAGGCGCCGTCATGAGCGGCACACCGCTGGTACTTATTGCTCGCCAGACC‑‑CTGGCGGGCGCGGGCTTGGTAGGTTTCGAgg                                          <  1:1252898/150‑1 (MQ=255)
                                                                                                             ggCCGATTGCCTTGCATGCCCAGGCACTCAAGGCCTTGGCTAAGGCCAACAATGTGTCGTTCGAATACGAAGGCGCCGTCATGAGCGGCACACCGGTGCTGCGTATGGCGCGCCAGACC‑‑CTGGCGGGCGCGGGCTTGGTAGGTTTCGAgg                                          <  2:1063493/150‑1 (MQ=255)
                                                                                                             ggCCGATTGCCTTGCATGCCCAGGCACTCAAGGCCTTGGCTAAGGCCAACAATGTGTCGTTCGAATACGAAGGCGCCGTCATGAGCGGCACACCGGTGCTGCGTATGGCGCGCCAGACC‑‑CTGGCGGGCGCGGGCTTGGTAGGTTTCGAgg                                          <  2:1452797/150‑1 (MQ=255)
                                                                                                             ggCCGATTGCCTTGCATGCCCAGGCACTCAAGGCCTTGGCTAAGGCCAACAATGTGTCGTTCGAATACGAAGGCGCCGTCATGAGCGGCACACCGGTGCTGCGTATGGCGCGCCAGACC‑‑CTGGCGGGCGCGGGCTTGGTAGGTTTCGAgg                                          <  1:1086501/150‑1 (MQ=255)
                                                                                                              gCCGAGTACCTGGCATGCCCAGGCACTCAAGGCCTTGGCTAAGGCCAACAATGTGTCGTTCGAATACGAAGGCGCCGTCATGAGCGGCACACCGGTGCTGCGTATGGCGCGCCAGACC‑‑CTGGCGGGCGCGGGCTTGGTAGGTTTCGAggg                                         >  2:1902734/1‑150 (MQ=255)
                                                                                                                cGATTGCCTTGCATGCCCAGGCACTCAAGGCCTTGGCTAAGGCCAACAATGTGTCGTTCGAATACGAAGGCGCCGTCATGAGCGGCACACCGGTGCTGCGTATGGCGCGCCAGACC‑‑CTGGCGGGCGCGGGCTTGGTAGGTTTCGAGGGCa                                       <  2:2128626/150‑1 (MQ=255)
                                                                                                                cGATTGCCTTGCATGCCCAGGCACTCAAGGCCTTGGCTAAGGCCAACAATGTGTCGTTCGAATACGAAGGCGCCGTCATGAGCGGCACACCGGTGCTGCGTATGGCGCGCCAGACC‑‑CTGGCGGGCGCGGGCTTGGTAGGTTTCGAGGGCa                                       >  2:1881921/1‑150 (MQ=255)
                                                                                                                 gATTGCCTTGCATGCCCAGGCACTCAAGGCCTTGGCTAAGGCCAACAATGTGTCGTTCGAATACGAAGGCGCCGTCATGAGCGGCACACCGGTGCTGCGTATGGCGCGCCAGACC‑‑CTGGCGGGCGCGGGCTTGGTAGGTTTCGAGGGCAt                                      >  2:1051941/1‑150 (MQ=255)
                                                                                                                 gATTGCCTTGCATGCCCAGGCACTCAAGGCCTTGGCTAAGGCCAACAATGTGTCGTTCGAATACGAAGGCGCCGTCATGAGCGGCACACCGGTGCTGCGTATGGCGCGCCAGACC‑‑CTGGCGGGCGCGGGCTTGGTAGGTTTCGAGGGCAt                                      >  2:2011963/1‑150 (MQ=255)
                                                                                                                   ttGCCTTGCATGCCCAGGCACTCAAGGCCTTGGCTAAGGCCAACAATGTGTCGTTCGAATACGAAGGCGCCGTCATGAGCGGCACACCGGTGCTGCGTATGGCGCGCCAGACC‑‑CTGGCGGGCGCGGGCTTGGTAGGTTTCGAGGGCATtc                                    <  1:885971/150‑1 (MQ=255)
                                                                                                                       cTTGCATGCCCAGGCACTCAAGGCCTTGGCTAAGGCCAACAATGTGTCGTTCGAATACGAAGGCGCCGTCATGAGCGGCACACCGGTGCTGCGTATGGCGCGCCAGACC‑‑CTGGCGGGCGCGGGCTTGGTAGGTTTCGAGGGCATTCTCaa                                <  1:1804587/150‑1 (MQ=255)
                                                                                                                       cTTGCATGCCCAGGCACTCAAGGCCTTGGCTAAGGCCAACAATGTGTCGTTCGAATACGAAGGCGCCGTCATGAGCGGCACACCGGTGCTGCGTATGGCGCGCCAGACC‑‑CTGGCGGGCGCGGGCTTGGTAGGTTTCGAGGGCATTCTCaa                                <  1:1439012/150‑1 (MQ=255)
                                                                                                                       cTTGCATGCCCAGGCACTCAAGGCCTTGACTAAGGCCAACAATGTGTCGTTCGAATACGAAGGCGCCGTCATGAGCGGCACACCGGTGCTGCGTATGGCGCGCCAGACC‑‑CTGGCGGGCGCGGGCTTGGTAGGTTTCGAGGGCATTCTCaa                                <  2:1804307/150‑1 (MQ=255)
                                                                                                                          gCATGCCCAGGCACTCAAGGCCTTGGCTAAGGCCAACAATGTGTCGTTCGAATACGAAGGCGCCGTCATGa                                                                                                               >  2:1853553/1‑71 (MQ=255)
                                                                                                                          gCATGCCCAGGCACTCAAGGCCTTGGCTAAGGCCAACAATGTGTCGTTCGAATACGAAGGCGCCGTCATGa                                                                                                               <  1:1853553/71‑1 (MQ=255)
                                                                                                                             tGCCCAGGCACTCAAGGCCTTGGCTAAGGCCAACAATGTGTCGTTCGAATACGAAGGCGCCGTCATGAGCGGCACACCGGTGCTGCGTATGGCGCGCCAGACC‑‑CTGGCGGGCGCGGGCTTGGTAGGTTTCGAGGGCATTCTCAACGGTAc                          <  1:1475815/150‑1 (MQ=255)
                                                                                                                              gCCCAGGCACTCAAGGCCTTGGCTAAGGCCAACAATGTGTCGTTCGAATACGAAGGCGCCGTCATGAGCGGCACACCGGTGCTGCGTATGGCGCGCCAGACC‑‑CTGGCGGGCGCGGGCTTGGTAGGTTTCGAGGGCATTCTCAACGGTAcc                         >  1:339019/1‑150 (MQ=255)
                                                                                                                                 cAGGCACTCAAGGCCTTGGCTAAGGCCAACAATGTGTCGTTCGAATACGAAGGCGCCGTCATGAGCGGCACACCGGTGCTGCGTATGGCGCGCCAGACC‑‑CTGGCGGGCGCGGGCTTGGTAGGTTTCGAgg                                          <  2:1680623/130‑1 (MQ=255)
                                                                                                                                 cAGGCACTCAAGGCCTTGGCTAAGGCCAACAATGTGTCGTTCGAATACGAAGGCGCCGTCATGAGCGGCACACCGGTGCTGCGTATGGCGCGCCAGACC‑‑CTGGCGGGCGCGGGCTTGGTAGGTTTCGAgg                                          >  1:1680623/1‑130 (MQ=255)
                                                                                                                                   ggCAGTAAAGGTCTTGGCTAAGGCCAACAATGTGTCGTTCGAATACGAAGGCGCCGTCATGAGCGGCACACCGGTGCTGCGTATGGCGCGCCAGACC‑‑CTGGCGGGCGCGGGCTTGGTAGGTTTCGAGGGCATTCTCAACGGTACCTCTaa                    <  2:1878093/150‑1 (MQ=255)
                                                                                                                                   ggCACTCAAGGCCTTGGCTAAGGCCAACAATGTGTCGTTCGAATACGAAGGCGCCGTCATGAGCGGCACACCGGTGCTGCGTATGGCGCGCCAGACC‑‑CTGGCGGGCGCGGGCTTGGTAGGTTTCGAGGGCATTCTCAACGGTACCTCTaa                    <  1:987369/150‑1 (MQ=255)
                                                                                                                                       cTCAAGGCCTTGGCTAAGGCCAACAATGTGTCGTTCGAATACGAAGGCGCCGTCATGAGCGGCACACCGGTGCTGCGTATg                                                                                        <  2:2135799/81‑1 (MQ=255)
                                                                                                                                       cTCAAGGCCTTGGCTAAGGCCAACAATGTGTCGTTCGAATACGAAGGCGCCGTCATGAGCGGCACACCGGTGCTGCGTATg                                                                                        >  1:2135799/1‑81 (MQ=255)
                                                                                                                                       cTCAAGGCCATGGATAAGGCCAACAAGGTGGCGTTCGAATACGAAGGCGCCGTCATGAGGGGCACACCGGTGCTGCGTATGGCGCGCCAGACC‑‑CTGGCGGGCGCGGGCTTGGTAGTTTTCGAGGGCATTCTCAACGGTACCTCTAACTAt                <  2:627617/150‑1 (MQ=255)
                                                                                                                                        tCAAGGCCTTGGCTAAGGCCAACAATGTGTCGTTCGAATACGAAGGCGCCGTCATGa                                                                                                               <  2:2187102/57‑1 (MQ=255)
                                                                                                                                        tCAAGGCCTTGGCTAAGGCCAACAATGTGTCGTTCGAATACGAAGGCGCCGTCATGa                                                                                                               >  1:2187102/1‑57 (MQ=255)
                                                                                                                                          ttggCCTTAGCCAAGGCCAACAATGTGTCGTTCGAATACGAAGGCGCCGTCATGAGCGGCACACCGGTGCTGCGTATGGCGCGCCAGACC‑‑CTGGCGGGCGCGGGCTTGGTAGGTTTCGAgg                                          <  2:1105155/119‑1 (MQ=255)
                                                                                                                                          ttggCCTTAGCCAAGGCCAACAATGTGTCGTTCGAATACGAAGGCGCCGTCATGAGCGGCACACCGGTGCTGCGTATGGCGCGCCAGACC‑‑CTGGCGGGCGCGGGCTTGGTAGGTTTCGAgg                                          >  1:1105155/3‑121 (MQ=255)
                                                                                                                                          aaGGCCTTGGCTAAGGCCAACAATGTGTCGTTCGAATACGAAGGCGCCGTCATGAGCGGCACACCGGTGCTGCGTATGGCGCGCCAGACC‑‑CTGGCGGGCGCGGGCTTGGTAGGTTTCGAGGGCATTCTCAACGGTACCTCTAACTAtgtg             >  2:1003354/1‑150 (MQ=255)
                                                                                                                                          aaGGCCTTGGCTAAGGCCAACAATGTGTCGTTCGAATACGAAGGCGCCGTCATGAGCGGCACACCGGTGCTGCGTATGGCGCGCCAGACC‑‑CTGGCGGGCGCGGGCTTGGTAGGTTTCGAGGGCATTCTCAACGGTACCTCTAACTAtgtg             >  2:826013/1‑150 (MQ=255)
                                                                                                                                              ccTTGGCTAAGGCCAACAATGTGTCGTTCGAATACGAAGGCGCCGTCATGAGCGGCACACCGGTGCTGCGTATGGCGCGCCAGACC‑‑CTGGCGGGCGCGGGCTTGGTAGGTTTCGAgg                                          <  2:1441052/117‑1 (MQ=255)
                                                                                                                                              ccTTGGCTAAGGCCAACAATGTGTCGTTCGAATACGAAGGCGCCGTCATGAGCGGCACACCGGTGCTGCGTATGGCGCGCCAGACC‑‑CTGGCGGGCGCGGGCTTGGTAGGTTTCGAgg                                          >  1:1441052/1‑117 (MQ=255)
                                                                                                                                              ccTTGGCTAAGGCCAACAATGTGTCGTTCGAATACGAAGGCGCCGTCATGAGCGGCACACCGGTGCTGCGTATGGCGCGCCAGACC‑‑CTGGCGGGCGCGGGCTTGGTAGGTTTCGAGGGCATTCTCAACGGTACCTCTAACTATGTGCTCa         >  1:1860978/1‑150 (MQ=255)
                                                                                                                                                ttGGCTAAGGCCAACAATGTGTCGTTCGAATACGAAGGCGCCGTCATGAGCGGCACACCGGTGCTGCGTATG‑‑GCGCAAGGCCTGCAGTCGGGCGCGGGCTTGGTAGGTTTAGAGCGGAACTCCAACCGTACCGTGAACTATGTGCTCAcg       >  1:1853676/1‑150 (MQ=255)
                                                                                                                                                    cTAAGGCCAACAATGTGTCGTTCGAATACGAAGGCGCCGTCATGAGCGGCACACCGGTGCTGCGTATGGCGCGCCAGACC‑‑CTGGCGGGCGCGGGCTTGGTAGGTTTCGAGGGCATTCTCAACGGTACCTCTAACTATGTGCTCACGCGCa   >  2:1203433/1‑150 (MQ=255)
                                                                                                                                                     tAAGGCCAACAATGTGTCGTTCGAATACGAAGGCGCCGTCATGAGCGGCACACCGGTGCTGCGTATGGCGCGCCAGACC‑‑CTGGCGGGCGCGGGCTTGGTAGGTTTCGAGGGCATTCTCAACGGTACCTCTAACTATGTGCTCACGCGCAt  >  2:1591929/1‑150 (MQ=255)
                                                                                                                                                     tAAGGCCAACAATGTGTCGTTCGAATACGAAGGCGCCGTCATGAGCGGCACACCGGTGCTGCGTATGGCGCGCCAGACC‑‑CTGGCGGGCGCGGGCTTGGTAGGTTTCGAGGGCATTCTCAACGGTACCTCTAACTATGTGCTCACGCGCAt  <  1:1051941/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                        
GACATCATCGCCGAGGCCACATTCACCAACCCCAAGGACGGTGAGCCGGCGGCAACCTTCTGCCGCTGGGCGCTGGAGGCGGGCAAGCACGTGGTCACCACCAACAAGGGGCCGATTGCCTTGCATGCCCAGGCACTCAAGGCCTTGGCTAAGGCCAACAATGTGTCGTTCGAATACGAAGGCGCCGTCATGAGCGGCACACCGGTGCTGCGTATGGCGCGCCAGACC‑‑CTGGCGGGCGCGGGCTTGGTAGGTTTCGAGGGCATTCTCAACGGTACCTCTAACTATGTGCTCACGCGCAT  >  NC_002947/772617‑772915

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