Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,807,091 T→C 16.2% intergenic (+17/‑41) proW → / → proX glycine betaine/proline ABC transporter permease/glycine betaine/proline ABC transporter periplasmic binding protein; cold shock protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,807,0910TC16.2% 210.2 / 33.0 99intergenic (+17/‑41)proW/proXglycine betaine/proline ABC transporter permease/glycine betaine/proline ABC transporter periplasmic binding protein; cold shock protein
Reads supporting (aligned to +/- strand):  ref base T (46/36);  new base C (5/11);  total (52/47)
Fisher's exact test for biased strand distribution p-value = 1.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.90e-01

GGCGCGACTCACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACACAAACT  >  NC_000913/2806992‑2807188
                                                                                                   |                                                                                                 
ggCGCGACTCACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCACCCATTCATTAAGTAACTCTGCACTTGCCCGGt                                                                                                   <  2:992787/100‑1 (MQ=255)
 gcgcGACTCACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTg                                                                                                  <  2:1076326/100‑1 (MQ=255)
  cgcgACTCACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGa                                                                                                 <  2:830600/100‑1 (MQ=255)
  cgcgACTCACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGa                                                                                                 <  1:283712/100‑1 (MQ=255)
  cgcgACTCACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGa                                                                                                 >  1:1808288/1‑100 (MQ=255)
     gACTCACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGc                                                                                              <  1:2513645/100‑1 (MQ=255)
        tCACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCgg                                                                                           <  2:1950071/100‑1 (MQ=255)
         cACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCggg                                                                                          >  1:91521/1‑100 (MQ=255)
         cACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCggg                                                                                          >  1:1010491/1‑100 (MQ=255)
         cACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCggg                                                                                          >  1:1847860/1‑100 (MQ=255)
          aCGCAGTCGCTGGAACCGTTGCTGGGACACCACTGGCCCTGTGGGTCCGCTGACCCGCCCATTCATTAAGTAACGCTGCACTTGCCCGGCGTCACCGGGc                                                                                         <  1:628349/100‑1 (MQ=255)
          aCGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGc                                                                                         <  1:359846/100‑1 (MQ=255)
          aCGCAGTCGCGGCAACCGTAGCTGGTACACCACCGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGc                                                                                         >  1:492641/1‑100 (MQ=255)
           cGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCa                                                                                        >  2:532221/1‑100 (MQ=255)
           cGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGCGTCACCGGGCa                                                                                        <  1:1335815/100‑1 (MQ=255)
             cAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGCGTCACCGGGCaag                                                                                      <  2:1974063/100‑3 (MQ=255)
             cAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGCGTCACCGGGCaag                                                                                      <  2:853284/100‑3 (MQ=255)
             cAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGCGTCACCGGGCaag                                                                                      <  2:870236/100‑3 (MQ=255)
             cAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGCGTCACCGGGCaag                                                                                      <  2:1652822/100‑3 (MQ=255)
              aGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGCGTCACCGGGCaagt                                                                                     <  2:1713378/100‑4 (MQ=255)
               gTCGCGGCAACCGTGGCTGGTACACCTCTTACCCTGTTGGGCTGCTGACGCGCCCATTCATTCAGCAACTCTGCACTTGCCCGGTGACTCCGTGCa                                                                                        >  1:2465941/1‑96 (MQ=255)
               gTCGCGGCAACCGTGGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTACCTCTGCACTTGCCCGGCGTCACCGGGCaagtg                                                                                    <  1:2459938/100‑5 (MQ=255)
               gTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGCGTCACCGGGCaagtg                                                                                    <  2:1449433/100‑5 (MQ=255)
               gTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGCGTCACCGGGCaagtg                                                                                    <  2:543556/100‑5 (MQ=255)
                tCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATc                                                                                   <  1:521625/100‑1 (MQ=255)
                tCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATc                                                                                   <  2:919754/100‑1 (MQ=255)
                tCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGCGTCACCGGGCaagtgc                                                                                   <  2:1078539/100‑6 (MQ=255)
                 cgcgGCAACCGGCGCTGGTACACCACTGGCCCTGTTGGCCTGCTGGCCCACCGATGCACTAAGAAAGAACACACTTTCCAGGGGAAGCCGGGCATAATc                                                                                   >  1:2429320/1‑99 (MQ=255)
                  gcgGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCAc                                                                                 <  1:669379/100‑1 (MQ=255)
                  gcgGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCAc                                                                                 >  1:138340/1‑100 (MQ=255)
                   cgGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTCACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCAcc                                                                                >  1:174816/1‑100 (MQ=255)
                   ccgCAACAGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGCCGCCGGGCATTATCa                                                                                  <  1:875694/96‑1 (MQ=255)
                    ggCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCAccc                                                                               >  1:1723342/1‑100 (MQ=255)
                    ggCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCACGGTGACGCCGGGCATTATCAccc                                                                               >  1:1100473/1‑100 (MQ=255)
                      cAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCt                                                                              >  1:1881767/1‑99 (MQ=255)
                      cAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTg                                                                             <  2:1556586/100‑1 (MQ=255)
                      cAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTg                                                                             <  1:338836/100‑1 (MQ=255)
                        aCCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGcc                                                                           <  2:1104949/100‑1 (MQ=255)
                          cGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCaa                                                                         >  2:2112916/1‑100 (MQ=255)
                                tGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAg                                                                   <  2:1062115/100‑1 (MQ=255)
                                tGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAg                                                                   <  1:1097169/100‑1 (MQ=255)
                                 ggTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAgg                                                                  <  2:1070582/100‑1 (MQ=255)
                                 ggTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAgg                                                                  <  1:1901931/100‑1 (MQ=255)
                                 ggTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAgg                                                                  >  2:2570202/1‑100 (MQ=255)
                                   tACACCACTGGCCCTGTTGGTCTGCTGAACCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGaa                                                                >  1:1959186/1‑100 (MQ=255)
                                     caccacTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCAAAAAACAAAAAAAAGGAATa                                                              <  1:102431/100‑1 (MQ=255)
                                      accacTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGCGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATaa                                                             >  1:2478240/1‑100 (MQ=255)
                                        cacTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACACTGCCAAAAAAAGGAATAACa                                                           >  1:202535/1‑100 (MQ=255)
                                            ggCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGc                                                       >  1:1774491/1‑100 (MQ=255)
                                            ggCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGc                                                       >  1:11837/1‑100 (MQ=255)
                                             gCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCg                                                      >  1:472231/1‑100 (MQ=255)
                                              cccTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGc                                                       <  1:592136/98‑1 (MQ=255)
                                              cccTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGa                                                     >  1:199772/1‑100 (MQ=255)
                                              cccTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGa                                                     >  1:2291004/1‑100 (MQ=255)
                                               ccTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCAGGCATTATCACCCTGCCAAAAAAAGGCCTAACAATGCGAc                                                    >  1:1592953/1‑100 (MQ=255)
                                                  gTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATa                                                 >  2:2361655/1‑100 (MQ=255)
                                                  gTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATa                                                 >  2:580535/1‑100 (MQ=255)
                                                   ttGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAATAAAAGGAATAACAATGCGACATAg                                                <  1:787850/100‑1 (MQ=255)
                                                   ttGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAg                                                <  1:1476224/100‑1 (MQ=255)
                                                   ttGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAg                                                <  1:47945/100‑1 (MQ=255)
                                                    tGGTCTGCTGACCCGCCCCTTCATTAAGTAACTCTGCATTTGCCCGGTGCCGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGAGACATAGc                                               >  1:459649/1‑100 (MQ=255)
                                                         tgctgACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACt                                          >  1:6840/1‑100 (MQ=255)
                                                          gctgACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATGATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACt                                          >  1:513783/1‑99 (MQ=255)
                                                               tccgatctttcattAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGc                                    <  2:888007/92‑1 (MQ=255)
                                                               cccGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGc                                    <  2:965469/100‑1 (MQ=255)
                                                                ccGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCg                                   >  2:193328/1‑100 (MQ=255)
                                                                  gCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGAc                                 >  1:369429/1‑100 (MQ=255)
                                                                    ccATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAg                               >  2:1421684/1‑100 (MQ=255)
                                                                      attcattAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCg                             >  1:883022/1‑100 (MQ=255)
                                                                       ttcattAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAACGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGt                            >  1:1722795/1‑100 (MQ=255)
                                                                       ttcattAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGcagatcgga                            >  1:888004/1‑92 (MQ=255)
                                                                         cattAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCAAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGttt                          >  1:1973437/1‑100 (MQ=255)
                                                                         cattAAGTAACTCTGCACTTGCCCGGTGACGCAGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGttt                          >  2:208685/1‑100 (MQ=255)
                                                                          attAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTg                         >  1:1789979/1‑100 (MQ=255)
                                                                           ttAAGTACCTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTg                         <  2:566042/99‑1 (MQ=255)
                                                                           ttAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGc                        <  2:1525280/100‑1 (MQ=255)
                                                                            tAAGTAACTCTGCACTTGCCCGGTGACGCCGAGCATTATCACACTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGcc                       >  2:1748144/1‑100 (MQ=255)
                                                                             aaGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCa                      >  2:1541573/1‑100 (MQ=255)
                                                                             aaGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCa                      >  1:1879457/1‑100 (MQ=255)
                                                                              aGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCAc                     >  2:589329/1‑100 (MQ=255)
                                                                              aGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGCTTGCCAc                     <  2:1353687/100‑1 (MQ=255)
                                                                                tAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGc                   <  2:847953/100‑1 (MQ=255)
                                                                                  aCTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCt                  <  1:1366819/99‑1 (MQ=255)
                                                                                   ctctGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTa                >  1:1773716/1‑100 (MQ=255)
                                                                                   ctctGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTa                >  2:812738/1‑100 (MQ=255)
                                                                                     ctGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTAtc              <  2:778927/100‑1 (MQ=255)
                                                                                       gataatGCCCGGCGTCACCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTAtctc            >  1:409638/6‑100 (MQ=255)
                                                                                       gCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTAtctc            <  2:1029294/100‑1 (MQ=255)
                                                                                       gCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTAtctc            >  2:1406853/1‑100 (MQ=255)
                                                                                       gCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTAtctc            >  1:338580/1‑100 (MQ=255)
                                                                                         taatGCCCGGCGTCACCGGGCATTATCACCCTGACAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTa          >  1:2536848/4‑100 (MQ=255)
                                                                                         aCTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTa          <  2:649060/100‑1 (MQ=255)
                                                                                         aCTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTa          <  1:2242725/100‑1 (MQ=255)
                                                                                          aatGCCCGGCGTCACCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTac         >  1:1167258/3‑100 (MQ=255)
                                                                                           atGCCCGGCGTCACCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTaca        >  1:1796266/2‑100 (MQ=255)
                                                                                            tGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTacac       <  2:2338801/100‑1 (MQ=255)
                                                                                               ccGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACACaaa    >  2:487873/1‑100 (MQ=255)
                                                                                               ccGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCATATTTTTTGTGACAGCGTTTGCCACGCTTATCTCTACACaaa    >  1:1817806/1‑100 (MQ=255)
                                                                                                 ggTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACACAAACt  <  2:367953/100‑1 (MQ=255)
                                                                                                 ggTGACGCCGGGCATTATCACCCTACCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACACAAACt  <  1:113111/100‑1 (MQ=255)
                                                                                                   |                                                                                                 
GGCGCGACTCACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACACAAACT  >  NC_000913/2806992‑2807188

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