Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,092,780 C→A 5.2% G15V (GGC→GTC)  frvA ← putative enzyme IIA component of PTS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,092,7800CA 5.2% 323.2 / 3.3 109G15V (GGC→GTC) frvAputative enzyme IIA component of PTS
Reads supporting (aligned to +/- strand):  ref base C (49/54);  new base A (4/2);  total (53/56)
Fisher's exact test for biased strand distribution p-value = 4.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.51e-01

Only 100 of 114 total aligned reads displayed.

CGTGCGGCACGGCGACACCGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAACTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGATTAAGTGATGAGAGGCGATATGTCGGTAAACAATGACAATTTGAGACAGAGTGAAAG  >  NC_000913/4092631‑4092929
                                                                                                                                                     |                                                                                                                                                     
cGTGCGGCACGGCGACACCGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAACTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGc                                                                                                                                                       <  2:1343972/150‑1 (MQ=255)
  tGCGGCACGGCTACACCGGAACCAAATCCAGTAGAGTGGATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAACTTGTGTGAGGCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGccc                                                                                                                                                     <  2:924423/150‑1 (MQ=255)
  tGCGGCACGGCGACACCGGAACCAAATCCAGTAGAGTGCCTGTTTTCGCGCAGCAGCAGGGTTTGCAGGAACTGGTGTGAGTCGGTGATAAAACCGTTTTGTACCGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGccc                                                                                                                                                     <  2:268570/150‑1 (MQ=255)
  tGCGGCACGGCGACACCGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAACTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGccc                                                                                                                                                     >  2:117705/1‑150 (MQ=255)
   gCGGCACGGCGACACCGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAACTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCt                                                                                                                                                    <  2:769970/150‑1 (MQ=255)
     ggCACGGCGACACCGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAACTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGa                                                                                                                                                  >  2:252746/1‑150 (MQ=255)
     ggCACGGCGACACCGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAACTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGa                                                                                                                                                  >  1:349284/1‑150 (MQ=255)
     ggCACGGCGACACCGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAACTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGa                                                                                                                                                  >  1:1330759/1‑150 (MQ=255)
           gCGACACCGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAACTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGACAACTGCTTCAGAACGGAATAAGCTCCATTGCCCTGAATATTc                                                                                                                                            >  1:1681573/1‑150 (MQ=255)
           gCGACACCGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAACTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTAATGGCGCCAACTGCTTCAGAACGGTATAAGCGCCATTGCCCTGAATATTc                                                                                                                                            <  2:277498/150‑1 (MQ=255)
            cGACACCGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAACTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCa                                                                                                                                           >  1:950049/1‑150 (MQ=255)
            cGACACCGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAACTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTCCCCTGAATATTCa                                                                                                                                           <  1:1239291/150‑1 (MQ=255)
              acacCGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAACTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTACAGAACGGAATAATCGCCATTGCCCTGAATATTAAgg                                                                                                                                         >  2:1351298/1‑150 (MQ=255)
              acacCGGAACCAAATCCAGTAGAGTGCAATTTTTCGTTCATCAGCAGGGTTTGGAGGAAATGTTGTTATTCGGTGATAAAACCGTTTTGTAACGCTATTGTAGCCAACTGCTTCAGAACGGAATAAGCGCCTTTGCCCTGAATATTCAgg                                                                                                                                         <  1:1133575/150‑1 (MQ=255)
               cacCGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAACTGGTGTGAGTCGGTGATAAAACCGTTTTTTAACGCTATTGTCGCCAACTGCTTCAGAACGTAATAAGCGCCAATGCCCTGAATATTCAGGt                                                                                                                                        >  2:311163/1‑150 (MQ=255)
                acCGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAACTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTc                                                                                                                                       >  2:1637995/1‑150 (MQ=255)
                 ccGGAACCAAATCCAGTAGTGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAACTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCa                                                                                                                                      <  2:1228625/150‑1 (MQ=255)
                  cGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAAATGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGACAACTGCTTCAGAAAGGAATAAGCGCCATTGCACTGAATATTCAGGTCaa                                                                                                                                     >  1:908166/1‑150 (MQ=255)
                       ccAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAACTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAg                                                                                                                                >  2:791444/1‑150 (MQ=255)
                       ccAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAACTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGACCTGAATATTCAGGGCAAtggcg                                                                                                                                <  1:21032/150‑5 (MQ=255)
                        cAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGTGTTTGCAGGAACTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGc                                                                                                                               <  1:1265044/150‑1 (MQ=255)
                        cAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAACTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGc                                                                                                                               >  2:816578/1‑150 (MQ=255)
                        cAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAACTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGc                                                                                                                               >  2:1668389/1‑150 (MQ=255)
                        cAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAACTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGc                                                                                                                               <  1:1271588/150‑1 (MQ=255)
                            tCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAACTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGc                                                                                                                           <  2:256321/150‑1 (MQ=255)
                            tCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGGTTGCAGGAACTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGc                                                                                                                           <  1:815013/150‑1 (MQ=255)
                             ccAGTAGAGTGCTTTTTTTCGCGCAGCAGCAGGGTTTGCAGGAACTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCa                                                                                                                          <  2:1169542/150‑1 (MQ=255)
                             ccAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAACTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCa                                                                                                                          <  1:214592/150‑1 (MQ=255)
                             ccAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAACTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCa                                                                                                                          >  2:1273848/1‑150 (MQ=255)
                                 tAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAACTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTaa                                                                                                                      <  2:1166925/150‑1 (MQ=255)
                                 tAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAACTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTaa                                                                                                                      <  1:80572/150‑1 (MQ=255)
                                  agagTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAACTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAag                                                                                                                     >  1:310030/1‑150 (MQ=255)
                                  agagTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAACTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAag                                                                                                                     <  1:1818948/150‑1 (MQ=255)
                                     gTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAACTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGc                                                                                                                  >  1:1395113/1‑150 (MQ=255)
                                      tGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAACTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACCCCTTGCAGTAAGAGCt                                                                                                                 >  2:1793003/1‑150 (MQ=255)
                                      tGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAACTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCt                                                                                                                 <  2:1279331/150‑1 (MQ=255)
                                         aTTTTTTCGCGCAGCAGCAGGGTTTGCAGGAACTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGcc                                                                                                              <  2:19602/150‑1 (MQ=255)
                                         aTTTTTTCGCGCAGCAGCAGGGTTTGCAGGAACTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGcc                                                                                                              <  1:842335/150‑1 (MQ=255)
                                          ttttttCGCGCAGCAGCAGGGTTTGCAGGAACTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTGCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCa                                                                                                             <  2:1331168/150‑1 (MQ=255)
                                           tttttCGCGCAGCAGCAGGGTTTGCAGGAACTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCAt                                                                                                            <  1:1171053/150‑1 (MQ=255)
                                            ttttCTCGCAGCAGCAGGGTTTGCAGGAACTGGTGTGAGTCGGTGATAAAACCTTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAATAGCTGCCAtt                                                                                                           <  1:115239/150‑1 (MQ=255)
                                            ttttCGCGCAGCAGCAGGGTTTGCAGGAACTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCAtt                                                                                                           <  2:1328452/150‑1 (MQ=255)
                                            ttttCGCGCAGCAGCAGGGTTTGCAGGAACTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCAtt                                                                                                           >  1:1789361/1‑150 (MQ=255)
                                                cgcgCAGCAGCAGGGTTTGCAGGAACTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATc                                                                                                       >  2:1080963/1‑150 (MQ=255)
                                                 gcgcAGCAGCAGGGTTTGCAGGAACTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCg                                                                                                      <  1:89829/150‑1 (MQ=255)
                                                 gcgcAGCAGCAGGGTTTGCAGGAACTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCg                                                                                                      <  1:1816359/150‑1 (MQ=255)
                                                  cgcAGCAGCAGGGTTTGCAGGAACTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCgg                                                                                                     <  1:384197/150‑1 (MQ=255)
                                                       cagcagGGTTTGCAGGAACTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTcc                                                                                                <  2:15677/150‑1 (MQ=255)
                                                        agcagGGTTTGCAGGAACTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGCGCTGCCATTAATCGGACTCCa                                                                                               >  2:1584238/1‑150 (MQ=255)
                                                             ggTTTGCAGGAACTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATaaaa                                                                                          <  1:390929/150‑1 (MQ=255)
                                                             ggTTTGCAGGAACTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATaaaa                                                                                          <  1:670277/150‑1 (MQ=255)
                                                              gTTTGCAGGAACTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAat                                                                                         >  2:1404161/1‑150 (MQ=255)
                                                                ttGCAGGAACTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAatat                                                                                       <  1:1583803/150‑1 (MQ=255)
                                                                  gCAGGAACTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAatatat                                                                                     >  2:626195/1‑150 (MQ=255)
                                                                     ggAACTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATAtttt                                                                                  <  2:1743156/150‑1 (MQ=255)
                                                                       aaatGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCGAACTGCTTCAGAACGGAATAAGCGCCATTGCCGTGAATATTCATGTTTATACAGCTTGGAGTTAGAGCTGCCATTAATCGTACTCCATAAAATATATTTTaa                                                                                <  1:993647/147‑1 (MQ=255)
                                                                           ggTGTGAGTCGGTGATAAAACCGTTTTGTAAGGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTCAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAATATATATTTTAATCaa                                                                            >  1:981753/1‑150 (MQ=255)
                                                                           ggTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCaa                                                                            >  2:1710118/1‑150 (MQ=255)
                                                                               tgAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGt                                                                        <  1:798032/150‑1 (MQ=255)
                                                                               tgAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGt                                                                        <  1:726895/150‑1 (MQ=255)
                                                                                gAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGGATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTa                                                                       >  2:645030/1‑150 (MQ=255)
                                                                                  gacGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGACAAAACAGCTTGAAGTAAGAGATGCCATTAATAGTACTCCATAAAATATATTTTAAACAATTGTaaa                                                                     >  2:1342717/3‑149 (MQ=255)
                                                                                    cGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAAtgg                                                                   <  1:432798/150‑1 (MQ=255)
                                                                                     ggTGATAATACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAAtggt                                                                  <  1:289583/150‑1 (MQ=255)
                                                                                     ggTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGATTATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGAGTCCATAAAATATATTTTAATCAATTGTAAtggt                                                                  >  2:1465785/1‑150 (MQ=255)
                                                                                        gATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGa                                                               >  2:838150/1‑150 (MQ=255)
                                                                                         aTAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTAAATCAATTGTAATGGTGGAg                                                              >  1:116197/1‑150 (MQ=255)
                                                                                          tAAAACCGTTTTGTAACGATATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATAATCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGatc                                                             >  1:1751819/1‑148 (MQ=255)
                                                                                           aaaaCCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCg                                                            <  2:196452/150‑1 (MQ=255)
                                                                                            aaaCCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTACTCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCgt                                                           >  2:880303/1‑149 (MQ=255)
                                                                                              aCCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGAtt                                                         >  2:1483603/1‑150 (MQ=255)
                                                                                               ccGTTTTGTATCGCTATTGTCGCAAACTACTTCAGAACGGAATAAGCTCCATTGCCCTGATTATTAAGGTCAATACATCTTGCAGTAAGAGCTGCCATAAATCGGACTCCATAATATATATTTTAATCAATTGTAATGGTGGAGCGAGTa                                                        <  2:534350/150‑1 (MQ=255)
                                                                                               ccGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGATTa                                                        >  1:1666687/1‑150 (MQ=255)
                                                                                                  ttttGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGATTAAGt                                                     <  1:738185/150‑1 (MQ=255)
                                                                                                   tttGTAACGCTATTGTAGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCTGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGATTAAGtg                                                    <  1:136779/150‑1 (MQ=255)
                                                                                                     tGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGATTAAGtgat                                                  >  2:852140/1‑150 (MQ=255)
                                                                                                         aCGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGATGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGATTAAGTGATgaga                                              >  2:1333947/1‑150 (MQ=255)
                                                                                                             tATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGATTAAGTGATGAGAGGCg                                          <  2:349284/150‑1 (MQ=255)
                                                                                                                tGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGATTAAGTGATGAGAGGCGata                                       <  2:1330759/150‑1 (MQ=255)
                                                                                                                 gTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCTATAAAATATATTTTAATCAATTGTAATGGTGGAGCGATTAAGTGATGAGAGGCGatat                                      >  1:1452529/1‑150 (MQ=255)
                                                                                                                    gCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTATTCGGCCTCCATAAAATATATTTTAATCAATTGTAAGGGTGGAGCGATTTAGTGATGAGAGGCGATATGTc                                   >  2:407796/1‑150 (MQ=255)
                                                                                                                    gCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGATTAAGTGATGAGAGGCGATATGTc                                   >  1:1568645/1‑150 (MQ=255)
                                                                                                                      cAACTGCTTCATAACGGAATAAGCGCCATTGCCCTGAATATTCAGGGCAATAAAGCCTGCAGTAAGCGCTGCCATTAATCGGCCTCCATAAAATATATTTTAATCAATTGTAATGGGGGAGCGTTTAAGTGATGAGAGGCGATATGTCgg                                 <  2:950049/150‑1 (MQ=255)
                                                                                                                           gCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGATTAAGTGATGAGAGTCGATATGTCGGTAAAc                            >  2:1432919/1‑150 (MQ=255)
                                                                                                                            cTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGATTAAGTGATGAGAGGCGATATGTCGGTAAACa                           <  1:572957/150‑1 (MQ=255)
                                                                                                                            cTTCAGAACGGAATAAGCGCCATTGCCCTGAACATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGATTAAGTGATGATAGGCGATATGTCGGTAAACa                           >  1:1083317/1‑150 (MQ=255)
                                                                                                                              tCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGATTAAGTGATGAGAGGCGATATGTCGGTAAACAAt                         <  1:863012/150‑1 (MQ=255)
                                                                                                                               cAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGATTAAGTGATGAGAGGCGATATGTCGGTAAACAATg                        <  1:816578/150‑1 (MQ=255)
                                                                                                                                aGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGATTAAGTGATGAGAGGCGATATGTCGGTAAACAATGa                       >  1:595432/1‑150 (MQ=255)
                                                                                                                                    cGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGATTAAGTGATGAGAGGCGATATGTCGGTAAACAATGACAAt                   <  1:1643263/150‑1 (MQ=255)
                                                                                                                                     tgAATAAGCGCCATTGCCCTGAATATTCAGTTCAATAAAGATTGAAGTAAGAGCTGCCATTAATAGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGATTAAGTGATGAGAGGCGATATGTCGGTAAACAATGACAAtt                  <  1:244749/149‑1 (MQ=255)
                                                                                                                                      gAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGATTAAGTGATGAGAGGCGATATGTCGGTAAACAATGACAAttt                 >  2:546322/1‑150 (MQ=255)
                                                                                                                                           agctgtaTTGACCTGAATATTCAGGGCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGATTAAGTGATGAGAGGCGATATGTCGGTAAACAATGACAATTTGAGAc            >  2:1502059/7‑150 (MQ=255)
                                                                                                                                           aGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGATTAAGTGATGAGAGGCGATATGTCGGTAAACAATGACAATTTGAGAc            >  1:864062/1‑150 (MQ=255)
                                                                                                                                           aGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGATTAAGTGATGAGAGGCGATATGTCGGTAAACAATGACAATTTGAGAc            >  2:1457375/1‑150 (MQ=255)
                                                                                                                                           aGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGATTAAGTGATGAGAGGCGATATGTCGGTAAACAATGACAATTTGAGAc            >  1:1148524/1‑150 (MQ=255)
                                                                                                                                              tgtaTTGACCTGAATATTCAGGGCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGATTAAGTGATGAGAGGCGATATGTCGGTAAACAATGACAATTTGAGACaga         >  1:1321010/4‑150 (MQ=255)
                                                                                                                                              gcCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGATTAAGTGATGAGAGGCGATATGTCGGTAAACAATGACAATTTGAGACaga         >  2:1512388/1‑150 (MQ=255)
                                                                                                                                                 aTTGACCTGAATATTCAGGGCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGATATTAAGTGATTAGATGCGATATGTAGTTAAAAAATGACATTTTGAGACAGAGTg      >  1:1198260/1‑150 (MQ=255)
                                                                                                                                                     cccTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGATTAAGTGATGAGAGGCGATATGTCGGTAAACAATGACAATTTGAGACAGAGTGAAAg  >  2:1191341/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CGTGCGGCACGGCGACACCGGAACCAAATCCAGTAGAGTGCATTTTTTCGCGCAGCAGCAGGGTTTGCAGGAACTGGTGTGAGTCGGTGATAAAACCGTTTTGTAACGCTATTGTCGCCAACTGCTTCAGAACGGAATAAGCGCCATTGCCCTGAATATTCAGGTCAATACAGCTTGCAGTAAGAGCTGCCATTAATCGGACTCCATAAAATATATTTTAATCAATTGTAATGGTGGAGCGATTAAGTGATGAGAGGCGATATGTCGGTAAACAATGACAATTTGAGACAGAGTGAAAG  >  NC_000913/4092631‑4092929

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