Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,261,326 C→T 6.2% intergenic (+394/‑310) asd → / → PP_1990 aspartate‑semialdehyde dehydrogenase/transposase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,261,3260CT 6.2% 889.5 / 31.8 290intergenic (+394/‑310)asd/PP_1990aspartate‑semialdehyde dehydrogenase/transposase
Reads supporting (aligned to +/- strand):  ref base C (140/132);  new base T (12/6);  total (152/138)
Fisher's exact test for biased strand distribution p-value = 2.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.76e-01

Only 100 of 296 total aligned reads displayed.

GTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGACCCGCCGTAAGGGCGGAAAGGTGACTATGCGCCACCTGCGTGAATGAATACCTACTCGGTGTGAATGCCCACATCCAAGAAGCAAGAAGCAAGAAGCAAGAAGCAAGAAGCAAATTCGTCTGTTTTGACTCAGC  >  NC_002947/2261179‑2261475
                                                                                                                                                   |                                                                                                                                                     
     gggCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACa                                                                                                                                                >  1:1864070/1‑150 (MQ=255)
         agggAGGCCCGCCACGCGCGGGGCCCCAGGTTGGGGCCTGCAGCGCGGTCCACCCGGTAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACaaaaa                                                                                                                                            <  2:4629799/148‑1 (MQ=17)
         aTGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTGGTTCCTTTTTGCCACGACACAAAg                                                                                                                                           >  2:4605014/1‑150 (MQ=255)
           ggATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAgg                                                                                                                                          >  2:1519723/1‑150 (MQ=255)
             aTGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGGAGCGCAGGGGCCGGATGATGGGAGCCGACGGTGTTTGGTTCCTTTTTGCCACGACAAAAAGGGa                                                                                                                                        >  2:5506888/1‑150 (MQ=255)
              ggCCCGTCACGCGGGGGCCCCCGGGTTGGGCCCGGCCGGGCGGCCCCCCCGGGACCACGCCCGGCGCCCGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGAc                                                                                                                                       <  1:1549490/149‑1 (MQ=14)
               gCCCGTCACGCGGTGGGCCCCAGGATGGGGCCTGCTGCGCGGTCCTCCGGGGAGCAAGCCCGGAGCGTGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGAcc                                                                                                                                      <  2:3244994/150‑1 (MQ=17)
                 ccGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGACccg                                                                                                                                    >  1:4473199/1‑150 (MQ=255)
                  cGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGACccgc                                                                                                                                   >  2:4071681/1‑150 (MQ=255)
                   gTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGGGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGACccgcc                                                                                                                                  >  2:5235584/1‑150 (MQ=255)
                   gTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGACccgcc                                                                                                                                  >  1:5959980/1‑150 (MQ=255)
                       aGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGTAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGACCCGCCGTaa                                                                                                                              >  2:4559901/1‑150 (MQ=255)
                          ggtGGGCCCCTGGTTGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGACCCGCCGTAAggg                                                                                                                           <  2:5768300/148‑1 (MQ=255)
                            ttgctcccAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGACCCGCCGTAAGGGCg                                                                                                                         >  1:1627515/6‑150 (MQ=255)
                             gggCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGACCCGCCGTAAGGGCgg                                                                                                                        <  2:3319927/150‑1 (MQ=255)
                              gcccccAGTATGGGGCCTGCTGGGCGGACCGCCCGGGATAATGCCCGGTGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGACCCGCCGTAAGGGCGGa                                                                                                                       <  1:3351695/148‑1 (MQ=17)
                               gCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGCGAGGGGAGCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGACCCGCCGTAAGGGCGGaa                                                                                                                      >  1:3727585/1‑150 (MQ=255)
                                ccccTGGATGGGGCCTGCTGGGGGGTCCTCCGGGGAGCACGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGACCCGCCGTAAGGGCGGaaa                                                                                                                     <  2:6451490/150‑1 (MQ=17)
                                ccccAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGACCCGCCGTAAGGGCGGaaa                                                                                                                     >  2:6294193/1‑150 (MQ=255)
                                 cccAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGACCCGCCCTTAcggcgggtcc                                                                                                                    <  2:5959882/150‑5 (MQ=255)
                                 cccAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGCGAGGGGACCCAGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGACCCGCCGTAAGGGCGGAAAg                                                                                                                    >  1:2064116/1‑150 (MQ=255)
                                   cAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGACCCGCCGTAAGGGCGGAAAGGt                                                                                                                  >  1:5866972/1‑150 (MQ=255)
                                   cAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGACCCGCCGTAAGGGCGGAAAGGt                                                                                                                  >  1:248241/1‑150 (MQ=255)
                                   cAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGACCCGCCGTAAGGGCGGAAAGGt                                                                                                                  >  1:2254054/1‑150 (MQ=255)
                                    aGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGACCCGCCGTAAGGGCGGAAAGGTg                                                                                                                 >  1:1317577/1‑150 (MQ=255)
                                     ggttGGGGCCTGCCGCGGGGGCCTCCGGGGAGCTAGCCCGGAGCGCGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGACCCGCCGTAAGGGCGGAAAGGTGa                                                                                                                <  2:1766113/147‑1 (MQ=17)
                                          gggCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGTCACGACAAAAAggg                                                                                                                                         <  1:3151969/120‑1 (MQ=255)
                                           ggCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGACCCGCCGTAAGGGCGGAAAGGTGACTATgc                                                                                                          >  2:1902170/1‑150 (MQ=255)
                                           ggCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGACCCGCCGTAAGGGCGGAAAGGTGACTATgc                                                                                                          >  2:2237764/1‑150 (MQ=255)
                                             ccTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGACCCGCCGTAAGGGCGGAAAGGTGACTATgcgc                                                                                                        >  1:1018803/1‑150 (MQ=255)
                                                  aGCGCGGTCCTCCCGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGACCCGCCGTAAGGGCGGAAAGGTGACTATGCGCCACCt                                                                                                   >  1:6216907/1‑150 (MQ=255)
                                                     gcgGTCCTCCCGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGACCCGCCGTAAGGGCGGAAAGGTGACTATGCGCCACCTGCg                                                                                                >  1:3554772/1‑150 (MQ=255)
                                                     gcgGTCCTCCCGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGACCCGCCGGAAGGGCGGTAAGGTGACTATTCGCCACCGGCg                                                                                                >  1:327657/1‑150 (MQ=255)
                                                      cgGTCCTCCCGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGACCCGCCGTAAGGGCGGAAAGGTGACTATGCGCCACCTGCGt                                                                                               <  2:4893597/150‑1 (MQ=255)
                                                       ggTCCTCCCGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGACCCGCCGTAAGGGCGGAAAGGTTACTATTCGCCACCAGCGtg                                                                                              >  1:1989156/1‑150 (MQ=255)
                                                       ggTCCTCCCGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGACCCGCCGTAAGGGCGGAAAGGTGACAATGCGCCACCTGCGtg                                                                                              >  1:3335771/1‑150 (MQ=255)
                                                        gccctccGGGGAGCATGCCGGTAGCGTGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGACCCGCCGTAAGGGCGGAAAGGTGACTATGCGCCACCTGCGtga                                                                                             <  1:5559678/148‑1 (MQ=18)
                                                            tccCGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGACCCGCCGTAAGGGCGGAAAGGTGACTATGCGCCACCTGCGtgaatga                                                                                         >  1:6765008/1‑150 (MQ=255)
                                                             cccGGGAGCAAGCCGGGAGCTAGGGGACCCCGGAGCGTAGCGCAGGGGCCTGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGACCCGCCGTAAGGGCGGAAAGGTGACTATGCGCCACCTGCGtgaatgaa                                                                                        <  2:5373129/150‑1 (MQ=255)
                                                              ccGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCTTAGCGCAGGGGCCGGATGATGGGAGCCGACGGATTTTGGTTCCTTTTTGCCACGACAAAAAGGGACCAGCCGTAAGGGCGGAAAGGTGACTAAGCGCCACCTGCGTgaatgaat                                                                                       >  1:4074943/1‑150 (MQ=255)
                                                              ccGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGACCCGCCGTAAGGGCGGAAAGGTGACTATGCGCCACCTGCGTgaatgaat                                                                                       >  2:880743/1‑150 (MQ=255)
                                                               cGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGACCCGCCGTAAGGGCGGAAAGGTGACTATGCGCCACCTGCGTGAATGAATa                                                                                      >  2:5328363/1‑150 (MQ=255)
                                                               cGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGACCCGCCGTAAGGGCGGAAAGGTGACGTGGCGCCACCGGCGtgagtgtctt                                                                                      >  2:4590719/1‑146 (MQ=17)
                                                               cGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGGTTTTGGTTCCTGTTTGCCCCGACACAAAGGGACCCGCCGTAAGGGCGGACAGGTGACTATGCGCCACCGGCGtgaatgactc                                                                                      >  2:1191426/1‑147 (MQ=18)
                                                                    gCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGACCCGCCGTAAGGGCGGAAAGGTGACTATGCGCCACCTGCGTGAATGAATACCTAc                                                                                 <  2:4473199/150‑1 (MQ=255)
                                                                      aaGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGACCCGCCGTAAGGGCGGAAAGGTGACTATGCGCCACCTGCGTGAATGAATACCTACTc                                                                               >  2:6065850/1‑150 (MQ=255)
                                                                        gCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGACCCGCCGTAAGGGCGGAAAGGTGACTATGCGCCACCGGCGTGAATGAATACCTACTCgg                                                                             >  1:5688590/1‑150 (MQ=255)
                                                                          ccGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGACCCGCCGTAAGGGCGGAAAGGTGACTATGCGCCACCTGCGTGAATGAATACCTACTCGgtg                                                                           >  2:4313279/1‑150 (MQ=255)
                                                                           cGGCGCGCGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGAGGGGAGCCGACGGATTTTGGGTCCTTTTTGCCACGACAAAAAGGGACCCGCCGTAAGGGCGGAAAGGTGACGATGCGCCACCTGCGTGAATGCATACCTACTCGgtgt                                                                          <  1:5996915/150‑1 (MQ=21)
                                                                           cGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGACCCGCCGTAAGGGCGGAAAGGTGACTATGCGCCACCTGCGTGAATGAATACCTACTCGgtgt                                                                          <  2:6414543/150‑1 (MQ=255)
                                                                             gcgggcggggACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGACCCGCCGTAAGGGCGGAAAGGTGACTATGCGCCACCTGCGTGAATGAATACCTACTCGGTGTGa                                                                        <  1:2397980/144‑1 (MQ=255)
                                                                              cgcgcGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGACCCGCCGTAAGGGCGGAAAGGTGACTATGCGCCACCTGCGTGAATGAATACCTACTCGGTGTGaa                                                                       <  1:4359185/149‑1 (MQ=255)
                                                                               gCGTGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGACCCGCCGTAAGGGCGGAAAGGTGACTATGCGCCACCTGCGTGAATGAATACCTACTCGGTGTGAAt                                                                      <  1:4831564/150‑1 (MQ=255)
                                                                                  aGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGACCCGCCGTAAGGGCGGAAAGGTGACTATGCGCCACCTGCGTGAATGAATACCTACTCGGTGTGAATGcc                                                                   >  2:4293542/1‑150 (MQ=255)
                                                                                   ggggACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGACCCGCCGTAAGGGCGGAAAGGTGACTATGCGCCACCTGCGTGAATGAATACCTACTCGGTGTGAATGccc                                                                  <  1:2948805/150‑1 (MQ=255)
                                                                                    gggACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGACCCGCCGTAAGGGCGGAAAGGTGACTATGCGCCACCTGCGTGAATGAATACCTACTCGGTGTGAATGCCca                                                                 <  2:6729342/150‑1 (MQ=255)
                                                                                     ggACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGATGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGACCCGCCGTAAGGGCGGAAAGGTGACTATGCGCCACCTGCGTGAATGAATACCTACTCGGTGTGAATGCCcac                                                                >  2:3095604/1‑150 (MQ=255)
                                                                                         cccGGAGCGCAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGACCCGCCGTAAGGGCGGAAAGGTGACTATGCGCCACCTGCGTGAATGAATACCTACTCGGTGTGAATGCCCACATcc                                                            >  2:5442651/1‑150 (MQ=255)
                                                                                           cGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGACCCGCCGTAAGGGCGGAAAGGTGACTATGCGCCACCTGCGTGAATGAATACCTACTCGGTGTGAATGCCCACATCcaa                                                          <  1:82096/150‑1 (MQ=255)
                                                                                           cGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGACCCGCCGTAAGGGCGGAAAGGTGACTATGCGCCACCTGCGTGAATGAATACCTACTCGGTGTGAATGCCCACATCcaa                                                          <  2:253600/150‑1 (MQ=255)
                                                                                              aGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGACCCGCCGTAAGGGCGGAAAGGTGACTATGCGCCACCTGCGTGAATGAATACCTACTCGGTGTGAATGCCCACATCcaagaa                                                       <  1:1106016/150‑1 (MQ=255)
                                                                                              aGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGACCCGCCGTAAGGGCGGAAAGGTGACTATGCGCCACCTGCGTGAATGAATACCTACTCGGTGTGAATGCCCACATCcaagaa                                                       <  2:1627515/150‑1 (MQ=255)
                                                                                              aGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGACCCGCCGTAAGGGCGGAAAGGTGACTATGCGCCACCTGCGTGAATGAATACCTACTCGGTGTGAATGCCCACATCcaagaa                                                       <  2:1927614/150‑1 (MQ=255)
                                                                                                cGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGACCCGCCGTAAGGGCGGAAAGGTGACTATGCGCCACCTGCGtgaat                                                                                           >  2:6524518/1‑112 (MQ=255)
                                                                                                 gTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGACCCGCCGTAAGGGCGGAAAGGTGACTATGCGCCACCTGCGTGAATGAATACCTACTCGGTGTGAATGCCCACATCcaagaagca                                                    >  1:5049526/1‑150 (MQ=255)
                                                                                                  tAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGACCCGCCGTAAGGGCGGAAAGGTGACTATGCGCCACCTGCGTGAATGAATACCTACTCGGTGTGAATGCCCACATCcaagaagcaa                                                   <  1:6476746/150‑1 (MQ=255)
                                                                                                    gcgcAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGACCCGCCGTAAGGGCGGAAAGGTGACTATGCGCCACCTGCGTGAATGAATACCTACTCGGTGTGAATGCCCACATCcaagaagcaaga                                                 <  1:4984982/150‑1 (MQ=255)
                                                                                                    gcgcAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGACCCGCCGTAAGGGCGGAAAGGTGACTAATCGCCCCCAGCGGGAATGAATAACTACTAGGTGTGAATGCCCACAACcaagaagcaaca                                                 >  1:2557162/1‑148 (MQ=18)
                                                                                                           ggCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGACCCGCCGTAAGGGCGGAAAGGTGACTATGCGCCACCTGCGTGAATGAATACCTACTCGGTGTGAATGCCCACATCcaagaagcaagaagcaaga                                          >  2:4465032/1‑150 (MQ=255)
                                                                                                             ccGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGACCCGCCGTAAGGGCGGAAAGGTGACTATGCGCCACCTGCGTGAATGAATACCTACTCGGTGTGAATGCCCACATCcaagaagcaagaagcaagaag                                        <  1:1808257/150‑1 (MQ=255)
                                                                                                                gatgatGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGACCCGCCGTAAGGGCGGAAAGGTGACTATGCGCCACCGGCGTGAATGAATACCTAATCGGTGTGAATGCCCACATCcaagaagcaagaagcaagaagcaa                                     >  1:6640873/1‑150 (MQ=255)
                                                                                                                 atgatgGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGACCCGCCGTAAGGGCGGAAAGGTGACTATGCGCCACCTGCGTGAATGAATACCTACTCGGTGTGAATGCCCACATCcaagaagcaagaagcaagaagcaag                                    <  1:2875609/150‑1 (MQ=255)
                                                                                                                  tgatgGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGACCCGCCGTAAGGGCGGAAAGGTGACTATGCGCCACCTGCGTGAATGAATACCTACTCGGTGTGAATGCCCACATCcaagaagcaagaagcaagaagcaaga                                   <  2:4415603/150‑1 (MQ=255)
                                                                                                                   gatgGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGACCCGCCGTAAGGGCGGAAAGGTGACTATGCGCCACCTGCGTGAATGAATACCTACTCGGTGTGAATGCCCACATCcaagaagcaagaagcaagaagcaagaa                                  <  2:1317577/150‑1 (MQ=255)
                                                                                                                        gAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGACCCGCCGTAAGGGCGGAAAGGTGACTATGCGCCACCTGCGTGAATGAATACCTACTCGGTGTGAATGCCCACATCcaagaagcaagaagcaagaagcaagaagcaag                             <  1:5759861/150‑1 (MQ=255)
                                                                                                                          gCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGACCCGCCGTAAGGGCGGAAAGGTGACTATGCGCCACCTGCGTGAATGAATACCTACTCGGTGTGAATGCCCACATCcaagaagcaagaagcaagaagcaagaagcaagaa                           <  1:5235535/150‑1 (MQ=255)
                                                                                                                          gCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGACCCGCCGTAAGGGCGGAAAGGTGACTATGCGCCACCTGCGTGAATGAATACCTACTCGGTGTGAATGCCCACATCcaagaagcaagaagcaagaagcaagaagcaagaa                           <  2:4959244/150‑1 (MQ=255)
                                                                                                                                  tttttGGTTCCTTTTTGCCACGACAATAAGGGACCCGCCGTAAGGGCGGAAAGGTGACTATGCGCCACCTGCGTGAATGAATACCTACTCGGTGTGAATGCCCACATCCAAGAAGCAAGAAGCAAGAAGCAAGAAGCAAGAAGCAAATTc                   >  2:6675468/1‑150 (MQ=255)
                                                                                                                                     ttGGTTCCTTTTTGCCACGACAAAAAGGGACCCGCCGTAAGGGCGGAAAGGTGACTATGCGCCACCTGCGTGAATGAATACCTACTCGGTGTGAATGCCCACATCCAAGAAGCAAGAAGCAAGAAGCAAGAAGCAAGAAGCAAATTCGTc                >  2:663249/1‑150 (MQ=255)
                                                                                                                                     ttGGTTCCTTTTTGCCACGACAAAAAGGGACCCGCCGTAAGGGCGGAAAGGTGACTATGCGCCACCTGCGTGAATGAATACCTACTCGGTGTGAATGCCCACATCCAAGAAGCAAGAAGCAAGAAGCAAGAAGCAAGAAGCAAATTCGTc                <  1:4034867/150‑1 (MQ=255)
                                                                                                                                      tGGTTCCTTTTTGCCACGACAAAAAGGGACCCGCCGTAAGGGCGGAAAGGTGACTATGCGCCACCTGCGTGAATGAATACCTACTCGGTGTGAATGCCCACATCCAAGAAGCAAGAAGCAAGAAGCAAGAAGCAAGAAGCAAATTCGTCt               <  2:5217456/150‑1 (MQ=255)
                                                                                                                                        gTTCCTTTTTGCCACGACAAAAAGGGACCCGCCGTAAGGGCGGAAAGGTGACTATGCGCCACCTGCGTGAATGAATACCTACTCGGTGTGAATGCCCACATCCAAGAAGCAAGAAGCAAGAAGCAAGAAGCAAGAATCAAATTCGTCTGt             >  2:3729909/1‑150 (MQ=255)
                                                                                                                                        gTTCCTTTTTGCCACGACAAAAAGGGACCCGCCGTAAGGGCGGAAAGGTGACTATGCGCCACCTGCGTGAATGAATACCTACTCGGTGTGAATGCCCACATCCAAGAAGCAAGAAGCAAGAAGCAAGAAGCAAGAAGCAAATTCGTCTGt             >  2:3677947/1‑150 (MQ=255)
                                                                                                                                         ttCCTTTTTGCCACGACAAAAAGGGACCCGCCGTAAGGGCGGAAAGGTGACTATGCGCCACCTGCGTGAATGAATACCTACTCGGTGTGAATGCCCACATCCAAGAAGCAAGAAGCAAGAAGCAAGAAGCAAGAAGCAAATTCGTCTGtt            >  2:3544637/1‑150 (MQ=255)
                                                                                                                                         ttCCTTTTTGCCACGACAAAAAGGGACCCGCCGTAAGGGCGGAAAGGTGACTATGCGCCACCTGCGTGAATGAATACCTACTCGGTGTGAAAGCCCACATCCAAGAAGCAAGAAGCAAGAAGCAAGAAGCAAGAAGCAAATTCGTCTGtt            >  2:2468991/1‑150 (MQ=255)
                                                                                                                                         ttCCTTTTTGCCACGACAAAAAGGGACCCGCCGTAAGGGCGGAAAGGTGACTAGGCGCCACCTGCGTGAATGAATACCTACTCGGTGTGAATGCCCACATCCAAGAAGCAAGAAGCAAGAAGCAAGAAGCAAGAAGCAAATTCGTCTGtt            <  1:4184433/150‑1 (MQ=255)
                                                                                                                                         tcccTTTTTGTCGTGGCAAAAAGGGACCCGCCGTAAGGGCGGAAAGGTGACTATGCGCCACCTGCGTGAATGAATACCTACTCGGTGTGAATGCCCACATCCAAGAAGCAAGAAGCAAGAAGCACGAAGCACGAAGCAAATTCGTCTGtt            >  2:2667141/3‑150 (MQ=14)
                                                                                                                                         tcccTTTTTGTCGTGGCAAAAAGGGACCCGCCGTAAGGGCGGAAAGGTGACTATGCGCCACCTGCGTGAATGAATACCTACTCGGTGTGAATGCCCACATCCAAGAAGCAAGAAGCAAGAAGCAAGAAGCAAGAAGCAAATTCGTCTGtt            >  2:789865/3‑150 (MQ=255)
                                                                                                                                         tcccTTTTTGTCGTGGCAAAAAGGGACCCGCCGTAAGGGCGGAAAGGTGACTATGCGCCACCTGCGTGAATGAATACCTACTCGGTGTGAATGCCCACATCCAAGAAGCAAGAAGCAAGAAGCAAGAAGCAAGAAGCAAATTCGTCTGtt            >  2:2195077/3‑150 (MQ=255)
                                                                                                                                         tcccTTTTTGTCGTGGCAAAAAGGGACCCGCCGTAAGGGCGGAAAGGTGACTATGCGCCACCTGCGTGAATGAATACCTACTCGGTGTGAATGCCCACATCCAAGAAGCAAGAAGCAAGAAGCAAGAAGCAAGAAGCAAATTCGTCTGtt            >  2:3901375/3‑150 (MQ=255)
                                                                                                                                         tcccTTTTTGTCGTGGCAAAAAGGAACCCGCCGTAAGGGCGGAAAGGTGACTATGCGCCACCTGCGTGAATGAATACCTACTCGGTGTGAATGCCCACATCCAAGAAGCAAGAAGCAAGAAGCAAGAAGCAAGAAGCAAAATCGTCTGtt            >  1:6330288/3‑150 (MQ=14)
                                                                                                                                          tCCTTTTTGCCACGACAAAAAGGGACCCGCCGTAAGGGCGGAAAGGTGACTATGCGCCACCTGCGTGAATGAATACCTACTCGGTGTGAATGCCCACATCCAAGAAGCAAGAAGCAAGAAGCAAGAAGCAAGAAGCAAATTCGTCTGttt           <  1:4317303/150‑1 (MQ=255)
                                                                                                                                           ccTTTTTGTCGTGGCAAAAAGGGACCCGCCGTAAGGGCGGAAAGGTGACTATGCGCCACCTGCGTGAATGAATACCTACTCGGTGTGAATGCCCACATCCAAGAAGCAAGAAGCAAGAAGCAAGAAGCAAGAAGCAAATTCGTCTGtttt          >  2:4156161/1‑150 (MQ=255)
                                                                                                                                           ccTTTTTGCCACGACAAAAAGGGACCCGCCGTAAGGGCGGAAAGGTGACTATGCGCCACCTGCGTGAATGTATACCTACTCGGTGTGAATGCCCACATCCAAGAAGCAAGAAGCAAGAAGCAAGAAGCAAGAAGCAAATTCGTCTGtttt          <  2:5603577/150‑1 (MQ=255)
                                                                                                                                            cTTTTTGCCACGACAAAAAGGGACCCGCCGTAAGGGCGGAAAGGTGACTATGCGCCACCTGCGTGAATGAATACCTACTCGGTGTGAATGCCCACATCCAAGAAGCAAGAAGCAAGAAGCAAGAAGCAAGAAGCAAATTCGTCTGTTTTg         >  2:5875062/1‑150 (MQ=255)
                                                                                                                                            cTTTTTGCCACGACAAAAAGGGACCCGCCGTAAGGGCGGAAAGGTGACTATGCGCCACCTGCGTGAATGAATACCTACTCGGTGGGAATGCCCACATCCAAGAAGCAAGAAGCAAGAAGCAAGAAGCAAGAAGCAAATTCGTCTGTTTTg         >  1:535626/1‑150 (MQ=255)
                                                                                                                                             tttttGCCACGACAAAAAGGGACCCGCCGTAAGGGCGGAAAGGTGACTATGCGCCACCTGCGTGAATGAATACCTACTCGGTGTGAATGCCCACATCCAAGAAGCAAGAAGCAAGAAGCAAGAAGCAAGAAGCAAATTCGTCTGTTTTGa        <  1:4559901/150‑1 (MQ=255)
                                                                                                                                             tttttGCCACGACAAAAAGGGACCCGCCGTAAGGGCGGAAAGGTGACTATGCGCCACCTGCGTGAATGAATACCTACTCGGTGTGAATGCCCACATCCAAGAAGCAAGAAGCAAGAAGCAAGAAGCAAGAAGCAAATTCGTCTGGTTTGa        >  1:6406099/1‑150 (MQ=255)
                                                                                                                                                 tGCCACGACAAAAAGGGACCCGCCGTAAGGGCGGAAAGGTGACTATGCGCCACCTGCGTGAATGAATACCTACTCGGTGTGAATGCCCACATCCAAGAAGCAAGAAGCAAGAAGCAAGAAGCAAGAAGCAAATTCGTCTGTTTTGACTCa    <  2:1018803/150‑1 (MQ=255)
                                                                                                                                                 tGCCACGACAAAAAGGGACCCGCCGTAAGGGCGGAAAGGTGACTATGCGCCACCTGCGTGAATGAATACCTACTCGGTGTGAATGCCCACATCCAAGAAGCAAGAAGCAAGAAGCAAGAAGCAAGAAGCAAATTCGTCTGTTTTGACTCa    <  1:2450036/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GTGACGGGCATGGATGCCCGTCAAGCGCTGGGCCCCAGGATGGGGCCTGCAGCGCGGTCCTCCCGGGAGCAAGCCCGGAGCGAGGGGACCCCGGAGCGTAGCGCAGGGGCCGGATGATGGGAGCCGACGGTTTTTGGTTCCTTTTTGCCACGACAAAAAGGGACCCGCCGTAAGGGCGGAAAGGTGACTATGCGCCACCTGCGTGAATGAATACCTACTCGGTGTGAATGCCCACATCCAAGAAGCAAGAAGCAAGAAGCAAGAAGCAAGAAGCAAATTCGTCTGTTTTGACTCAGC  >  NC_002947/2261179‑2261475

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