Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,082,621 G→C 5.4% G189A (GGC→GCC)  PP_3592 → RpiR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,082,6210GC 5.4% 299.7 / 6.8 92G189A (GGC→GCC) PP_3592RpiR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (46/41);  new base C (2/3);  total (48/44)
Fisher's exact test for biased strand distribution p-value = 6.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.67e-01

ACAACGTGGCGCAACGCCTGGCGAGCAAGCCGCGGGTGTTCATCGCCGGCTTCCAGACCGAGCGCGGCATCGCCATGTGCATGAGCCACCTGCTGCAATACCTGCGCGATGGCGTGCAACTGGTCGACATCAGCGCCGGGCACTTCGGCGAAGTGCTGCTGGGCCGTGCCGAAGACAGCGCCCTGGTGGTGTTCGAGGCGCGCCGTTATTCCCGCCATGCCCTGATGCTGTGCCAGAAGGCGCGCGAG‑GCTGGCATTGCGGTCACCCTGGTGACCGACACCTTCTGTGACTGGGCG  >  NC_002947/4082474‑4082769
                                                                                                                                                   |                                                                                                                                                     
acaACGTGGCGCAACGCCTGGCGAGCAAGCCGCGGGTGTTCATCGCCGGCTTCCAGACCGAGCGCGGCATCGCCATGTGCATGAGCCACCTGCTGCAATACCTGCGCGATGGCGTGCAACTGGTCGACATCAGCGCCGGGCACTTCGGc                                                                                                                                                      >  4:135448/1‑149 (MQ=255)
  aaCGTGGCGCAACGCCTGGCGAGCAAGCCGCGGGTGTTCATCGCCGGCTTCCAGACCGAGCGCGGCATCGCCATGTGCATGAGCCACCTGCTGCAATACCTGCGCGATGGCGTGCAACTGGTCGACATCAGCGCCGGGCACTTCGGCg                                                                                                                                                     >  3:425618/1‑148 (MQ=255)
   aCGTGGCGCAACGCCTGGCGAGCAAGCCGCGGGTGTTCATCGCCGGCTTCCAGACCGAGCGCGGCATCGCCATGTGCATGAGCCACCTGCTGCAATACCTGCGCGATGGCGTGCAACTGGTCGACATCAGCGCCGGGCACTTCGGCGaa                                                                                                                                                   <  6:269852/149‑1 (MQ=255)
   aCGTGGCGCAACGCCTGGCGAGCAAGCCGCGGGTGTTCATCGCCGGCTTCCAGACCGAGCGCGGCATCGCCATGTGCATGAGCCACCTGCTGCAATACCTGCGCGATGGCGTGCAACTGGTCGACATCAGCGCCGGGCACTTCGCCGaa                                                                                                                                                   <  6:121603/149‑1 (MQ=255)
   aCGTGGCGCAACGCCTGGCGAGCAAGCCGCGGGTGTTCATCGCCGGCTTCCAGACCGAGCGCGGCATCGCCATGTGCATGAGCCACCTGCTGCAATACCTGCGCGATGGCGTGCAACTGGTCGACATCAGCGCCGGGCACTTCGCCGaa                                                                                                                                                   <  6:127697/149‑1 (MQ=255)
    cGTGGCGCAACGCCTGGCGAGCAAGCCGCGGGTGTTCATCGCCGGCTTCCAGACCGAGCGCGGCATCGCCATGTGCATGAGCCACCTGCTGCAATACCTGCGCGATGGCGTGCAACTGGTCGACATCAGCGCCGGGCACTTCGGCGAAg                                                                                                                                                  <  4:130205/149‑1 (MQ=255)
    cGTGGCGCAACGCCTGGCGAGCAAGCCGCGGGTGTTCATCGCCGGCTTCCAGACCGAGCGCGGCATCGCCATGTGCATGAGCCACCTGCTGCAATACCTGCGCGATGGCGTGCAACTGGTCGACATCAGCGCCGGGCACTTCGGCGAAg                                                                                                                                                  >  8:146468/1‑149 (MQ=255)
     gTGGCGCAACGCCTGGCGAGCAAGCCGCGGGTGTTCATCGCCGGCTTCCAGACCGAGCGCGGCATCGCCATGTGCATGAGCCACCTGCTGCAATACCTGCGCGATGGCGTGCAACTGGTCGACATCAGCGCCGGGCACTTCGGCGAAGt                                                                                                                                                 >  3:218550/1‑149 (MQ=255)
     gTGGCGCAACGCCTGGCGAGCAAGCCGCGGGTGTTCATCGCCGGCTTCCAGACCGAGCGCGGCATCGCCATGTGCATGAGCCACCTGCTGCAATACCTGCGCGATGGCGTGCAACTGGTCGACATCAGCGCCGGGCACTTCGGCGAAGt                                                                                                                                                 >  6:372712/1‑149 (MQ=255)
             aCGCCTGGCGAGCAAGCCGCGGGTGTTCATCGCCGGCTTCCAGACCGAGCGCGGCATCGCCATGTGCATGAGCCACCTGCTGCAATACCTGCGCGATGGCGTGCAACTGGTCGACATCAGCGCCGGGCACTTCGGCGAAGTGCTGCTgg                                                                                                                                         <  4:86578/149‑1 (MQ=255)
              cGCCTGGCGAGCAAGCCGCGGGTGTTCATCGCCGGCTTCCAGACCGAGCGCGGCATCGCCATGTGCATGAGCCACCTGCTGCAATACCTGCGCGATGGCGTGCAACTGGTCGACATCAGCGCCGGGCACTTCGGCGAAGTGCTGCTggg                                                                                                                                        <  5:165178/149‑1 (MQ=255)
              cGCCTGGCGAGCAAGCCGCGGGTGTTCATCGCCGGCTTCCAGACCGAGCGCGGCATCGCCATGTGCATGAGCCACCTGCTGCAATACCTGCGCGATGGCGTGCAACTGGTCGACATCAGCGCCGGGCACTTCGGCGAAGTGCTGCTggg                                                                                                                                        >  8:232808/1‑149 (MQ=255)
              cGCCTGGCGAGCAAGCCGCGGGTGTTCATCGCCGGCTTCCAGACCGAGCGCGGCATCGCCATGTGCATGAGCCACCTGCTGCAATACCTGCGCGATGGCGTGCAACTGGTCGACATCAGCGCCGGGCACTTCGGCGAAGTGCTGCTggg                                                                                                                                        <  7:80870/149‑1 (MQ=255)
               gCCTGGCGAGCAAGCCGCGGGTGTTCATCGCCGGCTTCCAGACCGAGCGCGGCATCGCCATGTGCATGAGCCACCTGCTGCAATACCTGCGCGATGGCGTGCAACTGGTCGACATCAGCGCCGGGCACTTCGGCGAAGTGCTGCTGGGc                                                                                                                                       <  5:279050/149‑1 (MQ=255)
                    gCGAGCAAGCCGCGGGTGTTCATCGCCGGCTTCCAGACCGAGCGCGGCATCGCCATGTGCATGAGCCACCTGCTGCAATACCTGCGCGATGGCGTGCAACTGGTCGACATCAGCGCCGGGCACTTCGGCGAAGTGCTGCTGGGCCGTGc                                                                                                                                  >  8:412232/1‑149 (MQ=255)
                    gCGAGCAAGCCGCGGGTGTTCATCGCCGGCTTCCAGACCGAGCGCGGCATCGCCATGTACATGAGCCACCTGCTGCAATACCTGCGCGATGGCGTGCAACTGGTCGACATCAGCGCCGGGCACTTCGGCGAAGTGCTGCTGGGCCGTGc                                                                                                                                  >  7:15726/1‑149 (MQ=255)
                        gCAAGCCGCGGGTGTTCATCGCCGGCTTCCAGACCGAGCGCGGCATCGCCATGTGCATGAGCCACCTGCTGCAATACCTGCGCGATGGCGTGCAACTGGTCGACATCAGCGCCGGGCACTTCGGCGAAGTGCTGCTGGGCCGTGCCGaa                                                                                                                              <  1:71310/149‑1 (MQ=255)
                         cAAGCCGCGGGTGTTCATCGCCGGCTTCCAGACCGAGCGCGGCATCGCCATGTGCATGAGCCACCTGCTGCAATACCTGCGCGATGGCGTGCAACTGGTCGACATCAGCGCCGGGCACTTCGGCGAAGTGCTGCTGGGCCGTGCCGAAg                                                                                                                             <  2:125578/149‑1 (MQ=255)
                          aaGCCGCGGGTGTTCATCGCCGGCTTCCAGACCGAGCGCGGCATCGCCATGTGCATGAGCCACCTGCTGCAATACCTGCGCGATGGCGTGCAACTGGTCGACATCAGCGCCGGGCACTTCGGCGAAGTGCTGCTGGGCCGTGCCGAAGa                                                                                                                            >  6:109787/1‑149 (MQ=255)
                            gCCGCGGGTGTTCATCGCCGGCTTCCAGACCGAGCGCGGCATCGCCATGTGCATGAGCCACCTGCTGCAATACCTGCGCGATGGCGTGCAACTGGTCGACATCAGCGCCGGGCACTTCGGCGAAGTGCTGCTGGGCCGTGCCGAAGa                                                                                                                            >  2:226408/1‑147 (MQ=255)
                            gCCGCGGGTGTTCATCGCCGGCTTCCAGACCGAGCGCGGCATCGCCATGTGCATGAGCCACCTGCTGCAATACCTGCGCGATGGCGTGCAACTGGTCGACATCAGCGCCGGGCACTTCGGCGAAGTGCTGCTGGGCCGTGCCGAAGa                                                                                                                            >  4:60621/1‑147 (MQ=255)
                             ccGCGGGTGTTCATCGCCGGCTTCCAGACCGAGCGCGGCATCGCCATGTGCATGAGCCACCTGCTGCAATACCTGCGCGATGGCGTGCAACTGGTCGACATCAGCGCCGGGCACTTCGGCGAAGTGCTGCTGGGCCGTGCCGAAGACAg                                                                                                                         >  3:389562/1‑149 (MQ=255)
                               gcgGGTGTTCATCGCCGGCTTCCAGACCGAGCGCGGCATCGCCATGTGCATGAGCCACCTGCTGCAATACCTGCGCGATGGCGTGCAACTGGTCGACATCAGCGCCGGGAACTTCGGCGAAGTGCTGCTGGGCCGTGCCGAAGACAgcg                                                                                                                       >  5:439084/1‑149 (MQ=255)
                                      ttCATCGCCGGCTTCCAGACCGAGCGCGGCATCGCCATGTGCATGAGCCACCTGCTGCAATACCTGCGCGATGGCGTGCAACTGGTCGACATCAGCGCCGGGCACTTCGGCGAAGTGCTGCTGGGCCGTGCCGAAGACAGCGCCCtggt                                                                                                                >  2:193955/1‑149 (MQ=255)
                                      ttCATCGCCGGCTTCCAGACCGAGCGCGGCATCGCCATGTGCATGAGCCACCTGCTGCAATACCTGCGCGATGGCGTGCAACTGGTCGACATCAGCGCCGGGCACTTCGGCGAAGTGCTGCTGGGCCGTGCCGAAGACAGCGCCCtggt                                                                                                                <  1:165399/149‑1 (MQ=255)
                                      ttCATCGCCGGCTTCCAGACCGAGAGCGGCATCGCCATGTGCATGAGCCACCTGCTGCAATACATGCGCGATGGCGTGCAACTGCTCGACATCAGCGCCGGGCACTTCGGCTAAGGGCTGCTGGGCCGTGCCGAAGACAGCGCCCtggt                                                                                                                >  8:340075/1‑149 (MQ=255)
                                        cATCGCCGGCTTCCAGACCGAGCGCGGCATCGCCATGTGCATGAGCCACCTGCTGCAATACCTGCGCGATGGCGTGCAACTGGTCGACATCAGCGCCGGGCACTTCGGCGAAGTGCTGCTGGGCCGTGCCGAAGACAGCGCCCtggtgg                                                                                                              >  8:295808/1‑149 (MQ=255)
                                         aTCGCCGGCTTCCAGACCGAGCGCGGCATCGCCATGTGCATGAGCCACCTGCTGCAATACCTGCGCGATGGCGTGCAACTGGTCGACATCAGCGCCGGGCACTTCGGCGAAGTGCTGCTGGGCCGTGCCGAAGACAGCGCCCTggtggt                                                                                                             <  3:249046/149‑1 (MQ=255)
                                           cGCCGGCTTCCAGACCGAGCGCGGCATCGCCATGTGCATGAGCCACCTGCTGCAATACCTGCGCGATGGCGTGCAACTGGTCGACATCAGCGCCGGGCACTTCGGCGAAGTGCTGCTGGGCCGTGCCGAAGACAGCGCCCTGGTGgtgt                                                                                                           >  4:385068/1‑149 (MQ=255)
                                             ccGGCTTCCAGACCGAGCGCGGCATCGCCATGTGCATGAGCCACCTGCTGCAATACCTGCGCGATGGCGTGCAACTGGTCGACATCAGCGCCGGGCACTTCGGCGAAGTGCTGCTGGGCCGTGCCGAAGACAGCGCCCTGGTGGTGTTc                                                                                                         <  1:193955/149‑1 (MQ=255)
                                               ggCTTCCAGACCGAGCGCGGCATCGCCATGTGCATGAGCCACCTGCTGCAATACCTGCGCGATGGCGTGCAACTGGTCGACATCAGCGCCGGGCACTTCGGCGAAGTGCTGCTGGGCCGTGCCGAAGACAGCGCCCTGGTGGTGTTCGa                                                                                                       <  6:320931/149‑1 (MQ=255)
                                                gCTTCCAGACCGAGCGCGGCATCGCCATGTGCATGAGCCACCTGCTGCAATACCTGCGCGATGGCGTGCCACTGGTCGACATCAGCGCCGGGCACTTCGGCGCAGTGCTGCTGGGCCGGGCCGAAGACAGCGCCCTGGTGGGGTTCgcg                                                                                                      >  8:162825/1‑147 (MQ=255)
                                                gCTTCCAGACCGAGCGCGGCATCGCCATGTGCAGGAGCCACCTGCTGCAATACCTGCGCGATGGCGTGCAACTGGTCGACATCAGCGCCGGGCACTTCGGCGAAGTGCGGCTGGGCCGTGCCGAAGACAGCGCCCTGGTGGTGTTCGAg                                                                                                      <  8:134254/149‑1 (MQ=255)
                                                    ccAGACCGAGCGCGGCATCGCCATGTGCATGAGCCACCTGCTGCAATACCTGCGCGATGGCGTGCAACTGGTCGACATCAGCGCCGGGCACTTCGGCGAAGTGCTGCTGGGCCGTGCCGAAGACAGCGccc                                                                                                                    >  1:157765/1‑131 (MQ=255)
                                                    ccAGACCGAGCGCGGCATCGCCATGTGCATGAGCCACCTGCTGCAATACCTGCGCGATGGCGTGCAACTGGTCGACATCAGCGCCGGGCACTTCGGCGAAGTGCTGCTGGGCCGTGCCGAAGACAGCGccc                                                                                                                    <  2:157765/131‑1 (MQ=255)
                                                     cAGACCGAGCGCGGCATCGCCATGTGCATGAGCCACCTGCTGCAATACCTGCGCGATGGCGTGCAACTGGTCGACATCAGCGCCGGGCACTTCGGCGAAGTGCTGCTGGGCCGTGCCGAAGACAGCGACCTGGTGGTGTTCGAGgcgcg                                                                                                 >  8:405272/1‑149 (MQ=255)
                                                      aGACCGAGCGCGGCATCGCCATGTGCATGAGCCACCTGCTGCAATACCTGCGCGATGGCGTGCAACTGGTCGACATCAGCGCCGGGCACTTCGGCGAAGTGCTGCTGGGCCGTGCCGAAGACAGCGCCCTGGTGGTGTTCGAGgcgcgc                                                                                                <  4:65276/149‑1 (MQ=255)
                                                       gACCGAGCGCGGCATCGCCATGTGCATGAGCCACCTGCTGCAATACCTGCGCGATGGCGTGCAACTGGTCGACATCAGCGCCGGGCACTTCGGCGAAGTGCTGCTGGGCCGTGCCGAAGACAGCGCCCTGGTGGTGTTCGAGGCGCGcc                                                                                               >  8:201018/1‑149 (MQ=255)
                                                        aCCGAGCGCGGCATCGCCATGTGCATGAGCCACCTGCTGCAATACCTGCGCGATGGCGTGCAACTGGTCGACATCAGCGCCGGGCACTTCGGCGAAGTGCTGCTGGGCCGt                                                                                                                                    <  6:215920/111‑1 (MQ=255)
                                                        aCCGAGCGCGGCATCGCCATGTGCATGAGCCACCTGCTGCAATACCTGCGCGATGGCGTGCAACTGGTCGACATCAGCGCCGGGCACTTCGGCGAAGTGCTGCTGGGCCGt                                                                                                                                    >  5:215920/1‑111 (MQ=255)
                                                        aCCGAGCGCGGCATCGCCATGTGCATGAGCCACCTGCTGCAATACCTGCGCGATGGCGTGCAACTGGTCGACATCAGCGCCGGGCACTTCGGCGAAGTGCTGCTGGGCCGTGCCGAAGACAGCGCCCTGGTGGTGTTCGAGGCGCGCCg                                                                                              >  1:412498/1‑149 (MQ=255)
                                                                cgGCATCGCCATGTGCATGAGCCACCTGCTGCAATACCTGCGCGATGGCGTGCAACTGGTCGACATCAGCGCCGGGCACTTCGGCGAAGTGCTGCTGGGCCGTGCCGAAGACAGCGCCCTGGTGGTGTTCGAGGCGCGCCGTTATTccc                                                                                      <  4:218550/149‑1 (MQ=255)
                                                                 ggCATCGCCATGTGCATGAGCCACCTGCTGCAATACCTGCGCGATGGCGTGCAACTGGTCGACATCAGCGCCGGGCACTTCGGCGAAGTGCTGCTGGGCCGTGCCGAAGACAGCGCCCTGGTGGTGTTCGAGGCGCGCCGTTATTCCCg                                                                                     <  3:60621/149‑1 (MQ=255)
                                                                 ggCATCGCCATGTGCATGAGCCACCTGCTGCAATACCTGCGCGATGGCGTGCAACTGGTCGACATCAGCGCCGGGCACTTCGGCGAAGTGCTGCTGGGCCGTGCCGAAGACAGCGCCCTGGTGGTGTTCGAGGCGCGCCGTTATTCCCg                                                                                     <  4:425618/149‑1 (MQ=255)
                                                                   cATCGCCATGTGCATGAGCCACCTGCTGCAATACCTGCGCGATGGCGTGCAACTGGTCGACATCAGCGCCGGGCACTTCGGCGAAGTGCTGCTGGGCCGTGCCGAAGACAGCGCCCTGGTGGTGTTCGAGGCGCGCCGTTATTCCCGcc                                                                                   <  2:412498/149‑1 (MQ=255)
                                                                   cATCGCCATGTGCATGAGCCACCTGCTGCAATACCTGCGCGATGGCGTGCAACTGGTCGACATCAGCGCCGGGAACTTCGGCGAAGTGCTGCTGGGCCGTGCCGAAGACAGCGCCCTGGTGGTGTTCGAGGCGCGCCGTTATTCCCGc                                                                                    >  3:33284/1‑148 (MQ=255)
                                                                       gCCATGTACATGAGCCACCTGCTGCAATACCTGCGCGATGGCGTGCAACTGGTCGACATCAGCGCCGGGCACTTCGGCGAAGTGCTGCTGGGCCGTGCCGAAGACAGCGCCCTGGTGGTGTTCGAGGCGCGCCGTTATTCCCGCCATg                                                                                <  8:15726/148‑1 (MQ=255)
                                                                        ccATGTGCATGAGCCACCTGCTGCAATACCTGCGCGATGGCGTGCAACTGGTCGACATCAGCGCCGGGCACTTCGGCGAAGTGCTGCTGGGCCGTGCCGAAGACAGCGCCCtggtg                                                                                                               >  1:122176/1‑116 (MQ=255)
                                                                        ccATGTGCATGAGCCACCTGCTGCAATACCTGCGCGATGGCGTGCAACTGGTCGACATCAGCGCCGGGCACTTCGGCGAAGTGCTGCTGGGCCGTGCCGAAGACAGCGCCCtggtg                                                                                                               <  2:122176/116‑1 (MQ=255)
                                                                             tgCATGAGCCACCTGCTGCAATACCTGCGCGATGGCGTGCTACTGGTCGACATCAGCGCCGGGCACTTCGGCGAAGTGCTGCTGGGCCGTGCCGAAGACAGCGCCCTGGTGGTGTTCGAGGCGCGCCGTTATTCCCGCCATGCCCTGa                                                                          <  5:372712/148‑1 (MQ=255)
                                                                               cATGAGCCACCTGCTGCAATACCTGCGCGATGGCGTGCAACTGGTCGACATCAGCGCCGGGCACTTCGGCGAAGTGCTGCTGGGCCGTGCCGAAGACAGCGCCCTGGTGGTGTTCGAGGCGCGCCGTTATTCCCGCCATGCCCTGATGc                                                                       >  1:202400/1‑149 (MQ=255)
                                                                               cATGAGCCACCTGCTGCAATACCTGCGCGATGGCGTGCAACTGGTCGACATCAGCGCCGGGCACTTCGGCGAAGTGCTGCTGGGCCGTGCCGAAGACAGCGCCCTGGTGGTGTTAGAGGCGCGCCGTTATTCCCGCCATGCCCTGATGc                                                                       >  1:404253/1‑149 (MQ=255)
                                                                                aTGAGCCACCTGCTGCAATACGTGCGCGATGGCGTGCAACTGGTCGACATCAGCGCCGGGCACTTCGGCGAAGTGCTGCTGGGCCGTGCCGAAGACAGCGCCCTGGTGGTGTTCGAGGCGCGCCGTTATTCCCGCCATGCCCTGATGCt                                                                      <  8:326047/149‑1 (MQ=255)
                                                                                 tGAGCCACCTGCTGCAATACCTGCGCGATGGCGTGCAACTGGTCGACATCAGCGCCGGGCACTTCGGCGAAGTGCTGCTGGGCCGTGCCGAAGACAGCGCCCTGGTGGTGTTCGAGGCGCGCCGTTATTCCCGCCATGCCCTGATGCtg                                                                     >  3:41645/1‑149 (MQ=255)
                                                                                   aGCCACCTGCTGCAATACCTGCGCGATGGCGTGCAACTGGTCGACATCAGCGCCGGGCACTTCGGCGAAGTGCTGCTGGGCCGTGCCGAAGACAGCGCCCTGGTGGTGTTCGAGGCGCGCCGTTATTCCCGCCATGCCCTGATGCtgt                                                                    <  4:10312/148‑1 (MQ=255)
                                                                                   aGCCACCTGCTGCAATACCTGCGCGATGGCGTGCAACTGGTCGACATCAGCGCCGGGCACTTCGGCGAAGTGCTGCTGGGCCGTGCCGAAGACAGCGCCCTGGTGGTGTTCGAGGCGCGCCGTTATTCCCGCCATGCCCTGATGCtgt                                                                    <  7:295808/148‑1 (MQ=255)
                                                                                        ccTGCTGCAATACCTGCGCGATGGCGTGCAACTGGTCGACATCAGCGCCGGGCACTTCGGCGAAGTGCTGCTGGGCCGTGCCGAAGACAGCGCCCTGGTGGTGTTCGAGGCGCGCCGTTATTCCCGCCATGCCCTGATGCTGTGCCAg                                                               >  5:205439/1‑148 (MQ=255)
                                                                                            ctgcAATACCTGCGCGATGGCGTGCAACTGGTCGACATCAGCGCCGGGCACTTCGGCGAAGTGCTGCTGGGCCGTGCCGAAGACAGCGCCCTGGTGGTGTTCGAGGCGCGCCGTTATTCCCGCCATGCCCTGATGCTGTGCCAGAAg                                                            <  4:389562/147‑1 (MQ=255)
                                                                                               cAATACCTGCGCGATGGCGTGCAACTGGTCGACATCAGCGCCGGGCACTTCGGCGAAGTGCTGCTGGGCCGTGCCGAAGACAGCGCCCTGGTGGTGTTCGAGGCGCGCCGTTATTCCCGCCATGCCCTGATGCTGTGCCAGAAGgcgcg                                                       >  3:260303/1‑149 (MQ=255)
                                                                                                       gcgcGATGGCGTGCAACTGGTCGACATCAGCGCCGGGCACTTCGGCGAAGTGCTGCTGGGCCGTGCCGAAGACAGCGCCCTGGTGGTGTTCGAGGCGCGCCGTTATTCCCGCCATGCCCTGATGCTGTGCCAGAAGGCGCGCGAG‑GCTg                                              >  4:395364/1‑149 (MQ=255)
                                                                                                        cgcgATGGCGTGCAACTGGTCGACATCAGCGCCGGGCACTTCGGCGAAGTGCTGCTGGGCCGTGCCGAAGACAGCGCCCTGGTGGTGTTCGAGGCGCGCCGTTATTCCCGCCATGCCCTGATGCTGTGc                                                                  >  4:58663/1‑129 (MQ=255)
                                                                                                        cgcgATGGCGTGCAACTGGTCGACATCAGCGCCGGGCACTTCGGCGAAGTGCTGCTGGGCCGTGCCGAAGACAGCGCCCTGGTGGTGTTCGAGGCGCGCCGTTATTCCCGCCATGCCCTGATGCTGTGc                                                                  <  3:58663/129‑1 (MQ=255)
                                                                                                        cgcgATGGCGTGCAACTGGTCGACATCAGCGCCGGGCACTTCGGCGAAGTGCTGCTGGGCCGTGCCGAAGACAGCGCCCTGGTGGTGTTCGAGGCGCGCCGTTATTCCCGCCATGCCCTGATGCTGTGCCAGAAGGCGCGCGAG‑GCTgg                                             <  7:412232/149‑1 (MQ=255)
                                                                                                        cgcgATGGCGTGCAACTGGTCGACATCAGCGCCGGGCACTTCGGCGAAGTGCTGCTGGGCCGTGCCGAAGACAGCGCCCTGGTGGTGTTCGAGGCGCGCCGTTATTCCCGCCATGCCCTGATGCTGTGCCAGAAGGCGCGCGAG‑GCTg                                              <  6:214582/148‑1 (MQ=255)
                                                                                                              ggCGTGGAACTGGTGGACATCAGCGGGGGGCACTTCGGCGAAGTGCTGCTGGGCCGTGCCGAAGACAGCGCCCTGGTGGTGTTCGAGGCGCGCCGTTATTCCCGCCATGCCCTGATGGTGTGCCAGAATGCGCGCGAG‑GCTGGCATTg                                        <  4:33284/148‑1 (MQ=255)
                                                                                                               gCGTGCAACTGGTCGACATCAGCGCCGGGCACTTCGGCGAAGTGCTGCTGGGCCGTGCCGAAGACAGCGCCCTGGTGGTGTTCGAGGCGCGCCGTTATTCCCGCCATGCCCTGATGCTGTGCCAGAAGGCGCGCGAG‑GCTGGCATTGCg                                      <  2:253226/149‑1 (MQ=255)
                                                                                                                cGTGCAACTGGTCGACATCAGCGCCGGGCACTTCGGCGAAGTGCTGCTGGGCCGTGCCGAAGACAGCGCCCTGGTGGTGTTCGAGGCGCGCCGTTATTCCCGCCATGCCCTGATGCTGTGCCAGAAGGCGCGCGAG‑GCTGGCATTGCgg                                     >  3:162062/1‑149 (MQ=255)
                                                                                                                 gTGCAACTGGTCGACATCAGCGCCGGGCACTTCGGCGAAGTGCTGCTGGGCCGTGCCGAAGACAGCGCCCTGGTGGTGTTCGAGGCGCGCCGTTATTCCCGCCATGCCCTGATGCTGTGCCAGAAGGCGCGCGAG‑GCTGGCTTTGCGGt                                    >  3:271210/1‑149 (MQ=255)
                                                                                                                  tGCAACTGGTCGACATCAGCGCCGGGCACTTCGGCGAAGTGCTGCTGGGCCGTGCCGAAGACAGCGCCCTGGTGGTGTTCGAGGCGCGCCGTTATTCCCGCCATGCCCTGATGATGTGCCAGAAGGCGCGCGAG‑GCTGGCATTCCGGt                                    >  3:64854/1‑148 (MQ=255)
                                                                                                                    cAACTGGTCGACATCAGCGCCGGGCACTTCGGCGAAGTGCTGCTGGGCCGTGCCGAAGACAGCGCCCTGGT‑GTGTTCGAGGCGCGCCGTTATTCCCGCCATGCCCTGATGCTGTGCCAGAAGGCGCGCGAG‑GCTGGCATTGCGGTCAcc                                >  8:332576/1‑149 (MQ=255)
                                                                                                                     accTGGTGGACCTCAGCGCCGGGCACTTCGGGGAAGTGCTGCTGGGCCGTGCCGGAGACCGCGCGCTGGTGGTGTTGGAGGCGCGCCGTTATTCCCGCCATGGGGGGATGCTGGGCGGGAAGGCGCGCGAG‑GCTGGCCTTGCGGTCagc                                >  8:302951/3‑147 (MQ=255)
                                                                                                                       cTGGTCGACATCAGCGCCGGGCACTTCGGCGAAGTGCTGCTGGGCCGTGCCGAAGACAGCGCCCTGGTGGTGTTCGAGGCGCGCCGTTATTCCCGCCATGCCCTGATGCTGTGCCAGAAGGCGCGCGAG‑GCTGGCATTGCGGTCACCCt                              <  7:201018/149‑1 (MQ=255)
                                                                                                                        tGGTCGACATCAGCGCCGGGCACTTCGGCGAAGTGCTGCTGGGCCGTGCCGAAGACAGCGCCCTGGTGGTGTTCGAGGCGCGCCGTTATTCCCGCCATGCCCTGATGCTGTGCCAGAAGGCGCGCGAG‑GCTGGCATTGCGGTCACCCTg                             >  3:406978/1‑149 (MQ=255)
                                                                                                                        tGGTCGACATCAGCGCCGGGCACTTCGGCGAAGTGCTGCTGGGCCGTGCCGAAGACAGCGCCCTGGTGGTGTTCGAGGCGCGCCGTTATTCCCGCCATGCCCTGATGCTGTGCCAGAAGGCGCGCGAG‑GCTGGCATTGCGGTCACCCTg                             >  8:249495/1‑149 (MQ=255)
                                                                                                                          gTCGACATCAGCGCCGGGCACTTCGGCGAAGTGCTGCTGGGCCGTGCCGAAGACAGCGCCCTGGTGGTGTTCGAGGCGCGCCGTTATTCCCGCCATGCCCTGATGCTGTGCCAGAAGGCGCGCGAG‑GCTGGCATTGCGGTCACCCTGGt                           >  5:260500/1‑149 (MQ=255)
                                                                                                                             gACATCAGCGCCGGGCACTTCGGCGAAGTGCTGCTGGGCCGTGCCGAAGACAGCGCCCTGGTGGTGTTCGAGGCGCGCAGTTATTCCCGCCATGCCCTGATGCTGTGCCAGAAGGCGCGCGAG‑GCTGGCATTGCGGTCACCCTGGTGAc                        >  4:93049/1‑149 (MQ=255)
                                                                                                                             gACAGCAGCGCCGGGCACTTCGGCGAAGGGCTGCTGGGCCGGGCCGAAGACAGCGCCCTGGTGGTGTTCGCGGCGCGCCGGTATTCCCGCCATGCCCTGATGCTGTGCCAGAAGGCGCGCGAG‑GCTGGCATTGCGGGCACCCTGGTGAc                        <  7:162825/149‑1 (MQ=255)
                                                                                                                              aCATCAGCGCCGGGCACTTCGGCGAAGTGCTGCTGGGCCGTGCCGAAGACAGCGCCCTGGTGGTGTTCGAGGCGGGCCGTTATTCCCGCCATGCCCTGATGCTGTGCCAGAAGGCGCGCGAGCGCTGGCATTGCGGTCACCCTGGTGAc                        <  1:226408/149‑1 (MQ=255)
                                                                                                                              aCATCAGCGCCGGGCACTTCGGCGAAGTGCTGCTGGGCCGTGCCGAAGACAGCGCCCTGGTGGTGTTCGAGGCGCGCCGTTATTCCCGCCATGCCCTGATGCTGTGCCAGAAGGCGCGCGAG‑GCTGGCATTGCGGTCACCCTGGTGAcc                       <  2:370896/149‑1 (MQ=255)
                                                                                                                                aTCAGCGCCGGGCACTTCGGCGAAGTGCTGCTGGGCCGTGCCGAAGACAGCGCCCTGGTGGTGTTCGAGGCGCGCCGTTATTCCCGCCATGCCCTGATGCTGTGCCAGAAGGCGCGCGAG‑GCTGGCATTGCGGTCACCCTGGTGACCGa                     >  8:406771/1‑149 (MQ=255)
                                                                                                                                aTCAGCGCCGGGCACTTCGGCGAAGTGCTGCTGGGCCGTGCCGAAGACAGCGCCCTGGTGGTGTTCGAGGCGCGCCGTTATTCCCGCCATGCCCTGATGCTGTGCCAGAAGGCGCGCGAG‑GCTGGCATTGCGGTCACCCTGGTGACCGa                     >  3:13196/1‑149 (MQ=255)
                                                                                                                                  cAGCGCCGGGCACTTCGGCGAAGTGCTGCTGGGCCGTGCCGAAGACAGCGCCCTGGTGGTGTTCGAGGCGCGCCGTTATTCCCGCCATGCCCTGATGCTGTGCCAGAAGGCGCGCGAG‑GCTGGCATTGCGGTCACCCTg                             >  3:194525/1‑139 (MQ=255)
                                                                                                                                  cAGCGCCGGGCACTTCGGCGAAGTGCTGCTGGGCCGTGCCGAAGACAGCGCCCTGGTGGTGTTCGAGGCGCGCCGTTATTCCCGCCATGCCCTGATGCTGTGCCAGAAGGCGCGCGAG‑GCTGGCATTGCGGTCACCCTg                             <  4:194525/139‑1 (MQ=255)
                                                                                                                                    gcgcCGGGCACTTCGGCGAAGTGCTGCTGGGCCGTGCCGAAGACAGCGCCCTGGTGGTGTTCGAGGCGCGCCGTTATTCCCGCCATGCCCTGATGCTGTGCCAGAAGGCGCGCGAG‑GCTGGCATTGCGGTCACCCTGGTGACCGACAcc                 <  4:41645/149‑1 (MQ=255)
                                                                                                                                        cGGGCACTTCGGCGAAGTGCTGCTGGGCCGTGCCGAAGACAGCGCCCTGGTGGTGTTCGAGGCGCGCCGTTATTCCCGCCATGCCCTGATGCTGTGCCAGAAGGCGCGCGAG‑GCTGGCATTGCGGTCACCCTGGTGACCGACACCTTCt             <  4:13196/149‑1 (MQ=255)
                                                                                                                                         gggCACTTCGGCGAAGTGCTGCTGGGCCGTGCCGAAGACAGCGCCCTGGT‑GTGTTCGAGGCGCGCCGTTATTCCCGCCATGCCCTGATGCTGTGCCAGAAGGCGCGCGAG‑GCTGGCATTGCGGTCACCCTGGTGACCGACACCTTCtgt           <  7:332576/149‑1 (MQ=255)
                                                                                                                                         cagCACTTCGCCGAAGTGCTGCTGGGCCGTGCCGAAGACAGCGCCCTGGTGGTGTTCGAGGCGCGCCGTTATTCCCGCCATGCCCTGATGCTGTGCCAGAAGGCGCGCGAG‑GCTGGCATTGCGGTCACCCTGGTGACCGACACCTTCtg            >  7:89560/3‑149 (MQ=255)
                                                                                                                                           gCACTTCGCCGAAGTGCTGCTGGGCCGTGCCGAAGACAGCGCCCTGGTGGTGTTCGAGGCGCGCCGTTATTCCCGCCATGCCCTGATGCTGTGCCAGAAGGCGCGCGAG‑GCTGGCATTGCGGTCACCCTGGTGAcc                       <  4:25000/136‑1 (MQ=255)
                                                                                                                                           gCACTTCGCCGAAGTGCTGCTGGGCCGTGCCGAAGACAGCGCCCTGGTGGTGTTCGAGGCGCGCCGTTATTCCCGCCATGCCCTGATGCTGTGCCAGAAGGCGCGCGAG‑GCTGGCATTGCGGTCACCCTGGTGAcc                       >  3:25000/1‑136 (MQ=255)
                                                                                                                                              cTTCGGCGAAGTGCTGCTGGGCCGTGCCGAAGACAGCGCCCTGGTGGTGTTCGAGGCGCGCCGTTATTCCCGCCATGCCCTGATGCTGTGCCAGAAGGCGCGCGAG‑GCTGGCATTGCGGTCACCCTGGTGACCGACACCTTCTGTGACt       >  6:253776/1‑149 (MQ=255)
                                                                                                                                               ttCGGCGAAGTGCTGCTGGGCCGTGCCGAAGACAGCGCCCTGGTGGTGTTCGAGGCGCGCCGTTATTCCCGCCATGCCCTGATGCTGTGCCAGAAGGCGCGCGAG‑GCTGGCATTGCGGTCACCCTGGTGACCGACACCTTCTGTGACTg      >  6:433155/1‑149 (MQ=255)
                                                                                                                                                 cGGCGAAGTGCGTCTGGGCCGTGCCGAAGACAGCGGCCTGGTGGGGTTGGAGGCGCGCCGTTATTCCCGCCATGCCCTGATGCTGTGCCAGAAGGCGCGCGAG‑GCTGGCATTGCGGTCACCCTGGTGACCGACACCTTCTGTGACTgg     <  7:340075/148‑1 (MQ=255)
                                                                                                                                                   gCGAAGTGCTGCTGGGCCGTGCCGAAGACAGCGCCCTGGTGGTGTTCGAGGCGCGCCGTTATTCCCGCCATGCCCTGATGCTGTGCCAGAAGGCGCGCGAG‑GCTGGCATTGCGGTCACCCTGGTGACCGACACCTTCTGTGACTGGGCg  >  2:1643/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
ACAACGTGGCGCAACGCCTGGCGAGCAAGCCGCGGGTGTTCATCGCCGGCTTCCAGACCGAGCGCGGCATCGCCATGTGCATGAGCCACCTGCTGCAATACCTGCGCGATGGCGTGCAACTGGTCGACATCAGCGCCGGGCACTTCGGCGAAGTGCTGCTGGGCCGTGCCGAAGACAGCGCCCTGGTGGTGTTCGAGGCGCGCCGTTATTCCCGCCATGCCCTGATGCTGTGCCAGAAGGCGCGCGAG‑GCTGGCATTGCGGTCACCCTGGTGACCGACACCTTCTGTGACTGGGCG  >  NC_002947/4082474‑4082769

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: ‑