Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,042,366 T→G 8.2% L137R (CTC→CGC) ‡ PP_3563 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,042,3660TG 8.2% 191.4 / 8.4 75L137R (CTC→CGC) ‡PP_3563hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (33/34);  new base G (2/4);  total (37/38)
Fisher's exact test for biased strand distribution p-value = 6.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.43e-01

TCAGCCTGTCGGGGCCACGCGACGCACGCGCCGCCATGCATGCCATCGAGGAAAGCTTCTGCCAGGTGGTGCTCGACCCGCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACAGAAGGTGCTGCAGGAGGTCAACGCCGCGCGCAAGGTGCCACGGCAATGCGGCGGCCAGCCGTTTGCCGCCGCCCCGCCACTAAGCTGGAGCACAGTGCTGGCCGGGGTCGCCG  >  NC_002947/4042219‑4042515
                                                                                                                                                   |                                                                                                                                                     
ccAGCCTGTCGGGGCCACGCGACGCACGCGCCGCCATGCATGCCATCGAGGAAAGCTTCTGCCAGGTGGTGCTCGACCCGCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCa                                                                                                                                                     <  2:1226538/149‑1 (MQ=255)
  aGCCTGTCGGGGCCACGCGACGCACGCGCCGCCATGCATGCCATCGAGGAAAGCTTCTGCCAGGTGGTGCTCGACCCGCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGc                                                                                                                                                   >  2:1228491/1‑150 (MQ=255)
   gCCTGTCGGGGCCACGCGACGCACGCGCCGCCATGCATGCCATCGAGGAAAGCTTCTGCCAGGTGGTGATAGACCCGCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCg                                                                                                                                                  <  2:450427/150‑1 (MQ=255)
     cTGTCGGGGCCACGCGACGCACGCGCCGCCATGCATGCCATCGAGGAAAGCTTCTGCCAGGTGGTGCTCGACCCGCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGNGgc                                                                                                                                                <  2:1198803/150‑1 (MQ=255)
     cTGTCGGGGCCACGCGACGCACGCGCCGCCATGCATGCCATCGAGGAAAGCTTCTGCCAGGTGGTGCTCGACCCGCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGNGgc                                                                                                                                                <  2:1199324/150‑1 (MQ=255)
     cTGTCGGGGCCACGCGACGCACGCGCCGCCATGCATGCCATCGAGGAAAGCTTCTGCCAGGTGGTGCTCGACCCGCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGgc                                                                                                                                                <  1:1360926/150‑1 (MQ=255)
        tCGGGGCCACGCGACGCACGCGCCGCCATGCATGCCATCGAGGAAAGCTTCTGCCAGGTGGTGCTCGACCCGCGATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGgccgc                                                                                                                                             >  2:949692/1‑150 (MQ=255)
        tCGGGGCCACGCGACGCACGCGCCGCCATGCATGCCATCGAGGAAAGCTTCTGCCAGGTGGTGCTCGACCCGCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGgccgt                                                                                                                                             <  1:790977/150‑2 (MQ=255)
        tCGGGGCCACGCGACGCACGCGCCGCCATGCATGCCATCGAGGAAAGCTTCTGCCAGGTGGTGCTCGACCCGCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGAGCAGCGgccgt                                                                                                                                             <  2:1174124/150‑2 (MQ=255)
            ggCCACGCGACGCACGCGCCGCCATGCATGCCATCGAGGAAAGCTTCTGCCTGGTGGTGCTCGACCCGCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTgg                                                                                                                                         <  1:1681580/150‑1 (MQ=255)
            ggCCACGCGACGCACGCGCCGCCATGCATGCCATCGAGGAAAGCTTCTGCCAGGTGGTGCTCGACCCGCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTgg                                                                                                                                         <  1:1004011/150‑1 (MQ=255)
             gCCACGCGACGCACGCGCCGCCATGCATGCCATCGAGGAAAGCTTCTGCCAGGTGGTGCTCGACCCGCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTggg                                                                                                                                        >  2:518899/1‑150 (MQ=255)
             gCCACGCGACGCACGCGCCGCCATGCATGCCATCGAGGAAAGCTTCTGCCAGGTGGTGCTCGACCCGCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTggg                                                                                                                                        >  1:1320805/1‑150 (MQ=255)
             gCCACGCGACGCACGCGCCGCCATGCATGCCATCGAGGAAAGCTTCTGCCAGGTGGTGCTCGACCCGCAATTCGTCGATAGCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGAGCAGCGgccgtgccag                                                                                                                                        <  1:737174/150‑7 (MQ=255)
                aCGCGACGCACGCGCCGCCATGCATGCCATCGAGGAAAGCTTCTGCCAGGTGGTGCTCGACCCGCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGa                                                                                                                                     <  1:1751566/150‑1 (MQ=255)
                aCGCGACGCACGCGCCGCCATGCATGCCATCGAGGAAAGCTTCTGCCAGGTGGTGCTCGACCCGCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACAGACGGATACTCCGGGCGCAGCCGCTGCTCAGCGGCCGCCTGGGCGa                                                                                                                                     >  1:1374136/1‑150 (MQ=255)
                  gcgACGCACGCGCCGCCATGCATGCCATCCAGGAAAGCTTCTGCCAGGTGGTGCTCGACCCGCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACt                                                                                                                                   <  2:240272/150‑1 (MQ=255)
                    gACGCACGCGCCGCCATGCATGCCATCGAGGAAAGCTTCTGCCAGGTGGTGCTCGACCCGCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTgg                                                                                                                                 <  1:52350/150‑1 (MQ=255)
                         acgcGCCGCCATGCATGCCATCGAGGAAAGCTTCTGCCAGGTGGTGCTCGACCCGTAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGc                                                                                                                            >  2:1021871/1‑150 (MQ=255)
                          cgcgCCGCCATGCATGCCATCGAGGAAAGCTTCTGCCAGGTGGTGCTCGACCCGCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGcc                                                                                                                           <  2:62676/150‑1 (MQ=255)
                           gcgcCGCCATGCATGCCATCGAGGAAAGCTTCTGCCAGGTGGTGCTCGACCCGCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGAGCAGCGgccgt                                                                                                                                             >  1:977335/1‑130 (MQ=255)
                           gcgcCGCCATGCATGCCATCGAGGAAAGCTTCTGCCAGGTGGTGCTCGACCCGCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGAGCAGCGgccgt                                                                                                                                             <  2:977335/131‑2 (MQ=255)
                              ccgccATGCATGCCATCGAGGAAAGCTTCTGCCAGGTGGTGCTCGACCCGCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAgg                                                                                                                       <  2:149664/150‑1 (MQ=255)
                               cgccATGCATGCCATCGAGGAAAGCTTCTGCCAGGTGGTGCTCGACCCGCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAggg                                                                                                                      >  2:1722351/1‑150 (MQ=255)
                                gccATGCATGCCATCGAGGAAAGCTTCTGCCAGGTGGTGCTCGACCCGCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGa                                                                                                                     <  1:1722351/150‑1 (MQ=255)
                                    tgcatgcCATCGAGGAAAGCTTCTGCCAGGTGGTGCTCGACCCGCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACAGa                                                                                                                 >  2:559740/1‑150 (MQ=255)
                                     gcatgcCATCGAGGAAAGCTTCTGCCAGGTGGTGCTCGACCCGCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCt                                                                                                                                           <  2:164503/123‑1 (MQ=255)
                                     gcatgcCATCGAGGAAAGCTTCTGCCAGGTGGTGCTCGACCCGCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCt                                                                                                                                           >  1:164503/1‑123 (MQ=255)
                                       atgcCATCGAGGAAAGCTTCTGCCAGGTGGTGCTCGACCCGCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACAGAAgg                                                                                                              <  2:1250888/150‑1 (MQ=255)
                                           cATCGGGGAAAGCTTCTGCCAGGTGGTGCGCGACCCGCAATTCGTCGATAGCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCAGCGTGCGGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACAGAAGGtgct                                                                                                          <  1:519511/150‑1 (MQ=255)
                                           cATCGAGGAAAGCTTCTGCCAGGTGGTGCTCGACCCGCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTgg                                                                                                                                         >  1:1219236/1‑119 (MQ=255)
                                           cATCGAGGAAAGCTTCTGCCAGGTGGTGCTCGACCCGCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTgg                                                                                                                                         <  2:1219236/119‑1 (MQ=255)
                                            aTCGAGGAAAGCTTCTGCCAGGTGGTGCTCGACCCGCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACAGAAGGtgctg                                                                                                         <  1:1538053/150‑1 (MQ=255)
                                             tCGAGGAAAGCTTCTGCCAGGTGGTGCTCGACCCGCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACAGAAGGtgctgc                                                                                                        >  1:747478/1‑150 (MQ=255)
                                              cGAGGAAAGCTTCTGCCAGGTGGTGCTCGACCCGCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACAGAAGGTGCTGCa                                                                                                       <  2:1445114/150‑1 (MQ=255)
                                                 ggAAAGCTTCTGCCAGGTGGTGCTCGACCCGCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACAGAAGGTGCTGCagga                                                                                                    <  2:709226/150‑1 (MQ=255)
                                                  gAAAGCTTCTGCCAGGTGGTGCTCGACCCGCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACAGAAGGTGCTGCaggag                                                                                                   <  1:1282450/150‑1 (MQ=255)
                                                   aaaGCTTCTGCCAGGTGGTGCTCGACCCGCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCGCGGCCGCTGAGCAGCGgccgtgc                                                                                                                                           >  2:938271/1‑106 (MQ=255)
                                                   aaaGCTTCTGCCAGGTGGTGCTCGACCCGCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCGCGGCCGCTGAGCAGCGgccgtgc                                                                                                                                           <  1:938271/109‑4 (MQ=255)
                                                   aaaGCTTCTGCCAGGTGGTGCTCGACCCGCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACAGAAGGTGCTGCaggagg                                                                                                  >  2:790799/1‑150 (MQ=255)
                                                    aaGCTTCTGCCAGGTGGTGCTCGACCCGCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACAGAAGGTGCTGCaggagg                                                                                                  <  2:767160/149‑1 (MQ=255)
                                                    aaGCTTCTGCCAGGTGGTGCTCGACCCGCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCCCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACAGAAGGTGCTGCaggagg                                                                                                  >  1:767160/1‑149 (MQ=255)
                                                        ttCTGCCAGGTGGTGCTCGACCCGCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACAGAAGGTGCTGCAGGAGGTCAAc                                                                                             >  1:1426355/1‑150 (MQ=255)
                                                           tGCCAGGTGGTGCTCGACCCGCAATTCGTCGATATCGGAGCCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACAGAAGGTGCTGCAGGAGGTCAAcgcc                                                                                          <  1:17365/150‑1 (MQ=255)
                                                             ccAGGTGGTGCTCGACCCGCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACAGAAGGTGCTGCAGGAGGTCAAcgccgc                                                                                        >  1:1629955/1‑150 (MQ=255)
                                                               aGGTGGTGCTCGACCCGCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACAGAAGGTGCTGCAGGAGGTCAACGCcgcgc                                                                                      <  1:1473339/150‑1 (MQ=255)
                                                                      gCTCGCCCCGCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACAGAAGGTGCTGCAGGAGGTCAACGCCGCGCGCAAGGt                                                                               <  1:1662678/150‑1 (MQ=255)
                                                                        tCGACCCGCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACAGAAGGTGCTGCAGGAGGTCAACGCCGCGCGCAAGGTGc                                                                             >  2:859362/1‑150 (MQ=255)
                                                                             ccGCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACAGAAGGTGCTGCAGGAGGTCAACGCCGCGCGCAAGGTGcc                                                                            >  2:749016/1‑146 (MQ=255)
                                                                             ccGCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACAGAAGGTGCTGCAGGAGGTCAACGCCGCGCGCAAGGTGcc                                                                            <  1:749016/146‑1 (MQ=255)
                                                                                       tcgATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGTTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACAGAAGGTGCTGCAGGAGGTCAACGCCGCGCGCAAGGTGCCACGGCAATgcggcg                                                              >  2:1512723/1‑150 (MQ=255)
                                                                                       tcgATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACAGAAGGTGCTGCAGGAGGTCAACGCCGCGCGCAAGGTGCCACGGCAATgcggcg                                                              <  1:1409252/150‑1 (MQ=255)
                                                                                             tCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACAGAAGGTGCTGCAGGAGGTCAACGCCGCGCGCAAGGTGCCACGGCAATGCGGCGGCCAGc                                                        >  2:1525667/1‑150 (MQ=255)
                                                                                                gCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACAGAAGGTGCTGCAGGAGGTCAACGCCGCGCGCAAGGTGCCACGGCAATGCGGCGGCCAGCCGt                                                     >  2:40075/1‑150 (MQ=255)
                                                                                                 cGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACAGAAGGTGCTGCAGGAGGTCAACGCCGCGCGCAAGGTGCCACGGCAATGCGGCGGCCAGCCGtt                                                    <  2:112129/150‑1 (MQ=255)
                                                                                                  gTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCCCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACAGAAGGTGCTGCAGGAGGTCAACGCCGCGCGCAAGGTGCCACGGCAATGCGGCGGCCAGCCGttt                                                   >  2:1910159/1‑150 (MQ=255)
                                                                                                     aGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACAGAAGGTGCTGCAGGAGGTCAACGCCGCGCGCAAGGTGCCACGGCAATGCGGCGGCCAGCCGTTTgcc                                                >  2:1317822/1‑150 (MQ=255)
                                                                                                     aGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACAGAAGGTGCTGCAGGAGGTCAACGCCGCGCGCAAGGTGCCACGGCAATGCGGCGGCCAGCCGTTTgcc                                                >  2:1317554/1‑150 (MQ=255)
                                                                                                     aGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACAGAAGGTGCTGCAGGAGGTCAACGCCGCGCGCAAGGTGCCACGGCAATGCGGCGGCCAGCCGTTTgcc                                                >  1:542753/1‑150 (MQ=255)
                                                                                                        cAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACAGAAGGTGCTGCAGGAGGTCAACGCCGCGCGCAAGGTGCCACGGCAATGCGGCGGCCAGCCGTTTgccgcc                                             >  2:622412/1‑150 (MQ=255)
                                                                                                        cAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACAGAAGGTGCTGCAGGAGGTCAACGCCGCGCGCAAGGTGCCACGGCAATGCGGCGGCCAGCCGTTTgccgcc                                             <  1:949692/150‑1 (MQ=255)
                                                                                                          ggaggGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACAGAAGGTGCTGCAGGAGGTCAACGCCGCGCGCAAGGTGCCACGGCAATGCGGCGGCCAGCCGTTTgccgccgc                                           <  1:40075/150‑1 (MQ=255)
                                                                                                             gggCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACAGAAGGTGCTGCAGGAGGTCAACGCCGCGCGCAAGGTGCCACGGCAATGCGGCGGCCAGCCGTTTgccgccgccgc                                        >  1:969666/1‑148 (MQ=255)
                                                                                                                acgtgACTTGCGGATTGTTCTAGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCACGCCGAGGGACAGACGGTGCTGCAGGAGGTCAACGCCGCGCGCAAGGTGCCACGGCAATGCGGCGGCCAGCCGTTTGCCGccgccccgcc                                     >  2:25302/5‑150 (MQ=255)
                                                                                                                 cgcgACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACAGAAGGTGCTGCAGGAGGTCAACGCCGCGCGCAAGGTGCCACGGCAATGCGGCGGCCAGCCGTTTGCCGCCGCCCCGCCa                                    <  1:790799/150‑1 (MQ=255)
                                                                                                                         gcgcATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAgg                                                                                                                       <  2:1509061/59‑1 (MQ=255)
                                                                                                                         gcgcATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAgg                                                                                                                       >  1:1509061/1‑59 (MQ=255)
                                                                                                                           gcATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACAGAAGGTGCTGCa                                                                                                       <  1:1457272/73‑1 (MQ=255)
                                                                                                                           gcATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACAGAAGGTGCTGCa                                                                                                       >  2:1457272/1‑73 (MQ=255)
                                                                                                                            cATCGTGCTGGCACGGCCGCTGCTCCCTGCCCTCCTGGGCGACTGGCAGGCCGAGGGACAGAAGGTGCTGCAGGAGGTCAACGCCGCGCGCAAGGTGCCACGGCAATGCGGCGGCCAGCCGTTTGCCGCCGCCCCGCCACTAAGCTGGAg                         >  1:1483282/1‑150 (MQ=255)
                                                                                                                              tCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACAGAAGGTGCTGCAGGAGGTCAACGCCGCGCGCAAGGTGCCACGGCAATGCGGCGGCCAGCCGTTTGCCGCCGCCCCGCCACTAAGCTGGAGca                       >  1:1556353/1‑150 (MQ=255)
                                                                                                                                      gCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACAGAAGGTGCTGCAGGAGGTCAACGCCGCGCGCAAGGTGCCACGGCAATGCGGCGGCCAGCCGTTTGCCGCCGCCCCGCCACTAAGCTGGAGCACAGTGCTg               >  1:1839034/1‑150 (MQ=255)
                                                                                                                                      gCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACAGAAGGTGCTGCAGGAGGTCAACGCCGCGCGCAAGGTGCCACGGCAATGCGGCGGCCAGCCGTTTGCCGCCGCCCCGCCACTAAGCTGGAGCACAGTGCTg               >  1:1851724/1‑150 (MQ=255)
                                                                                                                                        aCGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACAGAAGGTGCTGCAGGAGGTCAACGCCGCGCGCAAGGTGCCACGGCAATGCGGCGGCCAGCCGTTTGCCGCCTCCCCGCCACTAAGCTGGAGCACAGTGCTGGc             >  1:727629/1‑150 (MQ=255)
                                                                                                                                        aCGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACAGAAGGTGCTGCAGGAGGTCAACGCCGCGCGCAAGGTGCCACGGCAATGCGGCGGCCAGCCGTTTGCCGCCGCCCCGCCACTAAGCTGGAGCACAGTGCTGGc             >  2:307825/1‑150 (MQ=255)
                                                                                                                                               ctgctCAGCGGCCGCCCGGTCGCCTGGCAGGCCGAGGGACAGAAGGTGCTGCAGGAGCTCAACGCCTCGCGCAAGGTGCCACGGCAATGCGGCGGCCAGCCGTTTGCCGGCGCGCCGCCACTAAGCTGCAGCACAGTGCTGGCCGGGGTc      >  2:53209/1‑150 (MQ=255)
                                                                                                                                                 gctCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACAGAAGGTGCTGCAGGAGGTCAACGCCGCGCGCAAGGTGCCACGGCAATGCGGCGGCCAGCCGTTTGCCGCCGCCCCGCCACTAAGCTGGAGCACAGTGCTGGCCGGGGTcgg    >  2:135344/1‑149 (MQ=255)
                                                                                                                                                  ctCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACAGAAGGTGCTGCAGGAGGTCAACGCCGCGCGCAAGGTGCCACGGCAATGCGGCGGCCGGCCGTTTGCCGCCGCCCCGCCACTAAGCTGGAGCACAGTGCTGGCCGGGGTcgcc   >  1:300096/1‑150 (MQ=255)
                                                                                                                                                   tCAGCGGCCGCCGGGGCGACTGGCAGGCCGAGGGACAGAAGGTGCTGCAGGAGGTCAACGCCGCGCGCAAGGTGCCACGGCAATGCGGCGGCCAGCCGTTTGCCGCCGCCCCGCCACTAAGCTGGAGCACAGTGCTGGCCGGGGTcgccg  >  2:1508693/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
TCAGCCTGTCGGGGCCACGCGACGCACGCGCCGCCATGCATGCCATCGAGGAAAGCTTCTGCCAGGTGGTGCTCGACCCGCAATTCGTCGATATCGGCGTCAGCCAGGAGGGCCGCGACTGGCGCATCGTGCTGGCACGGCCGCTGCTCAGCGGCCGCCTGGGCGACTGGCAGGCCGAGGGACAGAAGGTGCTGCAGGAGGTCAACGCCGCGCGCAAGGTGCCACGGCAATGCGGCGGCCAGCCGTTTGCCGCCGCCCCGCCACTAAGCTGGAGCACAGTGCTGGCCGGGGTCGCCG  >  NC_002947/4042219‑4042515

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: ‑