Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,193,008 A→G 6.7% intergenic (‑212/+51) 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,0080AG 6.7% 317.8 / 9.3 104intergenic (‑212/+51)yehE/mrpDUF2574 domain‑containing protein YehE/P‑loop NTPase family protein
Reads supporting (aligned to +/- strand):  ref base A (46/51);  new base G (5/2);  total (51/53)
Fisher's exact test for biased strand distribution p-value = 2.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01

Only 100 of 107 total aligned reads displayed.

TGTACATCCACTACACAGCCACGAAGGATGATAATGAAGCATTGCCTGTATGATCAATCGACTTTGTAGAGTTTCGGACGAAGGTCCGCAGAATATTCGCAGTATTAAATAAGTGTTCATTAATGCGGCAAAAACTAATACACCGCATCAATGTAACATCTCTGCGGTGTATGAAGAATAACCAGAATGGTTAATAGGCGTTAGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCACCCTGCCAGTAGAGCTGAGCTGCAACGCGGTCAGCCAGCTGGCGGTAGATAGCGGTAAA  >  NC_000913/2192859‑2193157
                                                                                                                                                     |                                                                                                                                                     
tGTACATCCACTACACAGCCACGAAGGATGATAATGAAGCATTGCCTGTATGATCAATCGACTTTGTAGAGTTTCGGACGAAGGTCCGCAGAATATTCGCAGTATTAAATAAGTGTTCATTAATGCGGCAAAAACTAATACACCGCATCa                                                                                                                                                       >  2:806325/1‑150 (MQ=255)
 gTACATCCACTACACAGCCACGAAGGATGATAATGAAGCATTGCCTGTATGATCAATCGACTTTGTAGAGTTTCGGACGAAGGTCCGCAGAATATTCGCAGTATTAAATAAGTGTTCATTAATGCGGCAAAAACTAATACACCGCATCaa                                                                                                                                                      >  2:563480/1‑150 (MQ=255)
      tCCACTACACAGCCACGAAGGATGATAATGAAGCATTGCCTGTATGATCAATCGACTTTGTAGAGTTTCGGACGAAGGTCCGCAGAATATTCGCAGTATTAAATAAGTGTTCATTAATGCGGCAAAAACTAATACACCGCATCAATGTaa                                                                                                                                                 >  1:296766/1‑150 (MQ=255)
       ccACTACACAGCCACGAAGGATGATAATGAAGCATTGCCTGTATGATCAATCGACTTTGTAGAGTTTCGGACGAAGGTCCGCAGAATATTCGCAGTATTAAATAAGTGTTCATTAATGCGGCAAAAACTAATACACCGCATCAATg                                                                                                                                                    >  1:361181/1‑146 (MQ=255)
       ccACTACACAGCCACGAAGGATGATAATGAAGCATTGCCTGTATGATCAATCGACTTTGTAGAGTTTCGGACGAAGGTCCGCAGAATATTCGCAGTATTAAATAAGTGTTCATTAATGCGGCAAAAACTAATACACCGCATCAATg                                                                                                                                                    <  2:361181/146‑1 (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)
                         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)
                                 aTGAAGCATTGCCTGTATGATCAATCGACTTTGTAGAGTTTCGGACGAAGGTCCGCAGAATATTCGCAGTATTAAATAAGTGTTCATTAATGCGGCAAAAACTAATACACCGCATCAATGTAACATCTCTGCGGTGTATGAAGAATAAcc                                                                                                                      >  1:189151/1‑150 (MQ=255)
                                 aTGAAGCATTGCCTGTATGATCAATCGACTTTGTAGAGTTTCGGACGAAGGTCCGCAGAATATTAGCAGTATTAAATAAGTGTTCATTAATGCGGCAAAAACTAATACACCGCATCAATGTAACATCTCTGCGGTGTATGAAGAATAAcc                                                                                                                      >  2:580740/1‑150 (MQ=255)
                                      gCATTGCCTGTATGATCAATCGACTTTGTAGAGTTTCGGACGAAGGTCCGCAGAATATTCGCAGTATTAAATAAGTGTTCATTAATGCGGCAAAAACTAATACACCGCATCAATGTAACATCTCTGCGGTGTATGAAGAATAACCAGAAt                                                                                                                 >  1:1545154/1‑150 (MQ=255)
                                       cATTGCCTGTATGATCAATCGACTTTGTAGAGTTTCGGACGAAGGTCCGCAGAATATTCGCAGTATTAAATAAGTGTTCATTAATGCGGCAAAAACTAATACACCGCATCAATGTAACATCTCTGCGGTGTATGAAGAATAACCAGAATg                                                                                                                >  2:1597936/1‑150 (MQ=255)
                                       cATTGCCTGTATGATCAATCGACTTTGTAGAGTTTCGGACGAAGGTCCGCAGAATATTCGCAGTATTAAATAAGTGTTCATTAATGCGGCAAAAACTAATACACCGCATCAATGTAACATCTCTGCGGTGTATGAAGAATAACCAGAATg                                                                                                                <  2:296766/150‑1 (MQ=255)
                                          tGCCTGTATGATCAATCGACTTTGTAGAGTTTCGGACGAAGGTCCGCAGAATATTCGCAGTATTAAATAAGTGTTCATTAATGCGGCAAAAACTAATACACCGCATCAATGTAACATCTCTGCGGTGTATGAAGAATAACCAGAATGGtt                                                                                                             <  1:1766475/150‑1 (MQ=255)
                                          tGCCTGTATGATCAATCGACTTTGTAGAGTTTCGGACGAAGGTCCGCAGAATATTCGCAGTATTAAATAAGTGTTCATTAATGCGGCAAAAACTAATACACCGCATCAATGTAACATCTCTGCGGTGTATGAAGAATAACCAGAATGGtt                                                                                                             <  2:684754/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                                                                          <  1:806325/150‑1 (MQ=255)
                                                                             aCGAAGGTCCGCAGAATATTCGCAGTATTAAATAAGTGTTCATTAATGCGGCAAAAACTAATACACCGCATCAATGTAACATCTCTGCGGTGTATGAAGAATAACCAGAATGGTTAATAGGCGTTAGACCGCGCGGAAGGAAATCTCGcc                                                                          <  2:802937/150‑1 (MQ=255)
                                                                              cGAAGGTCCGCAGAATATTCGCAGTATTAAATAAGTGTTCATTAATGCGGCAAAAACTAATACACCGCATCAATGTAACATCTCTGCGGTGTATGAAGAATAACCAGAATGGTTAATAGGCGTTAGACCGCGCGGAAGGAAATCTCGCCt                                                                         <  2:1459567/150‑1 (MQ=255)
                                                                              cGAAGGTCCGCAGAATATTCGCAGTATTAAATAAGTGTTCATTAATGCGGCAAAAACTAATACACCGCATCAATGTAACATCTCTGCGGTGTATGAAGAATAACCAGAATGGTTAATAGGCGTTAGACCGCGCGGAAGGAAATCTCGCCt                                                                         <  1:982952/150‑1 (MQ=255)
                                                                              cGAAGGTCCGCAGAATATTCGCAGTATTAAATAAGTGTTCATTAATGCGGCAAAAACTAATACACCGCATCAATGTAACATCTCTGCGGTGTATGAAGAATAACCAGAATGGTTAATAGGCGTTAGACCGCGCGGAAGGAAATCTCGCCt                                                                         >  1:400656/1‑150 (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:1548196/150‑1 (MQ=255)
                                                                                                 cGCAGTATTAAATAAGTGTTCATTAATGCGGCAAAAACTAATACACCGCATCAATGTAACATCTCTGCGGTGTATGAAGAATAACCAGAATGGTTAATAGGCGTTAGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCACCCTGcc                                                      <  2:1545154/150‑1 (MQ=255)
                                                                                                  gCAGTATTAAATAAGTGTTCATTAATGCGGCAAAAACTAATACACCGCATCAATGTAACATCTCTGCGGTGTATGAAGAATAACCAGAATGGTTAATAGGCGTTAGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCACCCTGCCa                                                     >  2:1624147/1‑150 (MQ=255)
                                                                                                  gCAGTATTAAATAAGTGTTCATTAATGCGGCAAAAACTAATACACCGCATCAATGTAACATCTCTGCGGTGTATGAAGAATAACCAGAATGGTTAATAGGCGTTAGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCACCCTGCCa                                                     <  2:1085844/150‑1 (MQ=255)
                                                                                                   cAGTATTAAATAAGTGTTCATTAATGCGGCAAAAACTAATACACCGCATCAATGTAACATCTCTGCGGTGTATGAAGAATAACCAGAATGGTTAATAGGCGTTAGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCACCCTGCCAg                                                    <  1:1436125/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)
                                                                                                                gtgtTCATTAATGCGGCAAAAACTAATACACCGCATCAATGTAACATCTCTGCGGTGTATGAAGAATAACCAGAATGGTTAATAGGCGTTAGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCACCCTGCCAGTAGagctgagctg                                       <  2:338095/150‑1 (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:1740107/1‑150 (MQ=255)
                                                                                                                                   aaaCTAATACACCGCATCAATGTAACATCTCTGCGGTGTATGAAGAATAACCAGAATGGTTAATAGGCGTTAGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCACCCTGCCAGTAGAGCTGAGCTGCAACGCGGTCAGCCAGCTg                    >  1:151190/1‑150 (MQ=255)
                                                                                                                                   aaaCTAATACACCGCATCAATGTAACATCTCTGCGGTGTATGAAGAATAACCAGAATGGTTAATAGGCGTTAGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCACCCTGCCAGTAGAGCTGAGCTGCAACGCGGTCAGCCAGCTg                    <  1:264289/150‑1 (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              <  2:400656/150‑1 (MQ=255)
                                                                                                                                         aTACACCGCATCAATGTAACATCTCTGCGGTGTATGAAGAATAACCAGAATGGTTAATAGGCGTTAGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCACCCTGCCAGTAGAGCTGAGCTGCAACGCGGTCAGCCAGCTGGCGGTa              <  1:1718102/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)
                                                                                                                                                     |                                                                                                                                                     
TGTACATCCACTACACAGCCACGAAGGATGATAATGAAGCATTGCCTGTATGATCAATCGACTTTGTAGAGTTTCGGACGAAGGTCCGCAGAATATTCGCAGTATTAAATAAGTGTTCATTAATGCGGCAAAAACTAATACACCGCATCAATGTAACATCTCTGCGGTGTATGAAGAATAACCAGAATGGTTAATAGGCGTTAGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCACCCTGCCAGTAGAGCTGAGCTGCAACGCGGTCAGCCAGCTGGCGGTAGATAGCGGTAAA  >  NC_000913/2192859‑2193157

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: