Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,522,227 G→T 7.0% intergenic (+71/+15) infA → / ← clpA translation initiation factor IF‑1/ATP‑dependent serine protease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,522,2270GT 7.0% 807.0 / 38.1 242intergenic (+71/+15)infA/clpAtranslation initiation factor IF‑1/ATP‑dependent serine protease
Reads supporting (aligned to +/- strand):  ref base G (132/92);  new base T (11/6);  total (143/99)
Fisher's exact test for biased strand distribution p-value = 7.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.84e-01

Only 100 of 248 total aligned reads displayed.

CTCACTGGCGACAAGGTCCGCGTCGAACTGACGCCGTACGACCTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCATCGCGCAGGTCGACATGCACCACGCCACCATGCTCGGCCAGCTCGCCAAACAGGATCTCCTCGGCCAGCGGCCGCTTGATCTTGTCCTG  >  NC_002947/4522079‑4522376
                                                                                                                                                    |                                                                                                                                                     
ctcACTGGCGACAAGGTCCGCGTCGAACTGACGCCGTACGACCTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAgg                                                                                                                                                      <  1:5566804/150‑1 (MQ=255)
  cACTGGCGACAAGGTCCGCGTCGAACTGACGCCGTACGACCTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAACAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCa                                                                                                                                                    >  2:5789858/1‑150 (MQ=255)
  cACTGGCGACAAGGTCCGCGTCGAACTGACGCCGTACGACCTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCa                                                                                                                                                    <  1:1667973/150‑1 (MQ=255)
  cACTGGCGACAAGGTCCGCGTCGAACTGACGCCGTACGACCTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCa                                                                                                                                                    >  2:2393167/1‑150 (MQ=255)
     tGGCGACAAGGTCCGCGTCGAACTGACGCCGTACGACCTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGcc                                                                                                                                                 >  1:2743656/1‑150 (MQ=255)
      ggCGACAAGGTCCGCGTCGAACTGACGCCGTACGACCTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGccc                                                                                                                                                <  2:3482518/150‑1 (MQ=255)
      ggCGACAAGGTCCGCGTCGAACTGACGCCGTACGACCTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGccc                                                                                                                                                <  1:5879303/150‑1 (MQ=255)
        cGACAAGGTCCGCGTCGAACTGACGCCGTACGACCTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCCTTGCAGCCCCt                                                                                                                                              <  1:5101805/150‑1 (MQ=255)
        cGACAAGGTCCGCGTCGAACTGACGCCGTACGACCTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCt                                                                                                                                              <  1:6602216/150‑1 (MQ=255)
        cGACAAGGTCCGCGTCGAACTGACGCCGTACGACCTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCt                                                                                                                                              <  1:2950914/150‑1 (MQ=255)
         gACAAGGTCCGCGTCGAACTGACGCCGTACGACCTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTg                                                                                                                                             <  2:5867202/150‑1 (MQ=255)
         gACAAGGTCCGCGTCGAACTGACGCCGTACGACCTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTg                                                                                                                                             >  2:4086852/1‑150 (MQ=255)
             aGGTCCGCGTCGAACTGACGCCGTACGACCTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGcaca                                                                                                                                         >  2:3074629/1‑150 (MQ=255)
               gTCCGCGTCGAACTGACGCCGTACGACCTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATc                                                                                                                                       <  1:4408847/150‑1 (MQ=255)
                   gcgTCGAACTGACGCCGTACGACCTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCAcgc                                                                                                                                   <  1:1973839/150‑1 (MQ=255)
                   gcgTCGAACTGACGCCGTACGACCTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCAcgc                                                                                                                                   <  1:7432599/150‑1 (MQ=255)
                   gcgTCGAACTGACGCCGTACGACCTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCAcgc                                                                                                                                   <  1:4782048/150‑1 (MQ=255)
                        gAACTGACGCCGTACGACCTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCt                                                                                                                              >  1:889982/1‑150 (MQ=255)
                         aaCTGACGCCGTACGACCTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTc                                                                                                                             >  2:3694686/1‑150 (MQ=255)
                            tGACGCCGTACGACCTGAGCAATGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGc                                                                                                                          >  1:552983/1‑150 (MQ=255)
                            tGACGCCGTACGACCTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGc                                                                                                                          >  2:3262169/1‑150 (MQ=255)
                             gACGCCGTACGACCTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCACGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGcc                                                                                                                         <  2:1658598/150‑1 (MQ=255)
                               cGCCGTACGACCTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGt                                                                                                                       <  2:3818486/150‑1 (MQ=255)
                                gCCGTACGACCTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTg                                                                                                                      >  1:6376817/1‑150 (MQ=255)
                                 ccGTACGACCTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTgg                                                                                                                     >  2:60499/1‑150 (MQ=255)
                                   gTACGACCTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGtt                                                                                                                   >  2:1340015/1‑150 (MQ=255)
                                        aCCTGAGCAAGGTCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGTCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCTACCTCAGCCGTGGTTTCGaa                                                                                                              <  2:692345/150‑1 (MQ=255)
                                        aCCTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGaa                                                                                                              <  1:3826034/150‑1 (MQ=255)
                                        aCCTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGaa                                                                                                              <  1:6178684/150‑1 (MQ=255)
                                        aCCTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGaa                                                                                                              <  2:2743656/150‑1 (MQ=255)
                                          cTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCCGCCCCTGCACATCACGCGACCTCAGCCGGGGGGTCGAAgg                                                                                                            >  2:4445824/1‑149 (MQ=255)
                                          cTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGt                                                                                                            >  2:731771/1‑150 (MQ=255)
                                           tGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTc                                                                                                           <  2:1229276/150‑1 (MQ=255)
                                            gAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCg                                                                                                          >  1:7162186/1‑150 (MQ=255)
                                              gCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGAGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGaa                                                                                                        >  1:2763970/1‑150 (MQ=255)
                                               cAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAac                                                                                                       <  2:1739521/150‑1 (MQ=255)
                                                aaGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAaca                                                                                                      <  1:6696173/150‑1 (MQ=255)
                                                 aGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAacac                                                                                                     >  2:4492493/1‑150 (MQ=255)
                                                  gggCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACAcc                                                                                                    >  2:560627/1‑150 (MQ=255)
                                                   ggCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCa                                                                                                   <  1:5971910/150‑1 (MQ=255)
                                                    gCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAg                                                                                                  >  1:4735242/1‑150 (MQ=255)
                                                     ccGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTGTTGTGCCTCAAGACATTGCAGGGGCTGCAAGGCCGCCCCTGCACATCACGCGACCTCAGCCGTGGTGTCGAAGTCGAACACCAGt                                                                                                 >  2:320410/1‑150 (MQ=255)
                                                            aCCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCa                                                                                                   >  2:6629893/1‑141 (MQ=255)
                                                             ccTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATGGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCATc                                                                                         <  2:1408898/150‑1 (MQ=255)
                                                                 ccGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCATCgcgc                                                                                     <  1:3489559/150‑1 (MQ=255)
                                                                    tGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCATCGCGCAgg                                                                                  >  2:2014420/1‑150 (MQ=255)
                                                                    tGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTGTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCATCGCGCAgg                                                                                  >  1:2248604/1‑150 (MQ=255)
                                                                        cgcTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCCTTGCAGCCCCTGc                                                                                                                                            >  2:2451406/1‑88 (MQ=39)
                                                                        cgcTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCCTTGCAGCCCCTGc                                                                                                                                            <  1:2451406/88‑1 (MQ=39)
                                                                           tAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCATCGCGCAGGTCGACAt                                                                           >  1:6545956/1‑150 (MQ=255)
                                                                           tAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCATCGCGCAGGTCGACAt                                                                           >  2:2868803/1‑150 (MQ=255)
                                                                                tCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCATCGCGCAGGTCGACATGcacc                                                                      >  1:3679749/1‑150 (MQ=255)
                                                                                tCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCATCGCGCAGGTCGACATGcacc                                                                      >  2:3587784/1‑150 (MQ=255)
                                                                                tCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACAACACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCATCGCGCAGGTCGACATGcacc                                                                      >  1:1047941/1‑150 (MQ=255)
                                                                                  aaGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCATCGCGCAGGTCGACATGcaccac                                                                    <  2:2280220/150‑1 (MQ=255)
                                                                                  aaGCCATGAAAAAGCCCGGCAAGGTGCCGGGCTTTTTTGTGCCTCAAGACATAGCAGGGGCTGCAAGGCAGCACCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCAGCGCGCAGGTCGACATGcaccac                                                                    >  2:5403650/1‑150 (MQ=255)
                                                                                        tGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCATCGCGCAGGTCGACATGCACCACGccacc                                                              >  1:5690741/1‑150 (MQ=255)
                                                                                         gAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCATCGCGCAGGTCGACATGCACCACGccacca                                                             >  1:6140355/1‑150 (MQ=255)
                                                                                         gAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGGGGTTTCGAAGTCGAACACCAGTTCGCCATCGCGCAGGTCGACATGCACCACGccacca                                                             >  2:855461/1‑150 (MQ=255)
                                                                                             aaGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAg                                                                                                                                                       >  2:6306069/1‑56 (MQ=255)
                                                                                             aaGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAg                                                                                                                                                       <  1:6306069/56‑1 (MQ=255)
                                                                                              aGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCATCGCGCAGGTCGACATGCACCACGCCACCATGCTc                                                        <  2:1744010/150‑1 (MQ=255)
                                                                                               gCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCATCGCGCAGGTCGACATGCACCACGCCACCATGCTCg                                                       >  1:4035899/1‑150 (MQ=255)
                                                                                                 ccGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCATCGCGCAGGTCGACATGCACCACGCCACCATGCTCGGc                                                     >  1:5828811/1‑150 (MQ=255)
                                                                                                  cGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCATCGCGCAGGTCGACATGCACCACGCCACCATGCTCGGcc                                                    <  1:3185593/150‑1 (MQ=255)
                                                                                                       aGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCATCGCGCAGGTCGACATGCACCACGCCACCATGCTCGGCCAGCTc                                               >  2:6525084/1‑150 (MQ=255)
                                                                                                       aGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCATCGCGCAGGTCGACATGCACCACGCCACCATGCTCGGCCAGCTc                                               >  2:2628873/1‑150 (MQ=255)
                                                                                                       aGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCATCGCGCAGGTCGACATGCACCACGCCACCATGCTCGGCCAGCTc                                               <  1:3285252/150‑1 (MQ=255)
                                                                                                                gCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCATCGCGCAGGTCGACATGCACCACGCCACCATGCTCGGCCAGCTCGCCAAACAg                                      >  2:1507118/1‑150 (MQ=255)
                                                                                                                 caatattGTGCCACAAGAAATTGCAGCGGCAGCAACGCAGCCCCAGCACAACACGCGACCACACCCGTCGTTTCGAAGTCCAACACCAGTTCGCCATCGCGCAGGTCGACAAGCACCACGCCACCATCCTCCCCCAGCTCGCCAAACAgg                                     <  2:4258793/145‑1 (MQ=255)
                                                                                                                 cTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCATCGCGCAGGTCGACATGCACCACGCCACCATGCTCGGCCAGCTCGCCAAACAgg                                     <  1:3636105/150‑1 (MQ=255)
                                                                                                                 cTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCATCGCGCAGGTCGACATGCACCACGCCACCATGCTCGGCCAGCTCGCCAAACAgg                                     <  2:718075/150‑1 (MQ=255)
                                                                                                                  ttttttGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCATCGCGCAGGTCGACATGCACCACGCCACCATGCTCGGCCAGCTCGCCAAACAGGa                                    <  2:5038982/150‑1 (MQ=255)
                                                                                                                       tgtgCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCATCGCGCAGGTCGACATGCACCACGCCACCATGCTCGGCCAGCTCGCCAAACAGGATctcc                               <  1:6770058/150‑1 (MQ=255)
                                                                                                                       tgtgCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCATCGCGCAGGTCGACATGCACCACGCCACCATCCTCGGCCAGCTCGCCAAACAGGATctcc                               >  1:6003368/1‑150 (MQ=255)
                                                                                                                          gCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCATCGCGCAGGTCGACATGCACCACGCCACCATGCTCGGCCAGCTCGCCAAACAGGATCTCCTCg                            >  1:1460745/1‑150 (MQ=255)
                                                                                                                           ccTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCATCGCGCAGGTCGACATGCACCACGCCACCATGCTCGGCCAGCTCGCCAAACAGGATCTCCTCgg                           >  1:1282400/1‑150 (MQ=255)
                                                                                                                            cTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCATCGCGCAGGTCGACATGCACCACGCCACCATGCTCGGCCAGCTCGCCAAACAGGATCTCCTCGGc                          >  1:4157848/1‑150 (MQ=255)
                                                                                                                            cTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCATCGCGCAGGTCGACATGCACCACGCCACCATGCTCGGCCAGCTCGCCAAACAGGATCTCCTCGGc                          >  2:595835/1‑150 (MQ=255)
                                                                                                                             tCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCATCGCGCAGGTCGACATGCACCACGCCACCATGCTCGGCCAGCTCGCCAAACAGGATCTCCTCGGcc                         >  1:3119200/1‑150 (MQ=255)
                                                                                                                             tCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCATCGCGCAGGTCGACATGCACCACGCCACCATGCTCGGCCAGCTCGCCAAACAGGATCTCCTCGGcc                         >  1:5847235/1‑150 (MQ=255)
                                                                                                                             tCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCATCGCGCAGGTCGACATGCACCACGCCACCATGCTCGGCCAGCTCGCCAAACAGGATCTCCTCGGcc                         <  1:2036005/150‑1 (MQ=255)
                                                                                                                               aaGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCATCGCGCAGGTCGACATGCACCACGCCACCATGCTCGGCCAGCTCGCCAAACAGGATCTCCTCGGCCAg                       >  2:6735564/1‑150 (MQ=255)
                                                                                                                               aaGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCATCGCGCAGGTCGACATGCACCACGCCACCATGCTCGGCCAGCTCGCCAAACAGGATCTCCTCGGCCAg                       >  1:7326432/1‑150 (MQ=255)
                                                                                                                                 gAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCATCGCGCAGGTCGACATGCACCACGCCACCATGCTCGGCCAGCTCGCCAAACAGGATCTCCTCGGCCAGCg                     >  2:4500171/1‑150 (MQ=255)
                                                                                                                                   aaTTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCATCGCGCAGGTCGACATGCACCACGCCACCATGCTCGGCCAGCTCGCCAAACAGGATCTCCTCGGCCAGCGGc                   <  1:1507118/150‑1 (MQ=255)
                                                                                                                                        cAGGGGCTGCCTTGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCATCGCGCAGGTCGACATGCACCACGCCACCATGCTCGGCCAGCTCGCCAAACAGGATCTCCTCGGCCAGCGGCCGCtt              >  2:892966/1‑150 (MQ=255)
                                                                                                                                        cAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCATCGCGCAGGTCGACATGCACCACGCCACCATGCTCGGCCAGCTCGCCAAACAGGATCTCCTCGGCCAGCGGCCGCtt              <  2:7279914/150‑1 (MQ=255)
                                                                                                                                          ggggCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCATCGCGCAGGTCGACATGCACCACGCCACCATGCTCGGCCAGCTCGCCAAACAGGATCTCCTCGGCCAGCGGCCGCTTGa            >  2:7213556/1‑150 (MQ=255)
                                                                                                                                          ggggCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCATCGCGCAGGTCGACATGCACCACGCCACCATGCTCGGCCAGCTCGCCAAACAGGATCTCCTCGGCCAGCGGCCGCTTGa            >  1:3365740/1‑150 (MQ=255)
                                                                                                                                           gggCTGCCTTGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCATCGCGCAGGTCGACATGCACCACGCCACCATGCTCGGCCAGCTCGCCAAACAGGATCTCCTCGGCCAGCGGCCGCTTGAt           >  1:516866/1‑150 (MQ=255)
                                                                                                                                           gggCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCATCGCGCAGGTCGACATGCACCACGCCACCATGCTCGGCCAGCTCGCCAAACAGGATCTCCTCGGCCAGCGGCCGCTTGag           >  2:4737339/1‑149 (MQ=255)
                                                                                                                                           gggCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCATCGCGCAGGTCGACATGCACCACGCCACCATGCTCGGCCAGCTCGCCAAACAGGATCTCCTCGGCCAGCGGCCGCTTGAt           >  2:2230763/1‑150 (MQ=255)
                                                                                                                                           gggCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCATCGCGCAGGTCGACATGCACCACGCCACCATGCTCGGCCAGCTCGCCAAACAGGATCTCCTCGGCCAGCGGCCGCTTGAt           <  2:7080970/150‑1 (MQ=255)
                                                                                                                                           gggCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCATCGCGCAGGTCGACATGCACCACGCCACCATGCTCGGCCAGCTCGCCAAACAGGATCTCCTCGGCCAGCGGCCGCTTGAt           <  1:4006298/150‑1 (MQ=255)
                                                                                                                                            ggCTGCAAGGCAGTCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCATCGCGCAGGTCGACATGCACCACGCCACCATGCTCGGCCAGCTCGCCAAACAGGATCTCCTCGGCCAGCGGCCGCTTGATc          <  2:608102/150‑1 (MQ=255)
                                                                                                                                               tGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCATCGCGCAGGTCGACATGCACCACGCCACCATGCTCGGCCAGCTCGCCAAACAGGATCTCCTCGGCCAGCGGCCGCTTGATCTTg       >  2:2164978/1‑150 (MQ=255)
                                                                                                                                               tGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCATCGCGCAGGTCGACATGCACCACGCCACCATGCTCGGCCAGCTCGCCAAACAGGATCTCCTCGGCCAGCGGCCGCTTGATCTTg       >  1:1841184/1‑150 (MQ=255)
                                                                                                                                               tGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCATCGCGCAGGTCGACATGCACCACGCCACCATGCTCGGCCAGCTCGCCAAACAGGATCTCCTCGGCCAGCGGCCGCTTGATCTTg       >  1:1129837/1‑150 (MQ=255)
                                                                                                                                                  aaGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCATCGCGCAGGTCGACATGCACCACGCCACCATGCTCGGCCAGCTCGCCAAACAGGATCTCCTCGGCCAGCGGCCGCTTGATCTTGTcc    >  2:6533254/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CTCACTGGCGACAAGGTCCGCGTCGAACTGACGCCGTACGACCTGAGCAAGGGCCGCATCACCTACCGTGCGCGCTAAGCTCAAGCCATGAAAAAGCCCGGCCAGGTGCCGGGCTTTTTTGTGCCTCAAGAAATTGCAGGGGCTGCAAGGCAGCCCCTGCACATCACGCGACCTCAGCCGTGGTTTCGAAGTCGAACACCAGTTCGCCATCGCGCAGGTCGACATGCACCACGCCACCATGCTCGGCCAGCTCGCCAAACAGGATCTCCTCGGCCAGCGGCCGCTTGATCTTGTCCTG  >  NC_002947/4522079‑4522376

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: