Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,038,463 T→G 5.1% L73R (CTG→CGG)  PP_3559 → glycine betaine ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,038,4630TG 5.1% 230.5 / 4.4 78L73R (CTG→CGG) PP_3559glycine betaine ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base T (40/34);  new base G (2/2);  total (42/36)
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.91e-01

GTTGGTATTGAACTACGGCGACCACCTGCGCAGTTTGTCTGACCAGTTGCTGCAACTGCTGGTCGGCCTGGAGAACCTGTTGCGCCTGTTGCCTTGGTGGCTGTTGTTGTTGCTGGTGGGTCTGCTGGCCTGGCACGCCAGCCGCAGCCTGCTGCGCAGTGCCGTGCTGGTCGCGCTGCTGGCCCTGATCGGTATGCTCGGGCTCTGGGACAAGTTGCTGCAGACCCTGGCCCTGGTGCTGGTCAGTACCGGTTTGTGCGTACTGGTGGGTGTACCACTGGGCATCCTGCTG  >  NC_002947/4038314‑4038605
                                                                                                                                                     |                                                                                                                                              
gttggtATTGAACTACGGCGACCACCTGCGCAGTTTGTCTGACCAGTTGCTGCAACTGCTGGTCGGCCTGGAGAACCTGTTGCGCCTGTTGCCTTGGTGGCTGTTGTTGTTGCTGGTGGGTCTGCTGGCCTGGCACGCCAGCCGCAGCct                                                                                                                                                >  1:592474/1‑150 (MQ=255)
        tGAACTACGGCGACCACCTGCGCAGTTTGTCTGACCAGTTGCTGCAACTGCTGGTCGGCCTGGAGAACCTGTTGCGCCTGTTGCCTTGGTGGCTGTTGTTGTTGCTGGTGGGTCTGCTGGCCTGGCACGCCAGCCGCAGCCTGCTGCGCa                                                                                                                                        <  2:592474/150‑1 (MQ=255)
         gAACTACGGCGACCACCTGCGCAGTTTGTCTGACCAGTTGCTGCAACTGCTGGTCGGCCTGGAGAACCTGTTGCGCCTGTTGCCTTGGTGGCTGTTGTTGTTGCTGGTGGGTCTGCTGGCCTGGCACGCCAGCCGCAGCCTGCTGCGCAg                                                                                                                                       >  1:1176291/1‑150 (MQ=255)
          aaCTACGGCGACCACCTGCGCAGTTTGTCTGACCAGTTGCTGCAACTGCTGGTCGGCCTGGAGAACCTGTTGCGCCTGTTGCCTTGGTGGCTGTTGTTGTTGCTGGTGGGTCTGCTGGCCTGGCACGCCAGCCGCAGCCTGCTGCGCact                                                                                                                                      >  1:378055/1‑148 (MQ=255)
             tACGGCGACCACCTGCGCAGTTTGTCTGACCAGTTGCTGCAACTGCTGGTCGGCCTGGAGAACCTGTTGCGCCTGTTGCCTTGGTGGCTGTTGTTGTTGCTGGTGGGTCTGCTGGCCTGGCACGCCAGCCGCAGCCTGCTGCGCAGTGcc                                                                                                                                   >  1:1057/1‑150 (MQ=255)
                      caccTGCGCAGTTTGTCTGACCCGTTCAGGCAACTGCTGGTCGGCCTGGAGAACCTGTTGCGCCTGTTGCCTTGGTGGCTGTTGTTGTTGCTGGTGGGTCTGCTGGCCTGGCACGCCAGCCGCAGCCTGCTGCGCAGTGCCGTGCTGGTc                                                                                                                          <  1:965577/150‑1 (MQ=255)
                      caccTGCGCAGTTTGTCTGACCAGTTGCTGCAACTGCTGGTCGGCCTGGAGAACCTGTTGCGCCTGTTGCCTTGGTGGCTGTTGTTGTTGCTGGTGGGTCTGCTGGCCTGGCACGCCAGCCGCAGCAGGCTgcggc                                                                                                                                        >  1:450748/1‑134 (MQ=255)
                      caccTGCGCAGTTTGTCTGACCAGTTGCTGCAACTGCTGGTCGGCCTGGAGAACCTGTTGCGCCTGTTGCCTTGGTGGCTGTTGTTGTTGCTGGTGGGTCTGCTGGCCTGGCACGCCAGCCGCAGCAGGCTgcggc                                                                                                                                        <  2:450748/136‑3 (MQ=255)
                          tGCGCAGTTTGTCTGACCAGTTGCTGCAACTGCTGGTCGGCCTGGAGAACCTGTTGCGCCTGTTGCCTTGGTGGCTGTTGTTGTTGCTGGTGGGTCTGCTGGCCTGGCACGCCAGCCGCAGCCTGCTGCGCAGTGCCGTGCTGGTCgcgc                                                                                                                      <  2:1176291/150‑1 (MQ=255)
                                  ttGTCTGACCAGTTGCTGCAACTGCTGGTCGGCCTGGAGAACCTGTTGCGCCTGTTGCCTTGGTGGCTGTTGTTGTTGCTGGTGGGTCTGCTGGCCTGGCACGCCAGCCGCAGCCTGCTGCGCAGTGCCGTGCTGGTCGCGCTGCTGGcc                                                                                                              <  2:1057/150‑1 (MQ=255)
                                       tGACCAGTTGCTGCAACTGCTGGTCGGCCTGGAGAACCTGTTGCGCCTGTTGCCTTGGTGGCTGTTGTTGTTGCTGGTGGGTCTGCTGGCCTGGCACGCCAGCCGCAGCCTGCTGCGCAGTGCCGTGCTGGTCGCGCTGCTGGCCCTGAt                                                                                                         >  1:734422/1‑150 (MQ=255)
                                           cAGTTGCTGCAACTGCTGGTCGGCCTGGAGAACCTGTTGCGCCTGTTGCCTTGGTGGCTGTTGTTGTTGCTGGTGGGTCTGCTGGCCTGGCACGCCAGCCGCAGCCTGCTGCGCAGTGCCGTGCTGGTCGCGCTGCTGGCCCTGATCGGt                                                                                                     >  1:73198/1‑150 (MQ=255)
                                                  tgcAACTGCTAGTCGGCCTGCAGAACCTGTTGCGCCTGTTGGCTTGGTGGCTGGTCTTGTTGCTGGTGGGGCTGCTGGCCTGGCACGCCAGCCGCAGCCTGCTGCGCAGTGCCGTGCTGGTCGCGCTGCTGGCCCTGATCGGTATGCTCg                                                                                              >  2:1277626/1‑150 (MQ=255)
                                                   gcAACTGCTGGTCGGCCTGGAGAACCTGTTGCGCCTGTTGCCTTGGTGGCTGTTGTTGTTGCTGGTGGGTCTGCTGGCCTGGCACGCCAGCCGCAGCCTGCTGCGCAGTGCCGTGCTGGTc                                                                                                                          >  2:544706/1‑121 (MQ=255)
                                                   gcAACTGCTGGTCGGCCTGGAGAACCTGTTGCGCCTGTTGCCTTGGTGGCTGTTGTTGTTGCTGGTGGGTCTGCTGGCCTGGCACGCCAGCCGCAGCCTGCTGCGCAGTGCCGTGCTGGTc                                                                                                                          <  1:544706/121‑1 (MQ=255)
                                                    cAACTGCTGGTCGGCCTGGAGAACCTGTTGCGCCTGTTGCCTTGGTGGCTGTTGTTGTTGCTGGTGGGTCTGCTGGCCTGGCACGCCAGCCGCAGCCTGCTGCGCAGTGCCGTGCTGGTCGCGCTGCTGGCCCTGATCGGTATGCTCggg                                                                                            >  1:109838/1‑150 (MQ=255)
                                                      aCTGCTGGTCGGCCTGGAGAACCTGTTGCGCCTGTTGCCTTGGTGGCTGTTGTTGTTGCTGGTGGGTCTGCTGGCCTGGCACGCCAGCCGCAGCCTGCTGCGCAGTGCCGTGCTGGTCGCGCTGCTGGccc                                                                                                             <  1:177717/131‑1 (MQ=255)
                                                      aCTGCTGGTCGGCCTGGAGAACCTGTTGCGCCTGTTGCCTTGGTGGCTGTTGTTGTTGCTGGTGGGTCTGCTGGCCTGGCACGCCAGCCGCAGCCTGCTGCGCAGTGCCGTGCTGGTCGCGCTGCTGGccc                                                                                                             >  2:177717/1‑131 (MQ=255)
                                                       ctgctgGTCGGCCTGGAGAACCTGTTGCGCCTGTTGCCTTGGTGGCTGTTGTTGTTGCTGGTGGGTCTGCTGGCCTGGCACGCCAGCCGCAGCCTGCTGCGCAGTGCCGTGCTGGTCGCGCTGCTGGCCCTGATCGGTATGCTCGGGctc                                                                                         >  1:397707/1‑150 (MQ=255)
                                                          ctgGTCGGCCTGGAGAACCTGTTGCGCCTGTTGCCTTGGTGGCTGTTGTTGTTGCTGGTGGGTCTGCTGGCCTGGCACGCCAGCCGCAGCCTGCTGCGCAGTGCCGTGCTGGTCGCGCTGCTGGCCCTGATCGGTATGCTCGGGCTCTgg                                                                                      <  2:15423/150‑1 (MQ=255)
                                                               cGGCCTGGAGAACCTGTTGCGCCTGTTGCCTTGGTGTCTCTTTTTGTTGCTGGTGGGTCTGCTGGCCTGGCACGCCAGCCGCAGCCTGCTGCGCAGTGCCGTGCTGGTCGCGCTGCTGGCCCTGATCGGTATGCTCGGGCTCTGGGACaa                                                                                 <  2:943538/150‑1 (MQ=255)
                                                                    tGGAGAACCTGTTGCGCCTGTTGCCTTGGTGGCTGTTGTTGTTGCTGGTGGGTCTGCTGGCCTGGCACGCCAGCCGCAGCCTGCTGCGCAGTGCCGTGCTGGTCGCGCTGCTGGCCCTGATCGGTATGCTCGGGCTCTGGGACAAGTtgc                                                                            >  2:975290/1‑150 (MQ=255)
                                                                         aaCCTGTTGCGCCTGTTGCCTTGGTGGCTGTTGTTGTTGCTGGTGGGTCTGCTGGCCTGGCACGCCAGCCGCAGCCTGCTGCGCAGTGCCGTGCTGGTCGCGCTGCTGGCCCTGATCGGTATGCTCGGGCTCTGGGACAAGTTGCTGCAg                                                                       >  1:373770/1‑150 (MQ=255)
                                                                          aCCTGTTGCGCCTGTTGCCTTGGTGGCTGTTGTTGTTGCTGGTGGGTCTGCTGGCCTGGCACGCCAGCCGCAGCCTGCTGCGCAGTGCCGTGCTGGTCGCGCTGCTGGCCCTGATCGGTATGCTCGGGCTCTGGGACAAGTTGCTGCAGa                                                                      >  2:893348/1‑150 (MQ=255)
                                                                            cTGTTGCGCCTGTTGCCTTGGTGGCTGTTGTTGTTGCTGGTGGGTCTGCTGGCCTGGCACGCCAGCCGCAGCAGGCTgcggct                                                                                                                                       <  2:642394/83‑4 (MQ=255)
                                                                            cTGTTGCGCCTGTTGCCTTGGTGGCTGTTGTTGTTGCTGGTGGGTCTGCTGGCCTGGCACGCCAGCCGCAGCAGGCTgcggct                                                                                                                                       >  1:642394/1‑80 (MQ=255)
                                                                             tGTTGCGCCTGTTGCCTTGGTGGCTGGTGTTGTTGCTGGTGGGTCTGCTGGCCTGGCACGCCAGCCGCAGCCTGCTGCGCAGTGCCGTGCTGGTCGCGCTGCTGGCCCTGATCGGTATGCTCGGGCTCTGGGACAAGTTGCTGCAGAccc                                                                   >  1:1137277/1‑150 (MQ=255)
                                                                              gTTGCGCCTGTTGCCTTGGTGGCTGTTGTTGTTGCTGGTGGGTCTGCTGGCCTGGCACGCCAGCCGCAGCCTGCTGCGCAGTGCCGTGCTGGTCGCGCTGCTGGCCCTGATCGGTATGCTCGGGCTCTgg                                                                                      >  1:1200265/1‑130 (MQ=255)
                                                                              gTTGCGCCTGTTGCCTTGGTGGCTGTTGTTGTTGCTGGTGGGTCTGCTGGCCTGGCACGCCAGCCGCAGCCTGCTGCGCAGTGCCGTGCTGGTCGCGCTGCTGGCCCTGATCGGTATGCTCGGGCTCTgg                                                                                      <  2:1200265/130‑1 (MQ=255)
                                                                               ttGCGCCTGTTGCCTTGGTGGCTGTTGTTGTTGCTGGTGGGTCTGCTGGCCTGGCACGCCAGCCGCAGCCTGCTGCGCAGTGCCGTGc                                                                                                                               <  1:442628/88‑1 (MQ=255)
                                                                               ttGCGCCTGTTGCCTTGGTGGCTGTTGTTGTTGCTGGTGGGTCTGCTGGCCTGGCACGCCAGCCGCAGCCTGCTGCGCAGTGCCGTGc                                                                                                                               >  2:442628/1‑88 (MQ=255)
                                                                               ttGCGCCTGTTGCCTTGGTGGCTGTTGTTGTTGCTGGTGGGTCTGCTGGCCTGGCACGCCAGCCGCAGCCTGCTGCGCAGTGCCGTGCTGGTCGCGCTGCTGGCCCTGATCGGTATGCTCGGGCTCTGGGACaa                                                                                 <  1:1231271/134‑1 (MQ=255)
                                                                               ttGCGCCTGTTGCCTTGGTGGCTGTTGTTGTTGCTGGTGGGTCTGCTGGCCTGGCACGCCAGCCGCAGCCTGCTGCGCAGTGCCGTGCTGGTCGCGCTGCTGGCCCTGATCGGTATGCTCGGGCTCTGGGACaa                                                                                 >  2:1231271/1‑134 (MQ=255)
                                                                                  cgcCTGTTGCCTTGGTGGCTGTTGTTGTTGCTGGTGGGTCTGCTGGCCTGGCACGCCAGCCGCAGCCTGCTGCGCAGTGCCGTGCTGGTAGCGCTGCTGGCCCTGATCGGTATGCTCGGGCTCTGGGACAAGTTGCTGCAGAccctggcc                                                              >  1:838064/1‑150 (MQ=255)
                                                                                     cTGTTGCCTTGGTGGCTGTTGTTGTTGCTGGTGGGTCTGCTGGCCTGGCACGCCAGCCGCAGCCTGCTGCGCAGTGCCGTGCTGGTCGCGCTGCTGGCCCTGATCGGTATGCTCGGGCTCTGGGACAAGTTGCTGCAGAccctggccctg                                                           <  2:397707/150‑1 (MQ=255)
                                                                                       gTTGCCTTGGTGGCTGTTGTTGTTGCTGGTGGGTCTGCTGGCCTGGCACGCCAGCCGCAGCCTGCTGCGCAGTGCCGTGCTGGTCGCGCTGCTGGCCCTGATCGGTATGCTCGGGCTCTGGGACAAGTTGCTGCAGACCCTGGCCCTGGt                                                         >  2:801487/1‑150 (MQ=255)
                                                                                           ccTTGGTGGCTGTTGTTGTTGCTGGTGGGTCTGCTGGCCTGGCACGCCAGCCGCAGCCTGCTGCGCAGTGCCGTGCTGGTCGCGCTGCTGGCCCTGATCGGTATGCTCGGGCTCTGGGACAAGTTGCTGCAGACCCTGGCCCTGGTGCTg                                                     >  2:252119/1‑150 (MQ=255)
                                                                                             ttGGTGGCTGTTGTTGTTGCTGGTGGGTCTGCTGGCCTGGCACGCCAGCCGCAGCCTGCTGCGCAGTGCCGTGCTGGTCGCGCTGCTGGCCCTGATCGGt                                                                                                     <  1:788984/100‑1 (MQ=255)
                                                                                             ttGGTGGCTGTTGTTGTTGCTGGTGGGTCTGCTGGCCTGGCACGCCAGCCGCAGCCTGCTGCGCAGTGCCGTGCTGGTCGCGCTGCTGGCCCTGATCGGt                                                                                                     >  2:788984/1‑100 (MQ=255)
                                                                                               ggtggCTGTTGTTGTTGCTGGTGGGTCTGCTGGCCTGGCACGCCAGCCGCAGCCTGCTGCGCAGTGCCGTGCTGGTCGCGCTGCTGGCCCTGATCGGTATGCTCGGGCTCTGGGACAAGTTGCTGCAGACCCTGGCCCTGGTGCTGGTCa                                                 <  1:318450/150‑1 (MQ=255)
                                                                                                    cTGTTGTTGTTGCTGGTGGGTCTGCTGGCCTGGCACGCCAGCCGCAGCCTGCTGCGCAGTGCCGTGCTGGTCGCGCTGCTGGCCCTGATCGGTATGCTCGGGCTCTGGGACAAGTTGCTGCAGACCCTGGCCCTGGTGCTGGTCAGTAcc                                            <  1:975290/150‑1 (MQ=255)
                                                                                                     tgttgttgtTGCTGGTGGGTCTGCTGGCCTGGCACGCCAGCCGCAGCCTGCTGCGCAGTGCCGTGCTGGTCGCGCTGCTGGCCCTGATCGGTATGCTCGGGCTCTGGGACAAGTTGCTGCAGACCCTGGCCCTGGTGCTGGTCAGTACCg                                           >  1:876634/1‑150 (MQ=255)
                                                                                                     tgttgttgtTGCTGGTGGGTCTGCTGGCCTGGCACGCCAGCCGCAGCCTGCTGCGCAGTGCCGTGCTGGTCGCGCTGCTGGCCCTGATCGGTATGCTCGGGCTCTGGGACAAGTTGCTGCAGACCCTGGCCCTGGTGCTGGTCAGTACCg                                           >  1:331016/1‑150 (MQ=255)
                                                                                                     tgttgttgtTGCTGGTGGGTCTGCTGGCCTGGCACGCCAGCCGCAGCCTGCTGCGCAGTGCCGTGCTGGTAGCGCTGCTGGCCCTGATCGGTATGCTCGGGCTCTGGGACAAGTTGCTGCAGACCCTGGCCCTGGTGCTGGTCAGTACCg                                           <  2:838064/150‑1 (MQ=255)
                                                                                                       ttgttgttgCTGGTGGGTCTGCTGGCCTGGCACGCCAGCCGCAGCCTGCTGCGCAGTGCCGTGCTGGTCGCGCTGCTGGCCCTGATCGGTATGCTCGGGCTCTGGGACAAGTTGCTGCAGACCCTGGCCCTGGTGCTGGTCAGTACCGGt                                         <  2:109838/150‑1 (MQ=255)
                                                                                                         gttgttgCTGGTGGGTCTGCTGGCCTGGCACGCCAGCCGCAGCCTGCTGCGCAGTGCCGTGCTGGTCGCGCTGCTGGCCCTGATCGGTATGCTCGGGCTCTGGGACAAGTTGCTGCAGACCCTGGCCCTGGTGCTGGTCAGTACCGGtt                                        >  2:286488/1‑149 (MQ=255)
                                                                                                         gttgttgCTGGTGGGTCTGCTGGCCTGGCACGCCAGCCGCAGCCTGCTGCGCAGTGCCGTGCTGGTCGCGCTGCTGGCCCTGATCGGTATGCTCGGGCTCTGGGACAAGTTGCTGCAGACCCTGGCCCTGGTGCTGGTCAGTACCGGtt                                        <  1:286488/149‑1 (MQ=255)
                                                                                                            gttgCTGGTGGGTCTGCTGGCCTGGCACGCCAGCCGCAGCCTGCTGCGCAGTGCCGTGCTGGTCGCGCTGCTGGCCCTGATCGGTATGCTCGGGCTCTGGGACAAGTTGCTGCAGACCCTGGCCCTGGTGCTGGTCAGTACCGGTTtgtg                                    >  1:868442/1‑150 (MQ=255)
                                                                                                             ttgCTGGTGGGTCTGCTGGCCTGGCACGCCAGCCGCAGCCTGCTGCGCAGTGCCGTGCTGGTCGCGCTGCTGGCCCTGATCGGTATGCTCGGGCTCTGGGACAAGTtgctg                                                                          >  2:889343/1‑111 (MQ=255)
                                                                                                             ttgCTGGTGGGTCTGCTGGCCTGGCACGCCAGCCGCAGCCTGCTGCGCAGTGCCGTGCTGGTCGCGCTGCTGGCCCTGATCGGTATGCTCGGGCTCTGGGACAAGTtgctg                                                                          <  1:889343/111‑1 (MQ=255)
                                                                                                                cTGGTGGGTCTGCTGGCCTGGCACGCCAGCCGCAGCCTGCTGCGCAGTGCCGTGCTGGTCGCGCTGCTGGCCCTGATCGGTATGCTCGGGCTCTGGGACAAGTTGCTCCAGACCCTGGCCCTGGTGCTGGTCAGTCCCGGCTTGCGCGTa                                <  2:373770/150‑1 (MQ=255)
                                                                                                                       gTCTGCTGGCCTGGCACGCCAGCCGCAGCCTGCTGCGCAGTGCCGTGCTGGTCGCGCTGCTGGCCCTGATCGGTATGCTCGGGCTCTGGGACAAGTTGCTGCAGACCCTGGCCCTGGTGCTGGTCa                                                 >  2:1261227/1‑126 (MQ=255)
                                                                                                                       gTCTGCTGGCCTGGCACGCCAGCCGCAGCCTGCTGCGCAGTGCCGTGCTGGTCGCGCTGCTGGCCCTGATCGGTATGCTCGGGCTCTGGGACAAGTTGCTGCAGACCCTGGCCCTGGTGCTGGTCa                                                 <  1:1261227/126‑1 (MQ=255)
                                                                                                                       gTCTGCTGGCCTGGCACGCCAGCCGCAGCCTGCTGCGCAGTGCCGTGCTGGTCGCGCTGCTGGCCCTGATCGGTATGCTCGGGCTCTGGGACAAGTTGCTGCAGACCCTGGCCCTGGTGCTGGTCAGTACCGGTTTGTGCGTACtggtgg                         <  2:734422/150‑1 (MQ=255)
                                                                                                                        tCTGCTGGCCTGGCACGCCAGCCGCAGCCTGCTGCGCAGTGCCGTGCTGGTCGCGCTGCTGGCCCTGATCGGTATGCTCGGGCTCTGGGACAAGTTGCTGCAGACCCTGGCCCTGGTGCTGGTCAGTACCGGTTTGTGCGTACTggtggg                        >  2:266952/1‑150 (MQ=255)
                                                                                                                          ggctgGCGTGCCAGGCCAGCCGCAGCCTGCTGCGCAGTGCCGTGCTGGTCGCGCTGCTGGCCCTGATCGGTATGCTCGGGCTCTGGGACAAGTTGCTGCAGACCCTGGCCCTGGTGCTGGTCAGTACCGGTTTGTGCGTACTggtggg                        <  2:445728/147‑1 (MQ=255)
                                                                                                                          ggctgGCGTGCCAGGCCAGCCGCAGCCTGCTGCGCAGTGCCGTGCTGGTCGCGCTGCTGGCCCTGATCGGTATGCTCGGGCTCTGGGACAAGTTGCTGCAGACCCTGGCCCTGGTGCTGGTCAGTACCGGTTTGTGCGTACTggtggg                        >  1:445728/2‑148 (MQ=255)
                                                                                                                               gcctggcACGCCAGCCGCAGCCTGCTGCGCAGTGCCGTGCTGGTCGCGCTGCTGGCCCTGATCGGTATGCTCGGGCTCTGGGACAAGTTGCTGCAGACCCTGGCCCTGGTGCTGGTCAGTACCGGTTTGTGCGTACTGGTGGGTGTACCa                 <  1:1269357/150‑1 (MQ=255)
                                                                                                                                 ctggcACGCCAGCCGCAGCCTGCTGCGCAGTGCCGTGCTGGTCGCGCTGCTGGCCCTGATCGGTATGCTCGGGCTCTGGGACAAGTTGCTGCAGACCCTGGCCCTGGTGCTGGTCAGTACCg                                           >  2:590884/1‑122 (MQ=255)
                                                                                                                                 ctggcACGCCAGCCGCAGCCTGCTGCGCAGTGCCGTGCTGGTCGCGCTGCTGGCCCTGATCGGTATGCTCGGGCTCTGGGACAAGTTGCTGCAGACCCTGGCCCTGGTGCTGGTCAGTACCg                                           <  1:590884/122‑1 (MQ=255)
                                                                                                                                 ctggcACGCCAGCCGCAGCCTGCTGCGCAGTGCCGTGCTGGTCGCGCTGCTGGCCCTGATCGGTATGCTCGGGCTCTGGGACAAGTTGCTGCAGACCCTGGCCCTGGTGCTGGTCAGTACCGGTTTGTGCGTACTGGTGGGTGTACCACt               >  1:502372/1‑150 (MQ=255)
                                                                                                                                 ctggcACGCCAGCCGCAGCCTGCTGCGCAGTGCCGTGCTGGTCGCGCTGCTGGCCCTGATCGGTATGCTCGGGCTCTGGGACAAGTTGCTGCAGACCCTGGCCCTGGTGCTGGTCAGTACCGGTTTGTGCGTACTGGTGGGTGTACCACt               <  2:1137277/150‑1 (MQ=255)
                                                                                                                                  tggcACGCCAGCCGCAGCCTGCTGCGCAGTGCCGTGCTGGTCGCGCTGCTGGCCCTGATCGGTATGCTCGGGCTCTGGGACAAGTTGCTGCAGACCCTGGCCCTGGTGCTGGTCAGTACCGGTTTGTGCGTACtggt                           >  2:195825/1‑137 (MQ=255)
                                                                                                                                  tggcACGCCAGCCGCAGCCTGCTGCGCAGTGCCGTGCTGGTCGCGCTGCTGGCCCTGATCGGTATGCTCGGGCTCTGGGACAAGTTGCTGCAGACCCTGGCCCTGGTGCTGGTCAGTACCGGTTTGTGCGTACtggt                           <  1:195825/137‑1 (MQ=255)
                                                                                                                                   ggcACGCCAGCCGCAGCCTGCTGCGCAGTGCCGTGCTGGTCGCGCTGCTGGCCCTGATCGGTATGCTCGGGCTCTGGGACAAGTTGCTGCAGACCCTGGCCCTGGTGCTGGTCAGTACCGGTTTGTGCGTACTGGTGGGTGTACCACTgg             <  2:868442/150‑1 (MQ=255)
                                                                                                                                      aCGCCAGCCGCAGCCTGCTGCGCAGTGCCGTGCTGGTCGCGCTGCTGGCCCTGATCGGTATGCTCGGGCTCTGGGACAAGTTGCTGCAGACCCTGGCCCTGGTGCTGGTCAGTACCGGTTtgt                                     <  2:9027/123‑1 (MQ=255)
                                                                                                                                      aCGCCAGCCGCAGCCTGCTGCGCAGTGCCGTGCTGGTCGCGCTGCTGGCCCTGATCGGTATGCTCGGGCTCTGGGACAAGTTGCTGCAGACCCTGGCCCTGGTGCTGGTCAGTACCGGTTtgt                                     >  1:9027/1‑123 (MQ=255)
                                                                                                                                      aCGCCAGCCGCAGCCTGCTGCGCAGTGCCGTGCTGGTCGCGCTGCTGGCCCTGATCGGTATGCTCGGGCTCTGGGACAAGTTGCTGCAGACCCTGGCCCTGGTGCTGGTCAGTACCGGTTTGTGCGTACTGGTGGGTGTACCACTgg             >  1:901337/1‑147 (MQ=255)
                                                                                                                                      aCGCCAGCCGCAGCCTGCTGCGCAGTGCCGTGCTGGTCGCGCTGCTGGCCCTGATCGGTATGCTCGGGCTCTGGGACAAGTTGCTGCAGACCCTGGCCCTGGTGCTGGTCAGTACCGGTTTGTGCGTACTGGTGGGTGTACCACTgg             <  2:901337/147‑1 (MQ=255)
                                                                                                                                        gCCAGCCGCAGCCTGCTGCGCAGTGCCGTGCTGGTCGCGCTGCTGGCCCTGATCGGTATGCTCGGGCTCTGGGACAAGTTGCTGCAGACCCTGGCCCTGGTGCTGGTCAGTACCGGTTTGTGCGTACTGGTGGGTGTACCACTGGGCATc        <  1:26966/150‑1 (MQ=255)
                                                                                                                                        gCCAGCCGCAGCCTGCTGCGCAGTGCCGTGCTGGTCGCGCTGCTGGCCCTGATCGGTATGCTCGGGCTCTGGGACAAGTTGCTGCAGACCCTGGCCCTGGTGCTGGTCAGTACCGGTTTGTGCGTACTGGTGGGTGTACCACTGGGCATc        >  2:26966/1‑150 (MQ=255)
                                                                                                                                          cAGCCGCAGCCTGCTGCGCAGTGCCGTGCTGGTCGCGCTGCTGGCCCTGATCGGTATGCTCGGGCTCTGGGACAAGTTGCTGCAGACCCTGGCCCTGGTGCTGGTCAGTACCGGTTTGTGCGTACTGGTGGGTGTa                    <  1:941146/136‑1 (MQ=255)
                                                                                                                                          cAGCCGCAGCCTGCTGCGCAGTGCCGTGCTGGTCGCGCTGCTGGCCCTGATCGGTATGCTCGGGCTCTGGGACAAGTTGCTGCAGACCCTGGCCCTGGTGCTGGTCAGTACCGGTTTGTGCGTACTGGTGGGTGTa                    >  2:941146/1‑136 (MQ=255)
                                                                                                                                            gCCGCAGCCTGCTGCGCAGTGCCGTGCTGGTCGCGCTGCTGGCCCTGATCGGTATGCTCGGGCTCTGGGACAAGTTGCTGCAGACCCTGGCCCTGGTGCTGGTCAGTACCGGTTTGTGCGTACTGGTGGGTGTACCACTGGGCATCctgc    >  1:1163058/1‑150 (MQ=255)
                                                                                                                                            gCCGCAGCCTGCTGCGCAGTGCCGTGCTGGTCGCGCTGCTGGCCCTGATCGGTATGCTCGGGCTCTGGGACAAGTTGCTGCAGACCCTGGCCCTGGTGCTGGTCAGTACCGGTTTGCGCGTACTGGTGGGTGTACCACTGGGCATCctgc    >  1:1163166/1‑150 (MQ=255)
                                                                                                                                             ccGCAGCCTGCTGCGCAGTGCCGTGCTGGTCGCGCTGCTGGCCCTGATCGGTATGCTCGGGCTCTGGGACAAGTTGCTGCAGACCCTGGCCCTGGTGCTGGTCAGTACCGGTTTGTGCGTACTGGTGGGTGTACCACTGGGCATCctgct   <  2:268995/150‑1 (MQ=255)
                                                                                                                                              cGCAGCCTGCTGCGCAGTGCCGTGCTGGTCGCGCTGCTGGCCCTGATCGGTATGCTCGGGCTCTGGGACAAGTTGCTGCAGACCCTGGCCCTGGTGCTGGTCAGTACCGGTTTGTGCGTACTGGTGGGTGTACCACTGGGCATCctgctg  >  1:273291/1‑150 (MQ=255)
                                                                                                                                              cGCAGCCTGCTGCGCAGTGCCGTGCTGGTCGCGCTGCTGGCCCTGATCGGTATGCTCGGGCTCTGGGACAAGTTGCTGCAGACCCTGGCCCTGGTGCTGGTCAGTACCGGTTTGTGCGTACTGGTGGGTGTACCACTGGGCATCctgctg  <  1:252119/150‑1 (MQ=255)
                                                                                                                                              cGCAGCCTGCTGCGCAGTGCCGTGCTGGTCGCGCTGCTGGCCCTGATCGGTATGCTCGGGCTCTGGGACAAGTTGCTGCAGACCCTGGCCCTGGTGCTGGTCAGTACCGGTCTGTGCGTACTGGTGGGTGTACCACTGGGCATTctgctg  >  1:654529/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                              
GTTGGTATTGAACTACGGCGACCACCTGCGCAGTTTGTCTGACCAGTTGCTGCAACTGCTGGTCGGCCTGGAGAACCTGTTGCGCCTGTTGCCTTGGTGGCTGTTGTTGTTGCTGGTGGGTCTGCTGGCCTGGCACGCCAGCCGCAGCCTGCTGCGCAGTGCCGTGCTGGTCGCGCTGCTGGCCCTGATCGGTATGCTCGGGCTCTGGGACAAGTTGCTGCAGACCCTGGCCCTGGTGCTGGTCAGTACCGGTTTGTGCGTACTGGTGGGTGTACCACTGGGCATCCTGCTG  >  NC_002947/4038314‑4038605

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