Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,389,586 Δ1 bp 5.6% coding (1199/1956 nt) PP_3865 ← tail protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,389,5860A. 5.6% 212.9 / 6.9 72coding (1199/1956 nt)PP_3865tail protein
Reads supporting (aligned to +/- strand):  ref base A (33/33);  new base . (2/2);  total (37/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.42e-01

CGGCCCGACCTTCTCCCAGTTGGCGACGATCAGACCGGCCGCCGCCGCGATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCTGATCCACGGCCGCAGTTGCTCGAGCACCGCATTGATGCCGGGCAGCAAAGCATTGCCGATAGCAACCGCTACCCT  >  NC_002947/4389438‑4389731
                                                                                                                                                    |                                                                                                                                                 
cGGCCCGACCTTCTCCCAGTTGGCGACGATCAGACCGGCCGCCGCCGCGATCCCCGCGGCGATCAGGCAAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGaa                                                                                                                                                  <  2:103554/150‑1 (MQ=255)
   cccGACCTTCTCCCAGTTGGCGACGATCAGACCGGCCGCCGCCGCGATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGgc                                                                                                                                               >  2:318196/1‑150 (MQ=255)
    ccGACCTTCTCCCAGTTGGCGACGATCAGACCGGCCGCCGCCGCGATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGaa                                                                                                                                                  >  2:319394/1‑146 (MQ=255)
    ccGACCTTCTCCCAGTTGGCGACGATCAGACCGGCCGCCGCCGCGATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGaa                                                                                                                                                  <  1:319394/146‑1 (MQ=255)
      gACCTTCTCCCAGTTGGCGACGATCAGACCGGCCGCCGCCGCGATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCCCCCGCCGCGCGCTGCAcgccgc                                                                                                                                            >  2:395749/1‑150 (MQ=255)
       aCCTTCTCCCAGTTGGCGACGATCAGACCGGCCGCCGCCGCGATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGgccgcc                                                                                                                                           >  1:1157908/1‑150 (MQ=255)
              cccAGTTGGCGACGATCAGACCGGCCGCCGCCGCGATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCaa                                                                                                                                    <  2:199350/150‑1 (MQ=255)
                 aGTTGGCGACGATCAGACCGGCCGCCGCCGCGATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGa                                                                                                                                 <  1:318196/150‑1 (MQ=255)
                  tatcccgacgaTACGACCGGCCGCCGCCGGGCTCCCCGCGGCGATACGGCCAATAGGCGTGGCCCAGACGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGAc                                                                                                                                <  2:588695/145‑1 (MQ=255)
                     ggCGACGATCAGACCGGCCGCCGCCGCGATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGAt                                                                                                                             >  1:850492/1‑150 (MQ=255)
                      gCGACGATCAGACCGGCCGCCGCCGCGATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGAGCGGTGAAGGCCGCCCCGGCAATGACGATg                                                                                                                            >  1:108959/1‑150 (MQ=255)
                      gCGACGATCAGACCGGCCGCCGCCGCGATCCCAGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATg                                                                                                                            >  1:887309/1‑150 (MQ=255)
                        gacgaTCAGACCGGCCGCCGCCGCGATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGagatc                                                                                                                          >  1:588695/1‑146 (MQ=255)
                             tCAGACCGGCCGCCGCCGCGATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCa                                                                                                                     <  2:850492/150‑1 (MQ=255)
                             tCAGACCGGCCGCCGCCGCGATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACACGCGTGGCCATTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCa                                                                                                                     >  1:721152/1‑150 (MQ=255)
                              cAGACCGGCCGCCGCCGCGATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCAGACCTAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCAc                                                                                                                    <  2:680322/150‑1 (MQ=255)
                                gACCGGCCGCCGCCGCGATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCa                                                                                                                  >  1:1214771/1‑150 (MQ=255)
                                   cGGCCGCCGCCGCGATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACt                                                                                                               >  2:532651/1‑150 (MQ=255)
                                       cgccgccgcGATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCgg                                                                                                           <  1:1023395/150‑1 (MQ=255)
                                            ccgcGATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGt                                                                                                      >  2:156916/1‑150 (MQ=255)
                                               cgATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGc                                                                                                   <  1:566460/150‑1 (MQ=255)
                                                gATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCt                                                                                                  >  2:559226/1‑150 (MQ=255)
                                                gATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCt                                                                                                  <  2:260394/150‑1 (MQ=255)
                                                     ccGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGatca                                                                                             >  2:1151696/1‑150 (MQ=255)
                                                     ccGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTG‑AGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGatcat                                                                                            >  1:1208706/1‑150 (MQ=255)
                                                                  gCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCacat                                                                                 <  2:553293/149‑1 (MQ=255)
                                                                  gCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCacat                                                                                 >  1:553293/1‑146 (MQ=255)
                                                                    cAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCCTGCAATGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGc                                                                              >  1:1054372/1‑150 (MQ=255)
                                                                     aaTCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCt                                                                             <  1:1151696/150‑1 (MQ=255)
                                                                       tCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCgg                                                                                                           >  1:515285/1‑118 (MQ=255)
                                                                       tCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCgg                                                                                                           <  2:515285/118‑1 (MQ=255)
                                                                       tCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCTGa                                                                           <  2:1214771/150‑1 (MQ=255)
                                                                       tCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCTGa                                                                           >  2:565546/1‑150 (MQ=255)
                                                                            gTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCTGATCCAc                                                                      <  1:395749/150‑1 (MQ=255)
                                                                              ggCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCTGATCCACg                                                                     <  2:992199/149‑1 (MQ=255)
                                                                              ggCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCTGATCCACg                                                                     >  1:992199/1‑149 (MQ=255)
                                                                                  gcgAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCGCCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCTGATCCACgg                                                                    >  1:639218/1‑146 (MQ=255)
                                                                                  gcgAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCTGATCCACgg                                                                    <  2:639218/146‑1 (MQ=255)
                                                                                        gCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCTGATCCACGGCCGCAGTTGc                                                          >  2:349964/1‑150 (MQ=255)
                                                                                          cacaCCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCTGATCCACGGCCGCAGTTGCTc                                                        >  1:249229/1‑150 (MQ=255)
                                                                                           acacCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGATGCAGGCCGCCCCGGCAAGGCGCATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCTGATCCACGGCCGCAGTTGCTCg                                                       >  1:156405/1‑150 (MQ=255)
                                                                                             acCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCTGATCCACGGCCGCAGTTGCTCGAg                                                     >  2:738114/1‑150 (MQ=255)
                                                                                             acCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCTGATCCACGGCCGCAGTTGCTCGAg                                                     <  1:349964/150‑1 (MQ=255)
                                                                                               cGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGc                                                                                                             <  2:575692/92‑1 (MQ=255)
                                                                                               cGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGc                                                                                                             >  1:575692/1‑92 (MQ=255)
                                                                                                   cACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTG‑AGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCTGATCCACGGCCGCAGTTGCTc                                                        >  1:1261816/1‑140 (MQ=255)
                                                                                                   cACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTG‑AGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCTGATCCACGGCCGCAGTTGCTc                                                        <  2:1261816/140‑1 (MQ=255)
                                                                                                    aCCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTG‑AGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCTGATCCACGGCCGCAGTTGCTCGAGCACCGCAt                                             <  2:1208706/150‑1 (MQ=255)
                                                                                                      ccGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGcc                                                                                     >  1:966213/1‑109 (MQ=255)
                                                                                                      ccGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGcc                                                                                     <  2:966213/109‑1 (MQ=255)
                                                                                                      ccGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCTGATCCACGGCCGCAGTTGCTCGAGCACCGCAtt                                            >  2:21285/1‑150 (MQ=255)
                                                                                                          gTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCGTATCATCTGCGCCATATGGCTGATCCACGGCCGCAGTTGCTCGAGCACCGCATTGATg                                        <  1:156916/150‑1 (MQ=255)
                                                                                                             gCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCTGATCCACGGCCGCAGTTGCTCGAGCACCGCATTGATGCCg                                     <  1:21285/150‑1 (MQ=255)
                                                                                                              ccACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCTGATCCACGGCCGCAGTTGCTCGAGCACCGCATTGATGCCgg                                    <  2:950302/150‑1 (MQ=255)
                                                                                                                aCTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCgg                                                                                                           >  2:850939/1‑77 (MQ=255)
                                                                                                                aCTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCgg                                                                                                           <  1:850939/77‑1 (MQ=255)
                                                                                                                 cTGTGGCCGCGAACACCGCCGACCGCAGCGCGGTGACGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGATTGTCGCTGATCATCTGCGCCATATGGCTGATCCACGGCCGCAGTTGCTCGAGCACCGCATTGATGCCGGgca                                 >  2:15828/1‑150 (MQ=255)
                                                                                                                  tgtgGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCTGATCCACGGCCGCAGTTGCTCGAGCACCGCATTGATGCCGGgcag                                >  2:407143/1‑150 (MQ=255)
                                                                                                                      gCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCTGATCCACGGCAGCAGTTGCTCGAGCACCGCATTGATGCCGGGCAGCaaa                            <  1:565546/150‑1 (MQ=255)
                                                                                                                        cgcgAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCTGATCCACGGCCGCAGTTGCTCGAGCACCGCATTGATGCCGGGCAGCAAAGc                          <  2:1092717/150‑1 (MQ=255)
                                                                                                                         gcgAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCTGATc                                                                         >  1:1058630/1‑102 (MQ=255)
                                                                                                                         gcgAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCTGATc                                                                         <  2:1058630/102‑1 (MQ=255)
                                                                                                                          cgAACACCGCCGCACGCAGAGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCTGATCCACGGCCGCAGTTGCTCGAGCACCGCATTGATGCCGGGCAGCAAAGCAt                        <  2:108959/150‑1 (MQ=255)
                                                                                                                           gAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTg                                                                                         >  2:949236/1‑84 (MQ=255)
                                                                                                                           gAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTg                                                                                         <  1:949236/84‑1 (MQ=255)
                                                                                                                              cacCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCTGATCCACGGCCGCAGTTGCTCGAGCACCGCATTGATGCCGGGCAGCAAAGCATTGcc                    >  2:1169405/1‑150 (MQ=255)
                                                                                                                              cacCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCTGATCCACGGCCGCAGTTGCTCGAGCACCGCATTGATGCCGGGCAGCAAAGCATTGcc                    <  1:624637/150‑1 (MQ=255)
                                                                                                                               acCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCTGATCCACGGCCGCAGTTGCTCGAGCACCGCATTGATGCCGGGCAGCAAAGCATTGCCg                   >  1:463585/1‑150 (MQ=255)
                                                                                                                                        cgcaGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCTGATCCACGGCCGCAGTTGCTTGAGAAACGCATTGGTGCGTGGCAGCAAAGCATTGCCGCTAGCCAcc          >  2:102628/1‑150 (MQ=255)
                                                                                                                                          caGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCTGATCCACGGCCGCAGTTGCTCGAGCACCGCATTGATGCCGGGCAGCAAAGCATTGCCGATAGCAACCGc        <  2:1100086/150‑1 (MQ=255)
                                                                                                                                          caGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCTGATCCACGGCCGCAGTTGCTCGAGCACCGCATTGATGCCGGGCAGCAAAGCATTGCCGATAGCAACCGc        <  2:1100051/150‑1 (MQ=255)
                                                                                                                                             cgcgGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCTGATCCACGGCCGCAGTTGCTCGAGCACCGCATTGATGCCGGGCAGCAAAGCATTGCCGATAGCAACCGCTAc     <  1:1169405/150‑1 (MQ=255)
                                                                                                                                             cgcgGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCTGATCCACGGCCGCAGTTGCTCGAGCACCGCATTGATGCCGGGCAGCAAAGCATTGCCGATAGCAACCGCTAc     <  1:559226/150‑1 (MQ=255)
                                                                                                                                              gcgGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCTGATCCACGGCCGCAGTTGCTCGAGCACCGCATTGATGCCGGGCAGCAAAGCATTGCCGATAGCAACCGCTAcc    <  1:532651/150‑1 (MQ=255)
                                                                                                                                                ggTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATAGCTGATCCACGGCCGCAGTTGCTCGAGCACCGCATTGATGCCCGGCAGCAAAGCATTGCCGATAGCAACCGGTAACCt  >  1:36404/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
CGGCCCGACCTTCTCCCAGTTGGCGACGATCAGACCGGCCGCCGCCGCGATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCTGATCCACGGCCGCAGTTGCTCGAGCACCGCATTGATGCCGGGCAGCAAAGCATTGCCGATAGCAACCGCTACCCT  >  NC_002947/4389438‑4389731

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