Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0084633,515,8760TA30.7% 543.1 / 71.6 282intergenic (+630/+254)PA14_RS16005/PA14_RS16010AAA family ATPase/DUF4354 family protein
Reads supporting (aligned to +/- strand):  ref base T (48/139);  new base A (82/1);  total (142/140)
Fisher's exact test for biased strand distribution p-value = 1.31e-33
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 289 total aligned reads displayed.

CCTAAGCCATAGCCCTTCCCCATACCTGAACGGCACACCTATTGCCACCCGCAAAACCACAACCTCCCCCGCTGCTAACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGG‑CCTTCACCCCAACCCTC‑TCCCGCACGCGGGAGAGGGGGCCGGGCGTGCCTTGTGGGTCGAGCGCCTCGCTCGGGCGCAGCC‑TGGAACGATGTCTACCCTCCCCTAGGCCCCGGGTATCTGGAATTCGGCCCCCAGCCCAGGCAACCAAAAGCACCGGAAACGCGC  >  NC_008463/3515729‑3516025
                                                                                                                                                    |                                                                                                                                                       
 ctAAGCCATAGCCCTTCCCCATACCTGAACGGCACACCTATTGCCACCCGCAAAACCACAACCTCCCCCGCTGCTAACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGG‑CCCTCACCCCAACCCtc‑tc                                                                                                                                                     <  1:10354059/150‑1 (MQ=255)
 ctAAGCCATAGCCCTTCCCCATACCTGAACGGCACACCTATTGCCACCCGCAAAACCACAACCTCCCCCGCTGCTAACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGG‑CCCTCACCCCAACCCtc‑tc                                                                                                                                                     <  2:8824197/150‑1 (MQ=255)
  tAAGCCATAGCCCTTCCCCATACCTGAACGGCACACCTATTGCCACCCGCAAAACCACAACCTCCCCCGCTGCTAACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGG‑CCCACCCCCCAACCCAC‑Tcc                                                                                                                                                    >  2:7010899/1‑150 (MQ=255)
  tAAGCCATAGCCCTTCCCCAGACCTGAACGGCACACCTATTGCCACCCGCAAAACCACAACCTCCCCCGCTGCTAACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGG‑CCCTCACCCCAACCCTC‑Tcc                                                                                                                                                    <  2:6917441/150‑1 (MQ=255)
       cATAGCCCTTCCCCATACCTGAACGGCACACCTATTGCCACCCGCAAAACCACAACCTCCCCCGCTGCTAACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGG‑CCCTCACCCCAACCCAC‑TACCCCAc                                                                                                                                               >  1:8775838/1‑150 (MQ=255)
        aTAGCCCTTCCCCATACCTGAACGGCACACCTATTGCCACCCGCAAAACCACAACCTCCCCCGCTGCTAACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTTG‑CCCACACCCCAAACCAC‑Acccccccg                                                                                                                                              >  2:786792/1‑145 (MQ=255)
         tAGCCCTTCCCCATACCTGAACGGCACACCTATTGCCACCCGCAAAACCTCAACCTCCCCCGCTGCTAACCGCACTAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGG‑GCCTCACCCCAACCCAC‑TCCCCCCcgc                                                                                                                                             >  1:3532672/1‑150 (MQ=255)
         tAGCCCTTCCCCATACCTGAACGGCACACCTATTGCCACCCGCAAAACCACAACCTCCCCCGCTGCTAACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGG‑GCCTCAACCCAAACCAA‑TACCGCAagc                                                                                                                                             >  2:3115677/1‑150 (MQ=255)
                  cccATCCCTGAACGGCACACCTATTGCCACCCGCAAAACCACCACCTCCCCCGCAGCTAACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGG‑CCCTCACCCCAACCCTC‑TCCCGCACGCGGGAGAggg                                                                                                                                    <  2:10186143/150‑1 (MQ=255)
                  cccATACCTGAACGGCACACCTATTGCCACCCGCAAAACCACAACCTCCCCCGCTGCTAACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGG‑CCCTCACCCCAACCCTC‑TCCCGCACGCGGGAGAggg                                                                                                                                    <  2:2658263/150‑1 (MQ=255)
                   ccATACCTGAACGGCACACCTATTGCCACCCGCAAAACCACAACCTCCCCCGCTGCTAACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGG‑GCCAAACCCCAACACGA‑ACCCCCACGCGGGAGAgggg                                                                                                                                   >  1:1197523/1‑150 (MQ=255)
                    cATACCTGAACGGCACACCTATTGCCACCCGCAAAACCACAACCTCCCCCGCTGCTAACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGG‑CCCTCACCCCCACCCAC‑ACCCGCACGCCGGCGAggggg                                                                                                                                  >  1:6588527/1‑150 (MQ=255)
                      tACCTGAACGGCACACCTATTGCCACCCGCATAACCACAACCTCCCCCGCTGCTAACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGG‑CCCTCACCCCAACCCTC‑TCCCGCACGCGGGAGAGGGGGcc                                                                                                                                <  2:5245443/150‑1 (MQ=255)
                      tACCTGAACGGCACACCTATTGCCACCCGCAAAACCACAACCTCCCCCGCTGCTAACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGG‑CCCTCACCCCAACCCTC‑TCCCGCACGCGGGAGAGGGGGcc                                                                                                                                <  1:4782863/150‑1 (MQ=255)
                      tACCTGAACGGCACACCTATTGCCACCCGCAAAACCACAACCTCCCCCGCTGCTAACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGGGTTGG‑CCCGCACCCCCACCCAA‑AACCGCACCCGGGAGAtgggggc                                                                                                                                >  2:3497273/1‑148 (MQ=255)
                               ggCACACCTATTGCCACCCGCAAAACCACAACCTCCCCCGCTGCTAACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGG‑CCCGCACCCCCAACCTC‑TCCCCCACGCGGGCGAGGGGGGCGGGGGTggc                                                                                                                       >  1:95307/1‑148 (MQ=255)
                                gCACACCTATTGCCACCCGCAAAACCACAACCTCCCCCGCTGCTAACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGG‑GCCAACCCCCCACCCAC‑ACCCCCACCCGGGGGAGGGGGCCggggggggcg                                                                                                                      >  1:6035374/1‑143 (MQ=255)
                                 cacaCCTATTGCCACCCGCAAAACCACAACCTCCCCCGCTGCTAACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGG‑GCCGCAACCCCACCCAC‑AACCCCCCCCGGGGGGGGGGggcgggcggggcgt                                                                                                                     >  1:1954137/1‑144 (MQ=255)
                                 cacaCCTATTGCCACCCGCAAAACCACAACCTCCCCCGCTGCTAACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGG‑CCCTCACCCCAACCCTC‑TCCCGCACGCGGGAGAGGGGGCCGGGCGTGCCtt                                                                                                                     <  1:773703/150‑1 (MQ=255)
                                  acacCTATTGCCACCCGCAAAACCACAACCTCCCCCGCTGCTAACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGGTGG‑CCCACACCCCAACCCAA‑GCCCCCCCCCCGGAGGGGGGGGCGGGCGTggcgag                                                                                                                    >  1:1073929/1‑145 (MQ=255)
                                     ccTATTGCCACCCGCACACCCCCACCCTCCCCCGCTGCTAACCGCAAGAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGG‑CCCTCACCCCACCCCTC‑TCCCGCACGCGGGAGAGGGGGCCGGGCGTGCCTTGTgg                                                                                                                 <  1:6798476/150‑1 (MQ=255)
                                       tATTGCCACCCGCAAAACCACACCCTCCCCCGCTGCTACCCGAATTAGCTGCGGCTTGTGTGTCGCGTTTCGGTTTTTGGATTGTGGTGTTGG‑CCCTCACCCCACCCCTC‑TCCCGCACGGGGTATTGGGGGCGGGGCGTGCCTTGTGGGt                                                                                                               <  2:2685339/150‑1 (MQ=255)
                                        cttGCCACCCGCAAAACCACAACCTCCCCCGCTGCTAACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTGTTGGAATGTGGTGTTGG‑CCCTCACCCCCTCCCTC‑GCCCGCACGCGGGGGCGGGGGCCGGGGGGGCCGGGTggggc                                                                                                              >  2:197241/2‑148 (MQ=255)
                                        aTTGCCACCCGCAAAACCACAACCTCCCCCGCTGCTCACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGG‑CCCTCACCCCAACCCTC‑TCCCGCACGCGGGAGAGGGGGCCGGGCGTGCCTTGTGGGTc                                                                                                              <  1:6109992/150‑1 (MQ=255)
                                        aTTGCCACCCGCAAAACCACAACCTCCCCCGCTGCTAACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGGCCCAACACCCAAACCC‑C‑ACCCCCACGCGGGAGAGGGGGGCCGGCGTGCCGTGTGGGTc                                                                                                              >  2:10702755/1‑150 (MQ=255)
                                         ttGCCACCCGCAAAACCACAACCTCCCCCGCTGCTAACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGG‑CCCTCACCCCAACCCTC‑TCCCGCACGCGGGAGAGGGGGCCGGGCGTGCCTTGTGGGTCg                                                                                                             <  1:8120628/150‑1 (MQ=255)
                                           gCCACCCGCAAAACCACAACCTCCCCCGCTGCTAACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGG‑GCCGCAACCCAAAACAC‑TACCCCCCGCGGGAGAGGGGGCCGGGCGTGGCGGGTGGGTCGAg                                                                                                           >  1:4211603/1‑150 (MQ=255)
                                            ccACCCGCAAAACCACAACCTCCCCCGCGGCTAACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGG‑CCCGCACCCCCAACCAC‑TCCCGCACGCGGGAGAGGGGGCCGGGCGTGGCGTGTGgggcgggc                                                                                                          >  2:6485805/1‑150 (MQ=255)
                                              aCCCGCAAAACCACAACCTCCCCCGCTGCTAACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGG‑CCCACACCCCAACCCAC‑TCCCCCACGCGGGGGAGGGGGCCGGGCGTGCCGTGTGGGTCGAgcgc                                                                                                        >  2:1467522/1‑150 (MQ=255)
                                                ccGCAAAACCACAACCTCCCCCGCTGCTAACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGG‑CCCTCACCCCCACCCTC‑TCCCGCACGCGGGAGAGGGGGCCGGGCGTGCCTTGTGGGTCGAGCGCct                                                                                                      <  2:325121/150‑1 (MQ=255)
                                                ccGCAAAACCACAACCGCCCCCGCTGATAACCGCAATAGCTGGGGATTGAGTGTCGCGGATCGTTTTTTGGAATGTGGTGTTGG‑CCCTCACCCCAACCCTA‑TCCCGCACGCGGGTGAGGGGGCCGGGCGTGCCTTGTGGGTCGAGCGCct                                                                                                      <  1:9859747/150‑1 (MQ=255)
                                                 cGCAAAACCACAACCTCCCCCGCTGCTAACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGG‑CCCTCACCCCAACCCAC‑TCCCGCCCGCGGGGGAGGGGGGCGGGCGGGGCGTGTGGGTCGAGCGccac                                                                                                     >  1:7919787/1‑148 (MQ=255)
                                                  gCAAAACCACAACCTCCCCCGCTGCTAACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGG‑CCCACACCCCCACCCAC‑GCCCCCCCCCCGGCCAGGGGGGCGGGCGGGGCGTGGGGGTAGGGCGCctcg                                                                                                    >  2:8309912/1‑150 (MQ=255)
                                                      aaCCACAACCTCCCCCGCTGCTAACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGG‑GCCCCACCCCAACCCAC‑ACCCCCACGCGGGGGAGGGGGGCGGGGGGGGCGTGTGGGGCGAGCGCcacgctcg                                                                                                >  1:10086398/1‑150 (MQ=255)
                                                        ccACAACCTCCCCCGCTGCTAACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGG‑CCCTCACCCCAACCCGC‑TACCGCACGCGGGTGAGGGGGCCGGGCGTGCCTTGTGGGTCGAGCGCCTCGCGCggg                                                                                              >  2:7375522/1‑150 (MQ=255)
                                                          acaACCTCCCCCGCTGCTAACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGG‑GCCGCACCCCCACCCAC‑AACCCCCCCCCGGAGTGGGGGCCGGGGGGGCCGTGTGGGTCGAGCGCCACGCTCggggg                                                                                            >  2:9160774/1‑148 (MQ=255)
                                                           caACCTCCCCCGCTGCTAACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGG‑GCCGCACCCCCAACCTC‑TACCGCACGCCGGTGTGGGGGGCGGGCGTGGCGGGTGGGGCGAGCGCCGCGCGCGGgggc                                                                                           >  2:9811986/1‑150 (MQ=255)
                                                           caACCTCCCCCGCTGCTAACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGG‑GCCCAACCCCAAACCAA‑ACCCCCGCCCGGGGGTGGGGGGGGGGGGTGGCGTGTTGGGGGGGCCCCACGCACGGgcgc                                                                                           >  1:8862919/1‑150 (MQ=255)
                                                           caACCTCCCCCGCTGCTAACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGG‑CCCTCCCCCCAACCCTC‑TCCCGCACGCGGGAGAGGGGGCCGGGCGTGCCTTGTGGGTCGAGCGCCTCGCTCGGgcgc                                                                                           <  1:2690406/150‑1 (MQ=255)
                                                           caACCTCCCCCGCTGCTAACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGG‑CCCACAACCCAACCCAC‑AACCCCACCCGGGAGGGGGGGGCGGGGGTTCGTTGTTGGTCGTTCGCCTCGCACgggggg                                                                                           >  1:3293661/1‑147 (MQ=255)
                                                            aaCCTCCCCCGCTGCTAACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGG‑CCCTCACCCCAAACCGC‑TCCCGCCCCCGGGAGAGGGGGCCGGGCGTGCCGTGTGGGTCGAGCGCCACGCTCGGGCGCa                                                                                          >  2:10271520/1‑150 (MQ=255)
                                                                tCCCCCGCTGCTAACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGG‑GGCGCACCCCCACCCAC‑GACCCCCCCCGGGGGGGGGGGGCGGGCGGGGCGGGTGGGGCGTGCGCCGCGCACGGGGGCacgca                                                                                      >  2:1823462/1‑146 (MQ=255)
                                                                 cccccGCTGCTAACCGCAATAGCTGCGGCGTGAGTGTCGCGGATCGGTTTGTGGAATGTGGTGTTGG‑GCCGACCCCCAACCCTC‑TCCCCCCAGCGGGAGAGGGGGGCGGGCGGGGCGGGGGGGTCGGGCGCCTCGCTCGGGCCGAGccgg                                                                                     >  2:10424409/1‑148 (MQ=255)
                                                                 cccccGCTGCTAACCGCAATAGATGCGGCTTGAGTGTCGAGGATCGGTTTTTGGAATGTGGTGTTGG‑CACTCAACCCAACACTC‑GACCGCAAGCCGGAGCTGGGGACAGGCGTGGCGTGTGGTGAGAGCGCCTCGCGAGGGCGCAGccag                                                                                     >  2:8257863/1‑148 (MQ=255)
                                                                  ccccGCTGCTAACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGG‑CCCTCACCCCAACCCTC‑TCCCGCACGCGGGAGAGGGGGCCGGGCGTGCCTTGTGGGTCGAGCGCCTCGCTCGGGCGCAGCC‑Tgg                                                                                   <  1:3944883/150‑1 (MQ=255)
                                                                   cccGCTGCTAACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGG‑CCCTCACCCCAACCCTC‑TCCCGCACGCGGGAGAGGGGGCCGGGCGTGCCTTGTGGGTCGAGCGCCTCGCTCGGGCGCAGCC‑TGGa                                                                                  <  2:1641462/150‑1 (MQ=255)
                                                                     cGCTGCTAACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGG‑CCCTCACCCCAACCCTC‑TCCCGCACGCGGGAGAGGGGGCCGGGCGTGCCTTGTGGGTCGAGCGCCTCGCTCGGGCGCAGCC‑TGGAAc                                                                                <  1:5612523/150‑1 (MQ=255)
                                                                      gctgctAACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGG‑GCCGAACCCCAACACAC‑ACCCGCACGCGGGGGAGGGGGCCGGGCGTGCCGTGTGGGGCGAGCGCCACGCTCGGGCGCCGCC‑TGGACCg                                                                               >  2:4829620/1‑150 (MQ=255)
                                                                       ctgctAACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGG‑CCCACCCCCCAAAACGA‑ACCCGCACCCGGGCGAGGGGGGCGGGCGGGGCGGGTGGGTCGGGCGCCGCGCTCGGGCGCAACC‑GGGGACgg                                                                              >  1:10417173/1‑149 (MQ=255)
                                                                        tgctAACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGG‑GCCGCACCCCCACCCAC‑AACCCCCCGCCGGGGAGGGGGCCGGGCGGGCCGGGTGGGGGGAGCGCCGCCCAACGGCGCAGCCAGGGAACGa                                                                              >  1:744294/1‑150 (MQ=255)
                                                                         gctAACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGG‑CCCTCACCCCAACCCTC‑TCCCGCACGCGGGAGAGGGGGCCGGGCGTGCCTTGTGGGTCGAGCGCCTCGCTCGGGCGCAGCC‑TGGAACGATg                                                                            >  2:7406362/1‑150 (MQ=255)
                                                                          ctAACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGG‑CCCACCCCCCCAACCAC‑GACCCCACCCCGGGGGTGGGGGCGGGCGGGGCGTGTGGGGCGGGGCCCCCGCAGGGGGGCACCC‑GGGaaacctgt                                                                           >  1:637702/1‑144 (MQ=255)
                                                                           tAACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGG‑CCCTCACCCCAACCCTC‑TCCCGCACGCGGGAGAGGGGGCCGGGCGTGCCTTGTGGGTCGAGCGCCTCGCTCGGGCGCAGCC‑TGGAACGATGTc                                                                          <  2:557335/150‑1 (MQ=255)
                                                                                 cAATCGCTGCGGCTTGCGTGTCGCCGATCGGATTATGGAATGTGGTGTTGG‑CCCACACCCCAACCCTC‑TCCCGCACGCGGGAGAGGGGGCCGGGCGTGCCGTGTGGGTCGAGCGCCTCGCTCGGGCGCAGCC‑TGGAACGATGTCTACCCt                                                                    <  2:8690130/150‑1 (MQ=255)
                                                                                   aTAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGG‑CCCTCACCCCCAACCAC‑TACCCCACGCCGGGGAGGGGGGCGGGCGTGCCTTGTGGGTGGAGCGCCGCCCGCGGGCGCCGCC‑CGGTACGGTGTATAaccacc                                                                  >  2:1676802/1‑150 (MQ=255)
                                                                                    tAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGG‑CCCTCACCCCAACCCTC‑TCCCGCACGCGGGAGAGGGGGCCGGGCGTGCCTTGTGGGTCGAGCGCCTCGCTCGGGCGCAGCC‑TGGAACGATGTCTACCCTccc                                                                 <  1:10241238/150‑1 (MQ=255)
                                                                                      gCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGG‑CCCTCACCCCAACCCTC‑TCCCGCACGCGGGAGAGGGGGCCGGGCGTGCCTTGTGGGTCGAGCGCCTCGCTCGGGCGCAGCC‑TGGAACGATGTCTACCCTCCCCt                                                               <  1:8558634/150‑1 (MQ=255)
                                                                                      gCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGG‑CCCTCACCCCAACCCTC‑TCCCGCACGCGGGAGAGGGGGCCGGGCGTGCCTTGTGGGTCGAGCGCCTCGCTCGGGCGCAGCC‑TGGAACGATGTCTACCCTCCCCt                                                               <  2:5882101/150‑1 (MQ=255)
                                                                                       cTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGG‑CCCTCACCCCAACCCTC‑TCCCGCACGCGGGAGAGGGGGCCGGGCGTGCCTTGTGGGTCGAGCGCCTCGCGCGGGCGCAGCC‑TGGAACGATGTCTACCCTCCCCTa                                                              <  2:397036/150‑1 (MQ=255)
                                                                                         gCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGG‑CCCTCACCCCAACCCTC‑TCCCGCACGCGGGAGAGGGGGCCGGGCGTGCCTTGTGGGTCGAGCGCCTCGCTCGGGCGCAGCC‑TGGAACGATGTCTACCCTCCCCTAgg                                                            <  2:683426/150‑1 (MQ=255)
                                                                                             cTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGG‑GCCTCACCCCAAACCTCAACCCCCCCACGGGTTGGGGGGCC‑GGCGTGGCATGTGGGTGGAGCGCCGCGCAAGGGCGCCGCC‑CGGGACAGTGTATAACCGCCCCAAGGcccc                                                        >  2:7276648/1‑150 (MQ=255)
                                                                                              ttGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGG‑CCCTCACCCCAACCCTC‑TCCCGCACGCGGGAGAGGGGGCCGGGCGTGCCTTGTGGGTCGAGCGCCTCGCTCGGGCGCAGCC‑TGGAACGATGTCTACCCTCCCCTAGGCCCCg                                                       <  2:9546680/150‑1 (MQ=255)
                                                                                              ttGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGG‑CCCACACCCCAACCCAC‑TCCCGCACCCGGGAGAGGGGGCCCGGCGTGCCGTGTGGGTCGAGCGCCACGCACGGGCGCAGCC‑GGGAACGATGTCTACCCACCCCCAGGCCCCg                                                       >  2:1338500/1‑150 (MQ=255)
                                                                                                  gtgtCGCGGATCGGTTTTTGGAATGTGGTGTTGG‑CCCTCACCCCAACCCTC‑TCCCGCACGCGGGAGAGGGGGCCGGGCGTGCCTTGTGGGTCGAGCGCCTCGCTCGGGCGCAGCC‑TGGAACGATGTCTACCCTCCCCTAGGCCCCGGGTa                                                   <  2:901618/150‑1 (MQ=255)
                                                                                                  gtgtCGCGGATCGGTTTTTGGAATGTGGTGTTGG‑CCCTCACCCCAACCCTC‑TCCCGCACGCGGGAGAGGGGGCCGGGCGTGCCTTGTGGGTCGAGCGCCTCGCTCGGGCGCAGCC‑TGGAACGATGTCTACCCTCCCCTAGGCCCCGGGTa                                                   <  2:5039666/150‑1 (MQ=255)
                                                                                                    gtCGCGGATCGGTTTTTGGAATGTGGTGTTGG‑CCCACAACCCAACCCAC‑TCCCCCACGCGGGGGAGGGGGCCGGGCGGGCCGTGTGGGTCGAGCGCCACGCACGGGCGCCGCC‑GGGGACGCTGTCTACCCACCCCAAGGCCCCCGGGATc                                                 >  2:2304953/1‑150 (MQ=255)
                                                                                                     tCGCGGATCGGTTTTTGGAATGTGGTGTTGG‑CCCACACCCCAACCCGC‑TCCCCCACGCCGGGGAGGGGGGCCGGGGGGCCATGTGGGTGGAGCGCCACGCTAGGGCGCACAC‑GGGAACCATGTCTAACCACCCCAAGGGCCCGGGTATCt                                                >  2:4349592/1‑150 (MQ=255)
                                                                                                          gATCGGTTTTTGGAATGTGGTGTTGG‑CCCTCACCCCAACCCTC‑TCCCGCACGCGGGAGAGGGGGCCGGGCGTGCCTTGTGGGTCGAGCGCCTCGCTCGGGCGCAGCC‑TGGAACGATGTCTACCCTCCCCTAGGCCCCGGGTATCTGGAAt                                           <  2:569876/150‑1 (MQ=255)
                                                                                                            tCGGTTTTTGGAATGTGGTGTTGG‑CCCTCACCCCACCCCTC‑TCCCGCACGCGGGAGAGGGGGCCGGGCGTGCCTTGTGGGTCGAGCGCCTCGCTCGGGCGCAGCC‑TGGAACGATGTCTACCCTCCCCTAGGCCCCGGGTATCTGGAATTc                                         <  1:2882886/150‑1 (MQ=255)
                                                                                                            tCGGTTTTTGGAATGTGGTGTTGG‑CCCTCACCCCAACCCTC‑TCCCGCACGCGGGAGAGGGGGCCGGGCGTGCCTTGTGGGTCGAGCGCCTCGCTCGGGCGCAGCC‑TGGAACGATGTCTACCCTCCCCTAGGCCCCGGGTATCTGGAATTc                                         <  2:8508169/150‑1 (MQ=255)
                                                                                                                tttttGGACTGTGGTGTTGG‑CCCTCACCCCAACCCTC‑TCCCGCACGCGGGAGAGGGGGCCGGGCGTGCCTTGTGGGTCGAGCGCCTCGCTCGGGCGCAGCC‑TGGAACGATGTCTACCCTCCCCTAGGCCCCGGGTATCTGGAATTCGGcc                                     <  1:7596861/150‑1 (MQ=255)
                                                                                                                   ttGGAATGTGGTGTTGG‑CCCTCACCCCAACCCTC‑TCCCGCACGCGGGAGAGGGGGCCGGGCGTGCCTTGTGGGTCGAGCGCCTCGCTCGGGCGCAGCC‑TGGAACGATGTCTACCCTCCCCTAGGCCCCGGGTATCTGGAATTCGGccccc                                  <  1:4439772/150‑1 (MQ=255)
                                                                                                                       aaTGTGGTGTTGG‑CCCTCACCCCAACCCTC‑TCCCGCACGCGGGAGAGGGGGCCGGGCGTGCCTTGTGGGTCGAGCGCCTCGCTCGGGCGCAGCC‑TGGAACGATGTCTACCCTCCCCTAGGCCCCGGGTATCTGGAATTCGGCCcccagcc                              <  1:8497370/150‑1 (MQ=255)
                                                                                                                         tgtgGTGTTGG‑CCCTCACCCCAACCCTC‑TCCCGCACGCGGGAGAGGGGGCCGGGCGTGCCTTGTGGGTCGAGCGCCTCGCTCGGGCGCAGCC‑TGGAACGATGTCTACCCTCCCCTAGGCCCCGGGTATCTGGAATTCGGCCcccagccca                            <  2:5219887/150‑1 (MQ=255)
                                                                                                                         tgtgGAGTTGG‑CCCGCACCCCCACCCTC‑TCCCGCACGCGGGAGAGGGGGCCGGGCGTGCCTTGTGGGTCGAGAGCCTCGCTCGGGCGCAGCC‑TGGAACGATGTCTACCCTCCCCTAGGCCCCGGGTATCTGGAATTCGGCCcccagccca                            <  2:7615222/150‑1 (MQ=255)
                                                                                                                            ggtgTTGG‑CCCTCACCCCACCCCTC‑TCCCGCACGCGGGAGAGGGGGCCGGGCGTGCCTTGTGGGTCGAGCGCCTCGCTCGGGCGCAGCC‑TGGAACGATGTCTACCCTCCCCTAGGCCCCGGGTATCTGGAATTCGGCCCCCAGCCCAGGc                         <  2:10621000/150‑1 (MQ=255)
                                                                                                                             gtgtTGG‑CCCTCACCCCAACCCTC‑ACCCGCACGCGGGAGAGGGGGCCCGGCGTGCCGTGGGGGTCGAGCGCCTCGCTCGGGCGCAGCC‑GGGAACGATGTCTACCCACCCCTAGGCCCCGGGGATCTGGAATTCGGCCCCCAGCCCACGcc                        >  2:448411/1‑149 (MQ=255)
                                                                                                                               gtTGG‑CCCTCACCCCAACCCTC‑TCCCGCACGCGGGAGAGGGGGCCGGGCGTGCCTTGTGGGTCGAGCGCCTCGCTCGGGCGCAGCC‑TGGAACGATGTCTACCCTCCCCTAGGCCCCGGGTATCTGGAATTCGGCCCCCAGCCCAGGCacc                      >  2:5366205/1‑148 (MQ=255)
                                                                                                                                ttGG‑CCCTCACCCCCACCCTC‑TCCCGCACGCGGGAGAGGGGGCCGGGCGTGCCTTGTGGGTCGAGCGCCTCGCTCGGGCGCAGCC‑TGGAACGATGTCTACCCTCCCCTAGGCCCCGGGTATCTGGAATTCGGCCCCCAGCCCAGGCAAcc                     <  1:252750/150‑1 (MQ=255)
                                                                                                                                  gg‑CCCTCACCCCAACCCTC‑TCCCGCACGCGGGAGAGGGGGCCGGGCGTGCCGTGTGGGTCGAGCGCCTCGCTCGGGCGCAGCC‑TGGAACGATGTCTACCCTCCCCTAGGCCCCGGGTATCTGGAATTCGGCCCCCAGCCCAGGCAACCaa                   >  2:10032015/1‑150 (MQ=255)
                                                                                                                                   g‑CCCTCACCCCAACCCTC‑TCCCGCACGCGGGAGAGGGGGCCGGGCGTGCCTTGTGGGTCGAGCGCCTCGCTCGGGCGCAGCC‑TGGAACGATGTCTACCCTCCCCTAGGCCCCGGGTATCTGGAATTCGGCCCCCAGCCCAGGCAACCaaa                  <  2:5752915/150‑1 (MQ=255)
                                                                                                                                   g‑CCCTCACCCCAACCCTC‑TCCCGCACGCGGGAGAGGGGGCCGGGCGTGCCTTGTGGGTCGAGCGCCTCGCTCGGGCGCAGCC‑TGGAACGATGTCTACCCTCCCCTAGGCCCCGGGTATCTGGAATTCGGCCCCCAGCCCAGGCAACCaaa                  <  2:844360/150‑1 (MQ=255)
                                                                                                                                     ccctCACCCCAACCCTC‑TCCCGCACGCGGGAGAGGGGGCCGGGCGTGCCTTGTGGGTCGAGCGCCTCGCTCGGGCGCAGCC‑TGGAACGATGTCTACCCTCCCCTAGGCCCCGGGTATCTGGAATTCGGCCCCCAGCCCAGGCAACCaaaa                 >  1:5267869/4‑150 (MQ=255)
                                                                                                                                     ccctCACCCCAACCCTC‑TCCCGCACGCGGGAGAGGGGGCCGGGCGTGCCTTGTGGGTCGAGCGCCTCGCTCGGGCGCAGCC‑TGGAACGATGTCTACCCTCCCCTAGGCCCCGGGTATCTGGAATTCGGCCCCCAGCCCAGGCAACCaaaa                 <  2:2163598/147‑1 (MQ=255)
                                                                                                                                      cctCACCCCAACCCTC‑TCCCGCACGCGGGAGAGGGGGCCGGGCGTGCCTTGTGGGTCGAGCGCCTCGCTCGGGCGCAGCC‑TGGAACGATGTCTACCCTCCCCTAGGCCCCGGGTATCTGGAATTCGGCCCCCAGCCCAGGCAACCAAAAg                <  2:7130288/148‑1 (MQ=255)
                                                                                                                                        tCACCCCAACCCTC‑TCCCGCACGCGGGAGAGGGGGCCGGGCGTGCCTTGTGGGTCGAGCGCCTCGCTCGGGCGCAGCC‑TGGAACGATGTCTACCCTCCCCTAGGCCCCGGGTATCTGGAATTCGGCCCCCAGCCCAGGCAACCAAAAGCa              <  1:369132/150‑1 (MQ=255)
                                                                                                                                          aCCCCAACCCTC‑TCCCGCACGCGGGAGAGGGGGCCGGGCGTGCCTTGTGGGTCGAGCGCCTCGCTCGGGCGCAGCC‑TGGAACGATGTCTACCCTCCCCTAGGCCCCGGGTATCTGGAATTCGGCCCCCAGCCCAGGCACCCAAAAGCAcc            >  2:6221302/1‑150 (MQ=255)
                                                                                                                                          aCCCCAACCCTC‑TCCCGCACGCGGGAGAGGGGGCCGGGCGTGCCTTGTGGGTCGAGCGCCTCGCTCGGGCGCAGCC‑TGGAACGATGTCTACCCTCCCCTAGGCCCCGGGTATCTGGAATTCGGCCCCCAGCCCAGGCAACCAAAAGCAcc            >  1:42929/1‑150 (MQ=255)
                                                                                                                                           ccccAACCCTC‑TCCCGCACGCGGGAGAGGGGGCCGGGCGTGCCTTGTGGGTCGAGCGCCTCGCTCGGGCGCAGCC‑TGGAACGATGTCTACCCTCCCCTAGGCCCCGGGTATCTGGAATTCGGCCCCCAGCCCAGGCAACCAAAAGCACCg           >  2:2961624/1‑150 (MQ=255)
                                                                                                                                             ccAACCCTC‑TCCCGCACGCGGGAGAGGGGGCCGGGCGTGCCTTGTGGGTCGAGCGCCTCGCTCGGGCGCAGCC‑TGGAACGATGTCTACCCTCCCCTAGGCCCCGGGTATCTGGAATTCGGCCCCCAGCCCAGGCACCCAAAAGCACCGGa         >  2:2786371/1‑150 (MQ=255)
                                                                                                                                             ccAACCCTC‑TCCCGCACGCGGGAGAGGGGGCCGGGCGTGCCTTGTGGGTCGAGCGCCTCGCTCGGGCGCAGCC‑TGGAACGATGTCTACCCTCCCCTAGGCCCCGGGTATCTGGAATTCGGCCCCCAGCCCAGGCAACCAAAAGCACCGGa         >  1:4715429/1‑150 (MQ=255)
                                                                                                                                              cAACCCTC‑TCCCGCACGCGGGAGAGGGGGCCGGGCGTGCCTTGTGGGTCGAGCGCCTCGCTCGGGCGCAGCC‑TGGAACGATGTCTACCCTCCCCTAGGCCCCGGGTATCTGGAATTCGGCCCCCAGCCCAGGCAACCAAAAGCACCGGaa        <  1:6924317/150‑1 (MQ=255)
                                                                                                                                                 cccTC‑TCCCGCACGCGGGAGAGGGGGCCGGGCGTGCCTTGTGGGTCGAGCGCCTCGCTCGGGCGCAGCC‑TGGAACGATGTCTACCCTCCCCTCGGCCCCGGGTATCTGGAATTCGGCCCCCAGCCCAGGCAACCAAAAGCACCGGAAAcg     <  2:10537767/150‑1 (MQ=255)
                                                                                                                                                 cccTC‑TCCCGCACGCGGGAGAGGGGGCCGGGCGTGCCTTGTGGGTCGAGCGCCTCGCTCGGGCGCAGCC‑TGGAACGATGTCTACCCTCCCCTACGCCCCGGGTATCTGGAATTCGGCCCCCAGCCCAGGCAACCAAAAGCACCGGAAAcg     <  2:2796836/150‑1 (MQ=255)
                                                                                                                                                 cccTC‑TCCCGCACGCGGGAGAGGGGGCCGGGCGTGCCTGGTGGGTCGAGCGCCTCGCTCGGGCGCAGCC‑TGGAACGATGTCTACCCTCCCCTAGGCCCCGGGTCTCTGGAATTCGGCCCCCAGCCCAGGCAACCCAAAGCTCCGGAAAcg     >  2:3175876/1‑150 (MQ=255)
                                                                                                                                                  ccTC‑TCCCGCACGCGGGAGAGGGGGCCGGGCGTGCCTTGTGGGTCGAGCGCCTCGCTCGGGCGCAGCC‑TGGAACGATGTCTACCCTCCCCTAGGCCCCGGGTATCTGGAATTCGGCCCCCAGCCCAGGCAACCAAAAGCACCGGAAAcgc    >  1:2184142/1‑150 (MQ=255)
                                                                                                                                                  ccTC‑TCCCGCACGCGGGAGAGGGGGCCGGGCGTGCCTTGTGGGTCGAGCGCCTCGCTCGGGCGCAGCC‑TGGAACGATGTCTACCCTCCCCTAGGCCCCGGGTATCTGGAATTCGGCCCCCAGCCCAGGCAACCAAAAGCACCGGAAAcgc    <  1:10488784/150‑1 (MQ=255)
                                                                                                                                                  ccTC‑TCCCGCACGCGGGAGAGGGGGCCGGGCGTGCCTTGTGGGTCGAGCGCCTCGCTCGGGCGCAGCC‑TGGAACGATGTCTACCCTCCCCTAGGCCCCGGGTATCTGGAATTCGGCCCCCAGCCCAGGCAACCAAAAGCACCGGAAAcgc    <  2:2576724/150‑1 (MQ=255)
                                                                                                                                                   ctc‑tCCCGCACGCGGGAGAGGGGGCCGGGCGTGCCTTGTGGGTCGAGCGCCTCGCTCGGGCGCAGCC‑TGGAACGATGTCTACCCTCCCCTAGGCCCCGGGTATCTGGAATTCGGCCCCCAGCCCAGGCAACCAAAAGCACCGGAAAcgcg   <  2:4231212/150‑1 (MQ=255)
                                                                                                                                                    tc‑tcCCGCACGCGGGAGAGGGGGCCGGGCGTGCCTTGTGGGTCGAGCGCCTCGCTCGGGCGCAGCC‑TGGAACGATGTCTACCCTCCCCTAGGCCCCGGGTATCTGGAATTCGGCCCCCAGCCCAGGCAACCAAAAGCACCGGAAACgcgc  <  1:895356/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                       
CCTAAGCCATAGCCCTTCCCCATACCTGAACGGCACACCTATTGCCACCCGCAAAACCACAACCTCCCCCGCTGCTAACCGCAATAGCTGCGGCTTGAGTGTCGCGGATCGGTTTTTGGAATGTGGTGTTGG‑CCTTCACCCCAACCCTC‑TCCCGCACGCGGGAGAGGGGGCCGGGCGTGCCTTGTGGGTCGAGCGCCTCGCTCGGGCGCAGCC‑TGGAACGATGTCTACCCTCCCCTAGGCCCCGGGTATCTGGAATTCGGCCCCCAGCCCAGGCAACCAAAAGCACCGGAAACGCGC  >  NC_008463/3515729‑3516025

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