Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,937,143 A→T 5.6% D289V (GAT→GTT)  rbsB → D‑ribose ABC transporter periplasmic binding protein; ribose chemotaxis receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,937,1430AT 5.6% 1294.3 / 43.6 356D289V (GAT→GTT) rbsBD‑ribose ABC transporter periplasmic binding protein; ribose chemotaxis receptor
Reads supporting (aligned to +/- strand):  ref base A (178/157);  new base T (5/15);  total (183/173)
Fisher's exact test for biased strand distribution p-value = 1.97e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.05e-01

Only 100 of 357 total aligned reads displayed.

GACGGTACACCGGATGGCGAAAAAGCGGTGAATGATGGCAAACTAGCAGCGACTATCGCTCAGCTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGACATCCCGAA  >  NC_000913/3936995‑3937292
                                                                                                                                                    |                                                                                                                                                     
 aCGGTACACCGGATGGCGAAAAAGCGGTGAATGATGGCAAACTAGCAGCGACTATCGCTCAGCTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATc                                                                                                                                                     <  2:4061064/150‑1 (MQ=255)
 aCGGTACACCGGATGGCGAAAAAGCGGTGAATGATGGCAAACTAGCAGCGACTATCGCTCAGCTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGATGATc                                                                                                                                                     >  1:2934892/1‑150 (MQ=255)
  cGGTACACCGGATGGCGAAAAAGCGGTGAATGATGGCAAACTAGCAGCGACTATCGCTCAGCTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCt                                                                                                                                                    >  1:2714233/1‑150 (MQ=255)
   ggTACACCGGATGGCGAAAAAGCGGTGAATGATGGCAAACTAGCAGCGACTATCGCTCAGCTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTg                                                                                                                                                   >  2:3134200/1‑150 (MQ=255)
   ggTACACCGGATGGCGAAAAAGCGGTGAATGATGGCAAACTAGCAGCGACTATCGCTCAGCTACCCGATCAGAATGGCGCGAAAGGCGACGAAACCGCAGATAAAGTGATGAAAGGCGAGAAAGATCAGGCTAAGAATCCGGTTGATCTg                                                                                                                                                   >  2:3870414/1‑150 (MQ=255)
      acacCGGATGGCGAAAAAGCGGTGAATGATGGCAAACTAGCAGCGACTATCGCTCAGCTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGTTCTGaaa                                                                                                                                                <  1:483450/150‑1 (MQ=255)
      acacCGGATGGCGAAAAAGCGGTGAATGATGGCAAACTAGCAGCGACTATCGCTCAGCTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGTTCTGaaa                                                                                                                                                <  1:1864425/150‑1 (MQ=255)
      acacCGGATGGCGAAAAAGCGGTGAATGATGGCAAACTAGCAGCGACTATCGCTCAGCTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGaaa                                                                                                                                                >  1:1703962/1‑150 (MQ=255)
         cctgATGGCGAATAAGCGGTGAATGATGGCAAACTAGCTGCGACTATCGCTGAGCTACCCGATCAGATTGGCGCGAAAGGCGTCGATACCGCAGATAAAGTGCTGAAAGGCGTGAATGTTCAGGCTAAGTATCCGGTTGATCTGAAACTg                                                                                                                                             <  2:972527/147‑1 (MQ=255)
         ccGGATGGCGAAAAAGCGGTGAATGATGGCAAACTAGCAGCGACTATCGCTCAGCTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGTTCTGAAACTg                                                                                                                                             <  2:703623/150‑1 (MQ=255)
          cGGATGGCGAAAAAGCGGTGAATGATGGCAAACTAGCAGCGACTATCGCTCAGCTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTgg                                                                                                                                            >  1:1624945/1‑150 (MQ=255)
                gCGAAAAAGCGGTGAATGATGGCAAACTAGCAGCGACTATCGCTCAGCTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTa                                                                                                                                      <  2:1284545/150‑1 (MQ=255)
                 cGAAAAAGCGGTGAATGATGGCAAACTAGCAGCGACTATCGCTCAGCTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTaa                                                                                                                                     <  1:4086804/150‑1 (MQ=255)
                 cGAAAAAGCGGTGAATGATGGCAAACTAGCAGCGACTATCGCTCAGCTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGAGAAAGTGCTGAAAGGCTAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTaa                                                                                                                                     >  2:1166594/1‑150 (MQ=255)
                 cGAAAAAGCGGTGAATGATGGCAAACTAGCAGCGACTATCGCTCAGCTACCCGATCAGATAGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTaa                                                                                                                                     >  2:2241181/1‑150 (MQ=255)
                  gAAAAAGCGGTGAATGATGGCAAACTAGCAGCGACTATCGCTCAGCTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAg                                                                                                                                    <  2:209499/150‑1 (MQ=255)
                    aaaaGCGGTGAATGATGGCAAACTAGCAGCGACTATCGCTCAGCTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATACAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCa                                                                                                                                  <  1:2237002/150‑1 (MQ=255)
                       aGCGGTGAATGATGGCAAACTAGCAGCGACTATCGCTCAGCTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCagta                                                                                                                               <  1:1446915/150‑1 (MQ=255)
                       aGCGGTGAATGATGGCAAACTAGCAGCGACTATCGCTCAGCTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCagta                                                                                                                               <  1:3535620/150‑1 (MQ=255)
                       aGCGGTGAATGATGGCAAACTAGCAGCGACTATCGCTCAGCTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCagta                                                                                                                               <  1:3287244/150‑1 (MQ=255)
                       aGCGGTGAATGATGGCAAACTAGCAGCGACTATCGCTCAGCTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCagta                                                                                                                               <  1:3018032/150‑1 (MQ=255)
                       aGCGGTGAATGATGGCAAACTAGCAGCGACTATCGCTCAGCTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCagta                                                                                                                               <  2:1128770/150‑1 (MQ=255)
                        gCGGTGAATGATGGCAAACTAGCAGCGACTATCGCTCAGCTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCagtag                                                                                                                              >  1:1637425/1‑150 (MQ=255)
                         cGGTGAATGATGGCAAACTAGCAGCGACTATCGCTCAGCTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCagtagt                                                                                                                             >  2:392656/1‑150 (MQ=255)
                           gTGAATGATGGCAAACTAGCAGCGACTATCGCTCAGCTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAATGTGCTGAAAGGCGAGAACGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGttt                                                                                                                           >  1:1235393/1‑150 (MQ=255)
                           gTGAATGATGGCAAACTAGCAGCGACTATCGCTCAGCTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGttt                                                                                                                           <  2:3614124/150‑1 (MQ=255)
                           gTGAATGATGGCAAACTAGCAGCGACTATCGCTCAGCTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGttt                                                                                                                           <  2:3260618/150‑1 (MQ=255)
                            tGAATGATGGCAAACTAGCAGCGACTATCGCTCAGCTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGtttt                                                                                                                          >  2:296211/1‑150 (MQ=255)
                            tGAATGATGGCAAACTAGCAGCGACTATCGCTCAGCTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGtttt                                                                                                                          >  2:292814/1‑150 (MQ=255)
                             gAATGATGGCAAACTAGCAGCGACTATCGCTCAGCTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTa                                                                                                                         <  2:1657657/150‑1 (MQ=255)
                              aaTGATGGCAAACTAGCAGCGACTATCGCTCAGCTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTaa                                                                                                                        >  2:3766582/1‑150 (MQ=255)
                               atgatgGCAAACTAGCAGCGACTATCGCTCAGCTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGTTCTGAAACTGGTTGTTAAGCAGTAGTTTTAAt                                                                                                                       <  2:3441196/150‑1 (MQ=255)
                               atgatgGCAAACTAGCAGCGACTATCGCTCAGCTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAAt                                                                                                                       >  1:4188152/1‑150 (MQ=255)
                                tgatgGCAAACTAGCAGCGACTATCGCTCAGCTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATc                                                                                                                      >  1:2757090/1‑150 (MQ=255)
                                tgatgGCAAACTAGCAGCGACTATCGCTCAGCTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATc                                                                                                                      >  2:281140/1‑150 (MQ=255)
                                tgatgGCAAACTAGCAGCGACTATCGCTCAGCTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATc                                                                                                                      >  1:1978895/1‑150 (MQ=255)
                                 gatgGCAAACTAGCAGCGACTATCGCTCAGCTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCa                                                                                                                     >  1:2134250/1‑150 (MQ=255)
                                  atgGCAAACTAGCAGCGACTATCGCTCAGCTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAg                                                                                                                    >  1:3768375/1‑150 (MQ=255)
                                    ggCAAACTAGCAGCGACTATCGCTCAGCTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGt                                                                                                                  <  2:2563506/150‑1 (MQ=255)
                                       aaaCTAGCAGCGACTATCGCTCAGCTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCCGTAGTTTTAATCAGGTTGt                                                                                                               >  1:3354603/1‑150 (MQ=255)
                                       aaaCTAGCAGCGACTATCGCTCAGCTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGt                                                                                                               >  2:571981/1‑150 (MQ=255)
                                         aCTAGCAGCGACTATCGCTCAGCTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTAt                                                                                                             >  1:3260633/1‑150 (MQ=255)
                                              cagcGACTATCGCTCAGCTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCt                                                                                                        >  1:3935424/1‑150 (MQ=255)
                                              cagcGACTATCGCTCAGCTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCt                                                                                                        >  2:3037394/1‑150 (MQ=255)
                                                gcGACTATCGCTCAGCTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGa                                                                                                      <  1:772118/150‑1 (MQ=255)
                                                gcGACTATCGCTCAGCTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGa                                                                                                      >  1:4108503/1‑150 (MQ=255)
                                                             aGCTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAAt                                                                                         >  1:3213317/1‑150 (MQ=255)
                                                              gCTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATa                                                                                        >  1:980592/1‑150 (MQ=255)
                                                              gCTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATa                                                                                        <  2:737594/150‑1 (MQ=255)
                                                               cTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATAc                                                                                       <  1:1541813/150‑1 (MQ=255)
                                                               cTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATAc                                                                                       >  2:826108/1‑150 (MQ=255)
                                                                 aCCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGt                                                                                     <  2:2719007/150‑1 (MQ=255)
                                                                  cccGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTc                                                                                    >  1:4317059/1‑150 (MQ=255)
                                                                  cccGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTc                                                                                    >  1:3234054/1‑150 (MQ=255)
                                                                         aGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAATTACGTCATCTACa                                                                             >  2:1688839/1‑150 (MQ=255)
                                                                         aGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACa                                                                             <  1:1315549/150‑1 (MQ=255)
                                                                         aGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACa                                                                             <  1:1017223/150‑1 (MQ=255)
                                                                         aGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACa                                                                             >  2:1769246/1‑150 (MQ=255)
                                                                         aGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGAGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATTAATACGTCATCGACa                                                                             <  2:1816787/150‑1 (MQ=255)
                                                                          gATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAg                                                                            <  1:517888/150‑1 (MQ=255)
                                                                          gATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAg                                                                            <  1:1354703/150‑1 (MQ=255)
                                                                               gcgcGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGaa                                                                       >  2:3469736/1‑150 (MQ=255)
                                                                                cgcgAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGTTCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAAc                                                                      >  2:428024/1‑150 (MQ=255)
                                                                                  cgAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAGCgt                                                                    >  1:3886146/1‑150 (MQ=255)
                                                                                   gAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAATGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACgtg                                                                   >  2:1227654/1‑150 (MQ=255)
                                                                                   gAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACgtg                                                                   >  2:1790670/1‑150 (MQ=255)
                                                                                    aaaGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACgtgt                                                                  <  1:1280885/150‑1 (MQ=255)
                                                                                       ggCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGTTCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAat                                                               <  1:3162198/150‑1 (MQ=255)
                                                                                          gTCGAAACCGCAGATAAAGTTCATAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGGTGATCTGAAACGGGTTGTTAAGCAGTAGTTTTAATCAGGTGGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGACCGTGTAAtata                                                            >  2:4020355/1‑150 (MQ=255)
                                                                                          gTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAAtata                                                            >  1:2500179/1‑150 (MQ=255)
                                                                                           tCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATaa                                                           >  2:2603479/1‑150 (MQ=255)
                                                                                                  cGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGGCATCGACAGATGAACGTGTAATATAAAGAAAAg                                                    >  2:3349004/1‑150 (MQ=255)
                                                                                                     aGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAg                                                 >  2:2895748/1‑150 (MQ=255)
                                                                                                     aGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAg                                                 >  1:3616571/1‑150 (MQ=255)
                                                                                                      taTAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCAGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATTGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAgg                                                <  2:1392383/149‑1 (MQ=255)
                                                                                                      gATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAgg                                                <  2:4101149/150‑1 (MQ=255)
                                                                                                       aTAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAggg                                               <  1:3021607/150‑1 (MQ=255)
                                                                                                             tGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGTTCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCAcgcg                                         >  2:3748173/1‑150 (MQ=255)
                                                                                                                tGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCa                                      >  2:1497391/1‑150 (MQ=255)
                                                                                                                           aaaGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACgg                           >  1:90367/1‑150 (MQ=255)
                                                                                                                           aaaGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACgg                           >  1:61092/1‑150 (MQ=255)
                                                                                                                             aGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTg                         <  2:2699423/150‑1 (MQ=255)
                                                                                                                             aGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTg                         <  2:2438190/150‑1 (MQ=255)
                                                                                                                               ttCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGgc                       >  1:1067585/1‑150 (MQ=255)
                                                                                                                               ttCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGgc                       >  1:3077715/1‑150 (MQ=255)
                                                                                                                                 cAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGgcgc                     >  1:70102/1‑150 (MQ=255)
                                                                                                                                  attcTAAGTATCCGGATGATCAGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTAAGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTCACACGGTGGCGCa                    >  2:658312/4‑150 (MQ=255)
                                                                                                                                  aGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCa                    <  1:200643/150‑1 (MQ=255)
                                                                                                                                   ggCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCATGCGCCACCGTGTTAGGGTGgcgcgt                   <  1:2065359/150‑3 (MQ=255)
                                                                                                                                   ggCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCAt                   >  1:1478957/1‑150 (MQ=255)
                                                                                                                                        aGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAACGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTAt              <  2:2775187/150‑1 (MQ=255)
                                                                                                                                         gTATCCGGTTGATCTGAAAGTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGATCGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATg             <  2:2369695/150‑1 (MQ=255)
                                                                                                                                         gTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATg             <  1:2862834/150‑1 (MQ=255)
                                                                                                                                         gTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATg             <  1:1174810/150‑1 (MQ=255)
                                                                                                                                           aTCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGa           >  2:1976500/1‑150 (MQ=255)
                                                                                                                                            tCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGAc          >  1:2283885/1‑150 (MQ=255)
                                                                                                                                            tCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGAc          >  1:2633765/1‑150 (MQ=255)
                                                                                                                                            tCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGAc          >  1:4234472/1‑150 (MQ=255)
                                                                                                                                                gTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGACATcc      >  1:1140911/1‑150 (MQ=255)
                                                                                                                                                  tGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGACATCCCg    >  2:1346678/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GACGGTACACCGGATGGCGAAAAAGCGGTGAATGATGGCAAACTAGCAGCGACTATCGCTCAGCTACCCGATCAGATTGGCGCGAAAGGCGTCGAAACCGCAGATAAAGTGCTGAAAGGCGAGAAAGTTCAGGCTAAGTATCCGGTTGATCTGAAACTGGTTGTTAAGCAGTAGTTTTAATCAGGTTGTATGACCTGATGGTGACATAAATACGTCATCGACAGATGAACGTGTAATATAAAGAAAAGCAGGGCACGCGCCACCCTAACACGGTGGCGCATTTTATGGACATCCCGAA  >  NC_000913/3936995‑3937292

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