Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,666,510 T→A 5.1% L359* (TTA→TAA)  csiE → stationary phase inducible protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,666,5100TA 5.1% 293.0 / 3.7 91L359* (TTA→TAA) csiEstationary phase inducible protein
Reads supporting (aligned to +/- strand):  ref base T (39/47);  new base A (2/3);  total (41/50)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TACGGCGTCCATCTTTCTGATGAGGAAAGTGGTCTGGTCGCGGTGATTTTCGGTGCCTGGCTAATGCAGGAAAACGATCTGCATGAAAAACAGATTATTCTTCTTACCGGGAATGATAGCGAGCGAGAAGCGCAGATTGAGCAGCAGTTACGCGAACTAACGTTACTGCCGCTCAACATTAAGCATATGTCGGTAAAGGCATTTTTGCAGACAGGCGCTCCGCGTGGCGCGGCACTGATTATTGCGCCTTATACCATGCCATTACCGCTCTTTTCACCACCGCTGATCTATACG  >  NC_000913/2666362‑2666655
                                                                                                                                                    |                                                                                                                                                 
tACGGCGTCCATCTTTCTGATGAGGAAAGTGGTCTGGTCGCGGTGATTTTCGGTGCCTGGCTAATGCAGGAAAACGATCTGCATGAAAAACAGATTATTCTTCTTACCGGGAATGATAGCGAGCGAGAAGCGCAGATTGAGCAGCAGtt                                                                                                                                                   <  7:298498/149‑1 (MQ=255)
       tCCATCTTTCTGATGAGGAAAGTGGTCTGGTCGCGGTGATTTTCGGTGCCTGGCTAATGCAGGAAAACGATCTGCATGAAAAACAGATTATTCTTCTTACCGGGAATGATAGCGAGCGAGAAGCGGAGATTGAGCAGCAGTTACGCGaa                                                                                                                                            <  6:349970/149‑1 (MQ=255)
       tCCATCTTTCTGATGAGGAAAGTGGTCTGGTCGCGGTGATTTTCGGTGCCTGGCTAATGCAGGAAAACGATCTGCATGAAAAACAGATTATTCTTCTTACCGGGAATGATAGCGAGCGAGAAGCGCAGATTGAGCAGCAGTTACGCGaa                                                                                                                                            <  1:97718/149‑1 (MQ=255)
       tCCATCTTTCTGATGAGGAAAGTGGTCTGGTCGCGGTGATTTTCGGTGCCTGGCTAATGCAGGAAAACGATCTGCATGAAAAACAGATTATTCTTCTTACCGGGAATGATAGCGAGCGAGAAGCGCAGATTGAGCAGCAGTTACGCGaa                                                                                                                                            <  1:336374/149‑1 (MQ=255)
        acATCTTTCTGATGAGGAAAGTGGTCTGGTCGCGGTGATTTTCGGTGCCTGGCTAATGCAGGAAAACGATCTGCATGAAAAACAGATTATTCTTCTTACCGGGAATGATAGCGAGCGAGAAGCGCAGATTGAGCAGCAGTTACGCGAAc                                                                                                                                           >  2:399844/2‑149 (MQ=255)
           tCTTTCTGATGAGGAAAGTGGTCTGGTCGCGGTGATTTTCGGTGCCTGGCTAATGCAGGAAAACGATCTGCATGAAAAACAGATTATTCTTCTTACCGGGAATGATAGCGAGCGAGAAGCGCAGATTGAGCAGCAGTTACGCGAACTaa                                                                                                                                        >  7:94957/1‑149 (MQ=255)
             tttCTGATGAGGAAAGTGGTCTGGTCGCGGTGATTTTCGGTGCCTGGCTAATGCAGGAAAACGATCTGCATGAAAAACAGATTATTCTTCTTACCGGGAATGATAGCGAGCGAGAAGCGCAGATTGAGCAGCAGTTACGCGAACTAACg                                                                                                                                      >  7:270491/1‑149 (MQ=255)
             tttCTGATGAGGAAAGTGGTCTGGTCGCGGTGATTTTCGGTGCCTGGCTAATGCAGGAAAACGATCTGCATGAAAAACAGATTATTCTTCTTACCGGGAATGATAGCGAGAGAGAAGCGCAGATTGAGCAGCAGTAACGCGAACTAACg                                                                                                                                      >  1:57910/1‑149 (MQ=255)
               tCTGATGAGGAAAGTGGTCTGGTCGCGGTGATTTTCGGTGCCTGGCTAATGCAGGAAAACGATCTGCATGAAAAACAGATTATTCTTCTTACCGGGAATGATAGCGAGCGAGAAGCGCAGATTGAGCAGCAGTTACGCGAACTAAc                                                                                                                                       <  8:257034/146‑1 (MQ=255)
               tCTGATGAGGAAAGTGGTCTGGTCGCGGTGATTTTCGGTGCCTGGCTAATGCAGGAAAACGATCTGCATGAAAAACAGATTATTCTTCTTACCGGGAATGATAGCGAGCGAGAAGCGCAGATTGAGCAGCAGTTACGCGAACTAACGt                                                                                                                                     <  7:328938/148‑1 (MQ=255)
                 tgatgaGGAAAGTGGTCTGGTCGCGGTGATTTTCGGTGCCTGGCTAATGCAGGAAAACGATCTGCATGAAAAACAGATTATTCTTCTTACCGGGAATGATAGCGAGCGAGAAGCGCAGATTGAGCAGCAGTTACGCGAACTAACGTTAc                                                                                                                                  >  3:183060/1‑149 (MQ=255)
                 tgatgaGGAAAGTGGTCTGGTCGCGGTGATTTTCGGTGCCTGGCTAATGCAGGAAAACGATCTGCATGAAAAACAGATTATTCTTCTTACCGGGAATGATAGCGAGCGAGAAGCGCAGATTGAGCAGCAGTTACGCGAACTAACGTTAc                                                                                                                                  >  2:161374/1‑149 (MQ=255)
                   atgaGGAAAGTGGTCTGGTGGCGGTGATTTTCGGGGCCTGGCTAATGCAGGAAAACGATCTGCATGAAAAACAGATTATTCTTGTTACCGGGAATGATAGCGAGCGAGTAGCGCAGATTGAGCAGCAGTTACGCGAACTAACGTTACt                                                                                                                                 <  8:290163/148‑1 (MQ=255)
                   atgaGGAAAGTGGTCTGGTCGCGGTGATTTTCGGTGCCTGGCTAATGCAGGAAAACGATCTGCATGAAAAACAGATTATTCTTCTTACCGGGAATGATAGCGAGCGAGAAGCGCAGATTGAGCAGCAGTTACGCGAACTAACGTTACt                                                                                                                                 <  1:115911/148‑1 (MQ=255)
                   atgaGGAAAGTGGTCTGGTCGCGGTGATTTTCGGTGCCTGGCTAATGCAGGAAAACGATCTGCATGAAAAACAGATTATTCTTCTTACCGGGAATGATAGCGAGCGAGAAGCGCAGATTGAGCAGCAGTTACGCGAACTAACGTTACTg                                                                                                                                >  1:236823/1‑149 (MQ=255)
                   aggcggAACGTGGTCTGGGGGCGGTGAGTTTCGGTGCCTGGCTAATGCAGGAAAACGATCTGCATGAAAAACAGATTATTCTTCTTACCGGGAATGATAGCGAGCGAGAAGCGCAGATTGAGCAGCAGTTACGCGAACTAACGTTACt                                                                                                                                 <  8:290167/144‑1 (MQ=255)
                     gaGGAAAGTGGTCTGGTCGCGGTGATTTTCGGTGCCTGGCTAATGCAGGAAAACGATCTGCATGAAAAACAGATTATTCTTCTTACCGGGAATGATAGCGAGCGAGAAGCGCAGATTGAGCAGCACTTACGCGAACTAACGTTACTGcc                                                                                                                              >  3:112746/1‑149 (MQ=255)
                      aGGAAAGGGGTCTGGGCGCGGTGATTTTCGGTGCCTGGCTAATGCAGGAAAACGATCTGCATGAAAAACAGATTATTCTTCTTACCGGGAATGATAGCGAGCGAGAAGCGCAGATTGAGCAGCAGTTACGCGAACTAACGTTACTGcc                                                                                                                              <  8:306546/148‑1 (MQ=255)
                       ggAAAGTGGTCTGGTCGCGGTGATTTTCGGTGCCTGGCTAATGCAGGAAAACGATCTGCATGAAAAACAGATTATTCTTCTTACCGGGAATGATAGCGAGCGAGAAGCGCAGATTGAGCAGCAGTTACGCGAACTAACGTTACTGCCGc                                                                                                                            <  5:224960/149‑1 (MQ=255)
                                tctggtcGCGGTGATTTTCGGTGCCTGGCTAATGCAGGAAAACGATCTGCATGAAAAACAGATTATTCTTCTTACCGGGAATGATAGCGAGCGAGAAGCGCAGATTGAGCGGCAGTAACGTTAGTTCGCGTAACTGCTGCTCaatctgc                                                                                                                   <  1:76640/149‑6 (MQ=37)
                                tctggtcGCGGTGATTTTCGGTGCCTGGCTAATGCAGGAAAACGATCTGCATGAAAAACAGATTATTCTTCTTACCGGGAATGATAGCGAGCGAGAAGCGCAGATTGAGCGGCAGTAACGTTAGTTCGCGTAACTGCTGCTCaatctgc                                                                                                                   <  7:189522/149‑6 (MQ=37)
                                    gtcGCGGTGATTTTCGGTGCCTGGCTAATGCAGGAAAACGATCTGCATGAAAAACAGATTATTCTTCTTACCGGGAATGATAGCGAGCGAGAAGCGCAGATTGAGCAGCAGTTACGCGAACTAACGTTACTGCCGCTCAACATTAAGc                                                                                                                <  4:183060/148‑1 (MQ=255)
                                    gtcGCGGTGATTTTCGGTGCCTGGCTAATGCAGGAAAACGATCTGCATGAAAAACAGATTATTCTTCTTACCGGGAATGATAGCGAGCGAGAAGCGCAGATTGAGCAGCAGTTACGCGAACTAACGTTACTGCCGCTCAACATTAAGc                                                                                                                <  1:95326/148‑1 (MQ=255)
                                       gcgGTGATTTTCGGTGCCTGGCTAATGCAGGAAAACGATCTGCATGAAAAACAGATTATTCTTCTTACCGGGAATGATAGCGAGCGAGAAGCGCAGATTGAGCAGCAGTTACGCGAACTAACGTTACTGCCGCTCAACATTAAGCatat                                                                                                            <  6:316724/149‑1 (MQ=255)
                                        cgGTGATTTTCGGTGCCTGGCTAATGCAGGAAAACGATCTGCATGAAAAACAGATTATTCTTCTTACCGGGAATGATAGCGAGCGAGAAGCGCAGATTGAGCAGCAGTTACGCGAACTAACGTTACTGCCGCTCAACATTAAGCATATg                                                                                                           >  3:369066/1‑149 (MQ=255)
                                         ggTGATTTTCGGTGCCTGGCTAATGCAGGAAAACGATCTGCATGAAAAACAGATTATTCTTCTTACCGGGAATGATAGCGAGCGAGAAGCGCAGATTGAGCAGCAGTTACGCGAACTAACGTTACTGCCGCTCAACATTAAGCATATGt                                                                                                          >  3:28829/1‑149 (MQ=255)
                                         ggTGATTTTCGGTGCCTGGCTAATGCAGGAAAACGATCTGCATGAAAAACAGATTATTCTTCTTACCGGGAATGATAGCGAGCGAGAAGCGCAGATTGAGCAGCAGTTACGCGAACTAACGTTACTGCCGCTCAACATTAAGCATATGt                                                                                                          >  8:115296/1‑149 (MQ=255)
                                          gTGATTTTCGGTGCCTGGCTAATGCAGGAAAACGATCTGCATGAAAAACAGATTATTCTTCTTACCGGGAATGATAGCGAGCGAGAAGCGCAGATTGAGCAGCAGTTACGCGAACTAACGTTACTGCCGCTCAACATTAAGCATATGTc                                                                                                         >  6:270290/1‑149 (MQ=255)
                                                  cGGTGCCTGGCTAATGCAGGAAAACGATCTGCATGAAAAACAGATTATTCTTCTTACCGGGAATGATAGCGAGCGAGAAGCGCAGATTGAGCAGCAGTTACGCGAACTAACGTTACTGCCGCTCAACATTAAGCATATGTCGGTAAAgg                                                                                                 >  8:278398/1‑149 (MQ=255)
                                                  cGGTGAATGGCTAATTCAGGAAAACGATCTGCATTAAAAACAGATTATTCTTGTTACCGGGAATGATTGCGAGCTAGTAGATCAGATTGAGCAGCAGTTTCGCGAACTAACGTTACTGCCGCTTAACATTAATCATATGTGGTTAAAg                                                                                                  <  2:57910/148‑1 (MQ=255)
                                                       ccTGGCTAATGCAGGAAAACGATCTGCATGAAAAACAGATTATTCTTCTTACCGGGAATGATAGCGAGCGAGAAGCGCAGATTGAGCAGCAGTTACGCGAACTAACGTTACTGCCGCTCAACATTAAGCATATGTCGGTAAAGGCAttt                                                                                            >  5:183507/1‑149 (MQ=255)
                                                        cTGGCTAATGCAGGAAAACGATCTGCATGAAAAACAGATTATTCTTCTTACCGGGAATGATAGCGAGCGAGAAGCGCAGATTGAGCGGCAGTAACGTTAGTTCGCGTAACTGCTGCTCaat                                                                                                                       <  3:122571/121‑2 (MQ=25)
                                                        cTGGCTAATGCAGGAAAACGATCTGCATGAAAAACAGATTATTCTTCTTACCGGGAATGATAGCGAGCGAGAAGCGCAGATTGAGCGGCAGTAACGTTAGTTCGCGTAACTGCTGCTCaat                                                                                                                       >  4:122571/1‑120 (MQ=25)
                                                        cTGGCTAATGCAGGAAAACGATCTGCATGAAAAACAGATTATTCTTCTTACCGGGAATGATAGCGAGCGAGAAGCGCAGATTGAGCAGCAGTTACGCGAACTAACGTTACTGCAGCTCAACATTAAGCATATGTCGGTAAAGGCAtttt                                                                                           >  4:135644/1‑149 (MQ=255)
                                                          ggCTAATGCAGGAAAACGATCTGCATGAAAAACAGATTATTCTTCTTACCGGGAATGATAGCGAGCGAGAAGCGCAGATTGAGCAGCAGTTACGCGAACTAACGTTACTGCCGCTCAACATTAAGCATATGTCGGTAAAGGCATTTTTg                                                                                         <  3:301719/149‑1 (MQ=255)
                                                           gCTAATGCAGGAAAACGATCTGCATGAAAAATAGATTATTCTTCTTACCGGGAATGATAGCGAGCGAGAAGCGCAGATTGAGCAGCAGTTACGCGAACTAACGTTACTGCCGCTCAACATTAAGCATATGTCGGTAAAGGCAttttt                                                                                          <  7:279423/147‑1 (MQ=255)
                                                            cTAATGCAGGAAAACGATCTGCATGAAAAACAGATTATTCTTCTTACCGGGAATGATAGCGAGCGAGAAGCGCAGATTGAGCAGCAGTTACGCGAACTAACGTTACTGCCGCTCAACATTAAGCATATGTCGGTAAAGGCATTTTTGCa                                                                                       <  7:278398/149‑1 (MQ=255)
                                                            cTAATGCAGGAAAACGATCTGCATGAAAAACAGATTATTCTTCTTACCGGGAATGATAGCGAGCGAGAAGCGCAGATTGAGCAGCAGTTACGCGAACTAACGTTACTGCCGCTCAACATTAAGCATATGTCGGTAAAGGCATTTTTGCa                                                                                       <  8:94957/149‑1 (MQ=255)
                                                                  cAGGAAAACGATCTGCATGAAAAACAGATTATTCTTCTTACCGGGAATGATAGCGAGCGAGAAGCGCAGATTGAGCAGCAGTTACGCGAACTAACGTTACTGCCGCTCAACATTAAGCATATGTCGGTAAAGGCATTTTTGCAGACAgg                                                                                 >  3:310683/1‑149 (MQ=255)
                                                                  cAGGAAAACGATCTGCATGAAAAACAGATTATTCTTCTTACCGGGAATGATAGCGAGCGAGAAGCGCAGATTGAGCAGCAGTTACGCGAACTAACGTTACTGCCGCTCAACATTAAGCATATGTCGGTAAAGGCATTTTTGCAGACAgg                                                                                 >  4:235921/1‑149 (MQ=255)
                                                                  cAGGAAAACGATCTGCATGAAAAACAGATTATTCTTCTTACCGGGAATGATAGCGAGCGAGAAGCGCAGATTGAGCAGCAGTTACGCGAACTAACGTTACTGCCGCTCAACATTAAGCATATGTCGGTAAAGGCATTTTTGCAGACAg                                                                                  <  8:377971/148‑1 (MQ=255)
                                                                    ggAAAACGATCTGCATGAAAAACAGATTATTCTTCTTACCGGGAATGATAGCGAGCGAGAAGCGCAGATTGAGCAGCAGTTACGCGAACTAACGTTACTGCCGCTCAACATTAAGCATATGTCGGTAAAGGCATTTTTGCAGACAGgcg                                                                               >  7:48885/1‑149 (MQ=255)
                                                                          cGATCTGCATGAAAAACAGATTATTCTTCTTACCGGGAATGATAGCGAGCGAGAAGCGCAGATTGAGCAGCAGTTACGCGAACTAACGTTACTGCCGCTCAACATTAAGCATATGTCGGTAAAGGCATTTTTGCAGACAGGCGCTCcgc                                                                         <  1:238030/149‑1 (MQ=255)
                                                                           gATCTGCATGAAAAACAGATTATTCTTCTTACCGGGAATGATAGCGAGCGAGAAGCGCAGATTGAGCAGCAGTTACGCGAACTAACGTTACTGCCGCTCAACATTAAGCATATGTCGGTAAAGGCATTTTTGCAGACAGGCGCTCcgcg                                                                        >  1:344970/1‑149 (MQ=255)
                                                                                 cATGAAAAACAGATTATTCTTCTTACCGGGAATGATAGCGAGCGAGAAGCGCAGATTGAGCAGCAGTTACGCGAACTAACGTTACTGCCGCTCAACATTAAGCATATGTCGGTAAAGGCAt                                                                                              >  6:281296/1‑121 (MQ=255)
                                                                                 cATGAAAAACAGATTATTCTTCTTACCGGGAATGATAGCGAGCGAGAAGCGCAGATTGAGCAGCAGTTACGCGAACTAACGTTACTGCCGCTCAACATTAAGCATATGTCGGTAAAGGCAt                                                                                              <  5:281296/121‑1 (MQ=255)
                                                                                 cATGAAAAACAGATTATTCTTCTTACCGGGAATGATAGCGAGCGAGAAGCGCAGATTGAGCAGCAGTTACGCGAACTAACGTTACTGCCGCTCAACATTAAGCATATGTCGGTAAAGGCATTTTTGCAGACAGGCGCTCCGCGTGgcgc                                                                  >  8:324984/1‑149 (MQ=255)
                                                                                  aTGAAAAACAGATTATTCTTCTTACCGGGAATGATAGCGAGCGAGCAGCGCAGATTGAGCAGCAGTTACGCGAACTAACGTTACTGCCGCTCAACATTAAGCATATGTCGGTAAAGGCATTTTTGCAGACAGGCGCTCCGCGTGGcgcg                                                                 <  6:247175/149‑1 (MQ=255)
                                                                                  aTGAAAAACAGATTATTCTTCTTACCGGGAATGATAGCGAGCGAGAAGCGCAGATTGAGCAGCAGTTACGCGAACTAACGTTACTGCCGCTCAACATTAAGCATATGTCGGTAAAGGCATTTTTGCAGACAg                                                                                  <  2:235938/132‑1 (MQ=255)
                                                                                  aTGAAAAACAGATTATTCTTCTTACCGGGAATGATAGCGAGCGAGAAGCGCAGATTGAGCAGCAGTTACGCGAACTAACGTTACTGCCGCTCAACATTAAGCATATGTCGGTAAAGGCATTTTTGCAGACAg                                                                                  >  1:235938/1‑132 (MQ=255)
                                                                                    gAAAAACAGATTATTCTTCTTACCGGGAATGATAGCGAGCGAGAAGCGCAGATTGAGCAGCAGTTACGCGAACTAACGTTACTGCCGCTCAACATTAAGCATATGTCGGTAAAGGCATTTTTGCAGACAGGCGCTCCGCGTGGCGCgg                                                                >  1:40724/1‑148 (MQ=255)
                                                                                    gAAAAACAGATTATTCTTCTTACCGGGAATGATAGCGAGCGAGAAGCGCAGATTGAGCAGCAGTTACGCGAACTAACGTTACTGCCGCTCAACATTAAGCATATGTCGGTAAAGGCATTTTTGCAGACAGGCGCTCCGCGTGGCGCgg                                                                >  4:167237/1‑148 (MQ=255)
                                                                                     aaaaaCAGATTATTCTTCTTACCGGGAATGATAGCGAGCGAGAAGCGCAGATTGAGCAGCAGTTACGCGAACTAACGTTACTGCCGCTCAACATTAAGCATATGTCGGTAAAGGCATTTTTGCAGACAGGCGCTCCGCGTGGCGCGGCa                                                              >  4:280207/1‑149 (MQ=255)
                                                                                      aacaCAGCTTCTTCTGCTTCCCGGGAATGCTCGCGAGCGAGACGCGCAGATTGATCAGCCTTTACGCGAACTAACGGTACTGCCGCTCAACATTACGCATATGTCGGTAAAGGAATTTTTGCAGACAGGCGCTCCGCGTGGCGCGGCAc                                                             <  8:358650/146‑1 (MQ=255)
                                                                                        aaCAGATTATTCTTCTTACCGGGAATGATAGCGAGCGAGAAGCGCAGATTGAGCAGCAGTTACGCGAACTAACGTTACTGCCGCTCAACATTAAGCATATGTCGGTAAAGGCATTTTTGCAGa                                                                                     <  6:394066/123‑1 (MQ=255)
                                                                                             attattCTTCTTACCGGGAATGATAGCGAGCGAGAAGCGCAGATTGAGCAGCAGTTACGCGAACTAACGTTACTGCCGCTCAACATTAAGCATATGTCGGTAAAGGCATTTTTGCAGACAGGCGCTCCGCGTGGCGCGGCACTGattat                                                      >  7:5632/1‑149 (MQ=255)
                                                                                             attattCTTCTTACCGGGAATGAGAGCGAGCGAGAAGCGCAGATTGAGCAGCAGTTACGCGAACTAACGTTACTGCCGCTCAACATTAATCATATGTCGGTAAAGGCATTTTTGCAGACAGGCGCTCCGCGTGGCGCGGCACTGattat                                                      <  4:106103/149‑1 (MQ=255)
                                                                                                attCTTCTTACCGGGAATGATAGCGAGCGAGAAGCGCAGATTGAGCAGCAGTTACGCGAACTAACGTTACTGCCGCTCAACATTAAGCATATGTCGGTAAAGGCATTTTTGCAGACAGGCGCTCCGCGTGGCGCGGCACTGATTATTgc                                                   >  6:15075/1‑149 (MQ=255)
                                                                                                 ttcttcTTACCGGGAATGATAGCGAGCGAGAAGCGCAGATTGAGCAGCAGTTACGCGAACTAACGTTACTGCCGCTCAACATTAAGCATATGTCGGTAAAGGCATTTTTGCAGACAGGCGCTCCGCGTGGCGCGGCACTGATTATTgcg                                                  >  5:131366/1‑149 (MQ=255)
                                                                                                 ttcttcTTACCGGGAATGATAGCGAGCGAGAAGCGCAGATTGAGCAGCAGTTACGCGAACTAACGTTACTGCCGCTCAACATTAAGCATATGTCGGTAAAGGCATTTTTGCAGACAGGCGCTCCGCGTGGCGCGGCACTGATTATTgc                                                   <  4:369066/148‑1 (MQ=255)
                                                                                                 ttcttcTTACCGGGAATGATAGCGAGCGAGAAGCGCAGATTGAGCAGCAGTTACGCGAACTAACGTTACTGCCGCTCAACATTAAGCATATGTCGGTAAAGGCATTTTTGCAGACAGGCGCTCCGCGTGGCGCGGCACTGATTATTgc                                                   <  6:183507/148‑1 (MQ=255)
                                                                                                   cttcttACCGGGAATGATAGCGAGCGAGAAGCGCAGATTGAGCAGCAGTTACGCGAACTAACGTTACTGCCGCTCAACATTAAGCATATGTCGGTAAAGGCATTTTTGCAGACAGGCGCTCCGCGTGGCGCGGCACTGATTATTgcgc                                                 <  4:28829/148‑1 (MQ=255)
                                                                                                    ttcttACCGGGAATGATAGCGAGCGAGAAGCGCAGATTGAGCAGCAGTTACGCGAACTAACGTTACTGCCGCTCAACATTAAGCATATGTCGGTAAAGGCATTTTTGCAGACAGGCGCTCCGCGTGGCGCGGCACTGATTATTGCGCCt                                               >  4:115664/1‑149 (MQ=255)
                                                                                                      cttACCGGGAATGATAGCGAGCGAGAAGCGCAGATTGAGCAGCAGTTACGCGAACTAACGTTACTGCCGCTCAACATTAAGCATATGTCGGTAAAGGCATTTTTGCAGACAGGCGCTCCGCGTGGCGCGGCACTGATTATTGCGcc                                                <  4:310683/146‑1 (MQ=255)
                                                                                                      cttACCGGGAATGATAGCGAGCGAGAAGCGCAGATTGAGCAGCAGTTACGCGAACTAACGTTACTGCCGCTCAACATTAAGCATATGTCGGTAAAGGCATTTTTGCAGACAGGCGCTCCGCGTGGCGCGGCACTGATTATTGCGCCtt                                              <  8:270491/148‑1 (MQ=255)
                                                                                                      cttACCGGGAATGATAGCGAGCGAGAAGCGCAGATTGAGCAGCAGTTACGCGAACTAACGTTACTGCCGCTCAACATTAAGCATATGTCGGTAAAGGCATTTTTGCAGACAGGCGCTCCGCGTGGCGCGGCACTGATTATTGCGCCtt                                              <  8:48885/148‑1 (MQ=255)
                                                                                                         aCCGGGAATGATAGCGAGCGAGAAGCGCAGATTGAGCAGCAGTTACGCGAACTAACGTTACTGCCGCTCAACATTAAGCATATGTCGGTAAAGGCATTTTTGCAGACAGGCGCTCCGCGTGGCGCGGCACTGATTATTGCGCCTta                                             <  2:236823/146‑1 (MQ=255)
                                                                                                          ccGGGAATGATAGCGAGCGAGAAGCGCAGATTGAGCAGCAGTTACGCGAACTAACGTTACTGCCGCTCAACATTAAGCATATGTCGGTAAAGGCATTTTTGCAGACAGGCGCTCCGCGTGGCGCGGCACTGATTATTGCGCCTTATAcc                                         <  5:270290/149‑1 (MQ=255)
                                                                                                           cGGGAATGATAGCGAGCGAGAAGCGCAGATTGAGCAGCAGTTACGCGAACTAACGTTACTGCCGCTCAACATTAAGCATATGTCGGTAAAGGCATTTTTGCAGACAGGCGCTCCGCGTGGCGCGGCACTGATTATTGCGCCTTATACCa                                        >  8:32285/1‑149 (MQ=255)
                                                                                                            gggAATGATAGCGAGCGAGAAGCGCAGATTGAGCAGCAGTTACGCGAACTAACGTTACTGCCGCTCAACATTAAGCATATGTCGGTAAAGGCATTTTTGCAGACAGGCGCTCCGCGTGGCGCGGCACTGATTATTGCGCCTTATACCAt                                       <  2:344970/149‑1 (MQ=255)
                                                                                                            gggAATGATAGCGAGCGAGAAGCGCAGATTGAGCAGCAGTTACGCGAACTAACGTTACTGCCGCTCAACATTAAGCATATGTCGGTAAAGGCATTTTTGCAGACAGGCGCTCCGCGTGGCGCGGCACTGATTATTGCGCCTTATACCAt                                       >  4:31193/1‑149 (MQ=255)
                                                                                                             ggAATGATAGCGAGCGAGAAGCGCAGATTGAGCACCAGTTACGCGAACTAACGTTACTGCCGCTCAACATTAAGCATATGTCGGTAAAGGCATTTTTGCAGACAGGCGCTCCGCGTGGCGCGGCACTGATTATTGCGCCTTATACCATg                                      >  4:274436/1‑149 (MQ=255)
                                                                                                                   aTAGCGAGCGAGAAGCGCAGATTGAGCAGCAGTTACGCGAACTAACGTTACTGCCGCTCAACATTAAGCATATGTCGGTAAAGGCATTTTTGCAGACAGGCGCTCCGCGTGGCGCGGCACTGATTATTGCGCCTTATACCATGCCATTa                                <  7:324984/149‑1 (MQ=255)
                                                                                                                    tAGCGAGCGAGAAGCGCAGATTGAGCAGCAGTTACGCGAACTAACGTTACTGCCGCTCAACATTAAGCATATGTCGGTAAAGGAATTTTTGCAGACAGGCGCTCCGCGTGGCGCGGCACTGATTATTGCGCCTTATACCATGCCATTAc                               >  5:149026/1‑149 (MQ=255)
                                                                                                                      gcgagcgaGAAGCGCAGATTGAGCAGCAGTTACGCGAACTAACGTTACTGCCGCTCAACATTAAGCATATGTCGGTAAAGGCATTTTTGCAGACAGGCGCTCCGCGTGGCGCGGCACTGATTATTGCGCCTTATACCATGCCATTACCg                             <  2:40724/149‑1 (MQ=255)
                                                                                                                         agcgagAAGCGCAGATTGAGCAGCAGTTACGCGAACTAACGTTACTGCCGCTCAACATTAAGCATATGTCGGTAAAGGCATTTTTGCAGACAGGCGCTCCGCGTGGCGCGGc                                                               <  7:243374/112‑1 (MQ=255)
                                                                                                                         agcgagAAGCGCAGATTGAGCAGCAGTTACGCGAACTAACGTTACTGCCGCTCAACATTAAGCATATGTCGGTAAAGGCATTTTTGCAGACAGGCGCTCCGCGTGGCGCGGc                                                               >  8:243374/1‑112 (MQ=255)
                                                                                                                           cgagAAGCGCAGATTGAGCAGCAGTTACGCGAACTAACGTTACTGCCGCTCAACATTAAGCATATGTCGGTAAAGGCATTTTTGCAGACAGGCGCTCCGCGTGGCGCGGCACTGATTATTGCGCCTTATACCATGCCATTACCGCTCtt                        <  3:31193/149‑1 (MQ=255)
                                                                                                                            gagaAGCGCAGATTGAGCAGCAGTTACGCGAACTAACGTTACTGCCGCTCAACATTAAGCATATGTTGGt                                                                                                      >  2:196066/1‑70 (MQ=255)
                                                                                                                            gagaAGCGCAGATTGAGCAGCAGTTACGCGAACTAACGTTACTGCCGCTCAACATTAAGCATATGTTGGt                                                                                                      <  1:196066/70‑1 (MQ=255)
                                                                                                                            gagaAGCGCAGATTGAGCAGCAGTTACGCGAACTAACGTTACTGCCGCTCAACATTAAGCATATGTCGGTAAAGGCATTTTTGCAGACAGGCGCTCCGCGTGGCGCGGCACTGATTATTGCGCCTTATACCATGCCATTACCGCTCttt                       >  4:131274/1‑149 (MQ=255)
                                                                                                                                aGCGCAGATTGAGCAGCAGTTACGCGAACTAACGTTACTGCCGCTCAACATTAAGCATATGTCGGTAAAGGCATTTTTGCAGACAGGCGCTCCGCGTGGCGCGGCACTGATTATTGCGCCTTATACCATGCCATTACCGCTCTTTTcac                   <  3:280207/149‑1 (MQ=255)
                                                                                                                                aGCGCAGATTGAGCAGCAGTTACGCGAACTAACGTTACTGCCGCTCAACATTAAGCATATGTCGGTAAAGGCATTTTTGCAGACAGGCGCTCCGCGTGGCGCGGCACTGATTATTGCGCCTTATACCATGCCATTACCGCTCTTTTcac                   <  1:399844/149‑1 (MQ=255)
                                                                                                                                   gcAGATTGAGGAGCAGTTACGAGAACTTACGTTACTGCATCTCAACATTAAGCATATGTCGGTAAAGGCATTTTTGCAGACAGGCGCGCCGAGTGGCGCGGCACTGATTATTTCGCCTTATACCATGGCATTACAGCTCTTTTcacca                 <  6:149026/148‑1 (MQ=255)
                                                                                                                                    cAGATTGAGCAGCAGTTACGCGAACTAACGTTACTGCCGCTCAACATTAAGCATATGTCGGTAAAGGCATTTTTGCAGACAGGCGCTCCGCGTGGCGCGGCACTGATTATTGCGCCTTATACCATGCCATTACCGCTCTTTTCaccacc               <  3:115664/149‑1 (MQ=255)
                                                                                                                                       aTTGAGCAGCAGTTACGCGAACTAACGTTACTGCCGCTCAACATTAAGCATATGTCGGTAAAGGCATTTTTGCAGACAGGCGCTCCGCGTGGCGCGGCACTGa                                                          >  3:129360/1‑103 (MQ=255)
                                                                                                                                       aTTGAGCAGCAGTTACGCGAACTAACGTTACTGCCGCTCAACATTAAGCATATGTCGGTAAAGGCATTTTTGCAGACAGGCGCTCCGCGTGGCGCGGCACTGa                                                          <  4:129360/103‑1 (MQ=255)
                                                                                                                                           agcagcAGTTACGCGAACTAACGTTACTGCCGCTCAACATTAAGCATATGTCGGTAAAGGCATTTTTGCAGACAGGCGCTCCGCGTGGCGCGGCACTGATTATTGCGCCTTATACCATGCCATTACCGCTCTTTTCACCACCGCTGAt         <  3:135644/148‑1 (MQ=255)
                                                                                                                                            gcagcaGTTACGCGAACTAACGTTACTGCCGCTCAACATTAAGCATATGTCGGTAAAGGCATTTTTGCAGACAGGCGCTCCGCGTGGCGCGGCACTGATTATTGCGCCTTATACCATGCCATTACCGCTCTTTTCACCACCGCTGATc        >  6:311969/1‑148 (MQ=255)
                                                                                                                                                 agTTACGCGAACTAACGTTACTGCCGCTCAACATTAAGCATATGTCGGTAAAGGCATTTTTGCAGACAGGCGCTCCGCGTGGCGCGGCACTGATTATTGCGCCTTATACCATGCCATTACCGCTCTTTTc                     >  2:309523/1‑130 (MQ=255)
                                                                                                                                                 agTTACGCGAACTAACGTTACTGCCGCTCAACATTAAGCATATGTCGGTAAAGGCATTTTTGCAGACAGGCGCTCCGCGTGGCGCGGCACTGATTATTGCGCCTTATACCATGCCATTACCGCTCTTTTc                     <  1:309523/130‑1 (MQ=255)
                                                                                                                                                  gTTACGCGAACTAACGTTACTGCCGCTCAACATTAAGCATATGTCGGTAAAGGCATTTTTGCAGACAGGCGCTCCGCGTGGCGCGGCACTGATTATTGCGCCTTATACCATGCCATTACCGCTCTTTTCACCACCGCTGATCTATACg  >  8:252314/1‑148 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
TACGGCGTCCATCTTTCTGATGAGGAAAGTGGTCTGGTCGCGGTGATTTTCGGTGCCTGGCTAATGCAGGAAAACGATCTGCATGAAAAACAGATTATTCTTCTTACCGGGAATGATAGCGAGCGAGAAGCGCAGATTGAGCAGCAGTTACGCGAACTAACGTTACTGCCGCTCAACATTAAGCATATGTCGGTAAAGGCATTTTTGCAGACAGGCGCTCCGCGTGGCGCGGCACTGATTATTGCGCCTTATACCATGCCATTACCGCTCTTTTCACCACCGCTGATCTATACG  >  NC_000913/2666362‑2666655

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