Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,041,854 C→G 6.4% intergenic (+148/‑59) PP_3562 → / → PP_3563 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,041,8540CG 6.4% 196.8 / 3.8 62intergenic (+148/‑59)PP_3562/PP_3563hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (29/29);  new base G (2/2);  total (31/31)
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.82e-01

ACGCAAACAGCTGACAACGCTGGCCCGATGCGAATGACTGGCACCTGACCGCACTTTTCCACTGGCCTGGGGTCTACTCCCTCTTTTGCCGCCAGTAACGATTGCCGGTGTTTTCATTGCCTTGCACGACCTCGGGCCAAAAGCCCCGGGCTTGTGTAGGATGAGCAGCGCACACCACGACGCTGCCAGTCATAGGGAGATCCACATGCGCGTCCTGACATCCGTTGTCGCCTTGTCACTGGGCCTGGTCGTTTCGACCGGGGCGCTGGCAGCCACCGGCGAAGAGGCGCAGTT  >  NC_002947/4041708‑4042001
                                                                                                                                                  |                                                                                                                                                   
acGCAAACAGCTGACAACGCTGGCCCGATGCGAATGACTGGCACCTGACCGCACTTTTCCACTGGCCTGGGGTCTACTCCCTCTTTTGCCGCCAGTAACGATTGCCGGTGTTTTCATTGCCTTGCACGACCTCGGGCCAAAAGCCCCgg                                                                                                                                                   >  7:128078/1‑149 (MQ=255)
      aCAGCTGACAACGCTGGCCCGATCCGAATGACTGGCACCTGACCGCACTTTTCCACTGGCCTGGGGTCTACTCCCTCTTTTGCCGCCAGTAACGATTGCCGGTGTTTTCATTGCCTTGCACGACCTCGGGCCAAAAGCCCCGGGCTtgt                                                                                                                                             <  3:262457/149‑1 (MQ=255)
         gCTGACAACGCTGGCCCGATGCGAATGACTGGCACCTGACCGCACTTTTCCACTGGCCTGGGGTCTACTCCCTCTTTTGCCGCCAGTAACGATTGCCGGTGTTTTCATTGCCTTGCACGACCTCGGGCCAAAAGCCCCGGGCTTGTGTa                                                                                                                                          <  4:37963/149‑1 (MQ=255)
          cTGACAAGGCTGGCCCGATGCGAATGACTGGCACCTGACCGCACTTTTCCACTGGCCTGGGGTCTACTCCCTCTTTTGCCGCCAGTAACGATTGCCGGTGTTTTCATTGCCTTGCACGACCTCGGGCCAAAAGCCCCGGGCTTGTGTAg                                                                                                                                         <  1:128886/149‑1 (MQ=255)
          cTGACAACGCTGGCCCGATGCGAATGACTGGCACCTGACCGCACTTTTCCACTGGCCTGGGGTCTACTCCCTCTTTTGCCGCCAGTAACGATTGCCGGTGTTTTCATTGCCTTGCACGACCTCGGGCCAAAAGCCCCGGGCTTGTGTAg                                                                                                                                         >  5:202267/1‑149 (MQ=255)
             aCAACCCTGGCCCGATGCGAATGACTGGCACCTGACCGCACTTTTCCACTGGCCTGGGGTCTACTCCCTCTTTTGCCGCCAGTAACGATTGCCGGTGTTTTCATTGCCTTGCACGACCTCGGGCCAAAAGCCCCGGGCTTGTGTAGGa                                                                                                                                       <  7:255869/148‑1 (MQ=255)
                 cGCTGGCCCGATGCGAATGACTGGCACCTGACCGCACTTTTCCACTGGCCTGGGGTCTACTCCCTCTTTTGCCGCCAGTAACGATTGCCGGTGTTTTCATTGCCTTGCACGACCTCGGGCCAAAAGCCCCGGGCTTGTGTAGGATGagc                                                                                                                                  <  2:307329/149‑1 (MQ=255)
                      gCCCGATGCGAATGACTGGCACCTGACCGCACTTTTCCACTGGCCTGGGGTCTACTCCCTCTTTTGCCGCCAGTAACGATTGCCGGTGTTTTCATTGCCTTGCACGACCTCGGGCCAAAAGCCCCGCGCTTGTGTAGGATGAGCAgcg                                                                                                                              <  7:51528/148‑1 (MQ=255)
                                  tGACTGGCAACTGACCGCACTTTTCCACTGGCCTGGGGTCTACTCCCTCTTTTGCCGCCAGTAACGATTGCCGGTGTTTTCATTGCCTTGCACGACCTCGGGCCAAAAGCCCCGGGCTTGTGTAGGATGAGCAGCGCACACCACGACGc                                                                                                                 <  7:201967/149‑1 (MQ=255)
                                     cTGGCACCTGACCGCACTTTTCCACTGGCCTGGGGTCTACTCCCTCTTTTGCCGCCAGTAACGATTGCCGGTGGTTTCATTGCCTTGCACGACCTCGGGCCAAAAGCCCCGGGCTTGTGTAGGATGAGCAGCGCACACCCCGACGCTGc                                                                                                              >  8:118745/1‑149 (MQ=255)
                                        gCACCTGACCGCACTTTTCCACTGGCCTGGGGTCTACTCCCTCTTTTGCCGCCAGTAACGATTGCCGGTGTTTTCATTGCCTTGCACGACCTCGGGCCAAAAGCCCCGGGCTTGTGTAGGATGAGCAGCGCACACCACGACGCTGCCAg                                                                                                           >  6:331871/1‑149 (MQ=255)
                                            cTGACCGCACTTTTCCACTGGCCTGGGGTCTACTCCCTCTTTTGCCGCCAGTAACGATTGCCGCTGTTTTCATTGCCTTGCACGACCTCGGGCCAAAAGCCCCGTGCTTGTGTAGGATGAGAAGCGCACACCACGACGCTGCCAGTCa                                                                                                        >  4:57497/1‑148 (MQ=255)
                                               aCCGCACTTTTCGACTGGCCTGGGGTCTACTCCCTCTTTTGCCGCCAGTAACGATTGCCGGTGTTTTCATTGCCTTGCACGACCTCGGGCCACAAGCCCGGGGCttt                                                                                                                                              <  8:292475/107‑2 (MQ=255)
                                               aCCGCACTTTTCCACTGGCCTGGGGTCTACTCCCTCTTTTGCCGCCAGTAACGATTGCCGGTGTTTTCATTGCCTTGCACGACCTCGGGCCACAAGCCCGGGGCttt                                                                                                                                              >  7:292475/1‑106 (MQ=255)
                                               aCCGCACTTTTCCACTGGCCTGGGGTCTACTCCCTCTTTTGCCGCCAGTAACGATTGCCGGTGTTTTCATTGCCTTGCACGACCTCGGGCCAAAAGCCCCGGGCTTGTGTAGGATGAGCAGCGCACACCACGACGCTGCCAGTCATAg                                                                                                     <  1:322274/148‑1 (MQ=255)
                                               aCCGCACTTTTCCACTGGCCTGGGGTCTACTCCCTCTTTTGCCGCCAGTAACGATTGCCGGTGTTTTCATTGCCTTGCACGACCTCGGGCCAAAAGCCCCGGGCTTGTGTAGGATGAGCAGCGCACACCACGACGCTGCCAGTCATAg                                                                                                     >  2:322274/1‑148 (MQ=255)
                                                 cGCACTTTTCCACTGGCCTGGGGTCTACTCCCTCTTTTGCCGCCAGTAACGATTGCCGGTGTTTTCATTGCCTTGCACGACCTCGGGCCAAAAGCCCCGGGCTTGTGTAGGATGAGCAGCGCACACCACGACGCTGCCAGTCATAggca                                                                                                  >  4:280764/1‑147 (MQ=255)
                                                 cGCACTTTTCCACTGGCCTGGGGTCTACTCCCTCTTTTGCCGCCAGTAACGATTGCCGGTGTTTTCATTGCCTTGCACGACCTCGGGCCAAAAGCCCCGGGCTTGTGTAGGATGAGCAGCGCACACCACGACGCTGCCAGTCATAgg                                                                                                    >  5:141392/1‑147 (MQ=255)
                                                      ttttCCACTGGCCTGGGGTCTACTCCCTCTTTTGCCGCCAGTAACGATTGCCGGTGTTTTCATTGCCTTGCACGACCTCGGGCCAAAAGCCCCGGGCTTGTGTAGGATGAGCAGCGCACACCACGACGCTGCCAGTCATAGGGAGATc                                                                                              >  7:297247/1‑148 (MQ=255)
                                                      ttttCCACTGGCCTGGGGTCTACTCCCTCTTTTGCCGCCAGTAACGATTGCCGGTGTTTTCATTGCCTTGCACGACCTCGGGCCAAAAGCCCCGGGCTTGTGTAGGATGAGCAGCGCACACCACGACGCTGCCAGTCATAGGGAGATc                                                                                              >  2:114039/1‑148 (MQ=255)
                                                        ttCCACTGGCCTGGGGTCTACTCCCTCTTTTGCCGCCAGTAACGATTGCCGGTGTTTTCATTGCCTTGCACGACCTCGGGCCAAAAGCCCCGGGCTTGTGTAGGATGAGCAGCGCACACCACGACGCTGCCAGTCATAGGGAGATcc                                                                                             <  8:297247/147‑1 (MQ=255)
                                                         tCCACTGGCCTGGGGTCTACTCCCTCTTTTGCCGCCAGTAACGATTGCCGGTGTTTTCATTGCCTTGCACGACCTCGGGCCAAAAGCCCCGGGCTTGTGTAGGATGAGCAGCGCACACCACGACGCTGCCAGTCATAGGGAGATCcac                                                                                           <  8:128078/148‑1 (MQ=255)
                                                                        tCTACTCCCTCTTTTGCCGCCAGTAACGATTGCCGGTGTTTTCATTGCCTTGCACGACCTCGGGCCAAAAGCCCCGGGCTTGTGTAGGATGAGCAGCGCACACCACGACGCTGCCAGTCATAGGGAGATCCACATgcgc                                                                                     >  3:113284/1‑139 (MQ=255)
                                                                        tCTACTCCCTCTTTTGCCGCCAGTAACGATTGCCGGTGTTTTCATTGCCTTGCACGACCTCGGGCCAAAAGCCCCGGGCTTGTGTAGGATGAGCAGCGCACACCACGACGCTGCCAGTCATAGGGAGATCCACATgcgc                                                                                     <  4:113284/139‑1 (MQ=255)
                                                                          tACTCCCTCTTTTGCCGCCAGTAACGATTGCCGGTGTTTTCATTGCCTTGCACGACCTCGGGCCAAAAGCCCCGGGCTTGTGTAGGATGAGCAGCGCACACCACGACGCTGCCAGTCATAGGGAGATCCACATGCGCGTCCTGACATcc                                                                         >  6:184016/1‑149 (MQ=255)
                                                                            cTCCCTCTTTTGCCGCCAGTAACGATTGCCGGTGTTTTCATTGCCTTGCACGACCTCGGGCCAAAAGCCCCGGGCTTGTGTAGGATGAGCAGCGCACACCACGACGCTGCCAGTCATAGGGAGATCCACATGCGCGTCCTGACATCCg                                                                        <  5:331871/148‑1 (MQ=255)
                                                                               ccTCTTTTGCCGCCAGTAACGATTGCCGGTGTTTTCATTGCCTTGCACGACCTCGGGCCAAAAGCCCCGGGCTTGTGTAGGATGAGCAGCGCACACCACGACGCTGCCAGTCATAGGGAGATCCACATGCGCGTCCTGACATCCGTTGt                                                                    >  1:118930/1‑149 (MQ=255)
                                                                                     ttGCCGCCAGTAACGATTGCCGGTGTTTTCATTGCCTTGCACGACCTCGGGCCAAAAGCCCCGGGCTTGTGTAGGATGAGCAGCGCACACCACGACGCTGCCAGTCATAGGGAGATCCACATGCGCGTCCTGACATCCGTTGTCGcc                                                                <  3:280764/147‑1 (MQ=255)
                                                                                     ttGCCGCCAGTAACGATTGCCGGTGTTTTCATTGCCTAGCACGACCTCGGGCCAAAAGCCCCGGGCTTGTGTAGGATGAGCAGCGCACACCACGACGCTGCCAGTCATAGGGAGATCCACATGCGCGTCCTGACATCCGTTGTCGc                                                                 <  8:255381/146‑1 (MQ=255)
                                                                                      tGCCGCCAGTAACGATTGCCGGTGTTTTCATTGCCTTGCACGACCTCGGGCCAAAAGCCCCGGGCTTGTGTAGGATGAGCAGCGCACACCACGACGCTGCCAGTCATAGGGAGATCCACATGCGCGTCCTGACATCCGTTGTCGCCtt                                                              <  8:148618/148‑1 (MQ=255)
                                                                                        ccgccAGTAACGATTGCCGGTGTTTTCATTGCCTTGCACGACCTCGGGCCAAAAGCCCCGGGCTTGTGTAGGATGAGCAGCGCACACCACGACGCTGCCAGTCATAGGGAGATCCACATGCGCGTCCTGACATCCGTTGTCGCCTTGt                                                            <  6:202267/148‑1 (MQ=255)
                                                                                           ccAGTAACGATTGCCGGTGTTTTCATTGCCTTGCACGACCTCGGGCCAAAAGCCCCGGGCTTGTGTAGGATGAGCAGCGCACACCACGACGCTGCCAGTCACAGGGAGATCCACATGCGCGTCCTGACATCCGTTGTCGCCTTGTCACt                                                        <  1:81754/149‑1 (MQ=255)
                                                                                              gTAACGATTGCCGGTGTTTTCATTGCCTTGCACGACCTCGGGCCAAAAGCCCCGGGCTTGTGTAGGATGAGCAGCGCACACCACGACGCTGCCAGTCATAGGGAGATCCACATGCGCGTCCTGACATCCGTTGTCGCCTTGTCACTggg                                                     >  4:51410/1‑149 (MQ=255)
                                                                                                  cGATTGCCGGTGTTTTCATTGCCTTGCACGACCTCGGGCCAAAAGCCCCGGGCTTGTGTAGGATGAGCAGCGCACACCACGACGCTGCCAGTCATAGGGAGATCCACATGCGCGTCCTGACATCCGTTGTCGCCTTGTCACTGGGCCTg                                                 >  4:30340/1‑149 (MQ=255)
                                                                                                    aTTGCCGGTGTTTTCATTGCCTTGCACGACCTCGGGCCAAAAGCCCCGGGCTTGTGTAGGATGAGCAGCGCACACCACGACGCTGCCAGTCATAGGGAGATCCACATGCGCGTCCTGACATCCGTTGTCGCCTTGTCACTGGGCCTGGt                                               <  3:67128/149‑1 (MQ=255)
                                                                                                        ccGGTGTTTTCATTGCCTTGCACGACCTCGGGCCAAAAGCCCCGGGCTTGTGTAGGATGAGCAGCGCACACCACGACGCTGCCAGTCATAGGGAGATCCACATGCGCGTCCTGACATCCGTTGTCGCCTTGTCACTGGGCCTGGTCGtt                                           >  4:61961/1‑149 (MQ=255)
                                                                                                        ccGGTGTTTTCATTGCCTTGCACGACCTCGGGCCAAAAGCCCCGGGCTTGTGTAGGATGAGCAGCGCACACCACGACGCTGCCAGTCATAGGGAGATCCACATGCGCGTCCTGACATCCGTTGTCGCCTTGTCACTGGGCCTGGTCGtt                                           >  3:129216/1‑149 (MQ=255)
                                                                                                        ccGGTGTTTTCATTGCCTTGCACGACCTCGGGCCAAAAGCCCCGGGCTTGTGTAGGATGAGCAGCGCACACCACGACGCTGCCAGTCATAGGGAGATCCACATGCGCGTCCTGACATCCGTTGTCGCCTTGTCACTGGGCCTGGTCGtt                                           >  6:102869/1‑149 (MQ=255)
                                                                                                         cGGTGTTTTCATTGCCTTGCACGACCTCGGGCCAAAAGCCCGGGGCTTGTGTAGGATGAGCAGCGCACACCACGACGCTGCCAGTCATAGGGAGATCCACATGCGCGTCCTGACATCCGtt                                                                      >  3:181050/1‑121 (MQ=255)
                                                                                                         cGGTGTTTTCATTGCCTTGCACGACCTCGGGCCAAAAGCCCGGGGCTTGTGTAGGATGAGCAGCGCACACCACGACGCTGCCAGTCATAGGGAGATCCACATGCGCGTCCTGACATCCGtt                                                                      <  4:181050/121‑1 (MQ=255)
                                                                                                                ttCATTGCCTTGCACGACCTCGGGCCAAAAGCCCCGGGCTTGTGTAGGATGAGCAGCGCACACCACGACGCTGCCAGTCATAggg                                                                                                   >  5:171409/1‑85 (MQ=255)
                                                                                                                ttCATTGCCTTGCACGACCTCGGGCCAAAAGCCCCGGGCTTGTGTAGGATGAGCAGCGCACACCACGACGCTGCCAGTCATAggg                                                                                                   <  6:171409/85‑1 (MQ=255)
                                                                                                                ttCATTGCCTTGCACGACCTCGGGCCAAAAGCCCCGGGCTTGTGTAGGATGAGCAGCGCACACCACGACGCTGCCAGTCATAGGGAGATCCACATGCGCGTCCTGACATCCGTTGTCGCCTTGTCACTGGGCCTGGTCGTTTCGACCgg                                   >  3:2740/1‑149 (MQ=255)
                                                                                                                   aTTGCCTTGCACGACCTCGGGCCATAAGCCCCGGGCTTGTGTAGGATGAGCAGCGCACACCACGACGCTGCCAGTCATAGGGAGATCCACATGCGCGTCCTGACATCCGTTGTCGCCTTGTCACTGGGCCTGGTCGTTTCGACCGGGgc                                <  3:21938/149‑1 (MQ=255)
                                                                                                                   aTTGCCTTGCACGACCTCGGGCCAAAAGCCCCGGGCTTGTGTAGGATGAGCAGCGCACACCACGACGCTGCCAGTCATAGGGAGATCCACATGCGCGTCCTGACATCCGTTGTCGCCTTGTCACTGGGCCTGGTCGTTTCGACCgggg                                 <  5:95656/148‑1 (MQ=255)
                                                                                                                       ccTTGCACGACCTCGGGCCAAAAGCCCCGGGCTTGTGTAGGATGAGCAGCGCACACCACGACGCTGCCAGTCATAGGGAGATCCACATGCGCGTCCTGACATCCGTTGTCGCCTTGTCACTGGGCCTGGTCGTTTCGACCGGGGCGCTg                            >  3:69498/1‑149 (MQ=255)
                                                                                                                        cTTGCACGACCTCGGGCCAAAAGCCCCGGGCTTGTGTAGGATGAGCAGCGCACACCacgacg                                                                                                                  <  2:265190/62‑1 (MQ=255)
                                                                                                                        cTTGCACGACCTCGGGCCAAAAGCCCCGGGCTTGTGTAGGATGAGCAGCGCACACCACGACGCTGCCAGTCATAGGGAGATCCACATGCGCGTCCTGACATCCGTTGTCGCCTTGTCACTGGGCCTGGTCGTTTCGACCGGGGCGCTgg                           >  4:62704/1‑149 (MQ=255)
                                                                                                                        cTTGCACGACCTCGGGCCAAAAGCCCCGGGCTTGTGTAGGATGAGCAGCGCACACCACGACGCTGCCAGTCATAGGGAGATCCACATGCGCGTCCTGACATCCGTTGTCGCCTTGTCACTGGGCCTGGTCGTTTCGACCGGGGCGCTgg                           >  3:243520/1‑149 (MQ=255)
                                                                                                                           gCACGACCTCGGGCCAAAAGCCCCGGGCTTGTGTAGGATGAGCAGCGCACACCACGACGCTGCCAGTCATAGGGAGATCCACATGCGCGTCCTGACATCCGTTGTCGCCTTGTCACTGGGCCTGGTCGTTTCGACCGGGGCGCTGGCa                         <  5:184016/148‑1 (MQ=255)
                                                                                                                           gCACGACCTCGGGCCAAAAGCCCCGGGCTTGTGTAGGATGAGCAGCGCACACCACGACGCTGCCAGTCATAGGGAGATCCACATGCGCGTCCTGACATCCGTTGTCGCCTTGTCACTGGGCCTGGTCGTTTCGACCGGGGCGCTGGCa                         <  5:102869/148‑1 (MQ=255)
                                                                                                                            cACGACCTCGGGCCAAAAGCCCCGGGCTTGTGTAGGATGAGCAGCGCACACCACGACGCTGCCAGTCATAGGGAGATCCACATGCGCGTCCTGACATCCGTTGTCGCCTTGTCACTGGGCCTGGTCGTTTCGACCGGGGCGCTGGCAg                        <  3:51410/148‑1 (MQ=255)
                                                                                                                            cACGACCTCGGGCCAAAAGCCCCGGGCTTGTGTAGGATGAGCAGCGCACACCACGACGCTGCCAGTCATAGGGAGATCCACATGCGCGTCCTGACATCCGTTGTCGCCTTGTCACTGGGCCTGGTCGTTTCGACCGGGGCGCTGGCAGc                       >  5:54792/1‑149 (MQ=255)
                                                                                                                                  cTCGGGCCAAAAGCCCCGGGCTTGTGTAGGATGAGCAGCGCACACCACGACGCTGCCAGTCATAGGGAGATCCACATGCGCGTCCTGACATCCGTTGTCGCCTTGTCACTGGGCCTGGTCGTTTCGACCGGGGCGCTGGCAGCCAc                    <  6:12060/146‑1 (MQ=255)
                                                                                                                                  cTCGGGCCAAAAGCCCCGGGCTTGTGTAGGATGAGCAGCGCACACCACGACGCTGCCAGTCATAGGGAGATCCACATGCGCGTCCTGACATCCGTTGTCGCCTTGTCACTGGGCCTGGTCGTTTCGACCGGGGCGCTGGCAGCCAc                    >  5:12060/1‑146 (MQ=255)
                                                                                                                                  cTCGGGCCAAAAGCCCCGGGCTTGTGTAGGATGAGCAGCGCACACCACGACGCTGCCAGTCATAGGGAGATCCACATGCGCGTCCTGACATCCGTTGTCGCCTTGTCACTGGGCCTGGTCGTTTCGACCGGGGCGCTGGCAGCCACCgg                 >  2:216692/1‑149 (MQ=255)
                                                                                                                                   tCGGGCCAAAAGCCCCGGGCTTGTGTAGGATGAGCAGCGCACACCACGACGCTGCCAGTCATAGGGAGATCCACATGCGCGTCCTGACATCCGTTGTCGCCTTGTCACTGGGCCTGGTCgaaa                                          <  7:174144/123‑4 (MQ=255)
                                                                                                                                   tCGGGCCAAAAGCCCCGGGCTTGTGTAGGATGAGCAGCGCACACCACGACGCTGCCAGTCATAGGGAGATCCACATGCGCGTCCTGACATCCGTTGTCGCCTTGTCACTGGGCCTGGTCgaaa                                          >  8:174144/1‑120 (MQ=255)
                                                                                                                                    cGGGCCAAAAGCCCCGGGCTTGTGTAGGATGAGCAGCGCACACCACGACGCTGCCAGTCATAGGGAGATCCACATGCGCGTCCTGACATCCGTTGTCGCCTTGTCACTGGGCCTGGTCGTTTCGACCGGGGCGCTGGCAGCCACCGGc                >  8:19464/1‑148 (MQ=255)
                                                                                                                                      ggCCAAAAGCCCCGGGCTTGTGTAGGATGAGCAGCGCACACCACGACGCTGCCAGTCATAGGGAGATCCACATGCGCGTCCTGACATCCGTTGTCGCCTTGTCACTGGGCCTGGTCGTTTCGACCGGGGCGCTGGCAGCCACCGGCGaa             <  4:243520/149‑1 (MQ=255)
                                                                                                                                       gCCAAAAGCCCCGGGCTTGTGTAGGATGAGCAGCGCACACCACGACGCTGCCAGTCATAGGGAGATCCACATGCGCGTCCTGACATCCGTTGTCGCCTTGTCACTGGGCCTGGTCGTTTCGACCGGGGCGCTGGCAGCCACCGGCg               >  2:180788/1‑146 (MQ=255)
                                                                                                                                              gCCCCGGGCTTGTGTAGGATGAGCAGCGCACACCACGACGCTGCCAGTCATAGGGAGATCCACATGCGCGTCCTGACATCCGTTGTCGCCTTGTCACTGGGCCTGGTCGTTTCGACCGGGGCGCTGGCAGCCACCGGCGAAGAGGCGCa     <  1:114039/149‑1 (MQ=255)
                                                                                                                                                cccGGGCTTGTGTAGGATGAGCAGCGCACACCACGACGCTGCCAGTCATAGGGAGATCCACATGCGCGTCCTGACATCCGTTGTCGCCTTGTCACTGGGCCTGGTCGTTTCGACCGGGGCGCTGGCAGCCACCGGCGAAGAGGCGCAGt   <  2:118930/149‑1 (MQ=255)
                                                                                                                                                 ccGGGCTTGTGTAGGATGAGCAGCGCACACCACGACGCTGCCAGTCATAGGGAGATCCACATGCGCGTCCTGACATCCGTTGTCGCCTTGTCACTGGGCCTGGTCGTTTCGACCGGGGCGCTGGCAGCCACCGGCGAAGAGGCGCAGtt  >  8:35391/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
ACGCAAACAGCTGACAACGCTGGCCCGATGCGAATGACTGGCACCTGACCGCACTTTTCCACTGGCCTGGGGTCTACTCCCTCTTTTGCCGCCAGTAACGATTGCCGGTGTTTTCATTGCCTTGCACGACCTCGGGCCAAAAGCCCCGGGCTTGTGTAGGATGAGCAGCGCACACCACGACGCTGCCAGTCATAGGGAGATCCACATGCGCGTCCTGACATCCGTTGTCGCCTTGTCACTGGGCCTGGTCGTTTCGACCGGGGCGCTGGCAGCCACCGGCGAAGAGGCGCAGTT  >  NC_002947/4041708‑4042001

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