Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,807,093 A→T 15.5% intergenic (+19/‑39) 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,0930AT15.5% 231.4 / 30.9 97intergenic (+19/‑39)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 A (48/32);  new base T (4/11);  total (53/44)
Fisher's exact test for biased strand distribution p-value = 2.36e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.72e-01

Only 100 of 105 total aligned reads displayed.

CGCGACTCACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACACAAACTTTTG  >  NC_000913/2806994‑2807192
                                                                                                   |                                                                                                   
cgcgACTCACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGa                                                                                                     >  1:1808288/1‑100 (MQ=255)
cgcgACTCACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGa                                                                                                     <  2:830600/100‑1 (MQ=255)
cgcgACTCACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGa                                                                                                     <  1:283712/100‑1 (MQ=255)
   gACTCACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGc                                                                                                  <  1:2513645/100‑1 (MQ=255)
      tCACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCgg                                                                                               <  2:1950071/100‑1 (MQ=255)
       cACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCggg                                                                                              >  1:1847860/1‑100 (MQ=255)
       cACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCggg                                                                                              >  1:1010491/1‑100 (MQ=255)
       cACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCggg                                                                                              >  1:91521/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:543556/100‑5 (MQ=255)
             gTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGCGTCACCGGGCaagtg                                                                                        <  2:1449433/100‑5 (MQ=255)
              tCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATc                                                                                       <  2:919754/100‑1 (MQ=255)
              tCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATc                                                                                       <  1:521625/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                                                                       <  1:1097169/100‑1 (MQ=255)
                               ggTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAgg                                                                      <  2:1070582/100‑1 (MQ=255)
                               ggTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAgg                                                                      >  2:2570202/1‑100 (MQ=255)
                               ggTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAgg                                                                      <  1:1901931/100‑1 (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:11837/1‑100 (MQ=255)
                                          ggCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGc                                                           >  1:1774491/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:580535/1‑100 (MQ=255)
                                                 ttGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAATAAAAGGAATAACAATGCGACATAg                                                    <  1:787850/100‑1 (MQ=255)
                                                 ttGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAg                                                    <  1:47945/100‑1 (MQ=255)
                                                 ttGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAg                                                    <  1:1476224/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)
                                                              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)
                                                                         ttAAGTACCTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTg                             <  2:566042/99‑1 (MQ=255)
                                                                         ttAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGc                            <  2:1525280/100‑1 (MQ=255)
                                                                          tAAGTAACTCTGCACTTGCCCGGTGACGCCGAGCATTATCACACTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGcc                           >  2:1748144/1‑100 (MQ=255)
                                                                           aaGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCa                          >  1:1879457/1‑100 (MQ=255)
                                                                           aaGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCa                          >  2:1541573/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)
                                                                                       taatGCCCGGCGTCACCGGGCATTATCACCCTGACAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTa              >  1:2536848/4‑100 (MQ=255)
                                                                                       aCTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTa              <  1:2242725/100‑1 (MQ=255)
                                                                                       aCTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTa              <  2:649060/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)
                                                                                                  gACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACACAAACtttt   >  1:1023028/1‑100 (MQ=255)
                                                                                                  gACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACACAAACtttt   >  1:1342765/1‑100 (MQ=255)
                                                                                                   aCGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGTGACAGCGTTTGCCACGCTTATCTCTACACAAACTTTtg  <  2:1272029/100‑1 (MQ=255)
                                                                                                   aCGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACACAAACTTTtg  <  2:669652/100‑1 (MQ=255)
                                                                                                   aCGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACACAAACTTTtg  <  2:622450/100‑1 (MQ=255)
                                                                                                   aCGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACACAAACTTTtg  >  2:334836/1‑100 (MQ=255)
                                                                                                   aCGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACACAAACTTTtg  >  1:1418303/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
CGCGACTCACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACACAAACTTTTG  >  NC_000913/2806994‑2807192

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: