Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,193,018 C→T 6.9% intergenic (‑222/+41) 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,0180CT 6.9% 310.8 / 8.7 101intergenic (‑222/+41)yehE/mrpDUF2574 domain‑containing protein YehE/P‑loop NTPase family protein
Reads supporting (aligned to +/- strand):  ref base C (42/51);  new base T (5/2);  total (48/53)
Fisher's exact test for biased strand distribution p-value = 2.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.14e-01

Only 100 of 104 total aligned reads displayed.

CACAGCCACGAAGGATGATAATGAAGCATTGCCTGTATGATCAATCGACTTTGTAGAGTTTCGGACGAAGGTCCGCAGAATATTCGCAGTATTAAATAAGTGTTCATTAATGCGGCAAAAACTAATACACCGCATCAATGTAACATCTCTGCGGTGTATGAAGAATAACCAGAATGGTTAATAGGCGTTAGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCACCCTGCCAGTAGAGCTGAGCTGCAACGCGGTCAGCCAGCTGGCGGTAGATAGCGGTAAATTCGCTCTCC  >  NC_000913/2192872‑2193167
                                                                                                                                                  |                                                                                                                                                     
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:1211990/1‑150 (MQ=255)
             gatgatAATGAAGCATTGCCTGTATGATCAATCGACTTTGTAGAGTTTCGGACGAAGGTCCGCAGAATATTCGCAGTATTAAATAAGTGTTCATTAATGCGGCAAAAACTAATACACCGCATCAATGTAACATCTCTGCGGTGTATgaag                                                                                                                                       <  2:205640/150‑1 (MQ=255)
              atgatAATGAAGCATTGCCTGTATGATCAATCGACTTTGTAGAGTTTCGGACGAAGGTCCGCAGAATATTCGCAGTATTAAATAAGTGTTCATTAATGCGGCAAAAACTAATACACCGCATCAATGTAACATCTCTGCGGTGTATgaaga                                                                                                                                      >  2:1928631/1‑150 (MQ=255)
                    aTGAAGCATTGCCTGTATGATCAATCGACTTTGTAGAGTTTCGGACGAAGGTCCGCAGAATATTCGCAGTATTAAATAAGTGTTCATTAATGCGGCAAAAACTAATACACCGCATCAATGTAACATCTCTGCGGTGTATGAAGAATAAcc                                                                                                                                >  1:189151/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)
                              gCCTGTATGATCAATCGACTTTGTAGAGTTTCGGACGAAGGTCCGCAGAATATTCGCAGTATTAAATAAGTGTTCATTAATGCGGCAAAAACTAATACACCGCATCAATGTAACATCTCTGCGGTGTATGAAGAATAACCAGAATGGTTa                                                                                                                      >  2:1772335/1‑150 (MQ=255)
                                 tGTATGATCAATCGACTTTGTAGAGTTTCGGACGAAGGTCCGCAGAATATTCGCAGTATTAAATAAGTGTTCATTAATGCGGCAAAAACTAATACACCGCATCAATGTAACATCTCTGCGGTGTATGAAGAATAACCAGAATGGTTAATa                                                                                                                   <  2:1349318/150‑1 (MQ=255)
                                    aTGATCAATCGACTTTGTAGAGTTTCGGACGAAGGTCCGCAGAATATTCGCAGTATTAAATAAGTGTTCATTAATGCGGCAAAAACTAATACACCGCATCAATGTAACATCTCTGCGgtgt                                                                                                                                             <  2:132140/121‑1 (MQ=255)
                                    aTGATCAATCGACTTTGTAGAGTTTCGGACGAAGGTCCGCAGAATATTCGCAGTATTAAATAAGTGTTCATTAATGCGGCAAAAACTAATACACCGCATCAATGTAACATCTCTGCGgtgt                                                                                                                                             >  1:132140/1‑121 (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:264289/1‑150 (MQ=255)
                                               aCTTTGTAGAGTTTCGGACGAAGGTCCGCAGAATATTCGCAGTATTAAATAAGTGTTCATTAATGCGGCAAAAACTAATACACCGCATCAATGTAACATCTCTGCGGTGTATGAAGAATAACCAGAATGGTTAATAGGCGTTAGACcgcg                                                                                                     <  2:962470/150‑1 (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                                                                           >  2:1428916/1‑147 (MQ=255)
                                                                            aGAATATTCGCAGTATTAAATAAGTGTTCATTAATGCGGCAAAAACTAATACACCGCATCAATGTAACATCTCTGCGGTGTATGAAGAATAACCAGAATGGTTAATAGGCGTTAGACCGCGCGGAAGGAAATCTCGCCTGGAATGAc                                                                           <  1:1428916/147‑1 (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:1085844/150‑1 (MQ=255)
                                                                                     gCAGTATTAAATAAGTGTTCATTAATGCGGCAAAAACTAATACACCGCATCAATGTAACATCTCTGCGGTGTATGAAGAATAACCAGAATGGTTAATAGGCGTTAGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCACCCTGCCa                                                               >  2:1624147/1‑150 (MQ=255)
                                                                                      cAGTATTAAATAAGTGTTCATTAATGCGGCAAAAACTAATACACCGCATCAATGTAACATCTCTGCGGTGTATGAAGAATAACCAGAATGGTTAATAGGCGTTAGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCACCCTGCCAg                                                              <  1:1436125/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)
                                                                                             aaaTAAGTGTTCATTAATGCGGCAAAAACTAATACACCGCATCAATGTAACATCTCTGCGGTGTATGAAGAATAACCAGAATGGTTAATAGGCGTTAGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCACCCTGCCAGTAgagct                                                       >  1:637164/1‑150 (MQ=255)
                                                                                              aaTAAGTGTTCATTAATGCGGCAAAAACTAATACACCGCATCAATGTAACATCTCTGCGGTGTATGAAGAATAACCAGAATGGTTAATAGGCGTTAGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCACCCTGCCAGTAgagctg                                                      <  2:422907/150‑1 (MQ=255)
                                                                                                tAAGTGTTCATTAATGCGGCAAAAACTAATACACCGCATCAATGTAACATCTCTGCGGTGTATGAAGAATAACCAGAATGGTTAATAGGCGTTAGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCACCCTGCCAGTAgagctgag                                                    <  2:94082/150‑1 (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: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                        <  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)
                                                                                                                                           gTAACATCTCTGCGGTGTATGAAGAATAACCAGAATGGTTAATAGGCGTTAGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCACCCTGCCAGTAGAGCTGAGCTGCAACGCGGTCATCCAGCTGGCGGTAGATAGCGGTAAATTc         <  2:1410534/150‑1 (MQ=255)
                                                                                                                                           gTAACATCTCTGCGGTGTATGAAGAATAACCAGAATGGTTAATAGGCGTTAGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCACCCTGCCAGTAGAGCTGAGCTGCAACGCGGTCAGCCAGCTGGCGGTAGATAGCGGTAAATTc         >  2:490311/1‑150 (MQ=255)
                                                                                                                                             aaCATCTCTGCGGTGTATGAAGAATAACCAGAATGGTTAATAGGCGTTAGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCACCCTGCCAGTAGAGCTGAGCTGCAACGCGGTCAGCCAGCTGGCGGTAGATAGCGGTAAATTCGc       <  1:1314810/150‑1 (MQ=255)
                                                                                                                                                aTCTCTGCGGTGTATGAAGAATAACCAGAATGGTTAATAGGCGTTAGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCACCCTGCCAGTAGAGCTGAGCTGCAACGCGGTCAGCCAGCTGGCGGTAGATAGCGGTAAATTCGctct    >  1:1974587/1‑150 (MQ=255)
                                                                                                                                                aTCTCTGCGGTGTATGAAGAATAACCAGAATGGTTAATAGGCGTTAGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCACCCTGCCAGTAGAGCTGAGCTGCAACGCGGTCAGCCAGCTGGCGGTAGATAGCGGTAAATTCGctct    >  2:617113/1‑150 (MQ=255)
                                                                                                                                                  ctctGCGGTGTATGAAGAATAACCAGAATGGTTAATAGGCGTTAGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCACCCTGCCAGTAGAGCTGAGCTGCAACGCGGTCAGCCAGCTGGCGGTAGATAGCGGTAAATTCGCTCTcc  <  1:1772335/150‑1 (MQ=255)
                                                                                                                                                  ctctGCGGTGTATGAAGAATAACCAGAATGGTTAATAGGCGTTAGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCACCCTGCCAGTAGAGCTGAGCTGCAACGCGGTCAGCCAGCTGGCGGTAGATAGCGGTAAATTCGCTCTcc  <  2:1504348/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
CACAGCCACGAAGGATGATAATGAAGCATTGCCTGTATGATCAATCGACTTTGTAGAGTTTCGGACGAAGGTCCGCAGAATATTCGCAGTATTAAATAAGTGTTCATTAATGCGGCAAAAACTAATACACCGCATCAATGTAACATCTCTGCGGTGTATGAAGAATAACCAGAATGGTTAATAGGCGTTAGACCGCGCGGAAGGAAATCTCGCCTGGAATGACTTCACCCTGCCAGTAGAGCTGAGCTGCAACGCGGTCAGCCAGCTGGCGGTAGATAGCGGTAAATTCGCTCTCC  >  NC_000913/2192872‑2193167

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: