Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,807,095 G→A 6.8% intergenic (+21/‑37) 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,0950GA 6.8% 332.8 / 10.0 103intergenic (+21/‑37)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 G (57/39);  new base A (2/5);  total (59/44)
Fisher's exact test for biased strand distribution p-value = 1.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01

Only 100 of 104 total aligned reads displayed.

TCCTCGCCATTATCCTCGATCGTCTGACGCAGGCCGTTGGGCGCGACTCACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCC‑AAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACACAAACTTTTGCTGCCGATCTGCCGGGCAAAGGCATTACTGTTAATCCAGTTCAGAGCACC  >  NC_000913/2806953‑2807242
                                                                                                                                              |                                                                                                                                                    
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)
  cTCGCCATTATCCTCGATCGTCTGACGCAGGCCGTTGGGCGCGACTCACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCAt                                                                                                                                              >  4:318911/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)
        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)
                          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)
                                        gcgcGACTCACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCgg                                                                                                                                                  >  4:305631/1‑107 (MQ=255)
                                        gcgcGACTCACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCgg                                                                                                                                                  <  3:305631/107‑1 (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: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)
                                                                                    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)
                                                                                                             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                  >  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             >  2:903946/1‑148 (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)
                                                                                                                                        gatGACGCCGGGCATTATCACCCTGCC‑AAAAAAACGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACACAAACTTTTGCTGCCGATCTGCCGGGCAAAGGCATTACTGTTAATCCAGTTCaga       <  3:739547/147‑1 (MQ=255)
                                                                                                                                              gCCGGGCATTATCACCCTGCC‑AAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACACAAACTTTTGCTGCCGATCTGCCGGGCAAAGGCATTACTGTTAATCCAGTTCAGAGCAcc  >  2:581992/1‑148 (MQ=255)
                                                                                                                                              |                                                                                                                                                    
TCCTCGCCATTATCCTCGATCGTCTGACGCAGGCCGTTGGGCGCGACTCACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCC‑AAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACACAAACTTTTGCTGCCGATCTGCCGGGCAAAGGCATTACTGTTAATCCAGTTCAGAGCACC  >  NC_000913/2806953‑2807242

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: ‑