Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,807,093 A→T 6.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,0930AT 6.5% 361.6 / 10.6 104intergenic (+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 (57/40);  new base T (2/5);  total (59/45)
Fisher's exact test for biased strand distribution p-value = 2.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

Only 100 of 105 total aligned reads displayed.

GATTGTGATCCTCGCCATTATCCTCGATCGTCTGACGCAGGCCGTTGGGCGCGACTCACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCC‑AAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACACAAACTTTTGCTGCCGATCTGCCGGGCAAAGGCATTACTGTTAATCCAGTTCAGA  >  NC_000913/2806945‑2807237
                                                                                                                                                    |                                                                                                                                                 
gATTGTGATCCTCGCCATTATCCTCGATCGTCTGACGCAGGCCGTTGGGCGCGACTCACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGa                                                                                                                                                   >  4:179496/1‑149 (MQ=255)
 aTTGTGATCCTCGCCATTATCCTCGATCGTCTGACGCAGGCCGTTGGGCGCGACTCACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGAc                                                                                                                                                  <  3:191449/149‑1 (MQ=255)
        tCCTCGCCATTATCCTCGATCGTCTGACGCAGGCCGTTGGGCGCGACTCACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCggg                                                                                                                                            <  2:597440/148‑1 (MQ=255)
         ccTCGCCATTATCCTCGATCGTCTGACGCAGGCCGTTGGGCGCGACTCACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGc                                                                                                                                           >  1:466353/1‑148 (MQ=255)
          cTCGCCATTATCCTCGATCGTCTGACGCAGGCCGTTGGGCGCGACTCACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCGTTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCAt                                                                                                                                         >  4:739547/1‑149 (MQ=255)
          cTCGCCATTATCCTCGATCGTCTGACGCAGGCCGTTGGGCGCGACTCACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCAt                                                                                                                                         >  4:425449/1‑149 (MQ=255)
          cTCGCCATTATCCTCGATCGTCTGACGCAGGCCGTTGGGCGCGACTCACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCAt                                                                                                                                         <  4:761790/149‑1 (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)
              ccATTATCCTCGATCGTCTGACGCAGGCCGTTGGGCGCGACTCACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATc                                                                                                                                     <  3:221231/149‑1 (MQ=255)
                aTTATCCTCGATCGTCTGACGCAGGCCGTTGGGCGCGACTCACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCa                                                                                                                                    >  2:880602/1‑148 (MQ=255)
                aTTATCCTCGATCGTCTGACGCAGGCCGTTGGGCGCGACTCACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCa                                                                                                                                    <  1:880602/148‑1 (MQ=255)
                aTTATCCTCGATCGTCTGACGCAGGCCGTTGGGCGCGACTCACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCa                                                                                                                                    <  2:154935/148‑1 (MQ=255)
                   aTCCTCGATCGTCTGACGCAGGCCGTTGGGCGCGACTCACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCAccc                                                                                                                                 <  2:212624/148‑1 (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:856497/148‑1 (MQ=255)
                              tCTGACGCAGGCCGTTGGGCGCGACTCACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCC‑aaaaaaa                                                                                                                     <  2:433378/148‑1 (MQ=255)
                              tCTGACGCAGGCCGTTGGGCGCGACTCACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCC‑AAAAAAAg                                                                                                                    <  3:348205/149‑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)
                                                gcgcGACTCACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCgg                                                                                                                                             >  4:305631/1‑107 (MQ=255)
                                                 cgcgACTCACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCC‑AAAAAAAGGAATAACGATg                                                                                                         >  1:209594/1‑141 (MQ=255)
                                                 cgcgACTCACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCC‑AAAAAAAGGAATAACGATg                                                                                                         <  2:209594/141‑1 (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:588449/148‑1 (MQ=255)
                                                           gCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCC‑AAAAAAAGGAATAACAATGCGACATAGCGTACtttt                                                                                        <  2:466353/148‑1 (MQ=255)
                                                            cAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCC‑AAAAAAAGGAATAACAATGCGACATAGCGTACttttt                                                                                       >  1:718771/1‑148 (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)
                                                                                               cTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCC‑AAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACACAAAc                                                    >  1:519553/1‑148 (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)
                                                                                                                          ttAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCC‑AAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACACAAACTTTTGCTGCCGATCTGCCGGGCAAAgg                         >  1:113062/1‑148 (MQ=255)
                                                                                                                          ttAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCC‑AAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACACAAACTTTTGCTGCCGATCTGCCGGGCAAAGGc                        >  4:345212/1‑149 (MQ=255)
                                                                                                                          ttAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCC‑AAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACACAAACTTTTGCTGCCGATCTGCCGGGCAAAGGc                        >  4:886867/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                       >  3:603634/1‑144 (MQ=255)
                                                                                                                                aaCTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCC‑AAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACACAAACTTTTGCTGCCGATCTGCCGGGCAAAGGCa                       <  4:603634/144‑1 (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             >  3:805605/1‑149 (MQ=255)
                                                                                                                                     tGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCC‑AAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACACAAACTTTTGCTGCCGATCTGCCGGGCAAAGGCATTACTGTTaa             <  4:5004/149‑1 (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)
                                                                                                                                           tGCCCGGTGACGCCGGGCATTATCACCCTGCC‑AAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACACAAACTTTTGCTGCCGATCTGCCGGGCAAAGGCATTACTGTTAATCCAg        >  2:903946/1‑148 (MQ=255)
                                                                                                                                           tGCCCGGCGTCACCGGGCATTATCACCCTGCC‑AAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACACAAACTTTTGCTGCCGATCTGCCGGGCAAAGGCATTACTGTTAATCCAGt       >  4:438227/1‑149 (MQ=255)
                                                                                                                                            gCCCGGTGACGCCGGGCATTATCACCCTGCC‑AAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACACAAACTTTTGCTGCCGATCTGCCGGGCAAAGGCATTACTGTTAATCCAGtt      >  4:615018/1‑149 (MQ=255)
                                                                                                                                              ccGGTGACGCCGGGCATTATCACCCTGCC‑AAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACACAAACTTTTGCTGCCGATCTGCCGGGCAAAGGCATTACTGTTAATCCAGTTCa    >  3:699705/1‑149 (MQ=255)
                                                                                                                                              ccGGTGACGCCGGGCATTATCACCCTGCC‑AAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACACAAACTTTTGCTGCCGATCTGCCGGGCAAAGGCATTACTGTTAATCCAGTTCa    >  4:383157/1‑149 (MQ=255)
                                                                                                                                               cGGTGACGCCGGGCATTATCACCCTGCC‑AAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACACAAACTTTTGCTGCCGATCTGCCGGGCAAAGGCATTACTGTTAATCCAGTTCag   >  4:818663/1‑149 (MQ=255)
                                                                                                                                                ggTGACGCCGGGCATTATCACCCTGCC‑AAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACACAAACTTTTGCTGCCGATCTGCCGGGCAAAGGCATTACTGTTAATCCAGTTCaga  <  3:425449/149‑1 (MQ=255)
                                                                                                                                                gatGACGCCGGGCATTATCACCCTGCC‑AAAAAAACGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACACAAACTTTTGCTGCCGATCTGCCGGGCAAAGGCATTACTGTTAATCCAGTTCaga  <  3:739547/147‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
GATTGTGATCCTCGCCATTATCCTCGATCGTCTGACGCAGGCCGTTGGGCGCGACTCACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCC‑AAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACACAAACTTTTGCTGCCGATCTGCCGGGCAAAGGCATTACTGTTAATCCAGTTCAGA  >  NC_000913/2806945‑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: ‑