Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,807,091 T→C 6.7% 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,0910TC 6.7% 340.2 / 10.9 104intergenic (+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 (58/39);  new base C (2/5);  total (60/44)
Fisher's exact test for biased strand distribution p-value = 1.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.33e-01

Only 100 of 109 total aligned reads displayed.

GGGATTGTGATCCTCGCCATTATCCTCGATCGTCTGACGCAGGCCGTTGGGCGCGACTCACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCC‑AAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACACAAACTTTTGCTGCCGATCTGCCGGGCAAAGGCATTACTGTTAATCCAGTTCAGA  >  NC_000913/2806943‑2807237
                                                                                                                                                    |                                                                                                                                                   
 ggATTGTGATCCTCGCCATTATCCTCGATCGTCTGACGCAGGCCGTTGGGCGCGACTCACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGt                                                                                                                                                     >  1:67278/1‑148 (MQ=255)
  gATTGTGATCCTCGCCATTATCCTCGATCGTCTGACGCAGGCCGTTGGGCGCGACTCACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGa                                                                                                                                                   >  4:179496/1‑149 (MQ=255)
   aTTGTGATCCTCGCCATTATCCTCGATCGTCTGACGCAGGCCGTTGGGCGCGACTCACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGAc                                                                                                                                                  <  3:191449/149‑1 (MQ=255)
           ccTCGCCATTATCCTCGATCGTCTGACGCAGGCCGTTGGGCGCGACTCACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGc                                                                                                                                           >  1:466353/1‑148 (MQ=255)
            cTCGCCATTATCCTCGATCGTCTGACGCAGGCCGTTGGGCGCGACTCACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCGTTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCAt                                                                                                                                         >  4:739547/1‑149 (MQ=255)
            cTCGCCATTATCCTCGATCGTCTGACGCAGGCCGTTGGGCGCGACTCACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCAt                                                                                                                                         >  4:318911/1‑149 (MQ=255)
            cTCGCCATTATCCTCGATCGTCTGACGCAGGCCGTTGGGCGCGACTCACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCAt                                                                                                                                         <  4:761790/149‑1 (MQ=255)
            cTCGCCATTATCCTCGATCGTCTGACGCAGGCCGTTGGGCGCGACTCACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCAt                                                                                                                                         >  4:425449/1‑149 (MQ=255)
            cTCGCCATTATCCTCGATCGTCTGACGCAGGCCGTTGGGCGCGACTCACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGCGTCACCGGGCaa                                                                                                                                         <  4:191763/149‑2 (MQ=255)
            cTCGCCATTATCCTCGATCGTCTGACGCAGGCCGTTGGGCGCGACTCACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGCGTCACCGGGCaa                                                                                                                                         <  4:478305/149‑2 (MQ=255)
            cTCGCCATTATCCTCGATCGTCTGACGCAGGCCGTTGGGCGCGACTCACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGCGTCACCGGGCa                                                                                                                                          <  1:304635/148‑1 (MQ=255)
            cTCGCCATTATCCTCGATCGTCTGACGCAGGCCGTTGGGCGCGACTCACGCAGTCGAGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGCGTCACCGGGCaa                                                                                                                                         <  3:555537/149‑2 (MQ=255)
             tCGCCATTATCCTCGATCGTCTGACGCAGGCCGTTGGGCGCGACTCACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCAtt                                                                                                                                        <  4:376497/149‑1 (MQ=255)
             tCGCCATTATCCTCGATCGTCTGACGCAGGCCGTTGGGCGCGACTCACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCAt                                                                                                                                         <  2:714485/148‑1 (MQ=255)
              cGCCATTATCCTCGATCGTCTGACGCAGGCCGTTGGGCGCGACTCACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCAtt                                                                                                                                        >  1:383510/1‑148 (MQ=255)
               gCCATTATCCTCGATCGTCTGACGCAGGCCGTTGGGCGCGACTCACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTa                                                                                                                                       >  1:905546/1‑148 (MQ=255)
                ccATTATCCTCGATCGTCTGACGCAGGCCGTTGGGCGCGACTCACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATc                                                                                                                                     <  3:221231/149‑1 (MQ=255)
                  aTTATCCTTGATCGTCTGACGCAGGCCGTTGGGCGCGACTCACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGt                                                                                                                                                     >  4:78355/1‑131 (MQ=255)
                  aTTATCCTCGATCGTCTGACGCAGGCCGTTGGGCGCGACTCACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGt                                                                                                                                                     <  3:78355/131‑1 (MQ=255)
                  aTTATCCTCGATCGTCTGACGCAGGCCGTTGGGCGCGACTCACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCa                                                                                                                                    <  1:880602/148‑1 (MQ=255)
                  aTTATCCTCGATCGTCTGACGCAGGCCGTTGGGCGCGACTCACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCa                                                                                                                                    <  2:154935/148‑1 (MQ=255)
                  aTTATCCTCGATCGTCTGACGCAGGCCGTTGGGCGCGACTCACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCa                                                                                                                                    >  2:880602/1‑148 (MQ=255)
                      tCCTCGATCGTCTGACGCAGGCCGTTGGGCGCGACTCACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCAcccag                                                                                                                               >  4:749439/1‑147 (MQ=255)
                         tCGATCGTCTGACGCAGGCCGTTGGGCGCGACTCACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGcc                                                                                                                             >  1:408328/1‑148 (MQ=255)
                          cGATCGTCTGACGCAGGCCGTTGGGCGCGACTCACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCC‑a                                                                                                                           >  1:541447/1‑148 (MQ=255)
                           gATCGTCTGACGCAGGCCGTTGGGCGCGACTCACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCC‑aaa                                                                                                                         >  4:24008/1‑149 (MQ=255)
                            aTCGTCTGACGCAGGCCGTTGGGCGCGACTCACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCC‑aaaa                                                                                                                        >  3:861360/1‑149 (MQ=255)
                               gTCTGACGCAGGCCGTTGGGCGCGACTCACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCC‑aaaaaaa                                                                                                                     <  3:55441/149‑1 (MQ=255)
                               gTCTGACGCAGGCCGTTGGGCGCGACTCACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCC‑aaaaaa                                                                                                                      <  2:541447/148‑1 (MQ=255)
                                tCTGACGCAGGCCGTTGGGCGCGACTCACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCC‑aaaaaaa                                                                                                                     <  2:433378/148‑1 (MQ=255)
                                tCTGACGCAGGCCGTTGGGCGCGACTCACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCC‑aaaaaaa                                                                                                                     <  2:856497/148‑1 (MQ=255)
                                   gACGCAGGCCGTTGGGCGCGACTCACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCC‑AAAAAAAGGa                                                                                                                  >  2:222254/1‑148 (MQ=255)
                                    aCGCAGGCCGTTGGGCGCGACTCACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCC‑AAAAAAAGGaa                                                                                                                 >  1:773684/1‑148 (MQ=255)
                                      gCAGGCCGTTGGGCGCGACTCACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCC‑AAAAAAAGGAATaa                                                                                                              >  3:264656/1‑149 (MQ=255)
                                        aGGCCGTTGGGCGCGACTCACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCC‑AAAAAAAGGAATAACa                                                                                                            >  4:811044/1‑149 (MQ=255)
                                         ggCCGTTGGGCGCGACTCACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCC‑AAAAAAAGGAATAACa                                                                                                            >  2:663970/1‑148 (MQ=255)
                                                 ggCGCGACTCACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCC‑AAAAAAAGGAATAACAATGCGACAt                                                                                                   >  4:130033/1‑149 (MQ=255)
                                                  gcgcGACTCACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCgg                                                                                                                                             <  3:305631/107‑1 (MQ=255)
                                                   cgcgACTCACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCC‑AAAAAAAGGAATAACGATg                                                                                                         <  2:209594/141‑1 (MQ=255)
                                                   cgcgACTCACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCC‑AAAAAAAGGAATAACGATg                                                                                                         >  1:209594/1‑141 (MQ=255)
                                                    gcgACTCACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCC‑AAAAAAAGGAATAACAATGCGACATAGc                                                                                                >  3:279956/1‑149 (MQ=255)
                                                       aCTCACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCC‑AAAAAAAGGAAAAACAATGCGACATAGCGTa                                                                                             <  3:716509/149‑1 (MQ=255)
                                                          cACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCC‑AAAAAAAGGAATAACAATGCGACATAGCGTACtt                                                                                          >  3:665935/1‑149 (MQ=255)
                                                          cACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCC‑AAAAAAAGGAATAACAATGCGACATAGCGTACt                                                                                           >  1:552455/1‑148 (MQ=255)
                                                           aCGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCC‑AAAAAAAGGAATAACAATGCGACATAGCGTACttt                                                                                         >  3:33353/1‑149 (MQ=255)
                                                            cGCAGACGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCC‑AAAAAAAGGAATAACAATGCGACATAGCGTACttt                                                                                         >  2:647830/1‑148 (MQ=255)
                                                             gCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCC‑AAAAAAAGGAATAACAATGCGACATAGCGTACtttt                                                                                        <  2:466353/148‑1 (MQ=255)
                                                             gCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCC‑AAAAAAAGGAATAACAATGCGACATAGCGTACtttt                                                                                        <  2:588449/148‑1 (MQ=255)
                                                              cAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCC‑AAAAAAAGGAATAACAATGCGACATAGCGTACttttt                                                                                       >  1:718771/1‑148 (MQ=255)
                                                               aGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACTCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCC‑AAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTg                                                                                      <  2:752085/148‑1 (MQ=255)
                                                                      gCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCC‑AAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCg                                                                              >  3:377555/1‑149 (MQ=255)
                                                                           cGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCC‑AAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTg                                                                          <  2:67278/148‑1 (MQ=255)
                                                                            gTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATGACCCTGCC‑AAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGcc                                                                        >  3:625814/1‑149 (MQ=255)
                                                                                      caccacTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCC‑AAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTAtc                                                               >  2:81367/1‑148 (MQ=255)
                                                                                       accacTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCC‑AAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTAtctc                                                             >  4:456264/1‑149 (MQ=255)
                                                                                        ccacTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCC‑AAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTAtctc                                                             >  2:498825/1‑148 (MQ=255)
                                                                                          acTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCC‑AAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTac                                                          <  3:179496/149‑1 (MQ=255)
                                                                                              gCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCC‑AAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACACaa                                                      <  4:33353/149‑1 (MQ=255)
                                                                                                ccTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCC‑AAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACACAAAc                                                    >  4:519316/1‑149 (MQ=255)
                                                                                                  tGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACACAAAc                                                    <  1:303899/148‑1 (MQ=255)
                                                                                                    ttGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCC‑AAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACACAAACtttt                                                <  3:24008/149‑1 (MQ=255)
                                                                                                    ttGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCC‑AAAAAAAGGAATAACAATGCGACATAGCGAACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACACAAACttt                                                 <  2:773684/148‑1 (MQ=255)
                                                                                                     tGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCC‑AAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACACAAACtttt                                                >  1:117705/1‑148 (MQ=255)
                                                                                                     tGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCC‑AAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACACAAACtttt                                                >  1:434208/1‑148 (MQ=255)
                                                                                                        tCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCC‑AAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACACAAACTTTtgct                                             <  2:886241/148‑1 (MQ=255)
                                                                                                            ctgACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCC‑AAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACACAAACTTTTGCTGCCGa                                        <  4:279956/149‑1 (MQ=255)
                                                                                                             tgACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCC‑AAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACACAAACTTTTGCTGCCGAt                                       <  3:9326/149‑1 (MQ=255)
                                                                                                                  cGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCC‑AAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACACAAACTTTTGCTGCCGATCTGc                                   >  1:799856/1‑148 (MQ=255)
                                                                                                                    cccATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCC‑AAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACACAAACTTTTGCTGCCGATCTGCCg                                 <  2:347808/148‑1 (MQ=255)
                                                                                                                     ccATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCC‑AAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTCGCCACGCTTATCTCTACACAAACTTTTGCTGCCGATCTGCCggg                               >  3:773869/1‑149 (MQ=255)
                                                                                                                     ccATTCATTAAATAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCC‑AAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACACAAACTTTTGCTGCCGATCTGCCggg                               >  3:783425/1‑149 (MQ=255)
                                                                                                                      cattcattAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCC‑AAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACACAAACTTTTGCTGCCGATCTGCCggg                               <  2:905546/148‑1 (MQ=255)
                                                                                                                       attcattAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCC‑AAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACACAAACTTTTGCTGCCGATCTGCCGGGCa                             >  3:5004/1‑149 (MQ=255)
                                                                                                                            ttAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCC‑AAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACACAAACTTTTGCTGCCGATCTGCCGGGCAAAgg                         >  1:113062/1‑148 (MQ=255)
                                                                                                                            ttAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCC‑AAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACACAAACTTTTGCTGCCGATCTGCCGGGCAAAGGc                        >  4:886867/1‑149 (MQ=255)
                                                                                                                            ttAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCC‑AAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACACAAACTTTTGCTGCCGATCTGCCGGGCAAAGGc                        >  4:345212/1‑149 (MQ=255)
                                                                                                                             tAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCC‑AAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACACAAACTTTTGCTGCCGATCTGCCGGGCAAAGGc                        >  2:216464/1‑148 (MQ=255)
                                                                                                                               aGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCC‑AAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACACAAACTTTTGCTGCCGATCTGCCGGGCAAAGGCAtt                     >  3:319719/1‑149 (MQ=255)
                                                                                                                                gTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCC‑AAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACACAAACTTTTGCTGCCGATCTGCCGGGCAAAGGCAtt                     <  1:216464/148‑1 (MQ=255)
                                                                                                                                gTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCC‑AAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACACAAACTTTTGCTGCCGATCTGCCGGGCAAAGGCAtt                     >  2:22535/1‑148 (MQ=255)
                                                                                                                                gTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCC‑AAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACACAAACTTTTGCTGCCGATCTGCCGGGCAAAGGCATTa                    >  4:793840/1‑149 (MQ=255)
                                                                                                                                  aaCTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCC‑AAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACACAAACTTTTGCTGCCGATCTGCCGGGCAAAGGCa                       <  4:603634/144‑1 (MQ=255)
                                                                                                                                  aaCTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCC‑AAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACACAAACTTTTGCTGCCGATCTGCCGGGCAAAGGCa                       >  3:603634/1‑144 (MQ=255)
                                                                                                                                  aaCTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCC‑AAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACACAAACTTTTGCTGCCGATCTGCCGGGCAAAGGCATTAc                   >  1:176491/1‑148 (MQ=255)
                                                                                                                                  aaCTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCC‑AAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACACAAACTTTTGCTGCCGATCTGCCGGGCAAAGGCATTACt                  <  3:130033/149‑1 (MQ=255)
                                                                                                                                    ctctGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCC‑AAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACACaaa                                                     <  2:531861/112‑1 (MQ=255)
                                                                                                                                    ctctGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCC‑AAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACACaaa                                                     >  1:531861/1‑112 (MQ=255)
                                                                                                                                     tctGCACTTGCCCGGTGACGCCGGGCATTATCACCCTTCC‑AAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACACAAACTTTTGCTGCCGATCTGCCGGGCAAAGGCATTACTGt                <  2:718771/148‑1 (MQ=255)
                                                                                                                                       tGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCC‑AAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACACAAACTTTTGCTGCCGATCTGCCGGGCAAAGGCATTACTGTTaa             <  4:5004/149‑1 (MQ=255)
                                                                                                                                       tGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCC‑AAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACACAAACTTTTGCTGCCGATCTGCCGGGCAAAGGCATTACTGTTaa             >  3:805605/1‑149 (MQ=255)
                                                                                                                                           cTTGCCCGGTGACGCCGGGCATTATCACCCTGCC‑AAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACACAAACTTTTGCTGCCGATCTGCCGGGCAAAGGCATTACTGTTAATcc          <  2:273571/148‑1 (MQ=255)
                                                                                                                                           aatGCCCGGCGTCACCGGGCATTATCACCCTGCC‑AAAAAAAGGAATAGCAATGCGACATAGCGTACTTTTTGCGACAGCGttt                                                                           >  3:616401/3‑83 (MQ=255)
                                                                                                                                           aatGCCCGGCGTCACCGGGCATTATCACCCTGCC‑AAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGttt                                                                           <  4:616401/81‑1 (MQ=255)
                                                                                                                                             tGCCCGGTGACGCCGGGCATTATCACCCTGCC‑AAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACACAAACTTTTGCTGCCGATCTGCCGGGCAAAGGCATTACTGTTAATCCAg        >  1:14188/1‑148 (MQ=255)
                                                                                                                                             tGCCCGGCGTCACCGGGCATTATCACCCTGCC‑AAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACACAAACTTTTGCTGCCGATCTGCCGGGCAAAGGCATTACTGTTAATCCAGt       >  4:438227/1‑149 (MQ=255)
                                                                                                                                              gCCCGGTGACGCCGGGCATTATCACCCTGCC‑AAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACACAAACTTTTGCTGCCGATCTGCCGGGCAAAGGCATTACTGTTAATCCAGtt      >  4:615018/1‑149 (MQ=255)
                                                                                                                                                ccGGTGACGCCGGGCATTATCACCCTGCC‑AAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACACAAACTTTTGCTGCCGATCTGCCGGGCAAAGGCATTACTGTTAATCCAGTTCa    >  4:383157/1‑149 (MQ=255)
                                                                                                                                                ccGGTGACGCCGGGCATTATCACCCTGCC‑AAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACACAAACTTTTGCTGCCGATCTGCCGGGCAAAGGCATTACTGTTAATCCAGTTCa    >  3:699705/1‑149 (MQ=255)
                                                                                                                                                 cGGTGACGCCGGGCATTATCACCCTGCC‑AAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACACAAACTTTTGCTGCCGATCTGCCGGGCAAAGGCATTACTGTTAATCCAGTTCag   >  4:818663/1‑149 (MQ=255)
                                                                                                                                                  ggTGACGCCGGGCATTATCACCCTGCC‑AAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACACAAACTTTTGCTGCCGATCTGCCGGGCAAAGGCATTACTGTTAATCCAGTTCaga  <  3:425449/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GGGATTGTGATCCTCGCCATTATCCTCGATCGTCTGACGCAGGCCGTTGGGCGCGACTCACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCC‑AAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACACAAACTTTTGCTGCCGATCTGCCGGGCAAAGGCATTACTGTTAATCCAGTTCAGA  >  NC_000913/2806943‑2807237

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑