Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,896,135 C→T 5.2% Q218Q (CAG→CAA PP_3439 ← AraC family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,896,1350CT 5.2% 342.3 / 5.8 115Q218Q (CAG→CAAPP_3439AraC family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (54/55);  new base T (3/3);  total (57/58)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.83e-01

Only 100 of 116 total aligned reads displayed.

TCCTGCAGGAGCGCCCAGGCTAGCGACGGCATTTCTGTAATGCCCCGGTGTCGCCGCCAGGTGCTGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATCAGGTAGGACCGAGACAGGCCGCAAGCCCCGCAGATGTCCTCGAGGCTGAGCGGGTCACTGCGATGGGC  >  NC_002947/3895988‑3896284
                                                                                                                                                    |                                                                                                                                                     
tCCTGCAGGAGCGCCCAGGCTAGCGACGGCATTTCTGTAATGCCCCGGTGTCGCCGCCAGGTGCTGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTg                                                                                                                                                     >  1:801935/1‑150 (MQ=255)
      aGGAGCGCCCAGGCTAGCGACGGCATTTCTGTAATGCCCCGGTGTCGCCGCCAGGTGCTGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCg                                                                                                                                               >  1:1839791/1‑150 (MQ=255)
       ggAGCGCCCAGGCTAGCGACGGCATTTCTGTGATGCCCCGGTGTCGCCGCCAGGTGCTGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCgg                                                                                                                                              <  1:319583/150‑1 (MQ=255)
        gAGCGCCCAGGCTAGCGACGGCATTTCTGTAATGCCCCGGTGTCGCCGCCAGGTGCTGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCggg                                                                                                                                             >  2:354139/1‑150 (MQ=255)
              ccAGGCTAGCGACGGCATTTCTGTAATGCCCCGGTGTCGCCGCCAGGTGCTGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCGGGCGCAGTTGCGccagg                                                                                                                                              >  1:2191928/1‑140 (MQ=255)
              ccAGGCTAGCGACGGCATTTCTGTAATGCCCCGGTGTCGCCGCCAGGTGCTGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCGGGCGCAGTTGCGCCAgg                                                                                                                                              <  2:2191928/143‑1 (MQ=255)
                  gCTAGCGACGGCATTTCTGTAATGCCCCGGTGTCGCCGCCAGGTGCTGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCAc                                                                                                                                   <  1:836059/150‑1 (MQ=255)
                        gACGGCATTTCTGTAATGCCCCGGTGTCGCCGCCAGGTGCTGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTg                                                                                                                             >  2:1878149/1‑150 (MQ=255)
                           ggCATTTCTGTAATGCCCCGGTGTCGCCGCCAGGTGCTGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCTGTCCAGGTGCACGCGCTGGTc                                                                                                                          >  2:2396811/1‑150 (MQ=255)
                           ggCATTTCTGTAATGCCCCGGTGTCGCCGCCAGGTGCTGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTc                                                                                                                          <  2:1648550/150‑1 (MQ=255)
                              aTTTCTGTAATGCCCCGGTGTCGCCGCCAGGTGCTGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCgg                                                                                                                                              >  2:360438/1‑127 (MQ=255)
                              aTTTCTGTAATGCCCCGGTGTCGCCGCCAGGTGCTGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCgg                                                                                                                                              <  1:360438/127‑1 (MQ=255)
                               tttCTGTAATGCCCCGGTGTCGCCGCCAGGTGCTGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGc                                                                                                                      <  1:2261364/150‑1 (MQ=255)
                                ttCTGTAATGCCCCGGTGTCGCCGCCAGGTGCTGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCa                                                                                                                     >  2:60333/1‑150 (MQ=255)
                                 tCTGTAATGCCCCGGTGTCGCCGCCAGGTGCTGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAg                                                                                                                    <  2:1050657/150‑1 (MQ=255)
                                   tGTAATGCCCCGGTGTCGCCGCCAGGTGCTGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGt                                                                                                                  >  1:1471668/1‑150 (MQ=255)
                                    gTAATGCCCCGGTGTCGCCGCCAGGTGCTGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTg                                                                                                                             <  2:2214617/138‑1 (MQ=255)
                                    gTAATGCCCCGGTGTCGCCGCCAGGTGCTGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTg                                                                                                                             >  1:2214617/1‑138 (MQ=255)
                                    gTAATGCCCCGGTGTCGCCGCCAGGTGCTGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTa                                                                                                                 <  2:1611641/150‑1 (MQ=255)
                                    gTAATGCCCCGGTGTCGCCGCCAGGTGCTGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTa                                                                                                                 <  1:1878149/150‑1 (MQ=255)
                                     tAATGCCCCGGTGTCGCCGCCAGGTGCTGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAg                                                                                                                <  1:788144/150‑1 (MQ=255)
                                                  tCGCCGCCAGGTGCTGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGtt                                                                                                                                       >  1:1949836/1‑114 (MQ=255)
                                                  tCGCCGCCAGGTGCTGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGtt                                                                                                                                       <  2:1949836/114‑1 (MQ=255)
                                                  tCGCCGCCAGGTGCTGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGc                                                                                                                      <  2:1646757/131‑1 (MQ=255)
                                                  tCGCCGCCAGGTGCTGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGc                                                                                                                      >  1:1646757/1‑131 (MQ=255)
                                                  tCGCCGCCAGGTGCTGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCa                                                                                                   >  2:1703842/1‑150 (MQ=255)
                                                   cgccgcCAGGTGCTGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAg                                                                                                  <  2:2207530/150‑1 (MQ=255)
                                                    gccgccAGGTGCTGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACgcgc                                                                                                                               <  2:432237/120‑1 (MQ=255)
                                                    gccgccAGGTGCTGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACgcgc                                                                                                                               >  1:432237/1‑120 (MQ=255)
                                                    gccgccAGGTGCTGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAgg                                                                                                 >  1:357864/1‑150 (MQ=255)
                                                     ccgccAGGTGCTGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCTCCCCACGCGCACCCCTGCATAATTCCGCCCTTTTGCGTCGTCTTCCGTGCCAGTTGCACTCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGc                                                                                                <  2:1463506/150‑1 (MQ=255)
                                                      cgccAGGTGCTGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGcc                                                                                               >  1:1976661/1‑150 (MQ=255)
                                                          aGGTGCTGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAg                                                                                           >  2:1587157/1‑150 (MQ=255)
                                                           ggTGCTGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGc                                                                                          >  2:1556414/1‑150 (MQ=255)
                                                           ggTGCTGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGc                                                                                          <  2:1839791/150‑1 (MQ=255)
                                                            gTGCTGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCg                                                                                         >  1:213728/1‑150 (MQ=255)
                                                                tGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTgg                                                                                     <  2:2400299/150‑1 (MQ=255)
                                                                 gTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGc                                                                                    <  2:591194/150‑1 (MQ=255)
                                                                 gTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGc                                                                                    >  1:1277173/1‑150 (MQ=255)
                                                                  tttGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGaa                                                                                            <  2:118987/141‑1 (MQ=255)
                                                                  tttGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGaa                                                                                            >  1:118987/1‑141 (MQ=255)
                                                                     gAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGaaa                                                                                <  1:1675214/150‑1 (MQ=255)
                                                                     gAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGaaa                                                                                <  1:1317757/150‑1 (MQ=255)
                                                                       aTGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAgg                                                                              >  2:571597/1‑150 (MQ=255)
                                                                        tGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGTAAGGc                                                                             >  2:1628358/1‑150 (MQ=255)
                                                                        tGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGc                                                                             <  1:621887/150‑1 (MQ=255)
                                                                        tGCCCGCTGCAGGTGTGCCTGATCAGCGAACCAGCCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGc                                                                             <  1:1830328/150‑1 (MQ=255)
                                                                               tGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTg                                                                                                                             >  2:2247064/1‑95 (MQ=255)
                                                                               tGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTg                                                                                                                             <  1:2247064/95‑1 (MQ=255)
                                                                               tGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATc                                                                      <  2:711959/150‑1 (MQ=255)
                                                                                 cAGGGGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATCAg                                                                    <  1:1431120/150‑1 (MQ=255)
                                                                                   ggTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGTAAGGCCCGAATCAGGt                                                                  <  1:1628358/150‑1 (MQ=255)
                                                                                        gCCTGATTAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATCAGGTAGGAc                                                             >  2:1359394/1‑150 (MQ=255)
                                                                                         ccTGATCAGTGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATgcg                                                                                                         >  1:687994/1‑105 (MQ=255)
                                                                                         ccTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAATCGTTGGCGAAAGGCACGAATCAGGTAGGAcc                                                            >  2:1144649/1‑150 (MQ=255)
                                                                                         ccTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATCAGGTAGGAcc                                                            >  1:1534594/1‑150 (MQ=255)
                                                                                           tGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGTGTTGGCGAAAGGCACGAATCAGGTAGGACCga                                                          <  2:801935/150‑1 (MQ=255)
                                                                                           tGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATCAGGTAGGACCga                                                          <  2:213728/150‑1 (MQ=255)
                                                                                           tGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATCAGGTAGGACCga                                                          >  1:774891/1‑150 (MQ=255)
                                                                                           tGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATCAGGTAGGACCga                                                          >  2:794589/1‑150 (MQ=255)
                                                                                           tGATCAGCGAACCCGGCCTCCTGTGCAACCTCAGAAAATAGGCAAAAATGGACCCACTAAGTCCAAGTCAGGAGGACGCGCTTGTCGATCAGGTGGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATCAGGTAGGACCga                                                          <  2:1463361/150‑1 (MQ=255)
                                                                                            gATCAGCGAACCCGGCCTCCTGTGAAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCGCGCGCTGGTCTAGCAGCGAGGTGTAGGGAAAGAGTCCGAACCGTTTGCGAAAGGCACGAATCAGGTAGGACCgag                                                         >  2:663216/1‑150 (MQ=255)
                                                                                                aGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGTATTTGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATCAGGTAGGACCGAGACAg                                                     >  1:166716/1‑150 (MQ=255)
                                                                                                 gCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATCAGGTAGGACCGAGACAgg                                                    <  1:1556414/150‑1 (MQ=255)
                                                                                                     aCCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATCAGGTAGGa                                                              >  2:633997/1‑136 (MQ=255)
                                                                                                     aCCCCGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATCAGGTAGGa                                                              <  1:633997/136‑1 (MQ=255)
                                                                                                       ccGGCCTCCTGTGCAACTTCAGCAATTGGGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATCAGGTAGGACCGAGACAGGCCGCaa                                              <  1:1173264/150‑1 (MQ=255)
                                                                                                       ccGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATCAGGTAGGACCGAGACAGGCCGCaa                                              >  1:1585163/1‑150 (MQ=255)
                                                                                                         ggCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATCAGGTAGGACCGAGACAGGCCGCAAGc                                            >  2:2253377/1‑150 (MQ=255)
                                                                                                         ggCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATCAGGTAGGACCGAGACAGGCCGCAAGc                                            <  1:354139/150‑1 (MQ=255)
                                                                                                             tcctGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGCAGCGTTGGCGAAAGGCCCGAATCAGGTAGGACCGAGACAGGTCGCCAGccccc                                        >  2:925734/1‑149 (MQ=255)
                                                                                                             tcctGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATCAGGTAGGACCGAGACAGGCCGCAAGCCCCg                                        >  1:472090/1‑150 (MQ=255)
                                                                                                                  tgCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCGGTTGCACGCGCTGGTCTAGCAGGTAGCCAGGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATCAGGTAGGACCGAGACAGGCCGCAAGCCCCGCAGAt                                   <  2:2416807/150‑1 (MQ=255)
                                                                                                                  tgCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCCGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATCAGGTAGGACCGAGACAGGCCGCAAGCCCCGCAGAt                                   >  1:1634718/1‑150 (MQ=255)
                                                                                                                         ttCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATCAGGTAGGACCGAGACAGGCCGCAAGCCCCGCAGATGTCCTc                             >  1:950515/1‑150 (MQ=255)
                                                                                                                           cAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATCAGGTAAGTCCGAGACAGGCCGCAAGCCCCGCAGATGTCCTCGa                           <  1:1144649/150‑1 (MQ=255)
                                                                                                                            aGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGACAGGCACGAATCAGGTAGGACCGAGACAAGCCGCAAGCCCCGCAGATGTCCTCCAg                          >  2:853601/1‑150 (MQ=255)
                                                                                                                            aGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATCAGGTAGGACCGAGACAGGCCGCAAGCCCCGCAGATGTCCTCGAg                          >  2:1976981/1‑150 (MQ=255)
                                                                                                                                aTTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCAc                                                                           >  1:1122563/1‑97 (MQ=255)
                                                                                                                                aTTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCAc                                                                           <  2:1122563/97‑1 (MQ=255)
                                                                                                                                aTTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATCAGGTAGGACCGAGACAGGCCGCAAGCCCCGCAGATGTCCTCGAGGCTg                      >  1:1081811/1‑150 (MQ=255)
                                                                                                                                   gAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATCAGGTAGGACCGAGACAGGCCGCAAGCCCCGCAGATGTCCTCGAGGCTGAGc                   <  2:1585163/150‑1 (MQ=255)
                                                                                                                                   aactggCCCGGGCGCAGTTGCGCCAGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTgg                                                                                     >  2:2320388/5‑84 (MQ=25)
                                                                                                                                   aactggCCCGGGCGCAGTTGCGCCAGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTgg                                                                                     <  1:2320388/80‑1 (MQ=25)
                                                                                                                                   aactggCCCGGGCGCAGTTGCGCCAGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATCAGGTAGGACCGAGACAGGCCGCAAGCCCCGCAGATGTCCTCGAGGCTGAGc                   >  2:1835804/5‑150 (MQ=255)
                                                                                                                                    aGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATCAGGTAGGACCGAGACAGGCCGCAAGCCCCGCAGATGTCCTCGAGGCTGAGCg                  <  1:330817/150‑1 (MQ=255)
                                                                                                                                     gCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATCAGGTAGGACCGAGACAGGCCGCAAGCCCCGCAGATGTCCTCGAGGCTGAGCgg                 <  2:1976661/150‑1 (MQ=255)
                                                                                                                                     gCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATCAGGTAGGACCGAGACAGGCCGCAAGCCCCGCAGATGTCCTCGAGGCTGAGCgg                 <  1:1587157/150‑1 (MQ=255)
                                                                                                                                       ggCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATCAGGTAGGACCGAGACAGGCCGCAAGCCCCGCAGATGTCCTCGAGGCTGAGCGGGt               <  1:794589/150‑1 (MQ=255)
                                                                                                                                          ccTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATCAGGTAGGACCGAGACAGGCCGCAAGCCCCGCAGATGTCCTCGAGGCTGAGCGGGTCAc            >  2:110165/1‑150 (MQ=255)
                                                                                                                                           cTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATCAGGTAGGACCGAGACAGGCCGCAAGCCCCGCAGATGTCCTCGAGGCTGAGCGGGTCACt           <  1:571597/150‑1 (MQ=255)
                                                                                                                                           cTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATCAGGTAGGACCGAGACAGGCCGCAAGCCCCGCAGATGTCCTCGAGGCTGAGCGGGTCACt           >  2:795374/1‑150 (MQ=255)
                                                                                                                                            tGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATCAGGTAGGACCGAGACAGGCCGCAAGCCCCGCAGATGTCCTCGAGGCTGAGCGGGTCACTg          >  1:1960804/1‑150 (MQ=255)
                                                                                                                                               cgcAACTGCGCCTGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATCAGGTAGGACCGAGACAGGCCGCAAGCCCCGCAGATGTCCTCGAGGCTGAGCGGGTCACTGCGa       <  2:1594732/150‑1 (MQ=255)
                                                                                                                                                gcAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATCAGGTAGGACCGAGACAGGCCGCAAGCCCCGCAGATGTCCTCGAGGCTGAGCGGGTCACTGc         <  1:1757672/147‑1 (MQ=255)
                                                                                                                                                gcAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATCAGGTAGGACCGAGACAGGCCGCAAGCCCCGCAGATGTCCTCGAGGCTGAGCGGGTCACTGc         >  2:1757672/1‑147 (MQ=255)
                                                                                                                                                 cAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATCAGGTAGGACCGAGACAGGCCGCAAGCCCCGCAGATGTCCTCGAGGCTGAGCGGGTCACTGCGATg     <  1:795374/150‑1 (MQ=255)
                                                                                                                                                  aaCTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATCAGGTAGGACCgaga                                                        >  1:2183777/1‑98 (MQ=255)
                                                                                                                                                  aaCTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATCAGGTAGGACCgaga                                                        <  2:2183777/98‑1 (MQ=255)
                                                                                                                                                    cTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATCAGGTAGGACCGAGACAGGCCGCAAGCCCCGCAGATGTCCTCGAGGCTGAGCGGGTCACTGCGATGGgc  <  1:568143/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TCCTGCAGGAGCGCCCAGGCTAGCGACGGCATTTCTGTAATGCCCCGGTGTCGCCGCCAGGTGCTGTTTGAATGCCCGCTGCAGGTGTGCCTGATCAGCGAACCCGGCCTCCTGTGCAACTTCAGCAATTGAGCGGCCCTGGCGCAACTGCGCCCGGGCCAGTTGCACGCGCTGGTCTAGCAGGTAGCCATGCGGCGTCAGGCCGAAGCGTTGGCGAAAGGCACGAATCAGGTAGGACCGAGACAGGCCGCAAGCCCCGCAGATGTCCTCGAGGCTGAGCGGGTCACTGCGATGGGC  >  NC_002947/3895988‑3896284

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