Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 660,945 A→C 5.6% intergenic (+74/‑93) speA → / → PP_0568 biosynthetic arginine decarboxylase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947660,9450AC 5.6% 606.6 / 17.1 178intergenic (+74/‑93)speA/PP_0568biosynthetic arginine decarboxylase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (75/93);  new base C (5/5);  total (80/98)
Fisher's exact test for biased strand distribution p-value = 7.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

Only 100 of 179 total aligned reads displayed.

ACTGCGCGTGAGCGGACCCAGTATCTGGATGCCTTGCGCCTGGGCTTGACCCGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAATGTTCAAAGGAGGACGATATGCAACGTCCAGGCTCCCATCGTTTACGTGGCGGTCGAGTGTCTGAATCCGGGCG  >  NC_002947/660797‑661093
                                                                                                                                                    |                                                                                                                                                    
aCTGCGCGTGAGCGGACCCAGTATCTGGATGCCTTGCGCCTGGGCTTGACCCGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAg                                                                                                                                                     <  2:2573080/150‑1 (MQ=255)
  tGCGCGTGAGCGGACCCAGTATCTGGATGCCTTGCGCCTGGGCTTGACCCGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGcc                                                                                                                                                   >  1:1627713/1‑150 (MQ=255)
       gTGAGCGGACCCAGTATCTGGATGCCTTGCGCCTGGGCTTGACCCGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCAc                                                                                                                                              <  1:3438222/150‑1 (MQ=255)
        tGAGCGGACCCAGTATCTGGATGCCTTGCGCCTGGGCTTGACCCGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCTTCGCCCTCata                                                                                                                                             <  1:3260019/150‑3 (MQ=255)
           gCGGACCCAGTATCTGGATGCCTTGCGCCTGGGCTTGACCCGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGc                                                                                                                                          <  1:4623172/150‑1 (MQ=255)
             ggACCCATTATCTGGATGCCTTGCGCCTGGGCTTGACCCGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGt                                                                                                                                        <  2:1982434/150‑1 (MQ=255)
                  cAGTATCTGGATGCCTTGCGCCTGGGCTTGACCCGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACaa                                                                                                                                   <  2:4829847/150‑1 (MQ=255)
                   aGTATCTGGATGCCTTGCGCCTGGGCTTGACCCGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAAc                                                                                                                                  >  1:2096605/1‑150 (MQ=255)
                   aGTATCTGGATGCCTTGCGCCTGGGCTTGACCCGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAAc                                                                                                                                  <  1:3802066/150‑1 (MQ=255)
                     tATCTGGATGCCTTGCGCCTGGGCTTGACCCGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCg                                                                                                                                >  2:311825/1‑150 (MQ=255)
                         tGGATGCCTTGCGCCTGGGCTTGACCCGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTaa                                                                                                                            <  2:3840393/150‑1 (MQ=255)
                          ggATGCCTTGCGCCTGGGCTTGACCCGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTaaa                                                                                                                           <  1:3508366/150‑1 (MQ=255)
                              gCCTTGCGCCTGGGCTTGACCCGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGt                                                                                                                              >  2:4951436/1‑143 (MQ=255)
                              gCCTTGCGCCTGGGCTTGACCCGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGt                                                                                                                              <  1:4951436/143‑1 (MQ=255)
                               ccTTGCGCCTGGGCTTGACCCGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATcc                                                                                                                      <  1:2473318/150‑1 (MQ=255)
                                cTTCCGCCTGGGCTTGACCCGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGACCTCACCAGCGTCACAACCGGTAAATATCCg                                                                                                                     >  1:4781889/1‑150 (MQ=255)
                                 ttGCGCCTGGGCTTGACCCGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCgg                                                                                                                    <  2:355910/150‑1 (MQ=255)
                                 ttGCGCCTGGGCTTGACCCGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCgg                                                                                                                    >  2:129927/1‑150 (MQ=255)
                                 ttGCGCCTGGGCTTGACCCGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCgg                                                                                                                    <  1:4295952/150‑1 (MQ=255)
                                      ccTGGGCTTGACCCGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTg                                                                                                               >  1:155234/1‑150 (MQ=255)
                                        tGGGCTTGACCCGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATGTGAGGGCTTCGCCCTCatatcgc                                                                                                                                         >  2:4985928/1‑116 (MQ=255)
                                        tGGGCTTGACCCGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATGTGAGGGCTTCGCCCTCatatcgc                                                                                                                                         <  1:4985928/122‑7 (MQ=255)
                                        tGGGCTTGACCCGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGt                                                                                                             >  2:3137895/1‑150 (MQ=255)
                                        tGGGCTTGACCCGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGt                                                                                                             >  2:119767/1‑150 (MQ=255)
                                            cTTGACCCGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCAt                                                                                                         >  2:2379768/1‑150 (MQ=255)
                                             ttGACCCGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATc                                                                                                        >  2:909845/1‑150 (MQ=255)
                                                aCCCGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCa                                                                                                     <  1:1026848/150‑1 (MQ=255)
                                                  ccGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCAcc                                                                                                   >  1:670356/1‑150 (MQ=255)
                                                  ccGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCAcc                                                                                                   >  1:4481164/1‑150 (MQ=255)
                                                  ccGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCAcc                                                                                                   >  1:1065354/1‑150 (MQ=255)
                                                        ttCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGcc                                                                                             <  1:5170877/150‑1 (MQ=255)
                                                         tCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCg                                                                                            >  1:2090748/1‑150 (MQ=255)
                                                         tCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCg                                                                                            <  2:1217241/150‑1 (MQ=255)
                                                            tACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTa                                                                                         <  2:3464971/150‑1 (MQ=255)
                                                              ccTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTAcc                                                                                       >  2:5103783/1‑150 (MQ=255)
                                                              ccTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTAcc                                                                                       >  2:2374621/1‑150 (MQ=255)
                                                              ccTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTAcc                                                                                       <  1:2459803/150‑1 (MQ=255)
                                                               cTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCt                                                                                      <  1:3128413/150‑1 (MQ=255)
                                                               cTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTACCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCt                                                                                      <  1:3356048/150‑1 (MQ=255)
                                                                tGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCtt                                                                                     >  2:4556986/1‑150 (MQ=255)
                                                                tGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCtt                                                                                     <  2:4445822/150‑1 (MQ=255)
                                                                tGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCtt                                                                                     >  1:251833/1‑150 (MQ=255)
                                                                 gtcgtcGTAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCttt                                                                                    <  2:1664374/150‑1 (MQ=255)
                                                                 gtcgtcGTAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCttt                                                                                    <  1:700001/150‑1 (MQ=255)
                                                                 gtcgtcGTAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCttt                                                                                    <  1:1524217/150‑1 (MQ=255)
                                                                     tcgtAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAAt                                                                                >  2:4398420/1‑150 (MQ=255)
                                                                     tcgtAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAAt                                                                                >  2:4605131/1‑150 (MQ=255)
                                                                      agtAATTGACGGGGCCTCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCTGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATg                                                                               <  2:2908494/149‑1 (MQ=255)
                                                                       gtAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGc                                                                              <  1:2374621/150‑1 (MQ=255)
                                                                            tGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAATGt                                                                         <  1:2801843/150‑1 (MQ=255)
                                                                                 gggCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAATGTTCaaa                                                                    <  1:4652684/150‑1 (MQ=255)
                                                                                 gggCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAATGTTCaaa                                                                    <  1:918595/150‑1 (MQ=255)
                                                                                 gggCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAATGTTCaaa                                                                    <  1:1148804/150‑1 (MQ=255)
                                                                                 gggCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAATGTTCaaa                                                                    <  2:2103666/150‑1 (MQ=255)
                                                                                 gggCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAATGTTCaaa                                                                    <  1:2998982/150‑1 (MQ=255)
                                                                                 gggCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAATGTTCaaa                                                                    <  1:3344533/150‑1 (MQ=255)
                                                                                     cGCTCCTGCAAAAGCCTGTGCTGCCTGGTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAATGTTCAAaggag                                                                <  2:3944096/150‑1 (MQ=255)
                                                                                      gCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAATGTTCAAaggagg                                                               >  1:5381195/1‑150 (MQ=255)
                                                                                      gCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAATGTTCAAaggagg                                                               >  2:5018288/1‑150 (MQ=255)
                                                                                        tCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAATGTTCAAAGGAGGAc                                                             <  1:4646908/150‑1 (MQ=255)
                                                                                           tGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAATGTTCAAAGGAGGACGat                                                          >  1:1098007/1‑150 (MQ=255)
                                                                                             cAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAATGTTCAAAGGAGGACGatat                                                        <  1:5451692/150‑1 (MQ=255)
                                                                                              aaaaGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAATGTTCAAAGGAGGACGATATg                                                       >  1:1767444/1‑150 (MQ=255)
                                                                                               aaaGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAATGTTCAAAGGAGGACGATATGc                                                      >  2:2908241/1‑150 (MQ=255)
                                                                                                 aGCCTGTGCTGCCTCGTGTGGGAGCGGCCTTGCGTCGCGATAGGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGACTGCAATGTTCAAAGGAGGACGATATGCaa                                                    <  2:280276/150‑1 (MQ=255)
                                                                                                  gCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCTCCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAATGTTCAAAGGAGGACGATATGCAAc                                                   <  2:1985757/150‑1 (MQ=255)
                                                                                                   ccTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAATGTTCAAAGGAGGACGATATGCAACg                                                  >  1:4925872/1‑150 (MQ=255)
                                                                                                         ctgcCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAATGTTCAAAGGAGGACGATATGCAACGTCCAgg                                            >  2:1569089/1‑150 (MQ=255)
                                                                                                         ctgcCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAATGTTCAAAGGAGGACGATATGCAACGTCCAgg                                            >  1:4192858/1‑150 (MQ=255)
                                                                                                         ctgcCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAATGTTCAAAGGAGGACGATATGCAACGTCCAgg                                            >  1:711016/1‑150 (MQ=255)
                                                                                                         ctgcCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAATGTTCAAAGGAGGACGATATGCAACGTCCAgg                                            >  1:5188650/1‑150 (MQ=255)
                                                                                                                   tgGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGt                                                                                                             <  1:1952794/75‑1 (MQ=255)
                                                                                                                      gAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAATGtt                                                                        <  1:1925280/109‑1 (MQ=255)
                                                                                                                      gAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAATGtt                                                                        >  2:1925280/1‑109 (MQ=255)
                                                                                                                      gAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAATGTTCAAAGGAGGACGATATGCAACGTCCAGGCTCCCATCGTTTa                               >  2:276921/1‑150 (MQ=255)
                                                                                                                       aGCGGCCTTGCGTCGCTATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAATGTTCAAAGGAGGACGATATGCAACGTCCAGGCTCCCATCGTTTAc                              <  1:5112655/150‑1 (MQ=255)
                                                                                                                       aGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAATGTTCAAAGGAGGACGATATGCAACGTCCAGGCTCCCATCGTTTAc                              <  1:1707246/150‑1 (MQ=255)
                                                                                                                        gCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAATGTTCAAAGGAGGACGATATGCAACGTCCAGGCTCCCATCGTTTACg                             >  2:2036283/1‑150 (MQ=255)
                                                                                                                           gCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAATGTTCAATGGAGGACGATATGCAACGTCCAGGCTCCCATCGTTTACGTgg                          >  1:952873/1‑150 (MQ=255)
                                                                                                                           gCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAATGTTCAAAGGAGGACGATATGCAACGTCCAGGCTCCCATCGTTTACGTgg                          >  1:4010746/1‑150 (MQ=255)
                                                                                                                             cTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAATGTTCAAAGGAGGACGATATGCAACGTCCAGGCTCCCATCGTTTACGTGGCg                        <  2:2708650/150‑1 (MQ=255)
                                                                                                                               tGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAATGTTCAAAGGAGGACGATATGCAACGTCCAGGCTCCCATCGTTTACGTGGCGGt                      <  2:2239953/150‑1 (MQ=255)
                                                                                                                                gCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAATGTTCAAAGGAGGACGATATGCAACGTCCAGGCTCCCATCGTTTACGTGGCGGTc                     <  1:4556986/150‑1 (MQ=255)
                                                                                                                                  gTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAATGTTCAAAGGAGGACGATATGCAACGTCCAGGCTCCCATCGTTTACGTGGCGGTCGa                   >  2:1402575/1‑150 (MQ=255)
                                                                                                                                  gTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAATGTTCAAAGGAGGACGATATGCAACGTCCAGGCTCCCATCGTTTACGTGGCGGTCGa                   <  1:2958145/150‑1 (MQ=255)
                                                                                                                                  gTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAATGTTCAAAGGAGGACGATATGCAACGTCCAGGCTCCCATCGTTTACGTGGCGGTCGa                   <  1:5090663/150‑1 (MQ=255)
                                                                                                                                  gTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAATGTTCAAAGGAGGACGATATGCAACGTCCAGGCTCCCATCGTTTACGTGGCGGTCGa                   <  1:2036283/150‑1 (MQ=255)
                                                                                                                                    cgcgATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAATGTTCAAAGGAGGACGATATGCAACGTCCAGGCTCCCATCGTTTACGTGGCGGTCGAgt                 <  2:1052253/150‑1 (MQ=255)
                                                                                                                                       gATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAATGTTCAAAGGAGGACGATATGCAACGTCCAGGCTCCCATCGTTTACGTGGCGGTCGAGTGTc              <  2:3721226/150‑1 (MQ=255)
                                                                                                                                       gATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAATGTTCAAAGGAGGACGATATGCAACGTCCAGGCTCCCATCGTTTACGTGGCGGTCGAGTGTc              <  2:1627713/150‑1 (MQ=255)
                                                                                                                                       gATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAATGTTCAAAGGAGGACGATATGCAACGTCCAGGCTCCCATCGTTTACGTGGCGGTCGAGTGTc              <  2:3132097/150‑1 (MQ=255)
                                                                                                                                        atatGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAATGTTCAAAGGAGGACGATATGCAACGTCCAGGCTCCCATCGTTTACGTGGCGGTCGAGTGTCt             <  1:2728230/150‑1 (MQ=255)
                                                                                                                                        atatGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAATGTTCAAAGGAGGACGATATGCAACGTCCAGGCTCCCATCGTTTACGTGGCGGTCGAGTGTCt             <  1:2311972/150‑1 (MQ=255)
                                                                                                                                           tGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAATGTTCAAAGGAGGACGATATGCAACGTCCAGGCTCCCATCGTTTACGTGGCGGTCGAGTGTCTGaa          >  1:1410164/1‑150 (MQ=255)
                                                                                                                                           tGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAATGTTCAAAGGAGGACGATATGCAACGTCCAGGCTCCCATCGTTTACGTGGCGGTCGAGTGTCTGaa          <  2:748073/150‑1 (MQ=255)
                                                                                                                                              gggCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAATGTTCAAAGGAGGACGATATGCAACGTCCAGGCTCCCATCGTTTACGTGGCGGTCGAGTGTCTGAATcc       <  1:1292977/150‑1 (MQ=255)
                                                                                                                                               ggCGAAGCCCTCACCAGCGTCACAACCGGTCAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCACTGTTCAAAGGAGGACGATATGCAACGTCCAGGCTCCCCTCGGTTACGTGGCGGTCGAGTGTCTGAATCCg      >  2:4010704/1‑150 (MQ=255)
                                                                                                                                                  gAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAATGTTCAAAGGAGGACGATATGTAACGTCCAGGCTCCCATCGTTTACGTGGCGGTCGAGTGTCTGAATCCGGgc   <  1:4205488/150‑1 (MQ=255)
                                                                                                                                                  gAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAATGTTCAAAGGAGGACGATATGCAACGTCCAGGCTCCCATCGTTTACGTGGCGGTCGAGTGTCTGAATCCGGgc   <  2:1267110/150‑1 (MQ=255)
                                                                                                                                                   aaGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAATGTTCAAAGGAGGACGATATGCAACGTCCAGGCTCCCATCGTTTACGTGGCGGTCGAGTGTCTGAATCCGGgcg  >  2:4894659/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
ACTGCGCGTGAGCGGACCCAGTATCTGGATGCCTTGCGCCTGGGCTTGACCCGGTCTTCGTACCTGTCGTCGTAATTGACGGGGCCGCTCCTGCAAAAGCCTGTGCTGCCTTTTGTGGGAGCGGCCTTGCGTCGCGATATGAGGGCGAAGCCCTCACCAGCGTCACAACCGGTAAATATCCGGCCCTGGTCCATCGCACCAATGCCGCTACCTTTGAATGCAATGTTCAAAGGAGGACGATATGCAACGTCCAGGCTCCCATCGTTTACGTGGCGGTCGAGTGTCTGAATCCGGGCG  >  NC_002947/660797‑661093

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: