Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 571,621 T→C 9.1% intergenic (+44/+47) ssb → / ← PP_0486 single‑stranded DNA‑binding protein/GntR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947571,6210TC 9.1% 381.3 / 19.9 131intergenic (+44/+47)ssb/PP_0486single‑stranded DNA‑binding protein/GntR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (55/64);  new base C (7/5);  total (62/69)
Fisher's exact test for biased strand distribution p-value = 5.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01

Only 100 of 135 total aligned reads displayed.

CAGGCCCAGCGCCCGCAACAGGCCCCTCAGCGCCCAGCGCCGCAACAGCCTGCGCCGCAGCCGGCGGCTGACTTCGACAGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGAATCCAGCAGGCTTGCCAGCTTGGTCAGGGCAAAGCCCAGTTCTTGGCGGCTGGGTGCCATCAGGGCGATGCGCACA  >  NC_002947/571473‑571770
                                                                                                                                                    |                                                                                                                                                     
cagGCCCAGCGCCCGCAACAGGCCCCTCAGCGCCCAGCGCCGCAACAGCCTGCGCCGCAGCCGGCGGCTGACTTCGACAGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGtt                                                                                                                                                      <  2:4861761/150‑1 (MQ=255)
cagGCCCAGCGCCCGCAACAGGCCCCTCAGCGCCCAGCGCCGCAACAGCCTGCGCCGCAGCCGGCGGCTGACTTCGACAGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGtt                                                                                                                                                      <  2:4663241/150‑1 (MQ=255)
 agGCCCAGCGCCCGCAACAGGCCCCTCAGCGCCCAGCGCCGCAACAGCCTGCGCCGCAGCCGGCGGCTGACTTCGACAGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTg                                                                                                                                                     >  1:108326/1‑150 (MQ=255)
    cccAGCGCCCGCAACAGGCCCCTCAGCGCCCAGCGCCGCAACAGCCTGCGCCGCAGCCGGCGGCTGACTTCGACAGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGccc                                                                                                                                                  <  2:921682/150‑1 (MQ=255)
    cccAGCGCCCGCAACAGGCCCCTCAGCGCCCAGCGCCGCAACAGCCTGCGCCGCAGCCGGCGGCTGACTTCGACAGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGccc                                                                                                                                                  <  2:2206723/150‑1 (MQ=255)
     ccAGCGCCCGCAACAGGCCCCTCAGCGCCCAGCGCCGCAACAGCCTGCGCCGCAGCCGGCGGCTGACTTCGACAGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCACGACCAGAAACCCAGGGCAGTTGCCCt                                                                                                                                                 >  1:5054043/1‑150 (MQ=255)
     ccAGCGCCCGCAACAGGCCCCTCAGCGCCCAGCGCCGCAACAGCCTGCGCCGCAGCCGGCGGCTGACTTCGACAGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCt                                                                                                                                                 >  2:4908909/1‑150 (MQ=255)
     ccAGCGCCCGCAACAGGCCCCTCAGCGCCCAGCGCCGCAACAGCCTGCGCCGCAGCCGGCGGCTGACTTCGACAGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCt                                                                                                                                                 >  1:3166813/1‑150 (MQ=255)
     ccAGCGCCCGCAACAGGCCCCTCAGCGCCCAGCGCCGCAACAGCCTGCGCCGCAGCCGGCGGCTGACTTCGACAGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCt                                                                                                                                                 >  2:3799703/1‑150 (MQ=255)
       aGCGCCCGCAACAGGCCCCTCAGCGCCCAGCGCCGCAACAGCCTGCGCCGCAGCCGGCGGCTGACTTCGACAGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTgg                                                                                                                                               <  1:1036057/150‑1 (MQ=255)
       aGCGCCCGCAACAGGCCCCTCAGCGCCCAGCGCCGCAACAGCCTGCGCCGCAGCCGGCGGCTGACTTCGACAGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTgg                                                                                                                                               <  2:192222/150‑1 (MQ=255)
       aGCGCCCGCAACAGGCCCCTCAGCGCCCAGCGCCGCAACAGCCTGCGCCGCAGCCGGCGGCTGACTTCGACAGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTgg                                                                                                                                               <  1:2440551/150‑1 (MQ=255)
        gcgcCCGCAACAGGCCCCTCAGCGCCCAGCGCCGCAACAGCCTGCGCCGCAGCCGGCGGCTGACTTCGACAGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAACTGCCCTggg                                                                                                                                              <  1:94199/150‑1 (MQ=255)
            ccGCAACAGGCCCCTCAGCGCCCAGCGCCGCAACAGCCTGCGCCGCAGCCGGCGGCTGACTTCGACAGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAATTGCCCTGGGtttt                                                                                                                                          >  1:3160606/1‑150 (MQ=255)
             cGCAACAGGCCCCTCAGCGCCCAGCGCCGCAACAGCCTGCGCCGCAGCCGGCGGCTGACTTCGACAGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGttttt                                                                                                                                         <  2:247427/150‑1 (MQ=255)
             cGCAACAGGCCCCTCAGCGCCCAGCGCCGCAACAGCCTGCGCCGCAGCCGGCGGCTGACTTCGACAGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAACTGCCCTGGGtttct                                                                                                                                         <  2:3642486/150‑3 (MQ=255)
             cGCAACAGGCCCCTCAGCGCCCAGCGCCGCAACAGCCTGCGCCGCAGCCGGCGGCTGACTTCGACAGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAACTGCCCTGGGtttct                                                                                                                                         <  2:3642415/150‑3 (MQ=255)
              gCAACAGGCCCCTCAGCGCCCAGCGCCGCAACAGCCTGCGCCGCAGCCGGCGGCTGACTTCGACAGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAACTGCCCTgg                                                                                                                                               <  2:259470/143‑1 (MQ=255)
               cAACAGGCCCCTCAGCGCCCAGCGCCGCAACAGCCTGCGCCGCAGCCGGCGGCTGACTTCGACAGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGttttttt                                                                                                                                       >  2:1201536/1‑150 (MQ=255)
                aaCAGGCCCCTCAGCGCCCAGCGCCGCAACAGCCTGCGCCGCAGCCGGCGGCTGACTTCGACAGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGTCCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTa                                                                                                                                      >  2:4364237/1‑150 (MQ=255)
                 aCAGGCCCCTCAGCGCCCAGCGCCGCAACAGCCTGCGCCGCAGCCGGCGGCTGACTTCGACAGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTAt                                                                                                                                     <  1:689069/150‑1 (MQ=255)
                  cAGGCCCCTCAGCGCCCAGCGCCGCAACAGCCTGCGCCGCAGCCGGCGGCTGACTTCGACAGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATg                                                                                                                                    >  1:4128002/1‑150 (MQ=255)
                    ggCCCCTCAGCGCCCAGCGCCGCAACAGCCTGCGCCGCAGCCGGCGGCTGACTTCGACAGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATggg                                                                                                                                  <  1:1538768/150‑1 (MQ=255)
                      ccccTCAGCGCCCAGCGCCGCAACAGCCTGCGCCGCAGCCGGCGGCTGACTTCGACAGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGgca                                                                                                                                <  2:3899220/150‑1 (MQ=255)
                          tCAGCGCCCAGCGCCGCAACAGCCTGCGCCGCAGCCGGCGGCTGACTTCGACAGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCAc                                                                                                                            >  1:2060287/1‑150 (MQ=255)
                          tCAGCGCCCAGCGCCGCAACAGCCTGCGCCGCAGCCGGCGGCTGACTTCGACAGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCAc                                                                                                                            >  1:2060136/1‑150 (MQ=255)
                                 ccagcgccGCAACAGCCTGCGCCGCAGCCGGCGGCTGACTTCGACAGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTAt                                                                                                                     >  1:5110050/1‑150 (MQ=255)
                                   agcgccGCAACAGCCTGCGCCGCAGCCGGCGGCTGACTTCGACAGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCa                                                                                                                   <  1:1119291/150‑1 (MQ=255)
                                   agcgccGCAACAGCCTGCGCCGCAGCCGGCGGCTGACTTCGACAGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCa                                                                                                                   <  1:2125127/150‑1 (MQ=255)
                                   agcgccGCAACAGCCTGCGCCGCAGCCGGCGGCTGACTTCGACAGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGATTTTTTATGGGCAGCACCGCGTATCa                                                                                                                   <  2:2077581/150‑1 (MQ=255)
                                     cgccgcAACAGCCTGCGCCGCAGCCGGCGGCTGACTTCGACAGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAgg                                                                                                                 <  2:2129823/150‑1 (MQ=255)
                                     cgccgcAACAGCCTGCGCCGCAGCCGGCGGCTGACTTCGACAGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAgg                                                                                                                 <  1:3333086/150‑1 (MQ=255)
                                         gcAACAGCCTGCGCCGCAGCCGGCGGCTGACTTCGACAGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATc                                                                                                             <  2:2354300/150‑1 (MQ=255)
                                           aaCAGCCTGCGCCGCAGCCGGCGGCTGACTTCGACAGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCGTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGc                                                                                                           >  1:4776555/1‑150 (MQ=255)
                                            aCAGCCTGCGCCGCAGCCGGCGGCTGACTTCGACAGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCt                                                                                                          <  1:3844446/150‑1 (MQ=255)
                                              aGCCTGCGCCGCAGCCGGCGGCTGACTTCGACAGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCttt                                                                                                        <  2:2412386/150‑1 (MQ=255)
                                               gCCTGCGCCGCAGCCGGCGGCTGACTTCGACAGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTc                                                                                                       >  1:3654555/1‑150 (MQ=255)
                                                ccTGCGCCGCAGCCGGCGGCTGACTTCGACAGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCa                                                                                                      >  1:2664478/1‑150 (MQ=255)
                                                ccTGCGCCGCAGCCGGCGGCTGACTTCGACAGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGATTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCa                                                                                                      <  1:533600/150‑1 (MQ=255)
                                                  tGCGCCGCAGCCGGCGGCTGACTTCGACAGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGt                                                                                                    >  2:1471299/1‑150 (MQ=255)
                                                   gcgcCGCAGCCGGCGGCTGACTTCGACAGCTTCGATGACGATAGTCCGTGCTGCTCAGAGCGGGAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTa                                                                                                   <  2:513141/150‑1 (MQ=255)
                                                        gcAGCCGGCGGCTGACTTCGACAGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGaacc                                                                                              >  1:4775338/1‑149 (MQ=255)
                                                              ggcggcTGACTTCGACAGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTACAGTAGAACAACCggg                                                                                        >  1:3749599/1‑150 (MQ=255)
                                                                    tGACTTCGACAGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCg                                                                                  >  2:3280373/1‑150 (MQ=255)
                                                                        ttCGACAGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGGACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGa                                                                              >  2:973083/1‑150 (MQ=255)
                                                                            aCAGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGAATcc                                                                          >  2:135047/1‑150 (MQ=255)
                                                                              aGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGAATCcag                                                                        <  2:1914667/150‑1 (MQ=255)
                                                                                cTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGAATCcagca                                                                      >  2:1524103/1‑150 (MQ=255)
                                                                                 ttCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGttt                                                                                                                                           <  2:3454683/80‑1 (MQ=255)
                                                                                     aTGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGAATCCAGCAGGCtt                                                                 >  2:4261117/1‑150 (MQ=255)
                                                                                       gACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCTGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGAATCCAGCAGGCTTGc                                                               <  1:232800/150‑1 (MQ=255)
                                                                                       gACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGAATCCAGCAGGCTTGc                                                               <  2:4180735/150‑1 (MQ=255)
                                                                                       gACGATATTCCGTTCTGATCAGAGCGGTAGGGCGGTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGAATCCAGCAGGCTTGc                                                               <  2:2934716/150‑1 (MQ=255)
                                                                                              ttCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGAATCCAGCAGGCTTGCCAGCTTg                                                        <  2:3595910/150‑1 (MQ=255)
                                                                                                  gTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGAATCCAGCAGGCTTGCCAGCTTGGTCa                                                    <  1:3636986/150‑1 (MQ=255)
                                                                                                  gTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGAATCCAGCAGGCTTGCCAGCTTGGTCa                                                    <  2:129689/150‑1 (MQ=255)
                                                                                                    tCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGAATCCAGCAGGCTTGCCAGCTTGGTCAgg                                                  <  1:3207838/150‑1 (MQ=255)
                                                                                                    tCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGAATCCAGCAGGCTTGCCAGCTTGGTCAgg                                                  <  1:3286525/150‑1 (MQ=255)
                                                                                                      tGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGAATCCAGCAGGCTTGCCAGCTTGGTCAGGGc                                                >  2:630799/1‑150 (MQ=255)
                                                                                                      tGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGAATCCAGCAGGCTTGCCAGCTTGGTCAGGGc                                                >  2:630801/1‑150 (MQ=255)
                                                                                                      tGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGAATCCAGCAGGCTTGCCAGCTTGGTCAGGGc                                                >  1:4513580/1‑150 (MQ=255)
                                                                                                         tCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGAATCCAGCAGGCTTGCCAGCTTGGTCAGGGCaaa                                             <  1:1607803/150‑1 (MQ=255)
                                                                                                           agagCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGAATCCAGCAGGCTTGCCAGCTTGGTCAGGGCAAAGc                                           >  1:1554353/1‑150 (MQ=255)
                                                                                                            gagCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGAATCCAGCAGGCTTGCCAGCTTGGTCAGGGCAAAGcc                                          <  2:2139027/150‑1 (MQ=255)
                                                                                                            gagCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGAATCCAGCAGGCTTGCCAGCTTGGTCAGGGCAAAGcc                                          <  1:4140700/150‑1 (MQ=255)
                                                                                                            gagCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGAATCCAGCAGGCTTGCCAGCTTGGTCAGGGCAAAGcc                                          <  2:3247056/150‑1 (MQ=255)
                                                                                                            gagCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGCATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGAATCCAGCAGGCTTGCCAGCTTGGTCAGGGCAAAGcc                                          >  1:2286548/1‑150 (MQ=255)
                                                                                                              gCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAATTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGAATCCAGCAGGCTTGCCAGCTTGGTCAGGGCAAAGCCCa                                        <  2:3160606/150‑1 (MQ=255)
                                                                                                               cGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGAATCCAGCAGGCTTGCCAGCTTGGTCAGGGCAAAGCCCAg                                       <  1:4264988/150‑1 (MQ=255)
                                                                                                                    gggCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGAATCCAGCAGGCTTGCCAGCTTGGTCAGGGCAAAGCCCAGTTCtt                                  <  1:2421570/150‑1 (MQ=255)
                                                                                                                    gggCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGAATCCAGCAGGCTTGCCAGCTTGGTCAGGGCAAAGCCCAGTTCtt                                  <  2:1310365/150‑1 (MQ=255)
                                                                                                                    gggCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTAGGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGAATCCAGCAGGCTTGCCAGCTTGGTCAGGGCAAAGCCCAGTTCtt                                  <  2:3046448/150‑1 (MQ=255)
                                                                                                                     ggCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGAATCCAGCAGGCTTGCCAGCTTGGTCAGGGCAAAGCCCAGTTCTTg                                 <  2:2374406/150‑1 (MQ=255)
                                                                                                                     ggCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGAATCCAGCAGGCTTGCCAGCTTGGTCAGGGCAAAGCCCAGTTCTTg                                 >  1:2180958/1‑150 (MQ=255)
                                                                                                                              cAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGAATCCAGCAGGCTTGCCAGCTTGGTCAGGGCAAAGCCCAGTTCTTGGCGGCTggg                        >  2:4720980/1‑150 (MQ=255)
                                                                                                                                 gACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGAATCCAGCAGGCTTGCCAGCTTGGTCAGGGCAAAGCCCAGTTCTTGGCGGCTGGGTGc                     >  1:744642/1‑150 (MQ=255)
                                                                                                                                 gACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGAATCCAGCAGGCTTGCCAGCTTGGTCAGGGCAAAGCCCAGTTCTTGGCGGCTGGGTGc                     >  1:4146757/1‑150 (MQ=255)
                                                                                                                                 gACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGAATCCAGCAGGCTTGCCAGCTTGGTCAGGGCAAAGCCCAGTTCTTGGCGGCTGGGTGc                     >  1:4146665/1‑150 (MQ=255)
                                                                                                                                 ataaaaaaaCCCAGGGCAACTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGAATCCAGCAGGCTTGCCAGCGGGGTCAGGGCAAAGCCCAGTTCTTGGCGGCTGGGTGc                     >  2:4480022/7‑150 (MQ=255)
                                                                                                                                    cAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGAATCCAGCAGGCTTGCCAGCTTGGTCAGGGCAAAGCCCAGTTCTTGGCGGCTGGGTGCCAt                  <  1:776096/150‑1 (MQ=255)
                                                                                                                                     aaaaaCCCAGGGCAACTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGAATCCAGCAGGCTTGCCAGCTTGGTCAGGGCAAAGCCCAGTTCTTGGCGGCTGGGTGCCATc                 >  1:1252475/3‑150 (MQ=255)
                                                                                                                                     aGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGAATCCAGCAGGCTTGCCAGCTTGGTCAGGGCAAAGc                                           <  1:3512060/124‑1 (MQ=255)
                                                                                                                                     aGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGAATCCAGCAGGCTTGCCAGCTTGGTCAGGGCAAAGc                                           >  2:3512060/1‑124 (MQ=255)
                                                                                                                                      gAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGAATCCAGCAGGCTTGCCAGCTTGGTCAGGGCAAAGCCCAGTTCTTGGCGGCTGGGTGCCATCa                >  1:3382980/1‑150 (MQ=255)
                                                                                                                                         aCCCAGGGCAACTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGAATCCAGCAGGCTTGCCAGCTTGGTCAGGGCAAAGCCCAGTTCTTGGCGGCTGGGTGCCATCAggg             >  1:2618564/1‑150 (MQ=255)
                                                                                                                                          cccAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACGACCGGGGCTCCGGCGAATCCAGCAGGCTTGCCAGCTTGGTCAGGGCAAAGCCCAGTTCTTGGCGGCTGGGTGCCATCAGGGc            <  1:1780627/150‑1 (MQ=255)
                                                                                                                                          cccAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGAATCCAGCAGGCTTGCCAGCTTGGTCAGGGCAAAGCCCAGTTCTTGGCGGCTGGGTGCCATCAGGGc            >  1:230686/1‑150 (MQ=255)
                                                                                                                                          cccAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGAATCCAGCAGGCTTGCCAGCTTGGTCAGGGCAAAGCCCAGTTCTTGGCGGCTGGGTGCCATCAGGGc            >  2:3795805/1‑150 (MQ=255)
                                                                                                                                          cccAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGAATCCAGCAGGCTTGCCAGCTTGGTCAGGGCAAAGCCCAGTTCTTGGCGGCTGGGTGCCATCAGGGc            >  1:4212064/1‑150 (MQ=255)
                                                                                                                                           ccAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGAATCCAGCAGGCTTGCCAGCTTGGTCAGGGCAAAGCCCAGTTCTTGGCGGCTGGGTGCCATCAGGGCg           >  1:1344562/1‑150 (MQ=255)
                                                                                                                                           ccAGGGCAACTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGAATCCAGCAGGCTTGCCAGCTTGGTCAGGGCAAAGCCCAGTTCTTGGCGGCTGGGt                       <  1:4310095/138‑1 (MQ=255)
                                                                                                                                           ccAGGGCAACTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGAATCCAGCAGGCTTGCCAGCTTGGTCAGGGCAAAGCCCAGTTCTTGGCGGCTGGGt                       >  2:4310095/1‑138 (MQ=255)
                                                                                                                                           caaGGGCAACTGCCCTGGGTTTTTTTATGGGAAGCAACGCTTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGAATACAGCAGGCTTGCCAGCTTGGTCAGGGCAAAGCCCAGTTCTTGGCGGCTGGGTGCCATCAGGGCg           >  2:2018797/3‑150 (MQ=255)
                                                                                                                                            cAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGAATCCAGCAGGCTTGCCAGCTTGGTCAGGGCAAAGCCCAGTTCTTGGCGGCTGGGTGCCATCAGGGCGa          >  1:3124195/1‑150 (MQ=255)
                                                                                                                                             aGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGAATCCAGCGGGCTTGCCAGCTTGGTCAGGGCAAAGCCCAGTTCTTGGCGGCTGGGTGCCATCAGGGCGAt         <  1:117852/150‑1 (MQ=255)
                                                                                                                                             aGGGCAGTTGCCCTGGGATTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGAATCCAGCAGGCTTGCCATCTTGGTCAGGGCAAAGCCCAGTTCTTGGCGGCTGGGTGCCATCAGGGCGAt         <  1:4364237/150‑1 (MQ=255)
                                                                                                                                             aGGGCAGTTGCCCTGGGATATTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGAATCCAGCAGGCTTGCCAGCTTGGTCAGGGCAAAGCCCAGTTCTTGGCGGCTGGGTGCCATCAGGGCGAt         <  2:108326/150‑1 (MQ=255)
                                                                                                                                              gggCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGAATCCAGCAGGCTTGCCAGCTTGGTCAGGGCAAAGCCCAGTTCTTGGCGGCTGGGTGCCATCAGGGCGATg        <  1:4872936/150‑1 (MQ=255)
                                                                                                                                               ggCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGAATCCAGCAGGCTTGCCAGCTTGGTCAGGGCAAAGCCCAGTTCTTGGCGGCTGGGTGCCATCAGGGCGATgc       >  2:2662994/1‑150 (MQ=255)
                                                                                                                                                  aGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGAATCCAGCAGGCTTGCCAGCTTGGTCAGGGCAAAGCCCAGTTCTTGGCGGCTGGGTGCCATCAGGGCGATGCGca    <  1:2187204/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CAGGCCCAGCGCCCGCAACAGGCCCCTCAGCGCCCAGCGCCGCAACAGCCTGCGCCGCAGCCGGCGGCTGACTTCGACAGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGAATCCAGCAGGCTTGCCAGCTTGGTCAGGGCAAAGCCCAGTTCTTGGCGGCTGGGTGCCATCAGGGCGATGCGCACA  >  NC_002947/571473‑571770

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