Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,236,226 T→C 5.7% intergenic (+78/‑36) yjbE → / → yjbF extracellular polysaccharide production threonine‑rich protein/extracellular polysaccharide production lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,236,2260TC 5.7% 938.5 / 30.4 279intergenic (+78/‑36)yjbE/yjbFextracellular polysaccharide production threonine‑rich protein/extracellular polysaccharide production lipoprotein
Reads supporting (aligned to +/- strand):  ref base T (104/159);  new base C (6/10);  total (110/169)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.49e-01

Only 100 of 283 total aligned reads displayed.

CTTGCGGCAACCGGTGGCGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGGTTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGTTATTACAGGCCTGTTCAGCCACGACTAAAGAGCTGGGCAATTCACTGTGGGACAGTCTGTTCGGCACGCCAGGCGTACA  >  NC_000913/4236077‑4236374
                                                                                                                                                     |                                                                                                                                                    
    cGGCAACCGGTGGCGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCgg                                                                                                                                                  >  2:2094153/1‑150 (MQ=255)
         aCCGGTGGCGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTACGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTgg                                                                                                                                             <  1:4001327/150‑1 (MQ=255)
            ggTGGCGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGGTTa                                                                                                                                          >  1:1799218/1‑150 (MQ=255)
             gTGGCGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACCGGAGTGTGGTTAt                                                                                                                                         <  1:2955435/150‑1 (MQ=255)
             gTGGCGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACCGGAGTGTGGTTAt                                                                                                                                         <  2:1351700/150‑1 (MQ=255)
             gTGGCGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACCGGAGTGTGGTTAt                                                                                                                                         <  1:1891930/150‑1 (MQ=255)
              tGGCGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACGCTCCGGTGTGGTTAtt                                                                                                                                        >  1:4242726/1‑150 (MQ=255)
               ggCGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGGTTATTc                                                                                                                                       <  2:3906332/150‑1 (MQ=255)
               ggCGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGGTTATTc                                                                                                                                       <  1:3263748/150‑1 (MQ=255)
                gCGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGAATGAAGGGATTAACCATAACCACACTCCGGTGTGGTTATTCt                                                                                                                                      >  1:1244710/1‑150 (MQ=255)
                         gTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGGTTATTCTGCCCCTCTg                                                                                                                             >  1:912954/1‑150 (MQ=255)
                         gTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGGTTATTCTGCCCCTCTg                                                                                                                             <  1:1469009/150‑1 (MQ=255)
                         gTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGGTTATTCTGCCCCTCTg                                                                                                                             >  1:3363263/1‑150 (MQ=255)
                           tCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGGTTATTCTGCCCctct                                                                                                                              >  1:3939019/1‑147 (MQ=255)
                           tCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGGTTATTCTGCCCctct                                                                                                                              <  2:3939025/147‑1 (MQ=255)
                           tCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGGTTATTCTGCCCCTCTGga                                                                                                                           >  1:2705162/1‑150 (MQ=255)
                              aaTACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGGTTATTCTGCCCCTCTGGagaa                                                                                                                        <  1:706816/150‑1 (MQ=255)
                              aaTACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGGTTATTCTGCCCCTCTGGagaa                                                                                                                        <  1:1080215/150‑1 (MQ=255)
                              aaTACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGGTTATTCTGCCCCTCTGGagaa                                                                                                                        <  1:1080382/150‑1 (MQ=255)
                              aaTACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGGTTATTCTGCCCCTCTGGagaa                                                                                                                        <  1:2431246/150‑1 (MQ=255)
                                tACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGGTTATTCTGCCCCTCTGGagaaga                                                                                                                      <  2:2281137/150‑1 (MQ=255)
                                  ccGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGGTTATTCTGCCCCTCTGGAGAAGAGt                                                                                                                    >  2:1458640/1‑150 (MQ=255)
                                     ggACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGGTTATTCTGCCCCTCTGGAGAAGAGTCGt                                                                                                                 >  1:2083465/1‑150 (MQ=255)
                                        ccaccaCAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGGTTATTCTGCCCCTCTGGAGAAGAGTCGTGaa                                                                                                              <  2:3688019/150‑1 (MQ=255)
                                        ccaccaCAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGGTTATTCTGCCCCTCTGGAGAAGAGTCGTGaa                                                                                                              <  2:2715683/150‑1 (MQ=255)
                                        ccaccaCAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGGTTATTCTGCCCCTCTGGAGAAGAGTCGTGaa                                                                                                              <  1:2249324/150‑1 (MQ=255)
                                        ccaccaCAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGGTTATTCTGCCCCTCTGGAGAAGAGTCGTGaa                                                                                                              <  2:1838361/150‑1 (MQ=255)
                                            cacaACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGTAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGGTTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGa                                                                                                          >  2:4250429/1‑150 (MQ=255)
                                            cacaACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGGTTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGa                                                                                                          <  1:236361/150‑1 (MQ=255)
                                             acaACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGGTTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGAc                                                                                                         <  2:3481014/150‑1 (MQ=255)
                                                aCCACGACCACCAGTACCCAGTAATAAAGTATCTATCCCCAAAATAATTCGAGTCAATGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGGTTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTg                                                                                                      <  2:2506782/150‑1 (MQ=255)
                                                 ccACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGGTTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGc                                                                                                     >  1:1706485/1‑150 (MQ=255)
                                                   aCGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGGTTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCAc                                                                                                   >  2:2567433/1‑150 (MQ=255)
                                                          ccaGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGGTTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCt                                                                                            >  2:1029375/1‑150 (MQ=255)
                                                           caGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGTTTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCtt                                                                                           <  1:1952142/150‑1 (MQ=255)
                                                           caGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGGTTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCtt                                                                                           <  2:4186361/150‑1 (MQ=255)
                                                             gTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGGTTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTAt                                                                                         <  2:1683945/150‑1 (MQ=255)
                                                              tACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCTTCTGTGGCTTGATGTATGAAGGTTTTAACCATAACCACACTCCGGTGTGGTTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATc                                                                                        <  1:3485608/150‑1 (MQ=255)
                                                                 ccAGTAATAAAGTATGTATCCCCAAATTAATTCGAGTCATTGCATCTGTGGCTAGAAGAATGAAGTGATTAACCATAACCACACTCCGGTGTGGTTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCTCTCATTCTTATCTGc                                                                                     <  2:3265406/150‑1 (MQ=255)
                                                                   aGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGGTTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCt                                                                                   >  2:855828/1‑150 (MQ=255)
                                                                   aGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGGTTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCt                                                                                   <  1:2142609/150‑1 (MQ=255)
                                                                    gTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGGTTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTg                                                                                  <  2:2789517/150‑1 (MQ=255)
                                                                       ataaAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAAAAACACTCCGGTGTGGTTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGtta                                                                               >  1:267180/1‑150 (MQ=255)
                                                                        taaAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGGTTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGttat                                                                              <  1:1439008/150‑1 (MQ=255)
                                                                        taaAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGGTTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGttat                                                                              >  2:3017624/1‑150 (MQ=255)
                                                                         aaaGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGGTTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGttatt                                                                             >  2:1636074/1‑150 (MQ=255)
                                                                            gtatgtatCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGGTTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGTTATTACa                                                                          <  1:497633/150‑1 (MQ=255)
                                                                              atgtatCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGGTTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGTTATTACAgg                                                                        >  2:3027607/1‑150 (MQ=255)
                                                                                     cccAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGGTTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGTTATTACAGGCCTGTTc                                                                 >  2:553741/1‑150 (MQ=255)
                                                                                     cccAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGGTTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGTTATTACAGGCCTGTTc                                                                 >  2:553926/1‑150 (MQ=255)
                                                                                      ccAAAATAATTCGAGTCATTGCATCTGTGTCTAGAAGTATGAAGGGATTAACCATAACCACTCTCCGGTGTTGTTATTCTGCCCCTCTGTAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGTTATTACAGGCCTGTTCa                                                                <  1:3116068/150‑1 (MQ=255)
                                                                                      ccAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGGTTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGTTATTACAGGCCTGTTCa                                                                <  1:3251753/150‑1 (MQ=255)
                                                                                      ccAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGGTTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGTTATTACAGGCCTGTTCa                                                                <  1:2617765/150‑1 (MQ=255)
                                                                                       cAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGGTTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGTTATTACAGGCCTGTTCAg                                                               >  2:3080001/1‑150 (MQ=255)
                                                                                       cAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGGTTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGTTATTACAGGCCTGTTCAg                                                               >  2:1945390/1‑150 (MQ=255)
                                                                                       cAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGGTTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGTTATTACAGGCCTGTTCAg                                                               >  2:675421/1‑150 (MQ=255)
                                                                                       cAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGGTTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGTTATTACAGGCCTGTTCAg                                                               <  2:478627/150‑1 (MQ=255)
                                                                                       cAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGGTTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGTTATTACAGGCCTGTTCAg                                                               >  1:3722057/1‑150 (MQ=255)
                                                                                        aaaaTAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGGTTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGTTATTACAGGCCTGTTCAGc                                                              >  2:2449847/1‑150 (MQ=255)
                                                                                          aatccttCGAGTCATCGCATCTGTGGCTAGATGTATGAAGGGATTAACCATAACCACACTCCGGTGTGGTTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGTTATTACAGGCCTGTTCAGCCa                                                            <  1:3642127/145‑1 (MQ=255)
                                                                                          aataatTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGGTTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGTTATTACAGGCCTGTTCAGCCa                                                            <  1:2222773/150‑1 (MQ=255)
                                                                                              atTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGGTTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGTTATTACAGGCCTGTTCAGCCACGAc                                                        <  1:1580315/150‑1 (MQ=255)
                                                                                              atTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGGTTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGTTATTACAGGCCTGTTCAGCCACGAc                                                        >  1:566891/1‑150 (MQ=255)
                                                                                                   aGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGGTTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGTTATTACAGGCCTGTTCAGCCACGACTAAag                                                   >  2:2835776/1‑150 (MQ=255)
                                                                                                    gTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGGTTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGTTATTACAGGCCTGTTCAGCCACGACTAAaga                                                  >  1:692082/1‑150 (MQ=255)
                                                                                                    gTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGGTTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGTTATTACAGGCCTGTTCAGCCACGACTAAaga                                                  <  2:2650118/150‑1 (MQ=255)
                                                                                                     tCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGGTTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGTTATTACAGGCCTGTTCAGCCACGACTAAagag                                                 <  2:337285/150‑1 (MQ=255)
                                                                                                     tCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGGTTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGTTATTACAGGCCTGTTCAGCCACGACTAAagag                                                 <  2:702672/150‑1 (MQ=255)
                                                                                                        ttGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGGTTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGTTATTACAGGCCTGTTCAGCCACGACTAAAGAGCTg                                              <  2:178165/150‑1 (MQ=255)
                                                                                                         tGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGGTTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGTTATTACAGGCCTGTTCAGCCACGACTAAAGAGCTgg                                             <  2:3990050/150‑1 (MQ=255)
                                                                                                          gCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGGTTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGTTATTACAGGCCTGTTCAGCCACGACTAAAGAGCTggg                                            >  1:1203373/1‑150 (MQ=255)
                                                                                                                     gaGAAGTATGACGGGATGAACCATAACCACACTCCGGTGTGGTTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGTTATTACAGGCCTGTTCAGCCACGACTAAAGAGCTGGGCAATTCACtgt                                 <  1:1892185/149‑1 (MQ=255)
                                                                                                                      aGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGGTTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGTTATTACAGGCCTGTTCAGCCACGACTAAAGAGCTGGGCAATTCACtgtg                                <  1:752873/150‑1 (MQ=255)
                                                                                                                           tATGAAGGGATTAACCATAACCACACTCCGGTGTGGTTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCCTTCTTATCTGCCTGTTATTACAGGCCTGTTCAGCCACGACTAAAGAGCTGGGCAATTCACTGTGGGACa                           >  2:1554146/1‑150 (MQ=255)
                                                                                                                           tATGAAGGGATTAACCATAACCACACTCCGGTGTGGTTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGTTATTACAGGCCTGTTCAGCCACGACTAAAGAGCTGGGCAATTCACTGTGGGACa                           >  1:2875789/1‑150 (MQ=255)
                                                                                                                           tATGAAGGGATTAACCATAACCACACTCCGGTGTGGTTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGTTATTACAGGCCTGTTCAGCCACGACTAAAGAGCTGGGCAATTCACTGTGGGACa                           <  2:3811174/150‑1 (MQ=255)
                                                                                                                             tGAAGGGATTAACCATAACCACACTCCGGTGTGGTTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGTTATTACAGGCCTGTTCAGCCACGACTAAAGAGCTGGGCAATTCACTGTGGGACAGt                         <  1:1750225/150‑1 (MQ=255)
                                                                                                                              gAAGGGATTAACCATAACCACACTCCGGTGTGGTTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGTTATTACAGGCCTGTTCAGCCACGACTAAAGAGCTGGGCAATTCACTGTGGGACAGTc                        >  2:2210324/1‑150 (MQ=255)
                                                                                                                              gAAGGGATTAACCATAACCACACTCCGGTGTGGTTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGTTATTACAGGCCTGTTCAGCCACGACTAAAGAGCTGGGCAATTCACTGTGGGACAGTc                        >  1:2306547/1‑150 (MQ=255)
                                                                                                                               aaGGGATTAACCATAACCACACTCCGGTGTGGTTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGTTATTACAGGCCTGTTCAGCCACGACTAAAGAGCTGGGCAATTCACTGTGGGACAGTCt                       >  2:2066324/1‑150 (MQ=255)
                                                                                                                               aaGGGATTAACCATAACCACACTCCGGTGTGGTTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGTTATTACAGGCCTGTTCAGCCACGACTAAAGAGCTGGGCAATTCACTGTGGGACAGTCt                       <  1:4165295/150‑1 (MQ=255)
                                                                                                                                  ggATTAACCATAACCACACTCCGGTGTGGTTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGTTATTACAGGCCTGTTCAGCCACGACTAAAGAGCTGGGCAATTCACTGTGGGACAGTCTGtt                    >  1:3575541/1‑150 (MQ=255)
                                                                                                                                   gATTAACCATAACCACACTCCGGTGTGGTTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGTTATTACAGGCCTGTTCAGCCACGACTAAAGAGCTGGGCAATTCACTGTGGGACAGTCTGTTc                   <  2:2704253/150‑1 (MQ=255)
                                                                                                                                     ttAACCATAACCACACTCCGGTGTGGTTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGTTATTACAGGCCTGTTCAGCCACGACTAAAGAGCTGGGCAATTCACTGTGGGACAGTCTGTTCgg                 >  2:1152030/1‑150 (MQ=255)
                                                                                                                                     ttAACCATAACCACACTCCGGTGTGGTTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCGGCACTCAGTCTGATCTGCCTGTTATTACAGGCCTGTTCAGCCACGACTAAAGAGCTGGGCAATTCACTGTGGGACAGTCTGTTCgg                 >  2:926786/1‑150 (MQ=255)
                                                                                                                                      gcagaataaccaCACCGGAGTGTGGTTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGTTATTACAGGCCTGTTCAGCCACGACTAAAGAGCTGGGCAATTCACTGTGGGACAGTCTGTTCGGc                >  2:3358318/6‑150 (MQ=255)
                                                                                                                                          cataaccaCACTCCGGTGTGGTTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGTTATTACAGGCCTGTTCAGCCACGACTAAAGAGCTGGGCAATTCACTGTGGGACAGTCTGTTCGGCACGc            <  2:3755032/150‑1 (MQ=255)
                                                                                                                                          aataaccaCACCGGAGTGTGGTTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGTTATTACAGGCCTGTTCAGCCACGACTAAAGAGCTGGGCAATTCACTGTGGGACAGTCTGTTCGGCACGc            >  2:2894352/2‑150 (MQ=255)
                                                                                                                                          aataaccaCACCGGAGTGTGGTTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGTTATTAAAGGCCTGTTCAGCCACGACTAAAGAGCTGGGCAATTCACTGTGGGACAGTCTGTTCGGCACGc            >  2:1458838/2‑150 (MQ=255)
                                                                                                                                            taaccaCACTCCGGTGTGGTTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGTTATTACAGGCCTGTTCAGCCACGACTAAAGAGCTGGGCAATTCACTGTGGGACAGTCTGTTCGGCACGCCa          >  1:4031379/1‑150 (MQ=255)
                                                                                                                                            taaccaCACTCCGGTGTGGTTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGTTATTACAGGCCTGTTCAGCCACGACTAAAGAGCTGGGCAATTCACTGTGGGACAGTCTGTTCGGCACGCCa          <  1:2548109/150‑1 (MQ=255)
                                                                                                                                            taaccaCACTCCGGTGTGGTTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGTTATTACAGGCCTGTTCAGCCACGACTAAAGAGCTGGGCAATTCACTGTGGGACAGTCTGTTCGGCACGCCa          <  1:1397940/150‑1 (MQ=255)
                                                                                                                                             aaccaCACTCCGGTGTGGTTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGTTATTACAGGCCTGTTCAGCCACGACTAAAGAGCTGGGCAATTCACTGTGGGACAGTCTGTTCGGCACGCCAg         >  1:3859098/1‑150 (MQ=255)
                                                                                                                                             aaccaCACTCCGGTGTGGTTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGTTATTACAGGCCTGTTCAGCCACGACTAAAGAGCTGGGCAATTCACTGTGGGACAGTCTGTTCGGCACGCCAg         >  1:3736380/1‑150 (MQ=255)
                                                                                                                                               ccaCTCTCCGGTGTGGTTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTAATTCTTATCTGCCTGTTATTACAGGCCTGTTCAGCCACGACTAAAGTGCTTGGCAATTCACTGTGGGACAGTCTGTTCGGCACGCCAGGc       <  2:206375/150‑1 (MQ=255)
                                                                                                                                               ccaCACTCCGGTGTGGTTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGTTATTACAGGCCTGTTCAGCCACGACTAAAGAGCTGGGCAATTCACTGTGGGACAGTCTGTTCGGCACGCCAGGc       <  2:3226444/150‑1 (MQ=255)
                                                                                                                                               ccaCACTCCGGTGTGGTTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGTTATTACAGGCCTGTTCAGCCACGACTAAAGAGCTGGGCAATTCACTGTGGGACAGTCTGTTCGGCACGCCAGGc       <  2:2390778/150‑1 (MQ=255)
                                                                                                                                               ccaCACTCCGGTGTGGTTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGTTATTACAGGCCTGTTCAGCCACGACTAAAGAGCTGGGCAATTCACTGTGGGACAGTCTGTTCGGCACGCCAGGc       <  2:4089375/150‑1 (MQ=255)
                                                                                                                                                 acacTCCGGTGTGGTTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGTTATTACAGGCCTGTTCAGCCACGACTAAAGAGCTGGGCAATTCACTGTGGGACAGTCTGTTCGGCACGCCAGGCGt     >  1:207244/1‑150 (MQ=255)
                                                                                                                                                  cacTCCGGTGTGGTTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGTTATTACAGGCCTGTTCAGCCACGACTAAAGAGCTGGGCAATTCACTGTGGGACAGTCTGTTCGGCACGCCAGGCGTa    <  2:18760/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CTTGCGGCAACCGGTGGCGGTGATGGTTCTAATACCGGGACCACCACAACCACGACCACCAGTACCCAGTAATAAAGTATGTATCCCCAAAATAATTCGAGTCATTGCATCTGTGGCTAGAAGTATGAAGGGATTAACCATAACCACACTCCGGTGTGGTTATTCTGCCCCTCTGGAGAAGAGTCGTGAAGCGACCTGCACTCATTCTTATCTGCCTGTTATTACAGGCCTGTTCAGCCACGACTAAAGAGCTGGGCAATTCACTGTGGGACAGTCTGTTCGGCACGCCAGGCGTACA  >  NC_000913/4236077‑4236374

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