Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,193,013 A→T 6.5% intergenic (‑217/+46) yehE ← / ← mrp DUF2574 domain‑containing protein YehE/P‑loop NTPase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,193,0130AT 6.5% 330.6 / 7.9 104intergenic (‑217/+46)yehE/mrpDUF2574 domain‑containing protein YehE/P‑loop NTPase family protein
Reads supporting (aligned to +/- strand):  ref base A (45/52);  new base T (5/2);  total (50/54)
Fisher's exact test for biased strand distribution p-value = 2.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 107 total aligned reads displayed.

TCCACTACACAGCCACGAAGGATGATAATGAAGCATTGCCTGTATGATCAATCGACTTTGTAGAGTTTCGGACGAAGGTCCGCAGAATATTCGCAGTATTAAATAAGTGTTCATTAATGCGGCAAAAACTAATACACCGCATCAATGTAACATCTCTGCGGTGTATGAAGAATAACCAGAATGGTTAATAGGCGTTAGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCACCCTGCCAGTAGAGCTGAGCTGCAACGCGGTCAGCCAGCTGGCGGTAGATAGCGGTAAATTCGC  >  NC_000913/2192865‑2193162
                                                                                                                                                    |                                                                                                                                                     
tCCACTACACAGCCACGAAGGATGATAATGAAGCATTGCCTGTATGATCAATCGACTTTGTAGAGTTTCGGACGAAGGTCCGCAGAATATTCGCAGTATTAAATAAGTGTTCATTAATGCGGCAAAAACTAATACACCGCATCAATGTaa                                                                                                                                                      >  1:296766/1‑150 (MQ=255)
tCCACTACACAGCCACGAAGGATGATAATGAAGCATTGCCTGTATGATCAATCGACTTTGTAGAGTTTCGGACGAAGGTCCGCAGAATATTCGCAGTATTAAATAAGTGTTCATTAATGCGGCAAAAACTAATACACCGCATCAATGTaa                                                                                                                                                      >  1:1444260/1‑150 (MQ=255)
 ccACTACACAGCCACGAAGGATGATAATGAAGCATTGCCTGTATGATCAATCGACTTTGTAGAGTTTCGGACGAAGGTCCGCAGAATATTCGCAGTATTAAATAAGTGTTCATTAATGCGGCAAAAACTAATACACCGCATCAATGTAAc                                                                                                                                                     <  2:1142944/150‑1 (MQ=255)
 ccACTACACAGCCACGAAGGATGATAATGAAGCATTGCCTGTATGATCAATCGACTTTGTAGAGTTTCGGACGAAGGTCCGCAGAATATTCGCAGTATTAAATAAGTGTTCATTAATGCGGCAAAAACTAATACACCGCATCAATGTAAc                                                                                                                                                     >  1:1837767/1‑150 (MQ=255)
 ccACTACACAGCCACGAAGGATGATAATGAAGCATTGCCTGTATGATCAATCGACTTTGTAGAGTTTCGGACGAAGGTCCGCAGAATATTCGCAGTATTAAATAAGTGTTCATTAATGCGGCAAAAACTAATACACCGCATCAATGTAAc                                                                                                                                                     <  2:285959/150‑1 (MQ=255)
 ccACTACACAGCCACGAAGGATGATAATGAAGCATTGCCTGTATGATCAATCGACTTTGTAGAGTTTCGGACGAAGGTCCGCAGAATATTCGCAGTATTAAATAAGAGTTCATAAATGCGGCAAAAACTAATACACCGCATCAATGTAAc                                                                                                                                                     >  2:1569630/1‑150 (MQ=255)
       cacaGCCACGAAGGATGATAATGAAGCATTGCCTGTATGATCAATCGACTTTGTAGAGTATCGGACGAAGGTCCGCAGAATATTCGCAGTATTAAATAAGTGTTCATTAATGCGGCAAAAACTAATACACCGCATAAATGTAACATctct                                                                                                                                               >  1:833094/1‑150 (MQ=255)
        acaTCCACGAAGGATGATAATGAAGCATTGCCTGTATGATCAATCGACTTTGTAGAGTTTCGGACGAAGGTCCGCAGAATATTCGCAGTATTAAATAAGTGTTCATTAATGCGGCAAAAACTAATACACCGCATCAATGTAACATCTCTg                                                                                                                                              <  2:1865356/150‑1 (MQ=255)
            ccACGAAGGATGATAATGAAGCATTGCCTGTATGATCAATCGACTTTGTAGAGTTTCGGACGAAGGTCCGCAGAATATTCGCAGTATTAAATAAGTGTTCATTAATGCGGCAAAAACTAATACACCGCATCAATGTAACATCTCTGCGgt                                                                                                                                          >  1:1337206/1‑150 (MQ=255)
             cACGAAGGATGATAATGAAGCATTGCCTGTATGATCAATCGACTTTGTAGAGTTTCGGACGAAGGTCCGCAGAATATTCGCAGTATTAAATAAGTGTTCATTAATGCGGCAAAAACTAATACACCGCATCAATGTAACATCTCTGCGgtg                                                                                                                                         <  1:1038213/150‑1 (MQ=255)
             cACGAAGGATGATAATGAAGCATTGCCTGTATGATCAATCGACTTTGTAGAGTTTCGGACGAAGGTCCGCAGAATATTCGCAGTATTAAATAAGTGTTCATTAATGCGGCAAAAACTAATACACCGCATCAATGTAACATCTCTGCGgtg                                                                                                                                         <  1:1073212/150‑1 (MQ=255)
                gAAGGATGATAATGAAGCATTGCCTGTATGATCAATCGACTTTGTAGAGTTTCGGACGAAGGTCCGCAGAATATTCGCAGTATTAAATAAGTGTTCATTAATGCGGCAAAAACTAATACACCGCATCAATGTAACATCTCTGCGGTGTAt                                                                                                                                      >  2:727731/1‑150 (MQ=255)
                 aaGGATGATAATGAAGCATTGCCTGTATGATCAATCGACTTTGTAGAGTTTCGGACGAAGGTCCGCAGAATATTCGCAGTATTAAATAAGTGTTCATTAATGCGGCAAAAACTAATACACCGCATCAATGTAACATCTCTGCGGTGTATg                                                                                                                                     <  2:389885/150‑1 (MQ=255)
                   ggATGATAATGAAGCATTGCCTGTATGATCAATCGACTTTGTAGAGTTTCGGACGAAGGTCCGCAGAATATTCGCAGTATTAAATAAGTGTTCATTAATGCGGCAAAAACTAATACACCGCATCAATGTAACATCTCTGCGGTGTATgaa                                                                                                                                   >  1:570817/1‑150 (MQ=255)
                    gatgatAATGAAGCATTGCCTGTATGATCAATCGACTTTGTAGAGTTTCGGACGAAGGTCCGCAGAATATTCGCAGTATTAAATAAGTGTTCATTAATGCGGCAAAAACTAATACACCGCATCAATGTAACATCTCTGCGGTGTATgaag                                                                                                                                  <  2:205640/150‑1 (MQ=255)
                    gatgatAATGAAGCATTGCCTGTATGATCAATCGACTTTGTAGAGTTTCGGACGAAGGTCCGCAGAATATTCGCAGTATTAAATAAGTGTTCATTAATGCGGCAAAAACTAATACACCGCATCAATGTAACATCTCTGCGGTGTATgaag                                                                                                                                  >  2:1211990/1‑150 (MQ=255)
                     atgatAATGAAGCATTGCCTGTATGATCAATCGACTTTGTAGAGTTTCGGACGAAGGTCCGCAGAATATTCGCAGTATTAAATAAGTGTTCATTAATGCGGCAAAAACTAATACACCGCATCAATGTAACATCTCTGCGGTGTATgaaga                                                                                                                                 >  2:1928631/1‑150 (MQ=255)
                           aTGAAGCATTGCCTGTATGATCAATCGACTTTGTAGAGTTTCGGACGAAGGTCCGCAGAATATTCGCAGTATTAAATAAGTGTTCATTAATGCGGCAAAAACTAATACACCGCATCAATGTAACATCTCTGCGGTGTATGAAGAATAAcc                                                                                                                           >  2:532946/1‑150 (MQ=255)
                           aTGAAGCATTGCCTGTATGATCAATCGACTTTGTAGAGTTTCGGACGAAGGTCCGCAGAATATTAGCAGTATTAAATAAGTGTTCATTAATGCGGCAAAAACTAATACACCGCATCAATGTAACATCTCTGCGGTGTATGAAGAATAAcc                                                                                                                           >  2:580740/1‑150 (MQ=255)
                                gCATTGCCTGTATGATCAATCGACTTTGTAGAGTTTCGGACGAAGGTCCGCAGAATATTCGCAGTATTAAATAAGTGTTCATTAATGCGGCAAAAACTAATACACCGCATCAATGTAACATCTCTGCGGTGTATGAAGAATAACCAGAAt                                                                                                                      >  1:1545154/1‑150 (MQ=255)
                                 cATTGCCTGTATGATCAATCGACTTTGTAGAGTTTCGGACGAAGGTCCGCAGAATATTCGCAGTATTAAATAAGTGTTCATTAATGCGGCAAAAACTAATACACCGCATCAATGTAACATCTCTGCGGTGTATGAAGAATAACCAGAATg                                                                                                                     <  2:296766/150‑1 (MQ=255)
                                 cATTGCCTGTATGATCAATCGACTTTGTAGAGTTTCGGACGAAGGTCCGCAGAATATTCGCAGTATTAAATAAGTGTTCATTAATGCGGCAAAAACTAATACACCGCATCAATGTAACATCTCTGCGGTGTATGAAGAATAACCAGAATg                                                                                                                     >  2:1597936/1‑150 (MQ=255)
                                   ttGCCTGTATGATCAATCGACTTTGTAGAGTTTCGGACGAAGGTCCGCAGAATATTCGCAGTATTAAATAAGTGTTCATTAATGCGGCAAAAACTAATACACCGCATCAATGTAACATCTCTGCGGTGTATGAAGAATAACCAGAATGGt                                                                                                                   <  1:1922820/150‑1 (MQ=255)
                                    tGCCTGTATGATCAATCGACTTTGTAGAGTTTCGGACGAAGGTCCGCAGAATATTCGCAGTATTAAATAAGTGTTCATTAATGCGGCAAAAACTAATACACCGCATCAATGTAACATCTCTGCGGTGTATGAAGAATAACCAGAATGGtt                                                                                                                  <  2:684754/150‑1 (MQ=255)
                                    tGCCTGTATGATCAATCGACTTTGTAGAGTTTCGGACGAAGGTCCGCAGAATATTCGCAGTATTAAATAAGTGTTCATTAATGCGGCAAAAACTAATACACCGCATCAATGTAACATCTCTGCGGTGTATGAAGAATAACCAGAATGGtt                                                                                                                  <  1:1766475/150‑1 (MQ=255)
                                     gCCTGTATGATCAATCGACTTTGTAGAGTTTCGGACGAAGGTCCGCAGAATATTCGCAGTATTAAATAAGTGTTCATTAATGCGGCAAAAACTAATACACCGCATCAATGTAACATCTCTGCGGTGTATGAAGAATAACCAGAATGGTTa                                                                                                                 >  2:1772335/1‑150 (MQ=255)
                                        tGTATGATCAATCGACTTTGTAGAGTTTCGGACGAAGGTCCGCAGAATATTCGCAGTATTAAATAAGTGTTCATTAATGCGGCAAAAACTAATACACCGCATCAATGTAACATCTCTGCGGTGTATGAAGAATAACCAGAATGGTTAATa                                                                                                              <  2:1349318/150‑1 (MQ=255)
                                           aTGATCAATCGACTTTGTAGAGTTTCGGACGAAGGTCCGCAGAATATTCGCAGTATTAAATAAGTGTTCATTAATGCGGCAAAAACTAATACACCGCATCAATGTAACATCTCTGCGgtgt                                                                                                                                        >  1:132140/1‑121 (MQ=255)
                                           aTGATCAATCGACTTTGTAGAGTTTCGGACGAAGGTCCGCAGAATATTCGCAGTATTAAATAAGTGTTCATTAATGCGGCAAAAACTAATACACCGCATCAATGTAACATCTCTGCGgtgt                                                                                                                                        <  2:132140/121‑1 (MQ=255)
                                           aTGATCAATCGACTTTGTAGAGTTTCGGACGAAGGTCCGCAGAATATTCGCAGTATTAAATAAGTGTTCATTAATGCGGCAAAAACTAATACACCGCATCAATGTAACATCTCTGCGGTGTATGAAGAATAACCATTCTGGTTattcttc                                                                                                           <  2:1444260/150‑7 (MQ=255)
                                            tGATCAATCGACTTTGTAGAGTTTCGGACGAAGGTCCGCAGAATATTCGCAGTATTAAATAAGTGTTCATTAATGCGGCAAAAACTAATACACCGCATCAATGTAACATCTCTGCGGTGTATGAAGAATAACCATTCTGGTTattcttca                                                                                                          <  1:1866547/150‑8 (MQ=255)
                                            tGATCAATCGACTTTGTAGAGTTTCGGACGAAGGTCCGCAGAATATTCGCAGTATTAAATAAGTGTTCATTAATGCGGCAAAAACTAATACACCGCATCAATGTAACATCTCTGCGGTGTATGAAGAATAACCAGAATGGTTAATAGGCg                                                                                                          >  2:1943697/1‑150 (MQ=255)
                                             gATCAATCGACTTTGTAGAGTTTCGGACGAAGGTCCGCAGAATATTCGCAGTATTAAATAAGTGTTCATTAATGCGGCAAAAACTAATACACCGCATCAATGTAACATCTCTGCGGTGTATGAAGAATAACCAGAATGGTTAATAGGCGt                                                                                                         >  1:1548196/1‑150 (MQ=255)
                                                    cGACTTTGTAGAGTTTCGGACGAAGGTCCGCAGAATATTCGCAGTATTAAATAAGTGTTCATTAATGCGGCAAAAACTAATACACCGCATCAATGTAACATCTCTGCGGTGTATGAAGAATAACCAGAATGGTTAATAGGCGTTAGACcg                                                                                                  >  1:1698062/1‑150 (MQ=255)
                                                      aCTTTGTAGAGTTTCGGACGAAGGTCCGCAGAATATTCGCAGTATTAAATAAGTGTTCATTAATGCGGCAAAAACTAATACACCGCATCAATGTAACATCTCTGCGGTGTATGAAGAATAACCAGAATGGTTAATAGGCGTTAGACcgcg                                                                                                <  2:962470/150‑1 (MQ=255)
                                                      aCTTTGTAGAGTTTCGGACGAAGGTCCGCAGAATATTCGCAGTATTAAATAAGTGTTCATTAATGCGGCAAAAACTAATACACCGCATCAATGTAACATCTCTGCGGTGTATGAAGAATAACCAGAATGGTTAATAGGCGTTAGACcgcg                                                                                                >  2:264289/1‑150 (MQ=255)
                                                       cTTTGTAGAGTTTCGGACGAAGGTCCGCAGAATATTCGCAGTATTAAATAAGTGTTCATTAATGCGGCAAAAACTAATACACCGCATCAATGTAACATCTCTGCGGTGTATGAAGAATAACCAGAATGGTTAATAGGCGTTAGACcgcgc                                                                                               <  2:1337206/150‑1 (MQ=255)
                                                        tttGTAGAGTTTCGGACGAAGGTCCGCAGAATATTCGCAGTATTAAATAAGTGTTCATTAATGCGGCAAAAACTAATACACCGCATCAATGTAACATCTCTGCGGTGTATGAAGAATAACCAGAATGGTTAATAGGCGTTAGACcgcgcg                                                                                              >  2:1908297/1‑150 (MQ=255)
                                                             agagTTTCGGACGAAGGTCCGCAGAATATTCGCAGTATTAAATAAGTGTTCATTAATGCGGCAAAAACTAATACACCGCATCAATGTAACATCTCTGCGGTGTATGAAGAATAACCAGAATGGTTAATAGGCGTTAGACCGCGCggaagg                                                                                         >  2:1314810/1‑150 (MQ=255)
                                                                gTTTCGGACGAAGGTCCGCAGAATATTCGCAGTATTAAATAAGTGTTCATTAATGCGGCAAAAACTAATACACCGCATCAATGTAACATCTCTGCGGTGTATGAAGAATAACCAGAATGGTTAATAGGCGTTAGACCGCGCGGAAGGaaa                                                                                      >  2:1146139/1‑150 (MQ=255)
                                                                  ttCGGACGAAGGTCCGCAGAATATTCGCAGTATTAAATAAGTGTTCATTAATGCGGCAAAAACTAATACACCGCATCAATGTAACATCTCTGCGGTGTATGAAGAATAACCAGAATGGTTAATAGGCGTTAGACCGCGCGGAAGGAAAtc                                                                                    >  1:1059091/1‑150 (MQ=255)
                                                                   tCGGACGAAGGTCCGCAGAATATTCGCAGTATTAAATAAGTGTTCATTAATGCGGCAAAAACTAATACACCGCATCAATGTAACATCTCTGCGGTGTATGAAGAATAACCAGAATGGTTAATAGGCGTTAGACCGCGCGGAAGGAAAtct                                                                                   >  1:1101136/1‑150 (MQ=255)
                                                                    cGGACGAAGGTCCGCAGAATATTCGCAGTATTAAATAAGTGTTCATTAATGCGGCAAAAACTAATACACCGCATCAATGTAACATCTCTGCGGTGTATGAAGAATAACCAGAATGGTTAATAGGCGTTAGACCGCGCGGAAGGAAAtctc                                                                                  >  2:398378/1‑150 (MQ=255)
                                                                     ggACGAAGGTCCGCAGAATATTCGCAGTATTAAATAAGTGTTCATTAATGCGGCAAAAACTAATACACCGCATCAATGTAACATCTCTGCGGTGTATGAAGAATAACCAGAATGGTTAATAGGCGTTAGACCGCGCGGAAGGAAATCTCg                                                                                 <  1:1146139/150‑1 (MQ=255)
                                                                      gACGAAGGTCCGCAGAATATTCGCAGTATTAAATAAGTGTTCATTAATGCGGCAAAAACTAATACACCGCATCAATGTAACATCTCTGCGGTGTATGAAGAATAACCAGAATGGTTAATAGGCGTTAGACCGCGCGGAAGGAAATCTCGc                                                                                <  1:1406376/150‑1 (MQ=255)
                                                                      gACGAAGGTCCGCAGAATATTCGCAGTATTAAATAAGTGTTCATTAATGCGGCAAAAACTAATACACCGCATCAATGTAACATCTCTGCGGTGTATGAAGAATAACCAGAATGGTTAATAGGCGTTAGACCGCGCGGAAGGAAATCTCGc                                                                                <  2:1101136/150‑1 (MQ=255)
                                                                       aCGAAGGTCCGCAGAATATTCGCAGTATTAAATAAGTGTTCATTAATGCGGCAAAAACTAATACACCGCATCAATGTAACATCTCTGCGGTGTATGAAGAATAACCAGAATGGTTAATAGGCGTTAGACCGCGCGGAAGGAAATCTCGcc                                                                               <  2:802937/150‑1 (MQ=255)
                                                                       aCGAAGGTCCGCAGAATATTCGCAGTATTAAATAAGTGTTCATTAATGCGGCAAAAACTAATACACCGCATCAATGTAACATCTCTGCGGTGTATGAAGAATAACCAGAATGGTTAATAGGCGTTAGACCGCGCGGAAGGAAATCTCGcc                                                                               <  1:806325/150‑1 (MQ=255)
                                                                        cGAAGGTCCGCAGAATATTCGCAGTATTAAATAAGTGTTCATTAATGCGGCAAAAACTAATACACCGCATCAATGTAACATCTCTGCGGTGTATGAAGAATAACCAGAATGGTTAATAGGCGTTAGACCGCGCGGAAGGAAATCTCGCCt                                                                              >  1:400656/1‑150 (MQ=255)
                                                                        cGAAGGTCCGCAGAATATTCGCAGTATTAAATAAGTGTTCATTAATGCGGCAAAAACTAATACACCGCATCAATGTAACATCTCTGCGGTGTATGAAGAATAACCAGAATGGTTAATAGGCGTTAGACCGCGCGGAAGGAAATCTCGCCt                                                                              <  2:1459567/150‑1 (MQ=255)
                                                                        cGAAGGTCCGCAGAATATTCGCAGTATTAAATAAGTGTTCATTAATGCGGCAAAAACTAATACACCGCATCAATGTAACATCTCTGCGGTGTATGAAGAATAACCAGAATGGTTAATAGGCGTTAGACCGCGCGGAAGGAAATCTCGCCt                                                                              <  1:982952/150‑1 (MQ=255)
                                                                         gAAGGTCCGCAGAATATTCGCAGTATTAAATAAGTGTTCATTAATGCGGCAAAAACTAATACACCGCATCAATGTAACATCTCTGCGGTGTATGAAGAATAACCAGAATGGTTAATAGGCGTTAGACCGCGCGGAAGGAAATCTCGCCTg                                                                             <  2:189151/150‑1 (MQ=255)
                                                                             gTCCGCAGAATATTCGCAGTATTAAATAAGTGTTCATTAATGCGGCAAAAACTAATACACCGCATCAATGTAACATCTCTGCGGTGTATGAAGAATAACCAGAATGGTTAATAGGCGTTAGACCGCGCGGAAGGAAATCTCGCCTGGAAt                                                                         <  2:1077590/150‑1 (MQ=255)
                                                                               ccGCAGAATATTCGCAGTATTAAATAAGTGTTCATTAATGCGGCAAAAACTAATACACCGCATCAATGTAACATCTCTGCGGTGTATGAAGAATAACCAGAATGGTTAATAGGCGTTAGACCGCGCGGAAGGAAATCTCGCCTGGAATGa                                                                       >  2:1095904/1‑150 (MQ=255)
                                                                                 gCAGAATATTCGCAGTATTAAATAAGTGTTCATTAATGCGGCAAAAACTAATACACCGCATCAATGTAACATCTCTGCGGTGTATGAAGAATAACCAGAATGGTTAATAGGCGTTAGACCGCGCGGAAGGAAATCTCGCCTGGAATGACt                                                                     >  1:338095/1‑150 (MQ=255)
                                                                                   aGAATATTCGCAGTATTAAATAAGTGTTCATTAATGCGGCAAAAACTAATACACCGCATCAATGTAACATCTCTGCGGTGTATGAAGAATAACCAGAATGGTTAATAGGCGTTAGACCGCGCGGAAGGAAATCTCGCCTGGAATGAc                                                                      <  1:1428916/147‑1 (MQ=255)
                                                                                   aGAATATTCGCAGTATTAAATAAGTGTTCATTAATGCGGCAAAAACTAATACACCGCATCAATGTAACATCTCTGCGGTGTATGAAGAATAACCAGAATGGTTAATAGGCGTTAGACCGCGCGGAAGGAAATCTCGCCTGGAATGAc                                                                      >  2:1428916/1‑147 (MQ=255)
                                                                                     aaTATTAGCAGTATTAAATAAGTGTTCATTAATGCGGCAAAAACTAATACACCGCATCAATGTAACATCTCTGCGGTGTATGAAGAATAACCAGAATGGTTAATAGGCGTTAGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCAc                                                                 <  1:580740/150‑1 (MQ=255)
                                                                                       tatTCGCAGTATTAAATAAGTGTTCATTAATGCGGCAAAAACTAATACACCGCATCAATGTAACATCTCTGCGGTGTATGAAGAATAACCAGAATGGTTAATAGGCGTTAGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCAccc                                                               >  1:1410534/1‑150 (MQ=255)
                                                                                           cGCAGTATTAAATAAGTGTTCATTAATGCGGCAAAAACTAATACACCGCATCAATGTAACATCTCTGCGGTGTATGAAGAATAACCAGAATGGTTAATAGGCGTTAGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCACCCTGcc                                                           <  2:1545154/150‑1 (MQ=255)
                                                                                           cGCAGTATTAAATAAGTGTTCATTAATGCGGCAAAAACTAATACACCGCATCAATGTAACATCTCTGCGGTGTATGAAGAATAACCAGAATGGTTAATAGGCGTTAGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCACCCTGcc                                                           <  2:1548196/150‑1 (MQ=255)
                                                                                            gCAGTATTAAATAAGTGTTCATTAATGCGGCAAAAACTAATACACCGCATCAATGTAACATCTCTGCGGTGTATGAAGAATAACCAGAATGGTTAATAGGCGTTAGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCACCCTGCCa                                                          <  2:1085844/150‑1 (MQ=255)
                                                                                             cAGTATTAAATAAGTGTTCATTAATGCGGCAAAAACTAATACACCGCATCAATGTAACATCTCTGCGGTGTATGAAGAATAACCAGAATGGTTAATAGGCGTTAGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCACCCTGCCAg                                                         <  2:175259/150‑1 (MQ=255)
                                                                                              aGTATTAAATAAGTGTTCATTAATGCGGCAAAAACTAATACACCGCATCAATGTAACATCTCTGCGGTGTATGAAGAATAACCAGAATGGTTAATAGGCGTTAGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCACCCTGCCAGt                                                        <  2:1297649/150‑1 (MQ=255)
                                                                                                   tAAATAAGTGTTCATTAATGCGGCAAAAACTAATACACCGCATCAATGTAACATCTCTGCGGTGTATGAAGAATAACCAGAATGGTTAATAGGCGTTAGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCACCCTGCCAGTAgagc                                                   <  1:1381388/150‑1 (MQ=255)
                                                                                                     aaTAAGTGTTCATTAATGCGGCAAAAACTAATACACCGCATCAATGTAACATCTCTGCGGTGTATGAAGAATAACCAGAATGGTTAATAGGCGTTAGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCACCCTGCCAGTAgagctg                                                 <  2:422907/150‑1 (MQ=255)
                                                                                                       tAAGTGTTCATTAATGCGGCAAAAACTAATACACCGCATCAATGTAACATCTCTGCGGTGTATGAAGAATAACCAGAATGGTTAATAGGCGTTAGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCACCCTGCCAGTAgagctgag                                               <  2:94082/150‑1 (MQ=255)
                                                                                                       tAAGTGTTCATTAATGCGGCAAAAACTAATACACCGCATCAATGTAACATCTCTGCGGTGTATGAAGAATAACCAGAATGGTTAATAGGCGTTAGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCACCCTGCCAGTAgagctgag                                               >  1:941630/1‑150 (MQ=255)
                                                                                                       tAAGTGTTCATTAATGCGGCAAAAACTAATACACCGCATCAATGTAACATCTCTGCGGTGTATGAAGAATAACCAGAATGGTTAATAGGCGTTAGACCGCGCGGAAGGAAATCTCGCATGGAATGACTTCACCCTGCCAGTAgagctgag                                               >  1:934787/1‑150 (MQ=255)
                                                                                                              tCATTAATGCGGCAAAAACTAATACACCGCATCAATGTAACATCTCTGCGGTGTATGAAGAATAACCAGAATGGTTAATAGGCGTTAGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCACCCTGCCAGTAGAGCTGAGCTGCAAc                                        >  1:722949/1‑150 (MQ=255)
                                                                                                                   aaTGCGGCAAAAACTAATACACCGCATCAATGTAACATCTCTGCGGTGTATGAAGAATAACCAGAATGGTTAATAGGCGTTAGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCACCCTGCCAGTAGAGCTGAGCTGCAACGCGGt                                   <  1:1597936/150‑1 (MQ=255)
                                                                                                                       cGGCAAAAACTAATACACCGCATCAATGTAACATCTCTGCGGTGTATGAAGTATAACCAGAATGGTTAATAGGCGTTAGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCACCCTGCCAGTAGAGCTGAGCTGCAACGCGGTcagc                               <  1:1943697/150‑1 (MQ=255)
                                                                                                                        ggCAAAAACTAATACACCGCATCAATGTAACATCTCTGCGGTGTATGAAGAATAACCAGAATGGTTAATAGGCGTTAGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCACCCTGCCAGTAGAGCTGAGCTGCAACGCGGTcagcc                              <  1:563480/150‑1 (MQ=255)
                                                                                                                         taaaaaaCTAATACACCGCATCAATGTAACATATCTGCGGTGTATGAAGAATAACCAGAATTGTTAATAGGCGTTAGACCGCGCGTAAGGAAATATCGACTGTAATGACTTCACCCTGCCAGTAGAGCTGAGATGCAACGCCGTcagccc                             >  1:445485/3‑149 (MQ=255)
                                                                                                                             aaaCTAATACACCGCATCAATGTAACATCTCTGCGGTGTATGAAGAATAACCAGAATGGTTAATAGGCGTTAGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCACCCTGCCAGTAGAGCTGAGCTGCAACGCGGTCAGCCAGCTg                         <  1:264289/150‑1 (MQ=255)
                                                                                                                             aaaCTAATACACCGCATCAATGTAACATCTCTGCGGTGTATGAAGAATAACCAGAATGGTTAATAGGCGTTAGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCACCCTGCCAGTAGAGCTGAGCTGCAACGCGGTCAGCCAGCTg                         >  1:1740107/1‑150 (MQ=255)
                                                                                                                             aaaCTAATACACCGCATCAATGTAACATCTCTGCGGTGTATGAAGAATAACCAGAATGGTTAATAGGCGTTAGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCACCCTGCCAGTAGAGCTGAGCTGCAACGCGGTCAGCCAGCTg                         >  1:151190/1‑150 (MQ=255)
                                                                                                                              tcttcaTACACCGCAGAGATGTTACATTGATGCGGTGTATGAAGAATAACCAGAATGGTTAATAGGCGTTAGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCACCCTGCCAGTAGAGCTGAGCTGCAACGCGGTCAGCCAGCTgg                        >  1:332495/6‑150 (MQ=38)
                                                                                                                               cttcaTACACCGCAGAGATGTTACATTGATGCGGTGTATGAAGAATAACCAGAATGGTTAATAGGCGTTAGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCACCCTGCCAGTAGAGCTGAGCTGCAACGCGGTCAGCCAGCTGGc                       <  2:332495/146‑1 (MQ=38)
                                                                                                                               aCTAATACACCGCATCAATGTAACATCTCTGCGGTGTATGAAGAATAACCAGAATGGTTAATAGGCGTTAGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCACCCTGCCAGTAGAGCTGAGCTGCAACGCGGTCAGCCAGCTGGc                       <  1:1420799/150‑1 (MQ=255)
                                                                                                                                cTAATACACCGCATCAATGTAACATCTCTGCGGTGTATGAAGAATAACCAGAATGGTTAATAGGCGTTAGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCACCCTGCCAGTAGAGCTGAGCTGCAACGCGGTCAGCCAGCTGGCg                      <  1:532946/150‑1 (MQ=255)
                                                                                                                                 tcaTACACCGCAGAGATGTTACATTGATGCGGTGTATGAAGAATAACCAGAATGGTTAATAGGCGTTAGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCACCCTGCCAGTAGAGCTGAGCTGCAACGCGGTCAGCCAGCTGGCgg                     >  2:306177/3‑150 (MQ=38)
                                                                                                                                 tAATACACCGCATCAATGTAACATCTCTGCGGTGTATGAAGAATAACCAGAATGGTTAATAGGCGTTAGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCACCCTGCCAGTAGAGCTGAGCTGCAACGCGGTCAGCCAGCTGGCgg                     <  1:727731/150‑1 (MQ=255)
                                                                                                                                  caTACACCGCAGAGATGTTACATTGATGCGGTGTATGAAGAATAACCAGAATGGTTAATAGGCGTTAGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCACCCTGCCAGTAGAGCTGAGCTGCAACGCGGTCAGCCAGCTGGCGGt                    >  2:1411648/2‑150 (MQ=39)
                                                                                                                                   aTACACCGCATCAATGTAACATCTCTGCGGTGTATGAAGAATAACCAGAATGGTTAATAGGCGTTAGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCACCCTGCCAGTAGAGCTGAGCTGCAACGCGGTCAGCCAGCTGGCGGTa                   <  1:1718102/150‑1 (MQ=255)
                                                                                                                                   aTACACCGCATCAATGTAACATCTCTGCGGTGTATGAAGAATAACCAGAATGGTTAATAGGCGTTAGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCACCCTGCCAGTAGAGCTGAGCTGCAACGCGGTCAGCCAGCTGGCGGTa                   <  2:400656/150‑1 (MQ=255)
                                                                                                                                   aTACACCGCATCAATGTAACATCTCTGCGGTGTATGAAGAATAACCAGAATGGTTAATAGGCGTTAGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCACCCTGCCAGTAGAGCTGAGCTGCAACGCGGTCAGCCAGCTGGCGGTa                   >  2:859859/1‑150 (MQ=255)
                                                                                                                                   aTACACCGCAGAGATGTTACATTGATGCGGTGTATGAAGAATAACCAGAATGGTTAATAGGCGTTAGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCACCCTGCCAGTAGAGCTGAGCTGCAACGCGGTCAGCCAGCTGGCGGTa                   <  1:1411648/150‑1 (MQ=255)
                                                                                                                                    tACACCGCATCAATGTAACATCTCTGCGGTGTATGAAGAATAACCAGAATGGTTAATAGGCGTTAGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCACCCTGCCAGTAGAGCTGAGCTGCAACGCGGTCAGCCAGCTGGCGGTAg                  <  2:705351/150‑1 (MQ=255)
                                                                                                                                      cacCGCATCAATGTAACATCTCTGCGGTGTATGAAGAATAACCAGAATGGTTAATAGGCGTTAGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCACCCTGCCAGTAGAGCTGAGCTGCAACGCGGTCAGCCAGCTGGCGGTAGAt                <  1:1479477/150‑1 (MQ=255)
                                                                                                                                      cacCGCAGAGATGTTACATTGATGCGGTGTATGAAGAATAACCAGAATGGTTAATAGGCGTTAGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCACCCTGCCAGTAGAGCTGAGCTGCAACGCGGTCAGCCAGCTGGCGGTAGAt                >  1:1280359/1‑150 (MQ=255)
                                                                                                                                      cacCGCAGAGATGTTACATTGATGCGGTGTATGAAGAATAACCAGAATGGTTAATAGGCGTTAGACCGCGCGGAAGGAAATCTCACCTGGAATGACTTCACCCTGCCAGTAGAGCTGAGCTGCAACGCGGTCAGCCAGCTGGCGGTAGAt                >  1:536526/1‑150 (MQ=255)
                                                                                                                                        ccGCATCAATGTAACATCTCTGCGGTGTATGAAGAATAACCAGAATGGTTAATAGGCGTTAGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCACCCTGCCAGTAGAGCTGAGCTGCAACGCGGTCAGCCAGCTGGCGGTAGATAg              <  1:859859/150‑1 (MQ=255)
                                                                                                                                         cGCATCAATGTAACATCTCTGCGGTGTATGAAGAATAACCAGAATGGTTAATAGGCGTTAGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCACCCTGCCAGTAGAGCTGAGCTGCAACGCGGTCAGCCAGCTGGCGGTAGATAGc             >  1:1766862/1‑150 (MQ=255)
                                                                                                                                           cATCAATGTAACATCTCTGCGGTGTATGAAGAATAACCAGAATGGTTAATAGGCGTTAGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCACCCTGCCAGTAGAGCTGAGCTGCAACGCGGTCAGCCAGCTGGCGGTAGATAGCgg           >  1:1504348/1‑150 (MQ=255)
                                                                                                                                               aaTGTAACATCTCTGCGGTGTATGAAGAATAACCAGAATGGTTAATAGGCGTTAGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCACCCTGCCAGTAGAGCTGAGCTGCAACGCGGTCAGCCAGCTGGCGGTAGATAGCGGTaaa       <  1:398378/150‑1 (MQ=255)
                                                                                                                                                aTGTAACATCTCTGCGGTGTATGAAGAATAACCAGAATGGTTAATAGGCGTTAGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCACCCTGCCAGTAGAGCTGAGCTGCAATGCGGTCAGCCAGCTGGCGGTAGATAGCGGTAAAt      >  2:1876663/1‑150 (MQ=255)
                                                                                                                                                  gTAACATCTCTGCGGTGTATGAAGAATAACCAGAATGGTTAATAGGCGTTAGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCACCCTGCCAGTAGAGCTGAGCTGCAACGCGGTCATCCAGCTGGCGGTAGATAGCGGTAAATTc    <  2:1410534/150‑1 (MQ=255)
                                                                                                                                                  gTAACATCTCTGCGGTGTATGAAGAATAACCAGAATGGTTAATAGGCGTTAGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCACCCTGCCAGTAGAGCTGAGCTGCAACGCGGTCAGCCAGCTGGCGGTAGATAGCGGTAAATTc    >  2:490311/1‑150 (MQ=255)
                                                                                                                                                    aaCATCTCTGCGGTGTATGAAGAATAACCAGAATGGTTAATAGGCGTTAGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCACCCTGCCAGTAGAGCTGAGCTGCAACGCGGTCAGCCAGCTGGCGGTAGATAGCGGTAAATTCGc  <  1:1314810/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TCCACTACACAGCCACGAAGGATGATAATGAAGCATTGCCTGTATGATCAATCGACTTTGTAGAGTTTCGGACGAAGGTCCGCAGAATATTCGCAGTATTAAATAAGTGTTCATTAATGCGGCAAAAACTAATACACCGCATCAATGTAACATCTCTGCGGTGTATGAAGAATAACCAGAATGGTTAATAGGCGTTAGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCACCCTGCCAGTAGAGCTGAGCTGCAACGCGGTCAGCCAGCTGGCGGTAGATAGCGGTAAATTCGC  >  NC_000913/2192865‑2193162

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