Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,201,878 C→T 7.5% A46V (GCC→GTC)  PP_4583 → putative Peptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,201,8780CT 7.5% 306.0 / 10.7 107A46V (GCC→GTC) PP_4583putative Peptidase
Reads supporting (aligned to +/- strand):  ref base C (47/50);  new base T (3/5);  total (51/56)
Fisher's exact test for biased strand distribution p-value = 7.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

Only 100 of 108 total aligned reads displayed.

AAGAGACCTTCAATGCCAACCAAGCCCCAACCCCCGATTGCCCAAGCCGACCAGGCCACCGACCCGTATGCCTGGCTGCAACAGCGCGACACCCCCGAGGTACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGG‑CCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTGGAGACCGACCTGTCATTGCCGGCGCCCTGGGGCCCCTACCTGTACTACACCCGTACCACCGCCGGCGACGAATATCCGCGTCACTACCGCTGC  >  NC_002947/5201730‑5202026
                                                                                                                                                     |                                                                                                                                                    
aaGAGACCTTCAATGCCAACCAAGCCCCAACCCCCGATTGCCCAAGCCGACCAGGCCACCGACCCGTATGCCTGGCTGCAACAGCGCGACACCCCCGAGGTACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGG‑cc                                                                                                                                                     >  2:1451938/1‑150 (MQ=255)
     aCCTTCAATGCCAACCAAGCCCCAACCCCCGATTGCCCAAGCCGACCAGGCCATCGACCCGTATGCCTGGCTGCAACAGCGCGACACCCCCGAGGTACTCGCTTACCTGTAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGG‑CCGACCa                                                                                                                                                <  1:2298915/150‑1 (MQ=255)
     aCCTTCAATGCCAACCAAGCCCCAACCCCCGATTGCCCAAGCCGACCAGGCCACCGACCCGTATGCCTGGCTGCAACAGCGCGACACCCCCGAGGTACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGG‑CCGACCa                                                                                                                                                >  2:928392/1‑150 (MQ=255)
       cTTCAATGCCAACCAAGCCCCAACCCCCGATTGCCCAAGCCGACCAGGCCACCGACCCGTATGCCTGGCTGCAACAGCGCGACACCCCCGAGGTACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGG‑CCGACCAgg                                                                                                                                              >  2:392988/1‑150 (MQ=255)
        ttCAATGCCAACCAAGCCCCAACCCCCGATTGCCCAAGCCGACCAGGCCACCGACCCGTATGCCTGGCTGCAACAGCGCGACACCCCCGAGGTACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGG‑TCGGCCAGGc                                                                                                                                             <  1:908121/150‑1 (MQ=255)
        ttCAATGCCAACCAAGCCCCAACCCCCGATTGCCCAAGCCGACCAGGCCACCGACCCGTATGCCTGGCTGCAACAGCGCGACACCCCCGAGGTACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGG‑CCGACCAGGc                                                                                                                                             >  1:1193864/1‑150 (MQ=255)
         tCAATGCCAACCAAGCCCCAACCCCCGATTGCCCAAGCCGACCAGGCCACCGACCCGTATGCCTGGCTGCAACAGCGCGACACCCCCGAGGTACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGG‑CCGACCAGGCa                                                                                                                                            <  1:1016443/150‑1 (MQ=255)
          cAATGCCAACCAAGCCCCAACCCCCGATTGCCCAAGCCGACCAGGCCACCGACCCGTATGCCTGGCTGCAACAGCGCGACACCCCCGAGGTACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGG‑TCGGCCAGGCag                                                                                                                                           <  2:1256818/150‑2 (MQ=255)
              gCCAACCAAGCCCCAACCCCCGATTGCCCAAGCCGACCAGGCCACCGTCCCGTATGCCTGGCTGCAACAGCGCGACACCCCCGAGGTACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGG‑CCGACCAGGCACCGtt                                                                                                                                       <  1:1012276/150‑1 (MQ=255)
               ccaaccaaGGCCCCACCCCCGATTGCCCACGCCGACCAGGCCACCGACCCGTATACCTGGCTGCAACAGCGCGACACCCCCGAGGTACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGG‑CCGACCAGGCACCGTTg                                                                                                                                      >  2:2060997/1‑150 (MQ=255)
                caaccaaGCCCCAACCCCCGATTGCCCAAGCCGACCAGGCCACCGACCCGTATGCCTGGCTGCAACAGCGCGACACCCCCGAGGTACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGG‑TCGGCCAGGCaggcttcc                                                                                                                                     >  1:2348926/1‑143 (MQ=255)
                   ccaaGCCCCAACCCCCGATTGCCCAAGCAGACCAGGCCACCGACCCGTATGCCTAGCTTCAACAGCGCGACACCCCCGAGGTACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGG‑CCGACCAGGCACCGTTGCGTg                                                                                                                                  >  1:1501379/1‑150 (MQ=255)
                    caaGCCCCAACCCCCGATTGCCCAAGCCGACCAGGCCACCGACCCGTATGCCTGGCTGCAACAGCGCGACACCCCCGAGGTACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGG‑CCGACCAGGCACCGTTGCGTGa                                                                                                                                 >  1:1883246/1‑150 (MQ=255)
                    caaGCCCCAACCCCCGATTGCCCAAGCCGACCAGGCCACCGACCCGTATGCCTGGCTGCAACAGCGCGACACCCCCGAGGTACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGG‑CCGACCAGGCACCGTTGCGTGa                                                                                                                                 >  1:681421/1‑150 (MQ=255)
                    caaGCCCCAACCCCCGATTGCCCAAGCCGACCAGGCCACCGACCCGTATGCCTGGCTGCAACAGCGCGACACCCCCGAGGTACTCGCTTACCTGCAAACTGAAAACGCCTACCAGGAAGCCTGCCTGG‑CCGACCAGGCACCGTTGCGTGa                                                                                                                                 >  2:804511/1‑150 (MQ=255)
                         cccAACCCCCGATTGCCCAAGCCGACCAGGCCACCGACCCGTATGCCTGGCTGCAACAGCGCGACACCCCCGAGGTACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGG‑CCGACCAGGCACCGTTGCGTGAGCGcc                                                                                                                            >  2:1807665/1‑150 (MQ=255)
                              cccccGATTGCCCAAGCCGACCAGGCCACCGACCCGTATGCCTGGCTGCAACAGCGCGACACCCCCGAGGTACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGG‑CCGACCAGGCACCGTTGCGTGAGCGCCTGTTc                                                                                                                       >  1:1200187/1‑150 (MQ=255)
                              cccccGATTGCCCAAGCCGACCAGGCCACCGACCCGTATGCCTGGCTGCAACAGCGCGACACCCCCGAGGTACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGG‑CCGACCAGGCACCGTTGCGTGAGCGCCTGTTc                                                                                                                       <  2:2025363/150‑1 (MQ=255)
                                cccGATTGCCCAAGCCGACCAGGCCACCGACCCGTATGCCTGGCTGCAACAGCGCGACACCCCCGAGGTACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGG‑CCGACCAGGCACCGTTGCGTGAGCGCCTGTTCga                                                                                                                     >  2:1425826/1‑150 (MQ=255)
                                 ccGATTGCCCAAGCCGACCAGGCCACCGACCCGTATGCCTGGCTGCAACAGCGCGACACCCCCGAGGTACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGG‑CCGACCAGGCACCGTTGCGTGAGCGCCTGTTCgaa                                                                                                                    <  1:392988/150‑1 (MQ=255)
                                     ttGCCCAAGCCGACCAGGCCACCGACCCGTATGCCTGGCTGCAACAGCGCGACACCCCCGAGGTACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGG‑CCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGaaa                                                                                                                <  2:681421/150‑1 (MQ=255)
                                        cccAAGCCGACCAGGCCACCGACCCGTATGTCTGCCTGCAACTGCGCGACACCCCCGAGGTACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGAAAGCCTGCCTGG‑CCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCa                                                                                                             <  1:376074/150‑1 (MQ=255)
                                         ccAAGCCGACCAGGCCACCGACCCGTATGCCTGGCTGCAACAGCGCGACACCCCCGAGGTACTCGCTTACCTGCAAGCTGAAAACGTCTCCCAGAAGGCCTCCCTGG‑CCGACCAGGCACCGTTGCGCGAGCGCCTGTTCGACAAAATCaa                                                                                                            >  1:1015184/1‑150 (MQ=255)
                                          cAAGCCGACCAGGCCACCGACCCGTATGCCTGGCTGCAACAGCGCGACACCCCCGAGGTACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGG‑CCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAg                                                                                                           <  1:2321038/150‑1 (MQ=255)
                                          cAAGCCGACCAGGCCACCGACCCGTATGCCTGGCTGCAACAGCGCGACACCCCCGAGGTACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGG‑CCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAg                                                                                                           <  1:501373/150‑1 (MQ=255)
                                           aaGCCGACCAGGCCACCGACCCGTATGCCTGGCTGCAACAGCGCGACACCCCCGAGGTACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGG‑CCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAgg                                                                                                          <  2:1200187/150‑1 (MQ=255)
                                             gCCGACCAGGCCACCGACCCGTATGCCTGGCTGCAACAGCGCGACACCCCCGAGGTACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGG‑CCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGc                                                                                                        >  2:1254260/1‑150 (MQ=255)
                                                gACCAGGCCACCGACCCGTATGCCTGGCTGCAACAGCGCGACACCCCCGAGGTACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGG‑CCGACCAGGCACCGTTGCGTGAGCGCCTGt                                                                                                                         >  1:1473015/1‑130 (MQ=255)
                                                gACCAGGCCACCGACCCGTATGCCTGGCTGCAACAGCGCGACACCCCCGAGGTACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGG‑CCGACCAGGCACCGTTGCGTGAGCGCCTGt                                                                                                                         <  2:1473015/130‑1 (MQ=255)
                                                 aCCAGGCCACCGACCCGTATGCCTGGCTGCAACAGCGCGACACCCCCGAGGTACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGG‑CCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCa                                                                                                    >  1:201509/1‑150 (MQ=255)
                                                  ccAGGCCACCGACCCGTATGCCTGGCTGCAACAGCGCGACACCCCCGAGGTACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGG‑CCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCAt                                                                                                   >  2:1837142/1‑150 (MQ=255)
                                                        cACCGACCCGTATGCCTGGCTGCAACAGCGCGACGCCCCCGAGGTACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGG‑CCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTGga                                                                                             <  1:1342322/150‑1 (MQ=255)
                                                        cACCGACCCGTATGCCTGGCTGCAACAGCGCGACACCCCCGAGGTACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGG‑CCGACCAgg                                                                                                                                              <  1:2002682/101‑1 (MQ=255)
                                                         aCCGACCCGTATGCCTGGCTGCAACAGCGCGACACCCCCGAGGTACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCATGAATGA‑ACTGCGCCACACCGGTGCGTTAGCGCATGTTCGAAGAAATCGAGGGCCGCATCGTGgag                                                                                            >  1:2089667/1‑150 (MQ=255)
                                                               ccGTATGCCTGGCTGCAACAGCGCGACACCCCCGAGGTACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGG‑CCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTGGAgaccgac                                                                                      <  2:1193864/150‑1 (MQ=255)
                                                                      ccTGGCTGCAACAGCGCGACACCCCCGAGGTACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGG‑CCGACCAGGCACCGTTGCGTGAGCGCCTGTTc                                                                                                                       <  1:1603223/110‑1 (MQ=255)
                                                                      ccTGGCTGCAACAGCGCGACACCCCCGAGGTACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGG‑CCGACCAGGCACCGTTGCGTGAGCGCCTGTTc                                                                                                                       >  2:1603223/1‑110 (MQ=255)
                                                                      ccTGGCTGCAACAGCGCGACACCCCCGAGGTACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGG‑CCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTGGAGACCGACCTGTCAt                                                                               >  2:1477645/1‑150 (MQ=255)
                                                                       cTGGCTGCAACAGCGCGACACCCCCGAGGTACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGG‑TCGGCCAGGCag                                                                                                                                           <  2:492132/89‑2 (MQ=255)
                                                                       cTGGCTGCAACAGCGCGACACCCCCGAGGTACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGG‑TCGGCCAGGCag                                                                                                                                           >  1:492132/1‑88 (MQ=255)
                                                                       cTGGCTGCAACAGCGCGACACCCCCGAGGTACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGG‑CCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTGGAGACCGACCTGTCAtt                                                                              >  2:1208191/1‑150 (MQ=255)
                                                                       cTGGCTGCAACAGCGCGACACCCCCGAGGTACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGG‑CCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTGGAGACCGACCTGTCAtt                                                                              <  2:525837/150‑1 (MQ=255)
                                                                       cTGGCTGCAACAGCGCGACACCCCCGAGGTACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGG‑CCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTGGAGACCGACCTGTCAtt                                                                              >  2:592276/1‑150 (MQ=255)
                                                                        tGGCTGCAACAGCGCGACACCCCCGAGGTACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGG‑CCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTGGAGACCGACCt                                                                                    >  1:2147448/1‑143 (MQ=255)
                                                                        tGGCTGCAACAGCGCGACACCCCCGAGGTACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGG‑CCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTGGAGACCGACCt                                                                                    <  2:2147448/143‑1 (MQ=255)
                                                                        tGGCTGCAACAGCGCGACACCCCCGAGGTACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGG‑CCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTGGAGACCGACCTGTCATTg                                                                             >  1:257668/1‑150 (MQ=255)
                                                                           cTGCAACAGCGCGACACCCCCGAGGTACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGG‑CCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTGGAGACCGACCTGTCATTGCCg                                                                          <  2:2064915/150‑1 (MQ=255)
                                                                                    cgcgACACCCCCGAGGTACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGG‑CCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTGGAGACCGACCTGTCATTGCCGGCGCCCTgg                                                                 <  1:2060997/150‑1 (MQ=255)
                                                                                    cgcgACACCCCCGAGGTACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGG‑CCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTGGAGACCGACCTGTCATTGCCGGCGCCCTgg                                                                 >  2:1147071/1‑150 (MQ=255)
                                                                                     gcgACACCCCCGAGGTACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGG‑CCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTGGAGACCGACCTGTCATTGCCGGCGCCCTggg                                                                <  1:592276/150‑1 (MQ=255)
                                                                                       gACACCCCCGAGGTACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGG‑CCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCNAGGGCCNCATCCTGGAGACNGACCTGTCATTGCCGGCGCCCTGGGGc                                                              <  1:1596369/150‑1 (MQ=255)
                                                                                       gACACCCCCGAGGTACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGG‑CCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTGGAGACCGACCTGTCATTGCCGGCGCCCTGGGGc                                                              <  1:1837142/150‑1 (MQ=255)
                                                                                       cacaaCCCCAAGGGTCTCGCCTACCGGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGGCCCCATCTGG‑ACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTGGAGACCGACCTGTCATTGCCGGCGCCCTGGGGc                                                              <  1:760994/149‑1 (MQ=255)
                                                                                            ccccGAGGTACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCGTGG‑CCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGCCCGCATCCTGGAGACCGACCGGTCATTGCCGGCGCCCTGGGGCCCCTa                                                         >  2:926666/1‑150 (MQ=255)
                                                                                            ccccGAGGTACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGG‑CCGACCAGGCACCGTTGCGt                                                                                                                                   >  1:980173/1‑76 (MQ=255)
                                                                                            ccccGAGGTACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGG‑CCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTGGAgac                                                                                          >  1:105436/1‑117 (MQ=255)
                                                                                            ccccGAGGTACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGG‑CCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTGGAgac                                                                                          <  2:105436/117‑1 (MQ=255)
                                                                                            ccccGAGGTACCCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGG‑CCGACCAGGCACCGTTGCGTGAgc                                                                                                                               <  2:742923/80‑1 (MQ=255)
                                                                                            ccccGAGGTACCCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGG‑CCGACCAGGCACCGTTGCGTGAgc                                                                                                                               >  1:742923/1‑80 (MQ=255)
                                                                                             cccGAGGTACTCGCTTACCTGCAATCTGAAAACGCCTACCAGGAAGCCTGCCTGG‑CCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTGGAGACCGACCTGTCATTGCCGGCGCCCTGGGGCCCCTAc                                                        <  2:201509/150‑1 (MQ=255)
                                                                                             cccGAGGTACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGG‑CCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTGGAGACCGACCTGTCATTGCCGGCGCCCTGGGGCCCCTAc                                                        <  1:329111/150‑1 (MQ=255)
                                                                                             cccGAGGTACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGG‑CCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTGGAGACCGACCTGTCATTGCCGGCGCCCTGGGGCCCCTAc                                                        >  1:1574718/1‑150 (MQ=255)
                                                                                               cGAGGTACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGG‑CCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTGGAGACCGACCTGTCATTGCCGGCGCCCTGGGGCCCCTACCt                                                      <  2:257668/150‑1 (MQ=255)
                                                                                                  ggTACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGG‑CCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTGGAGACCGACCTGTCATTGCCGGCGCCCTGGGGCCCCTACCTGta                                                   <  1:387573/150‑1 (MQ=255)
                                                                                                    tACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGG‑ACTTCCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTGGAGACCGACCTGTCATTGCCGGCGCCCTGGGGCCCCTACCTGtact                                                 <  1:1132840/150‑1 (MQ=255)
                                                                                                        cGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGG‑CCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTGGAGACCGACCTGTCATTGCCGGCGCCCTGGGGCCCCTACCTGTACTacac                                             <  1:1477645/150‑1 (MQ=255)
                                                                                                         gCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGG‑CCGACCAGGCACCGTTGc                                                                                                                                     <  1:2311837/61‑1 (MQ=255)
                                                                                                         gCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGG‑CCGACCAGGCACCGTTGc                                                                                                                                     >  2:2311837/1‑61 (MQ=255)
                                                                                                         gCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGG‑CCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTGGAGACCGACCTGTCATTGCCGGCGCCCTGGGGACCATAACTGTACTACAcc                                            >  2:2067560/1‑150 (MQ=255)
                                                                                                           ttACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGG‑CCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTGGAGACCGACCTGTCATTGCCGGCGCCCTGGGGCCCCTACCTGTACTACACCCg                                          >  2:1318121/1‑150 (MQ=255)
                                                                                                               cTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGG‑CCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTGGAGACCGACCTGTCATTGCCGGCGCCCTGGGGCCCCTACCTGTACTACACCCGTacc                                      >  1:792422/1‑150 (MQ=255)
                                                                                                               cTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGG‑CCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTGGAGACCGACCTGTCATTGCCGGCGCCCTGGGGCCCCTACCTGTACTACACCCGTacc                                      >  2:948305/1‑150 (MQ=255)
                                                                                                                 gCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGG‑CCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTGGAGACCGACCTGTCATTGCCGGCGCCCTGGGGCCCCTACCTGTACTACACCCGTaccac                                    >  1:1785619/1‑150 (MQ=255)
                                                                                                                   aaGCTGAAAACGCCTACCAGGAAGCCTGCCTGG‑CCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTGGAGACCGACCTGTCATTGCCGGCGCCCTGGGGCCCCTACCTGTACTACACCCGTACCAccg                                  >  2:2086678/1‑150 (MQ=255)
                                                                                                                     gCTGAAAACGCCTACCAGGAAGCCTGCCTGG‑CCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTGGAGACCGACCTGTCATTGCCGGCGCCCTGGGGCCCCTACCTGTACTACACCCGTACCAccgc                                 >  1:1570274/1‑149 (MQ=255)
                                                                                                                     gCTGAAAACGCCTACCAGGAAGCCTGCCTGG‑CCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTGGAGACCGACCTGTCATTGCCGGCGCCCTGGGGCCCCTACCTGTACTACACCCGTACCAccgc                                 <  2:1570274/149‑1 (MQ=255)
                                                                                                                          aaaCGCCTACCAGGAAGCCTGCCTGG‑CCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTGGAGACCGACCTg                                                                                   >  1:107586/1‑94 (MQ=255)
                                                                                                                          aaaCGCCTACCAGGAAGCCTGCCTGG‑CCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTGGAGACCGACCTg                                                                                   <  2:107586/94‑1 (MQ=255)
                                                                                                                          aaaCGCCTACCAGGAAGCCTGCCTGG‑CCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTGGAGACCGACCTGTCATTGCCGGCGCCCTGGGGCCCCTACCTGTACTACACCCGTACCACCGCCGGcga                           >  1:1379054/1‑150 (MQ=255)
                                                                                                                          aaaCGCCTACCAGGAAGCCTGCCTGG‑CCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTGGAGACCGACCTGTCATTGCCGGCGCCCTGGGGCCCCTACCTGTACTACACCCGTACCACCGCCGGcga                           <  1:1318121/150‑1 (MQ=255)
                                                                                                                          aaaCGCCTACCAGGAAGCCTGCCTGG‑CCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTGGAGACCGACCTGTCATTGCCGGCGCCCTGGGGCCCCTACCTGTACTACACCCGTACCACCGCCGGcga                           <  2:1501379/150‑1 (MQ=255)
                                                                                                                                cTACCAGGAAGCCTGCCTGG‑CCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTGGAGACCGACCTGTCATTGCCGGCGCCCTGAGGCCCCTACCTGTACTACACCCGTACCACCGCCGGCGACGAata                     >  1:870667/1‑150 (MQ=255)
                                                                                                                                 tACCAGGAAGCCTGCCTGG‑CCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTGGAgaccgac                                                                                      <  2:2204157/84‑1 (MQ=255)
                                                                                                                                 tACCAGGAAGCCTGCCTGG‑CCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTGGAgaccgac                                                                                      >  1:2204157/1‑84 (MQ=255)
                                                                                                                                   ccAGGAAGCCTGCCTGG‑CCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTGGAGACCGACCTGTCATTGCCGGCGCCCTGGGGCCCCTACCTGTACTACACCCGTACCACCGCCGGCGACGAATATcc                  >  2:812217/1‑150 (MQ=255)
                                                                                                                                    cAGGAAGCCTGCCTGG‑CCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTGGAGACCGACCTGTCATTGCCGGCGCCCTGGGGCCCCTACCTGTACTacac                                             <  1:1066700/122‑1 (MQ=255)
                                                                                                                                    cAGGAAGCCTGCCTGG‑CCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTGGAGACCGACCTGTCATTGCCGGCGCCCTGGGGCCCCTACCTGTACTacac                                             >  2:1066700/1‑122 (MQ=255)
                                                                                                                                     aGGAAGCCTGCCTGG‑CCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTGGAGACCGACCTGTCATTGCCGGCGCCCTGGGGCCCCTACCTGTACTACACCCGTACCACCGCCGGCGACGAATATCcgc                <  2:2089667/150‑1 (MQ=255)
                                                                                                                                     aGGAAGCCTGCCTGG‑CCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTGGAGACCGACCTGTCATTGCCGGCGCCCTGGGGCCCCTACCTGTACTACACCCGTACCACCGCCGGCGACGAATATCcgc                <  2:1883246/150‑1 (MQ=255)
                                                                                                                                     aGGAAGCCTGCCTGG‑CCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTGGAGACCGACCTGTCATTGCCGGCGCCCTGGGGCCCCTACCTGTACTACACCCGTACCACCGCCGGCGACGAATATCcgc                <  1:1147071/150‑1 (MQ=255)
                                                                                                                                     aGGAAGCCTGCCTGG‑CCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTGGAGACCGACCTGTCATTGCCGGCGCCCTGGGGCCCCTACCTGTACTACACCCGTACCACCGCCGGCGACGAATATCcgc                >  1:1456984/1‑150 (MQ=255)
                                                                                                                                     aGGAAGCCTGCCTGG‑CCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTGGAGACCGACCTGTCATTGCCGGCGCCCTGGGGCCCCTACCTGTACTACACCCGTACCACCGCCGGCGACGAATATCcgc                <  1:2086678/150‑1 (MQ=255)
                                                                                                                                      ggAAGCCTGCCTGG‑CCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTGGAGACCGACCTg                                                                                   <  2:2015810/82‑1 (MQ=255)
                                                                                                                                      ggAAGCCTGCCTGG‑CCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTGGAGACCGACCTg                                                                                   >  1:2015810/1‑82 (MQ=255)
                                                                                                                                      ggAAGCCTGCCTGG‑CCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTGGAGACCGACCTGTCATTGCCGGCGCCCTgg                                                                 <  1:2365732/100‑1 (MQ=255)
                                                                                                                                      ggAAGCCTGCCTGG‑CCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTGGAGACCGACCTGTCATTGCCGGCGCCCTgg                                                                 >  2:2365732/1‑100 (MQ=255)
                                                                                                                                            gtgcctgG‑TCGGCCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTGGAGACCGACCTGTCATTGCCGGCGCCCTGGGGCCCCTACCTGTACTACAccc                                           >  1:948634/2‑116 (MQ=255)
                                                                                                                                            gtgcctgG‑TCGGCCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTGGAGACCGACCTGTCATTGCCGGCGCCCTGGGGCCCCTACCTGTACTACAccc                                           <  2:948634/115‑1 (MQ=255)
                                                                                                                                              gcctgG‑CCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTGGAGACCGACCTGTCATTGCCGGCGCCCTGGGGCCCCTACCTGTACTACACCCGTACCACCGCCGGCGACGAATATCCGCGTCACTAcc       <  1:1451938/150‑1 (MQ=255)
                                                                                                                                                   g‑CCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTGGAGACCGACCTGTCATTGCCGGCGCCCTGGGGCCCCTACCTGTACTACACCCGTACCACCGCCGGCGACGAATATCCGCGTCACTACCGCTGc  <  2:1875863/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
AAGAGACCTTCAATGCCAACCAAGCCCCAACCCCCGATTGCCCAAGCCGACCAGGCCACCGACCCGTATGCCTGGCTGCAACAGCGCGACACCCCCGAGGTACTCGCTTACCTGCAAGCTGAAAACGCCTACCAGGAAGCCTGCCTGG‑CCGACCAGGCACCGTTGCGTGAGCGCCTGTTCGAAGAAATCAAGGGCCGCATCCTGGAGACCGACCTGTCATTGCCGGCGCCCTGGGGCCCCTACCTGTACTACACCCGTACCACCGCCGGCGACGAATATCCGCGTCACTACCGCTGC  >  NC_002947/5201730‑5202026

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