Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,041,195 C→G 9.3% intergenic (+173/‑489) PP_5680 → / → PP_4447 Cro/CI family transcriptional regulator/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,041,1950CG 9.3% 227.5 / 12.6 75intergenic (+173/‑489)PP_5680/PP_4447Cro/CI family transcriptional regulator/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (37/31);  new base G (5/2);  total (42/33)
Fisher's exact test for biased strand distribution p-value = 4.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.60e-01

AACACTTGTACTCCACACCCAGTCTCACATCTGACCCTGTCAGCGACGCTAATGACGGTGTTGATCCACAGAGGGTTAGGGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCTGGCCTGACCGCAGTAGAGC  >  NC_002947/5041047‑5041342
                                                                                                                                                     |                                                                                                                                                    
aaCACTTGTACTCCACACCCAGTCTCACATCTGACCCTGTCAGCGACGCTAATGACGGTGTTGATCCACAGAGGGTTAGGGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGgcgc                                                                                                                                                      >  2:346161/1‑149 (MQ=255)
  cacTTGTACTCCACACCCAGTCTCACATCTGACCCTGTCAGCGACGCTAATGACGGTGTTGATCCACAGAGGGTTAGGGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGgcgcgc                                                                                                                                                    <  1:393631/149‑1 (MQ=255)
     ttGTACTCCACACCCAGTCTCACATCTGACCCTGTCAGCGACGCTAATGACGGTGTTGATCCACAGAGGGTTAGGGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCt                                                                                                                                                  >  2:43672/1‑148 (MQ=255)
       gTACTCCACACCCAGTCTCACATCTGACGCTGTCAGCGACTCTAATGACGGTTTTGATCCAAAGAGGGTTAGGGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCtt                                                                                                                                                 <  3:159726/147‑1 (MQ=255)
          cTCCACACCCAGTCTCACATCTGACCCTGTCAGCGACGCTAATGACGGTGTTGATCCACAGAGGGTTAGGGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGcc                                                                                                                                            >  5:51526/1‑149 (MQ=255)
            ccACACCCAGTCTCACATCTGACCCTGGTCAGCGACGCTAATGACGGTGTTGATCCACAGAGGGTTAGGGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGccc                                                                                                                                           >  1:17665/1‑149 (MQ=255)
            ccACACCCAGTCTCACATCTGACCCTGTCAGCGACGCTAATGACGGTGTTGATCCACAGAGGGTTAGGGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCccg                                                                                                                                          >  3:191619/1‑149 (MQ=255)
              acacCCAGTCTCACATCTGACCCTGTCAGCGACGCTAATGACGGTGTTGATCCACAGAGGGTTAGGGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCccgcc                                                                                                                                        >  6:8605/1‑149 (MQ=255)
                acCCAGTCTCACATCTGACCCTGTCAGCGACGCTAATGACGGTGTTGATCCACAGAGGGTTAGGGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCccgcc                                                                                                                                        >  2:231315/1‑147 (MQ=255)
                 cccAGTCTCACATCTGACCCTGTCAGCGACGCTAATGACGGTGTTGATCCACAGAGGGTTAGGGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCTCGGCGGGCTTAAGGCGCGCCgatg                                                                                                                                    <  2:332621/149‑3 (MQ=255)
                  ccAGTCTCACATCTGACCCTGTCAGCGACGCTAATGACGGTGTTGATCCACAGAGGGTTAGGGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGGGCTTAAGGCGCGCCgat                                                                                                                                     >  7:171640/1‑147 (MQ=255)
                  ccAGTCTCACATCTGACCCTGTCAGCGACGCTAATGACGGTGTTGATCCACAGAGGGTTAGGGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGGGCTTAAGGCGCGCCgat                                                                                                                                     <  8:171640/148‑2 (MQ=255)
                  ccAGTCTCACATCTGACCCTGTCAGCGACGCTAATGACGGTGTTGATCCACAGAGGGTTAGGGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCgag                                                                                                                                    >  7:68216/1‑148 (MQ=255)
                  ccAGTCTCACATCTGACCCTGTCAGCGACGCTAATGACGGTGTTGATCCACAGAGGGTTAGGGATAGCCTTAAGGCGTAATTGTCTGAGGCCAGCACGTCTGCGATGAGGACATGACCTCTTCAGCGGCGCGCCGTAGGGCCGCCgag                                                                                                                                    >  6:338546/1‑148 (MQ=255)
                     gTCTCACATCTGACCCTGTCAGCGACGCTAATGACGGTGTTGATCCACAGAGGGTTAGGGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGagaaga                                                                                                                                <  5:8605/149‑1 (MQ=255)
                     gTCTCACATCTGACCCTGTCAGCGACGCTAATGACGGTGTTGATCCACAGAGGGTTAGGGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGagaaga                                                                                                                                <  1:346161/149‑1 (MQ=255)
                        tcACATCTGACCCTGTCAGCGACGCTAATGACGGTGTTGATCCACAGAGGGTTAGGGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGc                                                                                                                              >  8:76869/1‑148 (MQ=255)
                         cacaTCTGACCCTGTCAGCGACGCTAATGACGGTGTTGATCCACAGAGGGTTAGGGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGc                                                                                                                              >  7:159637/1‑147 (MQ=255)
                          acttCTGACCCTGTCAGCGACGCTAATGACGGTGTTGATCCACAGAGGGTTAGGGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGc                                                                                                                              <  7:66040/143‑1 (MQ=255)
                          acaTCTGACCCTGTCAGCGACGCTAATGACGGTGTTGATCCACAGAGGGTTAGGGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGc                                                                                                                              >  1:9823/1‑146 (MQ=255)
                                 aCCCTGTCAGCGACGCTAATGACGGTGTTGATCCACAGAGGGTTAGGGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGAcccc                                                                                                                    <  6:338770/149‑1 (MQ=255)
                                 aCCCTGTCAGCGACGCTAATGACGGTGTTGATCCACAGAGGGTTAGGGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGAc                                                                                                                       <  8:159637/146‑1 (MQ=255)
                                             gACGCTAATGACGGTGTTGATCCACAGAGGGTTAGGGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGAccc                                                                                                         <  1:43672/149‑1 (MQ=255)
                                              aCGCTAATGACGGTGTTGATCCACAGAGGGTTAGGGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGAcccc                                                                                                        >  6:185502/1‑149 (MQ=255)
                                                gCTAATGACGGTGTTGATCCACAGAGGGTTAGGGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCt                                                                                                       >  2:295320/1‑148 (MQ=255)
                                                  tAATGACGGTGTTGATCCACAGAGGGTTAGGGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGagaaga                                                                                                                                >  7:101091/1‑121 (MQ=255)
                                                  tAATGACGGTGTTGATCCACAGAGGGTTAGGGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGagaaga                                                                                                                                <  8:101091/121‑1 (MQ=255)
                                                             ttGATCCACAGAGGGTTAGGGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATgg                                                                                          <  2:17665/148‑1 (MQ=255)
                                                               gATCCACAGAGGGTTAGGGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCg                                                                                        <  1:231315/148‑1 (MQ=255)
                                                                   cacaGAGGGTTAGGGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTtc                                                                                   >  3:148854/1‑149 (MQ=255)
                                                                    acaGAGGGTTAGGGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTtc                                                                                   >  8:279582/1‑148 (MQ=255)
                                                                    acaGAGGGTTAGGGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTtc                                                                                   >  4:381026/1‑148 (MQ=255)
                                                                    acaGAGGGTTAGGGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCATAAGGCGGttagg                                                                                  <  6:51526/149‑4 (MQ=255)
                                                                      agagGGTTAGGGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTctct                                                                                <  5:14412/149‑1 (MQ=255)
                                                                          ggTTAGGGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTtctc                                                                                 >  3:57661/1‑144 (MQ=255)
                                                                          ggTTAGGGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTtctc                                                                                 <  4:57661/144‑1 (MQ=255)
                                                                          ggTTAGGGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCt                                                                            <  8:8261/149‑1 (MQ=255)
                                                                               gggATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCa                                                                        <  6:15971/148‑1 (MQ=255)
                                                                                ggATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGt                                                                      >  4:104252/1‑149 (MQ=255)
                                                                                  aTAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTct                                                                    <  7:76869/149‑1 (MQ=255)
                                                                                   tAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTctg                                                                   <  8:295038/149‑1 (MQ=255)
                                                                                      ccTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCg                                                                                        >  2:174416/1‑125 (MQ=255)
                                                                                        ttAAGGCGTAATTGTCTGATGCCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCtgctgc                                                               >  3:375534/1‑148 (MQ=255)
                                                                                           aGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAg                                                           <  2:9823/149‑1 (MQ=255)
                                                                                             gCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCa                                                                        >  7:307894/1‑134 (MQ=255)
                                                                                             gCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCa                                                                        <  8:307894/134‑1 (MQ=255)
                                                                                             gCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGc                                                         >  2:388435/1‑149 (MQ=255)
                                                                                                  aTTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATa                                                    <  7:369898/149‑1 (MQ=255)
                                                                                                   ttGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAg                                                   >  1:392942/1‑149 (MQ=255)
                                                                                                        tGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGCAGGCCGATAGCCGAg                                              >  6:271489/1‑149 (MQ=255)
                                                                                                          aTACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGCCCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTc                                            >  4:68197/1‑149 (MQ=255)
                                                                                                              cAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCTCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGcc                                                 <  4:218338/140‑1 (MQ=255)
                                                                                                              cAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCTCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGcc                                                 >  3:218338/1‑140 (MQ=255)
                                                                                                              cAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTccaac                                        <  7:201803/149‑1 (MQ=255)
                                                                                                                acacCTCTGCGATGAGGACATGACCTCTTCACCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTccaacc                                       >  3:19284/1‑148 (MQ=255)
                                                                                                                    ctctGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGCTAGCCGAGGCCCACCAATTc                                  >  1:182225/1‑149 (MQ=255)
                                                                                                                         cGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCTGATAGCCGAGTCCAACCAATTCTAt                               <  5:338546/147‑1 (MQ=255)
                                                                                                                         cGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTAtt                              <  3:104252/148‑1 (MQ=255)
                                                                                                                           aTGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTAttttt                           >  1:217594/1‑149 (MQ=255)
                                                                                                                           aTGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTAttttt                           >  3:72780/1‑149 (MQ=255)
                                                                                                                            tGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTc                          <  2:217594/149‑1 (MQ=255)
                                                                                                                            tGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTc                          >  3:243418/1‑149 (MQ=255)
                                                                                                                              aGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCg                         <  1:295320/148‑1 (MQ=255)
                                                                                                                              aGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCg                         <  5:185502/148‑1 (MQ=255)
                                                                                                                              aGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGa                        <  4:375534/149‑1 (MQ=255)
                                                                                                                                    tGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGa                                               >  4:335551/1‑120 (MQ=255)
                                                                                                                                    tGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCTGGc                  <  2:392942/149‑1 (MQ=255)
                                                                                                                                     tggctctTCTCGGCGGGCTTAAGGCGCGCCGATGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCTGGcc                 >  8:249695/4‑149 (MQ=38)
                                                                                                                                      aCCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCTg                    >  6:236529/1‑145 (MQ=255)
                                                                                                                                      aCCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCTg                    <  5:236529/145‑1 (MQ=255)
                                                                                                                                       ccTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCTGGCCt                >  4:103081/1‑148 (MQ=255)
                                                                                                                                        ctctTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCACCCAATAATATTTTTCGACCTGGCCTGa              >  6:30807/1‑149 (MQ=255)
                                                                                                                                         tctTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCTGGCCTGAc             >  3:361985/1‑149 (MQ=255)
                                                                                                                                            ttctCGGCGGGCTTAAGGCGCGCCGATGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCTGGCCTGACCg           >  7:360367/4‑149 (MQ=255)
                                                                                                                                             tctCGGCGGGCTTAAGGCGCGCCGATGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCTGGCCTGACCg           >  3:118277/3‑148 (MQ=255)
                                                                                                                                             tctCGGCGGGCTTAAGGCGCGCCGATGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCTGGCCTGACCg           >  3:112802/3‑148 (MQ=255)
                                                                                                                                               tCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCTGGCCTGACCGCAGt       <  8:68216/149‑1 (MQ=255)
                                                                                                                                                  gcgcgcCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCTGGCCTGACCGCAGTaga    <  3:381026/149‑1 (MQ=255)
                                                                                                                                                    gcgcCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCTGGCCTGACCGCAGTAGAGc  >  3:152855/1‑149 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
AACACTTGTACTCCACACCCAGTCTCACATCTGACCCTGTCAGCGACGCTAATGACGGTGTTGATCCACAGAGGGTTAGGGATAGCCTTAAGGCGTAATTGTCTGATACCAACACCTCTGCGATGAGGACATGACCTCTTCATCGGCGCGCCTTAAGCCCGCCGAGAAGAGCCAGAGACCCCGAGCAAGACCCCTAACCGCCTTATGGCGGTTTCTCTGCCTCCCAGTCTGCTGCGAAGGCCGATAGCCGAGTCCAACCAATTCTATTTTTCGACCTGGCCTGACCGCAGTAGAGC  >  NC_002947/5041047‑5041342

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