Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,521,819 A→G 8.8% intergenic (‑34/‑33) lip ← / → PP_4855 lipase/osmotically‑inducible lipoprotein OsmE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,521,8190AG 8.8% 187.9 / 7.5 68intergenic (‑34/‑33)lip/PP_4855lipase/osmotically‑inducible lipoprotein OsmE
Reads supporting (aligned to +/- strand):  ref base A (30/32);  new base G (2/4);  total (32/36)
Fisher's exact test for biased strand distribution p-value = 6.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01

CGAGCATGCCTGGCACCAGCACCAGCGGGTATCGCGTTGCCAAATCCTGCTGCATTGCCGCAGCTCCTTTGCGTGTACGACGGGCCAACACTACAGCGGCCATCATGACCATCGCAATCGAACTCCTGGCGCGCGCTGTGGTTCAAAACCACACAGACCCTTTGCAAGGAGCGGGACCATGTACAAGCAGACCCTGGCAATCCTTCTGGCAAGCGCCACCCTGGCCGCCTGCGGCAGCCGCCCGGAAAACCCGGTGGATTACGTCACCTACCGCGATGAGCCGCTGGTCAA  >  NC_002947/5521674‑5521964
                                                                                                                                                 |                                                                                                                                                 
cGAGCATGCGTGGCACCAGCACCAGCGGGTATCGCGTTGCCAAATCCTGCTGCATTGCCGCAGCTCCTTTTCGTGTACGACTGTCCTACACTACAGCGGCCATCATGACCATCGCAATCGAACTCCTGGCGCGCGCTGTGGTTCAAAAc                                                                                                                                                <  6:191780/149‑1 (MQ=255)
    aattccTGGCACCAGCACCATCGTGTATCGCGTTGCCAAATCCTGCTGCATTGCCGCTGCTCCTTTGCGTGTACGACGGGACAACACTACAGCGGCCATCATGACCATCGCAATCGAACTCCTGGCGCGCGCTGTGGTTTTGAACcaca                                                                                                                                            <  1:266817/145‑1 (MQ=255)
     aTGCCTGGCACCAGCACCAGCGGGTATCGCGTTGCCAAATCCTGCTGCATTGCCGCAGCTCCTTTGCGTGTACGACGGGCCAACACTACAGCGGCCATCATGACCATCGCAATCGAACTCCTGGCGCGCGCTGTGGTTCAAAACcaca                                                                                                                                            >  5:84113/1‑148 (MQ=255)
      tGCCTGGCACCAGCACCAGCGGGTATCGCGTTGCCAAATCCTGCTGCATTGCCGCAGCTCCTTTGCGTGTACGACGGGCCAACACTACAGCGGCCATCATGACCATCGCAATCGAACTCCTGGCGCGCGCTGTGGTTTTGAACcacagc                                                                                                                                          <  5:370331/149‑3 (MQ=255)
        ccTGGCACCAGCACCAGCGGGTATCGCGTTGCCAAATCCTGCTGCATTGCCGCAGCTCCTTTGCGTGTACGACGGGCCAACACTACAGCGGCCATCATGACCATCGCAATCGAACTCCTGGCGCGCGCTGTGGTTCAAAACCACACAGa                                                                                                                                        <  6:84113/149‑1 (MQ=255)
         cTGGCACCAGCACCAGCGGGTATCGCGTTGCCAAATCCTGCTGCATTGCCGCAGCTCCTTTGCGTGTACGACGGGCCAACACTACAGCGGCCATCACGACCATCGCAATCGAACTCCTGGCGCGCGCTGTGGTTCAAAACCACACAGAc                                                                                                                                       >  2:200763/1‑149 (MQ=255)
              accagcaccagcGGGTATCGCGTTGCCAAATCCTGCTGCATTGCCGCAGCTCCTTTGCGTGTACGACGGGCCAACACTACAGCGGCCATCATGACCATCGCAATCGAACTCCTGGCGCGCGCTGTGGTTCAAAACCACACAGACCCttt                                                                                                                                  <  7:50690/149‑1 (MQ=255)
              accagcaccagcGGGTATCGCGTTGCCAAATCCTGCTGCATTGCCGCAGCTCCTTTGCGTGTACGACGGGCCAACACTACAGCGGCCATCATGACCATCGCAATCGAACTCCTGGCGCGCGCTGTGGTTCAAAACCACACAGACCCttt                                                                                                                                  <  2:108955/149‑1 (MQ=255)
              accagcaccagcGGGTATCGCGTTGCCAAATCCTGCTGCATTGCCGCAGCTCCTTTGCGTGTACGACGGGCCAACACTACAGCGGCCATCATGACCATCGCAATCGAACTCCTGGCGCGCGCGGTGGTTCAAAACCACACAGACCCttt                                                                                                                                  <  4:9029/149‑1 (MQ=255)
                cagcaccagcGGGTATCGCGTTGCCAAATCCTGCTGCATTGCCGCAGCTCCTTTGCGTGTACGACGGGCCAACACTACAGCGGCCATCATGACCATCGCAATCGAACTCCTGGCGCGCGCTGTGGTTTTGAACcacagc                                                                                                                                          >  2:97730/1‑137 (MQ=255)
                cagcaccagcGGGTATCGCGTTGCCAAATCCTGCTGCATTGCCGCAGCTCCTTTGCGTGTACGACGGGCCAACACTACAGCGGCCATCATGACCATCGCAATCGAACTCCTGGCGCGCGCTGTGGTTTTGAACcacagc                                                                                                                                          <  1:97730/139‑3 (MQ=255)
                 agcaccagcGGGTATCGCGTTGCCAAATCCTGCTGCATTGCCGCAGCTCCTTTGCGTGTACGACGGGCCAACACTACAGCGGCCATCATGACCATCGCAATCGAACTCCTGGCGCGCGCTGTGGTTCAAAACCACACAGACCCTTTGCa                                                                                                                               <  7:346880/149‑1 (MQ=255)
                                 gcgTTGCCAAATCCTGCTGCATTGCCGCAGCTCCTTTGCGTGTACGACGGGCCAACACTACAGCGGCCATCATGACCATCGCAATCGAACTCCTGGCGCGCGCTGTGGTTCAAAACCACACAGACCCTTTGCAAGGAGCGGGACCATGt                                                                                                               <  6:166256/149‑1 (MQ=255)
                                 gcgTTGCCAAATCCTGCTGCATTGCCGCAGCTCCTTTGCGTGTACGACGGGCCAACACTACAGCGGCCATCATGACCATCGCAATCGAACTCCTGGCGCGCGCTGTGCTTCAAAACCACACAGACCCTTTGCAAGGAGCGGGACCATg                                                                                                                <  5:106080/148‑1 (MQ=255)
                                      gCCAAATCCTGCTGCATTGCCGCAGCTCCTTTGCGTGTACGACGGGCCAACACTACAGCGGCCATCATGACCATCGCAATCGAACTCCTGGCGCGCGCTGTGGTTCAAAACCACACAGACCCTTTGCAAGGAGCGGGACCATGTACaa                                                                                                           <  5:203363/148‑1 (MQ=255)
                                       ccAAATCCTGCTGCATTGCCGCAGCTCCTTTGCGTGTACGACGGGCCAACACTACAGCGGCCATCATGACCATCGCAATCGAACTCCTGGCGCGCGCTGTGGTTCAAAACCACACAGACCCTTTGCAAGGAGCGGGACCATGTACAAg                                                                                                          >  5:239575/1‑148 (MQ=255)
                                       ccAAATCCTGCTGCATTGCCGCAGCTCCTTTGCGTGTACGACGGGCCAACACTACAGCGGCCATCATGACCATCGCAATCGAACTCCTGGCGCGCGCTGTGGTTCAAAACCACACAGACCCTTTGCAAGGAGCGGGACCATGTACAAg                                                                                                          <  4:312340/148‑1 (MQ=255)
                                           aTCCTGCTGCATTGCCGCAGCTCCTTTGCGTGTACGACGGGCCAACACTACAGCGGCCATCATGACCATCGCAATCGAACTCCTGGCGCGCGCTGTGGTTCAAAACCACACAGACCCTTTGCAAGGAGCGGGACCATGt                                                                                                               >  5:210904/1‑139 (MQ=255)
                                           aTCCTGCTGCATTGCCGCAGCTCCTTTGCGTGTACGACGGGCCAACACTACAGCGGCCATCATGACCATCGCAATCGAACTCCTGGCGCGCGCTGTGGTTCAAAACCACACAGACCCTTTGCAAGGAGCGGGACCATGt                                                                                                               <  6:210904/139‑1 (MQ=255)
                                              ctgctgCATTGCCGCAGCTCCTTTGCGTGTACGACGGGCCAACACTACAGCGGCCATCATGACCATCGCAATCGAACTCCTGGCGCGCGCTGTGGTTCAAAACCACACAGACCCTTTGc                                                                                                                                <  1:168339/119‑1 (MQ=255)
                                              ctgctgCATTGCCGCAGCTCCTTTGCGTGTACGACGGGCCAACACTACAGCGGCCATCATGACCATCGCAATCGAACTCCTGGCGCGCGCTGTGGTTCAAAACCACACAGACCCTTTGc                                                                                                                                >  2:168339/1‑119 (MQ=255)
                                               tgctgcATTGCCGCAGCTCCTTTGCGTGTACGACGGGCCAACACTACAGCGGCCATCATGACCATCGCAATCGAACTCCTGGCGCGCGCTGTGGTTCAAAACCACACAGACCCTTTGCAAGGAGCGGGACCATGTACAAGCAGACCCTg                                                                                                 >  3:325671/1‑149 (MQ=255)
                                                 ctgcATTGCCGCAGCTCCTTTGCGTGTACGACGGGCCAACACTACAGCGGCCATCACGACCATCGCAATCGAACTCCTGGCGCGCGCTGTGGTTCAAAACCACACAGACCCTTTGCAAGGAGCGGGACCATGTACAAGCAGACCCTgg                                                                                                <  1:200763/148‑1 (MQ=255)
                                                   gcATTGCCGCAGCTCCTTTGCGTGTACGACGGGCCAACACTACAGCGGCCATCATGACCATCGCAATCGAACTCCTGGCGCGCGCTGTGGTTCAAAACCACACAGACCCTTTGCAAGGAGCGGGACCATGTACAAGCAGACCCTGGCaa                                                                                             >  7:20628/1‑149 (MQ=255)
                                                     aTTGCCGCAGCTCCTTTGCGTGTACGACGGGACAACACTACAGCGGCCATCATGACCATCGCAATCGAACTCCTGGCGCGCGCTGTGGTTCAAAACCACACAGACCCTTTGCAAGGAGCGGGACCATGTACAAGCAGACCCTGGCAATc                                                                                           >  1:40960/1‑149 (MQ=255)
                                                           gCAGCTCCTTTGCGTGTACGACGGGCCAACACTACAGCGGCCATCATGACCATCGCAATCGAACTCCTGGCGCGCGCTGTGGTTCAAAACCACACAGACCCTTTGCAAGGAGCGGGACCATGTACAAGCAGACCCTGGCAATCCTTCTg                                                                                     >  3:373917/1‑149 (MQ=255)
                                                                           tACGACGGGCCAACACTACAGCGGCCATCATGACCATCGCAATCGAACTCCTGGCGCGCGCTGTGGTTCAAAACCACACAGACCCTTTGCAAGGAGCGGGACCATGTACAAGCAGACCCTGGCAATCCTTCTGGCAAGCGCCACCCTgg                                                                     >  3:358200/1‑149 (MQ=255)
                                                                            acgacgGGCCAACACTACAGCGGCCATCATGACCATCGCAATCGAACTCCTGGCGCGCGCTGTGGTTCAAAACCACACAGACCCTTTGCAAGGAGCGGGACCATGTACAAGCAGACCCTGGCAATCCTTCTGGCAAGCGCCACCCTGgc                                                                    >  7:378741/1‑149 (MQ=255)
                                                                              gacgGGCCAACACTACAGCGGCCATCATGACCATCGCAATCGAACTCCTGGCGCGCGCTGTGGTTCAAAACCACACAGACCCTTTGCAAGGAGCGGGACCATGTACAAGCAGACCCTGGCAATCCTTCTGGCAAGCGCCACCCTGgccg                                                                  >  6:418046/1‑149 (MQ=255)
                                                                              gacgGGCCAACACTACAGCGGCCATCATGACCATCGCAATCGAACTCCTGGCGCGCGCTGTGGTTCAAAACCACACAGACCCTTTGCAAGGAGCGGGACCATGTACAAGCAGACCCTGGCAATCCTTCTGGCAAGCGCCACCCTGgccg                                                                  >  2:364493/1‑149 (MQ=255)
                                                                                   gCCAACACTACAGCGGCCATCATGACCATCGCAATCGAACTCCTGGCGCGCGCTGTGGTTCAAAACCACACAGACCCTTTGCAAGGAGCGGGACCATGTACAAGCAGACCCTGGCAATCCTTCTGGCAAGCGCCACCCTGGCCGCCTGc                                                             >  2:178912/1‑149 (MQ=255)
                                                                                    ccAACACTACAGCGGCCATCATGACCATCGCAATCGAACTCCTGGCGCGCGCTGTGGTTCAAAACCACACAGACCCTTTGCAAGGAGCGGGACCATGTACAAGCAGACCCTGGCAATCCTTCTGGCAAGCGCCACCCTGGCCGCCTGCg                                                            >  7:25835/1‑149 (MQ=255)
                                                                                     cAACACTACAGCGGCCATCATGACCATCGCAATCGAACTCCTGGCGCGCGCTGTGGTTCAAAACCACACAGACCCTTTGCAAGGAGCGGGACCATGTACAAGCAGACCCTGGCAATCCTTCTGGCAAGCGCCACCCTGGCCGCCTGCgg                                                           >  1:169330/1‑149 (MQ=255)
                                                                                              aGCGGCCATCATGACCATCGCAATCGAACTCCTGGCGCGCGCTGTGGTTCAAAACCACACAGACCCTTTGCAAGGAGCGGGACCATGTACAAGCAGACCCTGGCAATCCTTCTGGCAAGCGCCACCCTGGCCGCCTGCGGCAGCCGccc                                                  <  4:358200/149‑1 (MQ=255)
                                                                                               gCGGCCATCATGACCATCGCAATCGAACTCCTGGCGCGCGCTGTGGTTCAAAACCACACAGACCCTTTGCAAGGAGCGGGACCATGTACAAGCAGACCCTGGCAATCCTTCTGGCAAGCGCCACCCTGGCCGCCTGCGGCAGCCGCCCg                                                 <  4:373917/149‑1 (MQ=255)
                                                                                                    catcatGACCATCGCAATCGAACTCCTGGCGCGCGCTGTGGTTCAAAACCACACAGACCCTTTGCAAGGAGCGGGACCATGTACAAGCAGACCCTGGCAATCCTTCTGGCAAGCGCCACCCTGGCCGCCTGCGGCAGCCGCCCGGaaaa                                            >  8:308834/1‑149 (MQ=255)
                                                                                                     atcatGACCATCGCAATCGAACTCCTGGCGCGCGCTGTGGTTCAAAACCACACAGACCCTTTGCAAGGAGCGGGACCATGTACAAGCAGACCCTGGCAATCCTTCTGGCAAGCGCCACCCTGGCCGCCTGCGGCAGCCGCCCGGAAAAc                                           >  4:320057/1‑149 (MQ=255)
                                                                                                       catGACCATCGCAATCGAACTCCTGGCGCGCGCTGTGGTTCAAAACCACACAGACCCTTTGCAAGGAGCGGGACCATGTACAAGCAGACCCTGGCAATCCTTCTg                                                                                     <  1:108833/105‑1 (MQ=255)
                                                                                                       catGACCATCGCAATCGAACTCCTGGCGCGCGCTGTGGTTCAAAACCACACAGACCCTTTGCAAGGAGCGGGACCATGTACAAGCAGACCCTGGCAATCCTTCTg                                                                                     >  2:108833/1‑105 (MQ=255)
                                                                                                        atGACCATCGCAATCGAACTCCTGGCGCGCGCTGTGGTTCAAAACCACACAGACCCTTTGCAAGGAGCGGGACCATGTACAAGCAGACCCTGGCAATCCTTCTGGCAAGCGCCACCCTGGCCGCCTGCGGCAGCCGCCCGGAAAAccc                                         <  1:178912/148‑1 (MQ=255)
                                                                                                          gACCATCGCAATCGAACTCCTGGCGCGCGCTGTGGTTCAAAACCACACAGACCCTTTGCAAGGAGCGGGACCATGTACAAGCAGACCCTGGCAATCCTTCTGGCAAGCGCCACCCTGGCCGCCTGCGGCAGCCGCCCGGa                                               >  1:218042/1‑140 (MQ=255)
                                                                                                          gACCATCGCAATCGAACTCCTGGCGCGCGCTGTGGTTCAAAACCACACAGACCCTTTGCAAGGAGCGGGACCATGTACAAGCAGACCCTGGCAATCCTTCTGGCAAGCGCCACCCTGGCCGCCTGCGGCAGCCGCCCGGa                                               <  2:218042/140‑1 (MQ=255)
                                                                                                          gACCATCGCAATCGAACTCCTGGCGCGCGCTGTGGTTCAAAACCACACAGACCCTTTGCAAGGAGCGGGACCATGTACAAGCAGACCCTGGCAATCCTTCTGGCAAGCGCCACCCTGGCCGCCTGCGGCAGCCGCCCGGAAAACCCGGt                                      >  3:363742/1‑149 (MQ=255)
                                                                                                           aCCATCGCAATCGAACTCCTGGCGCGCGCTGTGGTTCAAAACCACACAGACCCTTTGCAAGGAGCGGGACCATGTACAAGCAGACCCTGGCAATCCTTCTGGCAAGCGCCACCCTGGCCGCCTGCGGCAGCCGCCCGGAAAACCCGGt                                      <  8:20628/148‑1 (MQ=255)
                                                                                                           aCCATCGCAATCGAACTCCTGGCGCGCGCTGTGGTTCAAAACCACACAGACCCTTTGCAAGGAGCGGGACCATGTACAAGCAGACCCTGGCAATCCTTCTGGCAAGCGCCACCCTGGCCGCCTGCGGCAGCCGCCCGGAAAACCCGGt                                      <  1:364493/148‑1 (MQ=255)
                                                                                                              aTCGCAATCGAACTCCTGGCGCGCGCTGTGGTTCAAAACCACACAGACCCTTTGCAAGGAGCGGGACCATGTACAAGCAGACCCTGGCAATCCTTCTGGCAAGCGCCACCCTGGCCGCCTGCGGCAGCCGCCCGGAAAACCCGGTGGAt                                  <  4:345583/149‑1 (MQ=255)
                                                                                                              aTCGCAATCGAACTCCTGGCGCGCGCTGTGGTTCAAAACCACACAGACCCTTTGCAAGGAGCGGGACCATGTACAAGCAGACCCTGGCAATCCTTCTGGCAAGCGCCACCCTGGCCGCCTGCGGCAGCCGCCCGGAAAACCCGGTGGAt                                  >  4:175971/1‑149 (MQ=255)
                                                                                                                   aaTCGAACTCCTGGCGCGCGCTGTGGTTCAAAACCACACAGACCCTTTGCAAGGAGCGGGACCATGTACAAGCAGACCCTGGCAATCCTTCTGGCAAGCGCCACCCTGGCCGCCTGCGGCAGCCGCCCGGAAAACCCGGTGGATTACGt                             >  6:143275/1‑149 (MQ=255)
                                                                                                                           tCCTGGCGCGCGCTGTGGTTCAAAACCACACAGACCCTTTGCAAGGAGCGGGACCATGTACAAGCAGACCCTGGCAATCCTTCTGGCAAGCGCCACCCTGGCCGCCTGCGGCAGCCGCCCGGAAAACCCg                                        <  5:347964/130‑1 (MQ=255)
                                                                                                                           tCCTGGCGCGCGCTGTGGTTCAAAACCACACAGACCCTTTGCAAGGAGCGGGACCATGTACAAGCAGACCCTGGCAATCCTTCTGGCAAGCGCCACCCTGGCCGCCTGCGGCAGCCGCCCGGAAAACCCg                                        >  6:347964/1‑130 (MQ=255)
                                                                                                                           tCCTGGCGCGCGCTGTGGTTCAAAACCACACAGACCCTTTGCAAGGAGCGGGACCATGTACAAGCAGACCCTGGCAATCCTTCTGGCAAGCGCCACCCTGGCCGCCTGCGGCAGCCGCCCGGAAAACCCGGTGGATTACGTCACCTAcc                     >  6:130712/1‑149 (MQ=255)
                                                                                                                           tCCTGGCGCGCGCTGTGGTTCAAAACCACACAGACCCTTTGCAAGGAGCGGGACCATGTACAAGCAGACCCTGGCAATCCTTCTGGCAAGCGCCACCCTGGCCGCCTGCGGCAGCCGCCCGGAAAACCCGGTGGATTACGTCACCTAcc                     >  4:289202/1‑149 (MQ=255)
                                                                                                                             cTGGCGCGCGCTGTGGTTCAAAACCACACAGACCCTTTGCAAGGAGCGGGACCATGTACAAGCAGACCCTGGCAATCCTTCTGGCAAGCGCCACCCTGGCCGCCTGCGGCAGCCGCCCGGAAAACCCGGTGGATTACGTCACCTACcgc                   >  4:225253/1‑149 (MQ=255)
                                                                                                                                gcgcgcgcTGTGGTTCAAACCCACACAGACCCTTTGCAAGGAGCGGGACCATGTACAAGCAGACCCTGGCAATCCTTCTGGCAAGCGCCACCCTGGCCGCCTGCGGCAGCCGCCCGGAAAACCCGGTGGATTACGTCACCTACCGCGAt                <  3:289202/149‑1 (MQ=255)
                                                                                                                                gcgcgcgcTGTGGTTCAAAACCACACAGACCCTTTGCAAGGAGCGGGACCATGTACAAGCAGACCCTGGCAATCCTTCTGGCAAGCGCCACCCTGGCCGCCTGCGGCAGCCGCCCGGAAAACCCgg                                       <  6:356285/126‑1 (MQ=255)
                                                                                                                                gcgcgcgcTGTGGTTCAAAACCACACAGACCCTTTGCAAGGAGCGGGACCATGTACAAGCAGACCCTGGCAATCCTTCTGGCAAGCGCCACCCTGGCCGCCTGCGGCAGCCGCCCGGAAAACCCgg                                       >  5:356285/1‑126 (MQ=255)
                                                                                                                                gcgcgcgcTGTGGTTCAAAACCACACAGACCCTTTGCAAGGAGCGGGACCATGTACAAGCAGACCCTGGCAATCCTTCTGGCAAGCGCCACCCTGGCCGCCTGCGGCAGCCGCCCGGAAAACCCGGTGGATTACGTCACCTACCGCGa                 <  3:175971/148‑1 (MQ=255)
                                                                                                                                gcgcgcgcTGTGGGTCAAAACCACACAGACCCTTTGCAAGGAGCGGGACCATGTACAAGCAGACCCTGGCAATCCTTCTGGCAAGCGCCACCCTGGCCGCCTGCGGCAGCCGCCCGGAAAACCCGGTGGATTACGTCACCTACCGCGAt                <  7:153522/149‑1 (MQ=255)
                                                                                                                                 cgcgcgcTGTGGTTCAAAACCACACAGACCCTTTGCAAGGAGCGGGACCATGTACAAGCAGACCCTGGCAATCCTTCTGGCAAGCGCCACCCTGGCCGCCTGCGGCAGCCGCCCGGAAAACCCGGTGGATTACGTCACCTACCGCGATg               <  5:143275/149‑1 (MQ=255)
                                                                                                                                      tgtgtgGTTTTGAACCACACAGACCCTTTGCAAGGAGCGGGACCATGTACAAGCAGACCCTGGCAATCCTTCTGGCAAGCGCCACCCTGGCCGCCTGCGGCAGCCGCCCGGAAAACCCGGTGGATTa                                <  7:161716/125‑1 (MQ=255)
                                                                                                                                      tgtgtgGTTTTGAACCACACAGACCCTTTGCAAGGAGCGGGACCATGTACAAGCAGACCCTGGCAATCCTTCTGGCAAGCGCCACCCTGGCCGCCTGCGGCAGCCGCCCGGAAAACCCGGTGGATTa                                >  8:161716/3‑127 (MQ=255)
                                                                                                                                       cTGTGGTTCAAAACCACACAGACCCTTTGCAAGGAGCGGGACCATGTACAAGCAGACCCTGGCAATCCTTCTGGCAAGCGCCACGCTGGCCGCCTGCGGCAGCCGCCCGGAAAACCCGGTGGATTACGTCACCTACCGCGATGAGCCGc         <  4:363742/149‑1 (MQ=255)
                                                                                                                                          tgGTTCAAAACCACACAGACCCTTTGCAAGGAGCGGGACCATGTACAAGCAGACCCTGGCAATCCTTCTGGCAAGCGCCACCCTGGCCGCCTGCGGCAGCCGCCCGGAAAACCCGGTGGATTACGTCACCTACCGCGATGAGCCGCTg       <  2:169330/148‑1 (MQ=255)
                                                                                                                                            gTTCAAAACCACACAGACCCTTTGCAAGGAGCGGGACCATGTACAAGCAGACCCTGGCAATCCTTCTGGCAAGCGCCACCCTGGCCGCCTGCGGCAGCCGCCCGGAAAACCCGGTGGATTACGTCACCTACCGCGATGAGCCGCTGGt     >  3:308498/1‑148 (MQ=255)
                                                                                                                                            gTTCAAAACCACACAGACCCTTTGCAAGGAGCGGGACCATGTACAAGCAGACCCTGGCAATCCTTCTGGCAAGCGCCACCCTGGCCGCCTGCGGCAGCCGCCCGGAAAACCCGGTGGATTACGTCACCTACCGCGATGAGCCGCTGGTc    <  5:418046/149‑1 (MQ=255)
                                                                                                                                             ttCAAAACCACACAGACCCTTTGCAAGGAGCGGGACCATGTACAAGCAGACCCTGGCAATCCTTCTGGCAAGCGCCACCCTGGCCGCCTGCGGCAGCCGCCCGGAAAACCCGGTGGATTACGTCACCTACCGCGAt                >  7:108946/1‑136 (MQ=255)
                                                                                                                                             ttCAAAACCACACAGACCCTTTGCAAGGAGCGGGACCATGTACAAGCAGACCCTGGCAATCCTTCTGGCAAGCGCCACCCTGGCCGCCTGCGGCAGCCGCCCGGAAAACCCGGTGGATTACGTCACCTACCGCGAt                <  8:108946/136‑1 (MQ=255)
                                                                                                                                               cAAAACCACACAGACCCTTTGCAAGGAGCGGGACCATGTACAAGCAGACCCTGGCAATCCTTCTGGCAAGCGCCACCCTGGCCGCCTGCGGCAGCCGCCCGGAAAACCCGGTGGATTACGTCACCTACCGCGATGAGCCGCTGGTCaa  <  4:308498/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
CGAGCATGCCTGGCACCAGCACCAGCGGGTATCGCGTTGCCAAATCCTGCTGCATTGCCGCAGCTCCTTTGCGTGTACGACGGGCCAACACTACAGCGGCCATCATGACCATCGCAATCGAACTCCTGGCGCGCGCTGTGGTTCAAAACCACACAGACCCTTTGCAAGGAGCGGGACCATGTACAAGCAGACCCTGGCAATCCTTCTGGCAAGCGCCACCCTGGCCGCCTGCGGCAGCCGCCCGGAAAACCCGGTGGATTACGTCACCTACCGCGATGAGCCGCTGGTCAA  >  NC_002947/5521674‑5521964

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: