Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,613,962 A→C 7.6% E99A (GAA→GCA) ‡ nikA → nickel/heme ABC transporter periplasmic binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,613,9620AC 7.6% 503.7 / 27.6 145E99A (GAA→GCA) ‡nikAnickel/heme ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base A (74/60);  new base C (9/2);  total (83/62)
Fisher's exact test for biased strand distribution p-value = 1.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

Only 100 of 147 total aligned reads displayed.

TCGCCCAGAGCATGGTTTATGAACCATTGGTGAAATATCAGGCAGACGGTTCGGTGATCCCGTGGCTGGCAAAAAGCTGGACTCATTCAGAAGATGGTAAAACCTGGACCTTCACCCTGCGTGATGACGTGAAATTCTCCAACGGTGAACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGATTGTTGATGTTAAAGCACTCAGTAAAACAGAGCTGCAAATTACCCTGAAAAGCGCCTACT  >  NC_000913/3613815‑3614111
                                                                                                                                                   |                                                                                                                                                     
tCGCCCAGAGCATGGTTTATGAACCATTGGTGAAATATCAGGCAGACGGTTCGGTGATCCCGTGGCTGGCAAAAAGCTGGACTCATTCAGAAGATGGTAAAACCTGGACCTTCACCCTGCGTGATGACGTGAAATTCTCCAACGGTGAAc                                                                                                                                                     >  1:2665832/1‑150 (MQ=255)
   cccAGAGCATGGTTTATGAACCATTGGTGAAATATCAGGCAGACGGTTCGGTGATCCCGTGGCTGGCAAAAAGCTGGACTCATTCAGAAGATGGTAAAACCTGGACCTTCACCCTGCGTGATGACGTGAAATTCTCCAACGGTGAACCGt                                                                                                                                                  >  1:2612664/1‑150 (MQ=255)
   cccAGAGCATGGTTTATGAACCATTGGTGAAATATCAGGCAGACGGTTCGGTGATCCCGTGGCTGGCAAAAAGCTGGACTCATTCAGAAGATGGTAAAACCTGGACCTTCACCCTGCGTGATGACGTGAAATTCTCCAACGGTGAACCGt                                                                                                                                                  <  2:2562058/150‑1 (MQ=255)
      agagCATGGTTTATGAACCATTGGTGAAATATCAGGCAGACGGTTCGGTGATCCCGTGGCTGGCAAAAAGCTGGACTCATTCAGAAGATGGTAAAACCTGGACCTTCACCCTGCGTGATGACGTGAAATTCTCCAACGGTTCACCGttgg                                                                                                                                               <  1:1441478/150‑3 (MQ=255)
                ttATGAACCATTGGTGAAATATCAGGCAGACGGTTCGGTGATCCCGTGGCTGGCAAAAAGCTGGACTCATTCAGAAGATGGTAAAACCTGGACCTTCACCCTGCGTGATGACGTGAAATTCTCCAACGGTGAACCGTTCGATGCCGAggc                                                                                                                                     >  1:870947/1‑150 (MQ=255)
                 tATGAACCATTGGTGAAATATCAGGCAGACGGTTCGGTGATCCCGTGGCTGGCAAAAAGCTGGACTCATTCAGAAGATGGTAAAACCTGGACCTTCACCCTGCGTGATGACGTGAAATTCTCCAACGGTGAACCGTTCGATGCCGAggcg                                                                                                                                    >  1:1371272/1‑150 (MQ=255)
                     aaCCATTGGTGAAATATCAGGCAGACGGTTCGGTGATCCCGTGGCTGGCAAAAAGCTGGACTCATTCAGAAGATGGTAAAACCTGGACCTTCACCCTGCGTGATGACGTGAAATTCTCCAACGGTGAACCGTTCGATGCCGAggcggcgg                                                                                                                                <  1:2250174/150‑1 (MQ=255)
                     aaCCATTGGTGAAATATCAGGCAGACGGTTCGGTGATCCCGTGGCTGGCAAAAAGCTGGACTCATTCAGAAGATGGTAAAACCTGGACCTTCACCCTGCGTGATGACGTGAAATTCTCCAACGGTGAACCGTTCGATGCCGAggcggcgg                                                                                                                                <  1:273408/150‑1 (MQ=255)
                         aTTGGTGAAATATCAGGCAGACGGTTCGGTGATCCCGTGGCTGGCAAAAAGCTGGACTCATTCAGAAGATGGTAAAACCTGGACCTTCACCCTGCGTGATGACGTGAAATTCTCCAACGGTGAACCGTTCGATGCCGAGGCGGCGGCAGa                                                                                                                            >  2:2498160/1‑150 (MQ=255)
                          ttGGTGAAATATCAGGCAGACGGTTCGGTGATCCCGTGGCTGGCAAAAAGCTGGACTCATTCAGAAGATGGTAAAACCTGGACCTTCACCCTGCGTGATGACGTGAAATTCTCCAACGGTGAACCGTTCGATGCCGAGGCGGCGGCAGaa                                                                                                                           <  1:2005246/150‑1 (MQ=255)
                            ggTGAAATATCAGGCAGACGGTTCGGTGATCCCGTGGCTGGCAAAAAGCTGGACTCATTCAGAAGATGGTAAAACCTGGACCTTCACCCTGCGTGATGACGTGAAATTCTCCAACGGTGAACCGTTCGATGCCGAGGCGGCGGCAGaaaa                                                                                                                         <  1:1147944/150‑1 (MQ=255)
                                    atCAGGCAGACGGTTCGGTGATCCCGTGGCTGGCAAAAAGCTGGACTCATTCAGAAGATGGTAAAACATGGACCTTCACCCTGCGTGATGACGTGAAATTCTCCAACGGTGAACCGTTCGATGCCGAGGCGGCGGCAGACAACTTCcgcg                                                                                                                 >  1:2264849/1‑150 (MQ=255)
                                      cAGGCAGACGGTTCGGTGATCCCGTGGCTGGCAAAAAGCTGGACTCATTCAGAAGATGGTAAAACCTGGACCTTCACCCTGCGTGATGACGTGAAATTCTCCAACGGTGAACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCa                                                                                                               <  2:1804906/150‑1 (MQ=255)
                                         gCAGACGTTTCGGTGATCCCGTGGCTGGCAAAATGCTGGACTCATTCAGAAGTTGGTAAAACCTGGACCTTCACCCTGCTTGATGACGTTAAATTCTCCAACGGTTAACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTg                                                                                                            <  1:2074751/150‑1 (MQ=255)
                                         gCAGACGGTTCGGTGATCCCGTGGCTGGCAAAAAGCTGGACTCATTCAGAAGATGGTAAAACCTGGACCTTCACCCTGCGTGATGACGTGAAATTCTCCAACGGTGAACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTg                                                                                                            >  2:778708/1‑150 (MQ=255)
                                         gCAGACGGTTCGGTGATCCCGTGGCTGGCAAAAAGCTGGACTCATTCAGAAGATGGTAAAACCTGGACCTTCACCCTGCGTGATGACGTGAAATTCTCCAACGGTGAACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTg                                                                                                            >  2:2737383/1‑150 (MQ=255)
                                         gCAGACGGTTCGGTGATCCCGTGGCTGGCAAAAAGCTGGACTCATTCAGAAGATGGTAAAACCTGGACCTTCACCCTGCGTGATGACGTGAAATTCTCCAACGGTGAACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTg                                                                                                            <  2:1463435/150‑1 (MQ=255)
                                            gACGGTTCGGTGATCCCGTGGCTGGCAAAAAGCTGGACTCATTCAGAAGATGGTAAAACCTGGACCTTCACCCTGCGTGATGACGTGAAATTCTCCAACGGTGAACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTc                                                                                                         <  2:362482/150‑1 (MQ=255)
                                            gACGGTTCGGTGATCCCGTGGCTGGCAAAAAGCTGGACTCATTCAGAAGATGGTAAAACCTGGACCTTCACCCTGCGTGATGACGTGAAATTCTCCAACGGTGAACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTc                                                                                                         <  2:1034951/150‑1 (MQ=255)
                                            gACGGTTCGGTGATCCCGTGGCTGGCAAAAAGCTGGACTCATTCAGAAGATGGTAAAACCTGGACCTTCACCCTGCGTGATGACGTGAAATTCTCCAACGGTGAACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTc                                                                                                         <  2:594430/150‑1 (MQ=255)
                                                gTTCGGTGATCCCGTGGCTGGCAAAAAGCTGGACTCATTCAGAAGATGGTAAAACCTGGACCTTCACCCTGCGTGATGACGTGAAATTCTCCAACGGTGAACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATa                                                                                                     <  2:1697140/150‑1 (MQ=255)
                                                gTTCGGTGATCCCGTGGCTGGCAAAAAGCTGGACTCATTCAGAAGATGGTAAAACCTGGACCTTCACCCTGCGTGATGACGTGAAATTCTCCAACGGTGAACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATa                                                                                                     <  2:1593040/150‑1 (MQ=255)
                                                   cGGTGATCCCGTGGCTGGCAAAAAGCTGGACTCATTCAGAAGATGGTAAAACCTGGACCTTCACCCTGCGTGATGACGGGAAATTCTCCAACGGTGAACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAAcc                                                                                                  <  1:2501122/150‑1 (MQ=255)
                                                    ggTGATCCCGTGGCTGGCAAAAAGCTGGACTCATTCAGAAGATGGTAAAACCTGGACCTTCACCCTGCGTGATGACGTGAAATTCTCCAACGGTGAACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCg                                                                                                 >  2:1555444/1‑150 (MQ=255)
                                                     gTGATCCCGTGGCTGGCAAAAAGCTGGACTCATTCAGAAGATGGTAAAACCTGGACCTTCACCCTGCGTGATGACGTGAAATTCTCCAACGGTGAACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGt                                                                                                >  1:340314/1‑150 (MQ=255)
                                                      tGATCCCGTGGCTGGCAAAAAGCTGGACTCATTCAGAAGATGGTAAAACCTGGACCTTCACCCTGCGTGATGACGTGAAATTCTCCAACGGTGAACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTc                                                                                               >  1:693191/1‑150 (MQ=255)
                                                      tGATCCCGTGGCTGGCAAAAAGCTGGACTCATTCAGAAGATGGTAAAACCTGGACCTTCACCCTGCGTGATGACGTGAAATTCTCCAACGGTGAACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTc                                                                                               >  1:1723026/1‑150 (MQ=255)
                                                      tGATCCCGTGGCTGGCAAAAAGCTGGACTCATTCAGAAGATGGTAAAACCTGGACCTTCACCCTGCGTGATGACGTGAAATTCTCCAACGGTGAACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTc                                                                                               >  1:692933/1‑150 (MQ=255)
                                                         tCCCGTGGCTGGCAAAAAGCTGGACTCATTCAGAAGATGGTAAAACCTGGACCTTCACCCTGCGTGATGACGTGAAATTCTCCAACGGTGAACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAAc                                                                                            <  2:1950146/150‑1 (MQ=255)
                                                              tggctggcAAAAAGCTGGACTCACTCAGAAGATGGTAAAACCTGGACCTTCACCCTGCGTGATGACGTGAAATTCTCCAACGGTGAACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCAc                                                                                       >  2:264227/1‑150 (MQ=255)
                                                              tggctggcAAAAAGCTGGACTCACTCAGAAGATGGTAAAACCTGGACCTTCACCCTGCGTGATGACGTGAAATTCTCCAACGGTGAACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCAc                                                                                       >  2:2714531/1‑150 (MQ=255)
                                                                      aaaaaGCTGGACTCATTCAGAAGATGGTAAAACCTGGACCTTCACCCTGCGTGATGACGTGAAATTCTCCAACGGTGAACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCctggct                                                                               >  1:1880000/1‑150 (MQ=255)
                                                                          aGCTGGACTCATTCAGAAGATGGTAAAACCTGGACCTTCACCCTGCGTGATGACGTGAAATTCTCCAACGGTGAACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAg                                                                           <  1:562195/150‑1 (MQ=255)
                                                                             tGGACTCATTCAGAAGATGGTAAAACCTGGACCTTCACCCTGCGTGATGACGTGAAATTCTCCAACGGTGAACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTg                                                                        <  1:2167303/150‑1 (MQ=255)
                                                                              ggACTCATTCAGAAGATGGTAAAACCTGGACCTTCACCCTGCGTGATGACGTGAAATTCTCCAACGGTGAACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTgg                                                                       >  1:2457247/1‑150 (MQ=255)
                                                                              ggACTCATTCAGAAGATGGTAAAACCTGGACCTTCACCCTGCGTGATGACGTGAAATTCTCCAACGGTGAACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTgg                                                                       <  1:2383058/150‑1 (MQ=255)
                                                                               gACTCATTCAGAAGATGGTAAAACCTGGACCTTCACCCTGCGTGATGACGTGAAATTCTCCAACGGTGAACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGc                                                                      >  2:2430499/1‑150 (MQ=255)
                                                                               gACTCATTCAGAAGATGGTAAAACCTGGACCTTCACCCTGCGTGATGACGTGAAATTCTCCAACGGTGAACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGc                                                                      >  2:731766/1‑150 (MQ=255)
                                                                                  tCATTCAGAAGATGGTAAAACCTGGACCTTCACCCTGCGTGATGACGTGAAATTCTCCAACGGTGAACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCaaa                                                                   >  1:466713/1‑150 (MQ=255)
                                                                                  tCATTCAGAAGATGGTAAAACCTGGACCTTCACCCTGCGTGATGACGTGAAATTCTCCAACGGTGAACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCaaa                                                                   >  2:2565672/1‑150 (MQ=255)
                                                                                  tCATTCAGAAGATGGTAAAACCTGGACCTTCACCCTGCGTGATGACGTGAAATTCTCCAACGGTGAACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCaaa                                                                   >  2:1384988/1‑150 (MQ=255)
                                                                                  tCATTCAGAAGATGGTAAAACCTGGACCTTCACCCTGCGTGATGACGTGAAATTCTCCAACGGTGAACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCaaa                                                                   >  2:206532/1‑150 (MQ=255)
                                                                                  tCATTCAGAAGATGGTAAAACCTGGACCTTCACCCTGCGTGATGACGTGAAATTCTCCAACGGTGAACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCaaa                                                                   >  1:790788/1‑150 (MQ=255)
                                                                                     ttCAGAAGATGGTAAAACCTGGACCTTCACCCTGCGTGATGACGTGAAATTCTCCAACGGTGAACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCaaa                                                                   <  2:1162337/147‑1 (MQ=255)
                                                                                       cAGAAGATGGTAAAACCTGGACCTTCACCCTGCGTGATGACGTGAAATTCTCCAACGGTGAACCGTTCTATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGa                                                              <  2:117079/150‑1 (MQ=255)
                                                                                       cAGAAGATGGTAAAACCTGGACCTTCACCCTGCGTGATGACGTGAAATTCTCCAACGGTGAACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCata                                                              >  2:2741319/1‑148 (MQ=255)
                                                                                       cAGAAGATGGTAAAACCTGGACCTTCACCCTGCGTGATGACGTGAAATTCTCCAACGGTGAACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGa                                                              <  2:456503/150‑1 (MQ=255)
                                                                                       cAGAAGATGGTAAAACCTGGACCTTCACCCTGCGTGATGACGTGAAATTCTCCAACGGTGAACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGa                                                              <  1:2302877/150‑1 (MQ=255)
                                                                                          aagaTGGTAAAACCTGGACCTTCACCCTGCGTGATGACGTGAAATTCTCCAACGGTGAACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGAttg                                                           >  2:422571/1‑150 (MQ=255)
                                                                                          aagaTGGTAAAACCTGGACCTTCACCCTGCGTGATGACGTGAAATTCTCCAACGGTGAACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGAgtg                                                           >  2:1565832/1‑150 (MQ=255)
                                                                                            gaTGGTAAAACCTGGACCTTCACCCTGCGTGATGACGTGAAATTCTCCAACGGTGAACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGAttgtt                                                         <  2:914138/150‑1 (MQ=255)
                                                                                            gaTGGTAAAACCTGGACCTTCACCCTGCGTGATGACGTGAAATTCTCCAACGGTGAACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGAttgtt                                                         <  2:2039423/150‑1 (MQ=255)
                                                                                               ggTAAAACCTGGACCTTCACCCTGCGTGATGACGTGAAATTCTCCAACGGTGAACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGATTGTTGAt                                                      <  1:602277/150‑1 (MQ=255)
                                                                                               ggTAAAACCTGGACCTTCACCCTGCGTGATGACGTGAAATTCTCCAACGGTGAACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGATTGTTGAt                                                      <  1:1094520/150‑1 (MQ=255)
                                                                                                 tAAAACCTGGACCTTCACCCTGCGTGATGACGTGAAATTCTCCAACGGTGAACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGATTGTTGATGt                                                    >  2:2382063/1‑150 (MQ=255)
                                                                                                 tAAAACCTGGACCTTCACCCTGCGTGATGACGTGAAATTCTCCAACGGTGAACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGATTGTTGATGt                                                    >  2:1063250/1‑150 (MQ=255)
                                                                                                   aaaCCTGGACCTTCACCCTGCGTGATGACGTGAAATTCTCCAACGGTGAACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGATTGTTGATGTTa                                                  >  1:2475733/1‑150 (MQ=255)
                                                                                                      ccTGGACCTTCACCCTGCGTGATGACGTGAAATTCTCCAACGGTGAACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGATTGTTGATGTTAAAg                                               <  2:1359901/150‑1 (MQ=255)
                                                                                                       cTGGACTTTCACCCTGCGTGATGACGTGAAATTCTCCAACGGTGAACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGATTGTTGATGTTAAAGc                                              >  1:1991528/1‑150 (MQ=255)
                                                                                                       cTGGACCTTCACCCTGCGTGATGACGTGAAATTCTCCAACGGTGAACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGATTGTTGATGTTAAAGc                                              >  1:2564918/1‑150 (MQ=255)
                                                                                                       cTGGACCTTCACCCTGCGTGATGACGTGAAATTCTCCAACGGTGAACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGATTGTTGATGTTAAAGc                                              >  2:558067/1‑150 (MQ=255)
                                                                                                        tGGACCTTCACCCTGCGTGATGACGTGAAATTCTCCAACGGTGAACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGATTGTTGATGTTAAAGCa                                             >  1:2105360/1‑150 (MQ=255)
                                                                                                        tGGACCTTCACCCTGCGTGATGACGTGAAATTCTCCAACGGTGAACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGATTGTTGATGTTAAAGCa                                             <  1:931838/150‑1 (MQ=255)
                                                                                                              ttCACCCTGCGTGATGACGTGAAATTCTCCAACGGTGAACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGATTGTTGATGTTAAAGCACTCAGt                                       <  2:1373237/150‑1 (MQ=255)
                                                                                                               tCACCCTGCGTGATGACGTGAAATTCTCCAACGGTGAACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGATTGTTGATGTTAAAGCACTCAGTa                                      <  2:315024/150‑1 (MQ=255)
                                                                                                                  cccTGCGTGATGACGTGAAATTCTCCAACGGTGAACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGATTGTTGATGTTAAAGCACTCAGTaaaa                                   <  2:1253095/150‑1 (MQ=255)
                                                                                                                    cTGCGTGATGACGTGAAATTCTCCAACGGTGAACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGATTGTTGATGTTAAAGCACTCAGTAAAACa                                 <  1:1555444/150‑1 (MQ=255)
                                                                                                                    cTGCGTGATGACGTGAAATTCTCCAACGGTGAACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGATTGTTGATGTTAAAGCACTCAGTAAAACa                                 <  2:573669/150‑1 (MQ=255)
                                                                                                                      gCGTGATGACGTGAAATTCTCCAACGGTGAACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGATTGTTGATGTTAAAGCACTCAGTAAAACaga                               <  2:870947/150‑1 (MQ=255)
                                                                                                                      gCGTGATGACGTGAAATTCTCCAACGGTGAACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGATTGTTGATGTTAAAGCACTCAGTAAAACaga                               <  1:886884/150‑1 (MQ=255)
                                                                                                                      gCGTGATGACGTGAAATTCTCCAACGGTGAACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGATTGTTGATGTTAAAGCACTCAGTAAAACaga                               <  1:744715/150‑1 (MQ=255)
                                                                                                                      gCGTGATGACGTGAAATTCTCCAACGGTGAACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGATTGTTGATGTTAAAGCACTCAGTAAAACaga                               <  2:2347102/150‑1 (MQ=255)
                                                                                                                           atgaCGTGAAATTCTCCAACGGTGAACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGATTGTTGATGTTAAAGCACTCAGTAAAACAGAGCTGc                          >  1:1665851/1‑150 (MQ=255)
                                                                                                                              aCGTGAAATTCTCCAACGGTGAACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGATTGTTGATGTTAAAGCACTCAGTAAAACAGAGCTGCaaa                       >  2:639332/1‑150 (MQ=255)
                                                                                                                              aCGTGAAATTCTCCAACGGTGAACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGATTGTTGATGTTAAAGCACTCAGTAAAACAGAGCTGCaaa                       >  1:1145608/1‑150 (MQ=255)
                                                                                                                              aCGTGAAATTCTCCAACGGTGAACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGATTGTTGATGTTAAAGCACTCAGTAAAACAGAGCTGCaaa                       >  1:21594/1‑150 (MQ=255)
                                                                                                                              aCGTGAAATTCTCCAACGGTGAACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGATTGTTGATGTTAAAGCACTCAGTAAAACAGAGCTGCaaa                       >  1:2601037/1‑150 (MQ=255)
                                                                                                                                      ttCTCCAACGGTGAACCGTTCGATGCCGAGGCGGCGGCAGGAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGATTGTTGATGTTAAAGCACTCAGTAAAACAGAGCTGCAAATTACCCTg               >  2:1836728/1‑150 (MQ=255)
                                                                                                                                      ttCTCCAACGGTGAACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGATTGTTGATGTTAAAGCACTCAGTAAAACAGAGCTGCAAATTACCCTg               >  2:950192/1‑150 (MQ=255)
                                                                                                                                      ttCTCCAACGGTGAACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGATTGTTGATGTTAAAGCACTCAGTAAAACAGAGCTGCAAATTACCCTg               >  2:1624608/1‑150 (MQ=255)
                                                                                                                                           cAACGGTGAACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGATTGTTGATGTTAAAGCACTCAGTAAAACAGAGCTGCAAATTACCCTGAAAAg          >  1:30118/1‑150 (MQ=255)
                                                                                                                                           cAACGGTGAACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGATTGTTGATGTTAAAGCACTCAGTAAAACAGAGCTGCAAATTACCCTGAAAAg          >  1:886614/1‑150 (MQ=255)
                                                                                                                                            aaCGGTTCACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGATTGTTGATGTTAAAGCACTCAGTAAAACAGAGCTGCAAATTACCCTGAAAAgc         >  1:1232154/1‑150 (MQ=255)
                                                                                                                                            aaCGGTTCACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGATTGTTGATGTTAAAGCACTCAGTAAAACAGAGCTGCAAATTACCCTGAAAAgc         >  1:879136/1‑150 (MQ=255)
                                                                                                                                            aaCGGTTCACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGATTGTTGATGTTAAAGCACTCAGTAAAACAGAGCTGCAAATTACCCTGAAAAgc         >  1:2187593/1‑150 (MQ=255)
                                                                                                                                            aaCGGTTCACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGATTGTTGATGTTAAAGCACTCAGTAAAACAGAGCTGCAAATTACCCTGAAAAgc         >  1:2592670/1‑150 (MQ=255)
                                                                                                                                            aaCGGTGAACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGATTGTTGATGTTAAAGCACTCAGTAAAACAGAGCTGCAAATTACCCTGAAAAgc         <  1:110206/150‑1 (MQ=255)
                                                                                                                                             aCGGTGAACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGATTGTTGATGTTAAAGCACTCAGTAAAACAGAGCTGCAAATTACCCTGAAAAgcg        <  1:1095892/150‑1 (MQ=255)
                                                                                                                                             aCGGTGAACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGATTGTTGATGTTAAAGCACTCAGTAAAACAGAGCTGCAAATTACCCTGAAAAgcg        <  1:1274657/150‑1 (MQ=255)
                                                                                                                                             aCGGTGAACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGATTGTTGATGTTAAAGCACTCAGTAAAACAGAGCTGCAAATTACCCTGAAAAgcg        <  1:2271655/150‑1 (MQ=255)
                                                                                                                                               ggTGAACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGATTGTTGATGTTAAAGCACTCAGTAAAACAGAGCTGCAAATTACCCTGAAAAGCGcc      >  2:1467467/1‑150 (MQ=255)
                                                                                                                                                gTGAACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGATTGTTGATGTTAAAGCACTCAGTAAAACAGAGCTGCAAATTACCCTGAAAAGCGCct     >  1:1981555/1‑150 (MQ=255)
                                                                                                                                                gTGAACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGATTGTTGATGTTAAAGCACTCAGTAAAACAGAGCTGCAAATTACCCTGAAAAGCGCct     >  1:1800648/1‑150 (MQ=255)
                                                                                                                                                gTGAACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGATTGTTGATGTTAAAGCACTCAGTAAAACAGAGCTGCAAATTACCCTGAAAAGCGCct     >  2:1026146/1‑150 (MQ=255)
                                                                                                                                                gTGAACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGATTGTTGATGTTAAAGCACTCAGTAAAACAGAGCTGCAAAATACCCTGAAAAGCGCct     >  2:1499278/1‑150 (MQ=255)
                                                                                                                                                gTGAACAGTTCGATGCCTAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGATTGTTGATTTTAAAGCACTCAGTAAAACAGAGCTGCAAATTACCCTGAAAAGCGCct     >  2:1026148/1‑150 (MQ=255)
                                                                                                                                                 tGAACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGATTGTTGATGTTAAAGCACTCAGTAAAACAGAGCTGCAAATTACCCTGAAAAGCGCcta    >  2:523130/1‑150 (MQ=255)
                                                                                                                                                 tGAACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGATTGTTGATGTTAAAGCACTCAGTAAAACAGAGCTGCAAATTACCCTGAAAAGCGCcta    >  1:65675/1‑150 (MQ=255)
                                                                                                                                                 tGAACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGATTGTTGATGTTAAAGCACTCAGTAAAACAGAGCTGCAAATTACCCTGAAAAGCGCcta    >  1:65682/1‑150 (MQ=255)
                                                                                                                                                   aaCCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGATTGTTGATGTTAAAGCACTCAGTAAAACAGAGCTGCAAATTACCCTGAAAAGCGCctact  <  2:232230/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
TCGCCCAGAGCATGGTTTATGAACCATTGGTGAAATATCAGGCAGACGGTTCGGTGATCCCGTGGCTGGCAAAAAGCTGGACTCATTCAGAAGATGGTAAAACCTGGACCTTCACCCTGCGTGATGACGTGAAATTCTCCAACGGTGAACCGTTCGATGCCGAGGCGGCGGCAGAAAACTTCCGCGCAGTGCTCGATAACCGTCAACGTCACGCCTGGCTGGAGCTGGCAAACCAGATTGTTGATGTTAAAGCACTCAGTAAAACAGAGCTGCAAATTACCCTGAAAAGCGCCTACT  >  NC_000913/3613815‑3614111

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