Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 14,132 C→T 7.1% intergenic (+53/‑36) dnaK → / → dnaJ chaperone Hsp70, with co‑chaperone DnaJ/chaperone Hsp40, DnaK co‑chaperone

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091314,1320CT 7.1% 477.2 / 17.2 155intergenic (+53/‑36)dnaK/dnaJchaperone Hsp70, with co‑chaperone DnaJ/chaperone Hsp40, DnaK co‑chaperone
Reads supporting (aligned to +/- strand):  ref base C (70/74);  new base T (3/8);  total (73/82)
Fisher's exact test for biased strand distribution p-value = 2.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.50e-01

Only 100 of 158 total aligned reads displayed.

AGCAGCAGACTGCCGGTGCTGATGCTTCTGCAAACAACGCGAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCGTGAAATCAGAAAGGCCTACAAACGCCTGGCCATGAAATACCACCCGGACCGTAACCAG  >  NC_000913/13985‑14281
                                                                                                                                                   |                                                                                                                                                     
agcagcagACTGCCGCTGCTGATGCTTCTGCAAACAACGCGAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCg                                                                                                                                                     <  1:243989/150‑1 (MQ=255)
  cagcagACTGCCGGTGCTGATGCTTCTGCAAACAACGCGAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGGAGAGGAAATTCCCCTTCGcc                                                                                                                                                   <  2:151939/150‑1 (MQ=255)
  cagcagACTGCCGGTGCTGATGCTTCTGCAAACAACGCGAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGcc                                                                                                                                                   <  2:1866768/150‑1 (MQ=255)
  aagcagACTGGCTGTGCTGATGCTTCTGCAAACAACGAGAAATATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGCCACGGGCGGAGAGGAAATTCCCCTTCGcc                                                                                                                                                   <  2:142073/149‑1 (MQ=255)
   agcagACTGCCGGTGCTGATGCTTCTGCAAACAACGCGAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGccc                                                                                                                                                  >  2:701/1‑150 (MQ=255)
    gcagACTTCCGGTGCTGATGCTTCTGCAAACAACGCGAAAGATGACGATGTAGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATAATGACACGGGCGAAGGGGAATTTCCTCTCCGCCCg                                                                                                                                                 <  2:1873253/150‑1 (MQ=255)
    gcagACTGCCGGTGCTGATGCTTCTGCAAACAACGCGAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGGAGAGGAAATTCCCCTTCGCCCg                                                                                                                                                 <  2:1547257/150‑1 (MQ=255)
    gcagACTGCCGGTGCTGATGCTTCTGCAAACAACGCGAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGGAGAGGAAATTCCCCTTCGCCCg                                                                                                                                                 <  2:256899/150‑1 (MQ=255)
     cagACTGCCGGTGCTGATGCTTCTGCAAACAACGCGAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGGAGAGGAAATTCCCCTTCGCCCGt                                                                                                                                                <  2:245619/150‑1 (MQ=255)
       gACTGCCGGTGCTGATGCTTCTGCAAACAACGCGAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGc                                                                                                                                              <  1:701921/150‑1 (MQ=255)
       gACTGCCGGTGCTGATGCTTCTGCAAACAACGCGAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGc                                                                                                                                              <  2:1618994/150‑1 (MQ=255)
         cTGCCGGTGCTGATGCTTCTGCAAACAACGCGAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCAt                                                                                                                                            >  1:1966381/1‑150 (MQ=255)
          tGCCGGTGCTGATGCTTCTGCAAACAACGCGAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCAtt                                                                                                                                           <  1:278200/150‑1 (MQ=255)
            ccGGTGCTGATGCTTCTGCAAACAACGCGAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCa                                                                                                                                         >  2:1035533/1‑150 (MQ=255)
              ggTGCTGATGCTTCTGCAAACAACGCGAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATc                                                                                                                                       <  2:1478083/150‑1 (MQ=255)
              ggTGCTGATGCTTCTGCAAACAACGCGAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATc                                                                                                                                       >  2:515595/1‑150 (MQ=255)
                tGCTGATGCTTCTGCAAACAACGCGAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTa                                                                                                                                     >  1:1447854/1‑150 (MQ=255)
                 gcagATGCTTCTGCAAACAACGCGAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAg                                                                                                                                    >  2:653166/4‑150 (MQ=255)
                  cTGATGCTTCTGCAAACAACGCGAAAGCTGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGTCACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAgg                                                                                                                                   <  1:297789/150‑1 (MQ=255)
                  cTGATGCTTCTGCAAACAACGCGAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCTTGCATTCATCTAgg                                                                                                                                   >  1:2144994/1‑150 (MQ=255)
                  cTGATGCTTCTGCAAACAACGCGAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAgg                                                                                                                                   >  2:1717958/1‑150 (MQ=255)
                   tGATGCTTCTGCAAACAACGCGAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAggg                                                                                                                                  >  1:1206750/1‑150 (MQ=255)
                    gATGCTTCTGCAAACAACGCGAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAgggg                                                                                                                                 >  1:1006644/1‑150 (MQ=255)
                       gCTTCTGCAAACAACGCGAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCaa                                                                                                                              >  1:828365/1‑150 (MQ=255)
                       gCTTCTGCAAACAACGCGAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCaa                                                                                                                              >  1:2112354/1‑150 (MQ=255)
                         ttCTGCAAACAACGCGAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAAtt                                                                                                                            <  2:772528/150‑1 (MQ=255)
                          tCTGCAAACAACGCGACAGATGACGATGTTGGCGACACTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCACTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATGTCCTCTCCGCCCGTGCATTCATCTAGGGGCAAttt                                                                                                                           <  2:383160/150‑1 (MQ=255)
                           cTGCAAACAACGCGAAAGATGACGATGTTGTCGAAGCTGAATTTCAAGAAGTCAAAGACAAAAAATAAACGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCACCGCCCGTCCATTCATAAAGGCGCAATTTa                                                                                                                          >  2:106359/1‑150 (MQ=255)
                            tGCAAACAACGCGAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTaa                                                                                                                         <  2:628705/150‑1 (MQ=255)
                            tGCAAACAACGCGAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTaa                                                                                                                         >  1:255340/1‑150 (MQ=255)
                             gCAAACAACGCGAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTGATCTAGGGGCAATTTaaa                                                                                                                        >  2:575534/1‑150 (MQ=255)
                             gCAAACAACGCGAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTaaa                                                                                                                        >  1:1700215/1‑150 (MQ=255)
                             gCAAACAACGCGAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCGATAAGCGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATATAGGGGCAATTTaaa                                                                                                                        >  2:425216/1‑150 (MQ=255)
                                 acaacGCGAAAGATGACGATGTTGTCGACACTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAg                                                                                                                    >  2:471616/1‑150 (MQ=255)
                                     cgcgAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAGAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATgg                                                                                                                <  2:1772216/150‑1 (MQ=255)
                                     cgcgAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATgg                                                                                                                >  2:2211457/1‑150 (MQ=255)
                                             aTGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCaa                                                                                                        <  1:1231864/150‑1 (MQ=255)
                                                  gATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGatta                                                                                                   <  2:289864/150‑1 (MQ=255)
                                                     gttgtCGACTCTGAATTTGAATAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGAttatta                                                                                                <  1:853715/150‑1 (MQ=255)
                                                      ttgtCGACGCTGATTTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATTCTGACACGGTCGAATTGGAATTTCCTCTCAGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGGTAAGCAAGATTATTAc                                                                                               <  2:1941945/150‑1 (MQ=255)
                                                      ttgtCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTAc                                                                                               >  2:1741267/1‑150 (MQ=255)
                                                       tgtCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACg                                                                                              >  1:1774183/1‑150 (MQ=255)
                                                         tCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACgag                                                                                            <  1:515595/150‑1 (MQ=255)
                                                          cGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACgaga                                                                                           >  2:540908/1‑150 (MQ=255)
                                                           gACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGTCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGAt                                                                                          <  1:1095570/150‑1 (MQ=255)
                                                            aCGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGAtt                                                                                         <  1:1567945/150‑1 (MQ=255)
                                                               cTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTa                                                                                      >  1:1404959/1‑150 (MQ=255)
                                                                tGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAg                                                                                     >  2:1503898/1‑150 (MQ=255)
                                                                 gAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAgg                                                                                    >  2:862335/1‑150 (MQ=255)
                                                                 gAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAgg                                                                                    <  1:1374262/150‑1 (MQ=255)
                                                                 gAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAgg                                                                                    <  1:1430283/150‑1 (MQ=255)
                                                                     ttGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGCCACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGtt                                                                                >  1:599117/1‑150 (MQ=255)
                                                                       gaagaaGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTc                                                                              >  2:1649782/1‑150 (MQ=255)
                                                                        aagaagTCAAACACAAAAAATAATCGCCCTATACACGGGTAATTATACTGACACGGGCGAAGGGGAATATCCTCTCCGCCCGTGCAATCATCAAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTcc                                                                             <  2:300959/150‑1 (MQ=255)
                                                                         agaagTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTCGGGGCAATTTAAAAAAGATGGCTAACCAAGATTATTACGAGATTTTATGCGTTTCCa                                                                            >  2:549147/1‑150 (MQ=255)
                                                                         agaagTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCa                                                                            >  1:1506535/1‑150 (MQ=255)
                                                                             gTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTTGGCGTTTCCAAAAc                                                                        >  2:146192/1‑150 (MQ=255)
                                                                             gTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAAc                                                                        <  2:1321115/150‑1 (MQ=255)
                                                                             gTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAAc                                                                        >  1:19259/1‑150 (MQ=255)
                                                                                 aaGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCg                                                                    <  1:2014031/150‑1 (MQ=255)
                                                                                     cAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAag                                                                <  2:238973/150‑1 (MQ=255)
                                                                                     cAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTACAAAACAGCGGaac                                                                >  1:1759314/1‑149 (MQ=255)
                                                                                       aaaaaTAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAagag                                                              <  1:430802/150‑1 (MQ=255)
                                                                                         aaaTAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTACGCGTTTCCAAAACAGCGGAAGAGCg                                                            <  2:117195/150‑1 (MQ=255)
                                                                                          aataatCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCGt                                                           <  1:862335/150‑1 (MQ=255)
                                                                                               tCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTGTCCGCCCGTGCATTCAACGAGGGGCAATTTAAAAAAGATGGCTAAGCTAGAGTATTACGAGATTTTAGGCGTTTACAAAACAGCGGACGAGCGGGAAAt                                                      >  2:1175794/1‑150 (MQ=255)
                                                                                                 gCCCTATAAACGTGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCGTGAAATCa                                                    <  2:825160/150‑1 (MQ=255)
                                                                                                 gCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCGTGAAATCa                                                    <  1:1750606/150‑1 (MQ=255)
                                                                                                    cTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCGTGAAATCAGaa                                                 <  1:1524600/150‑1 (MQ=255)
                                                                                                      ataAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCGTGAAATCAGAAAg                                               <  2:1700215/150‑1 (MQ=255)
                                                                                                       taAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCGTGAAATCAGAAAgg                                              >  2:1265695/1‑150 (MQ=255)
                                                                                                       taAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCGTGAAATCAGAAAgg                                              <  1:600037/150‑1 (MQ=255)
                                                                                                         aaCGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCGTGAAATCAGAAAGGcc                                            >  2:775975/1‑150 (MQ=255)
                                                                                                         aaCGGGTAAATATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTAATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTAGAAAACAGCGGAAGAGCGTGAAATCAGAAAGGcc                                            >  1:663438/1‑150 (MQ=255)
                                                                                                            gggTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGCTTATTACGAGATTTTAGGCGTTTCCAAACCAGCGGAAGAGCGTGAATTCAGAAAGGCCTAc                                         >  1:1310450/1‑150 (MQ=255)
                                                                                                            gggTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCGTGAAATCAGAAAGGCCTAc                                         <  2:828365/150‑1 (MQ=255)
                                                                                                             ggTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCGTGAAATCAGAAAGGCCTACa                                        >  1:720897/1‑150 (MQ=255)
                                                                                                             ggTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCGTGAAATCAGAAAGGCCTACa                                        <  1:471616/150‑1 (MQ=255)
                                                                                                             ggTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATATTAGGCGTTTCCAAAACAGCGGAAGAGCGTGAAATCAGAAAGGCCTACa                                        <  1:1550794/150‑1 (MQ=255)
                                                                                                               tAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCGTGAAATCAGAAAGGCCTACaaa                                      <  1:1781796/150‑1 (MQ=255)
                                                                                                               tAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCGTGAAATCAGAAAGGCCTACaaa                                      <  2:1006644/150‑1 (MQ=255)
                                                                                                                  ttATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCGTGAAATCAGAAAGGCCTACAAACGc                                   >  2:474890/1‑150 (MQ=255)
                                                                                                                     taCTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCGTCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGAGTTTCCAAAACAGCGGAAGAGCGTGAAATCAGAAAGGCATACAAACGCCTg                                >  2:1438651/1‑150 (MQ=255)
                                                                                                                      aCTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCGTGAAATCAGAAAGGCCTACAAACGCCTgg                               <  1:1028234/150‑1 (MQ=255)
                                                                                                                        tGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCGTGAAATCAGAAAGGCCTACAAACGCCTGGcc                             <  1:515861/150‑1 (MQ=255)
                                                                                                                         tgcacGGGCGGAGAGGAAATTCCCCTTCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCGTGAAATCAGAAAGGCCTACAAACGCCTGGCCa                            >  1:900920/3‑150 (MQ=255)
                                                                                                                           cacGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCCAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGGGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCGTGAAATCAGAAAGGCCTACAAACGCCTGGCCATg                          >  1:1682618/1‑150 (MQ=255)
                                                                                                                               ggCGGAGAGGAAATTCCCCTTCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCGTGAAATCAGAAAGGCCTACAAACGCCTGGCCATg                          <  2:367724/146‑1 (MQ=255)
                                                                                                                               ggCGGAGAGGAAATTCCCCTTCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCGTGAAATCAGAAAGGCCTACAAACGCCTGGCCATg                          >  1:367724/1‑146 (MQ=255)
                                                                                                                               ggCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCGTGAAATCAGAAAGGCCTACAAACGCCTGGCCATGAAAt                      >  1:1651906/1‑150 (MQ=255)
                                                                                                                                       ggAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCGTGAAATCAGAAAGGCCTACAAACGCCTGGCCATGAAATACCACCCg              >  2:1801933/1‑150 (MQ=255)
                                                                                                                                        gAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCGTGAAATCAGAAAGGCCTACAAACGCCTGGCCATGAAATACCACCCgg             <  2:2112354/150‑1 (MQ=255)
                                                                                                                                            ttCCTCTCCGCGCGTCAATTCATCTTGGGTCAATTTAAAAATGATGGCTATGGAAGATTATTACGAGATTTTCGGCGTTTCCAAAACAGCGGAAGAGCGTGATAGCAGAATGGCCTACAAACGCCTTGCCATGAAATACCACCCGGACCg         <  1:106359/150‑1 (MQ=255)
                                                                                                                                              ccTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGTAAGAGCGTGAAATCAGAAAGGCCTACAAACGCCTGGCCATGAAATACCACCCGGACCGTa       >  1:1482927/1‑150 (MQ=255)
                                                                                                                                              ccTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCGTGAAATCAGAAAGGCCTACAAACGCCTGGCCATGAAATACCACCCGGACCGTa       >  2:1264331/1‑150 (MQ=255)
                                                                                                                                              ccTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCGTGAAATCAGAAAGGCCTACAAACGCCTGGCCATGAAATACCACCCGGACCGTa       <  2:1116685/150‑1 (MQ=255)
                                                                                                                                               ctctCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCGTGAAATCAGAAAGGCCTACAAACGCATGGCCATGAAATACCACCCGGACCGTaa      >  1:569271/1‑150 (MQ=255)
                                                                                                                                                   ccGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCGTGAAATCAGAAAGGCCTACAAACGCCTGGCCATGAAATACCACCCGGACCGTAACCAg  >  2:1993951/1‑150 (MQ=255)
                                                                                                                                                   ccGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCGTGAAATCAGAAAGGCCTACAAACGCCTGGCCATGAAATACCACCCGGACCGTAACCAg  >  1:9963/1‑150 (MQ=255)
                                                                                                                                                   ccGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCGTGAAATCAGAAAGGCCTACAAACGCCTGGCCATGAAATACCACCCGGAACGAAAACAg  >  2:58876/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
AGCAGCAGACTGCCGGTGCTGATGCTTCTGCAAACAACGCGAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCGTGAAATCAGAAAGGCCTACAAACGCCTGGCCATGAAATACCACCCGGACCGTAACCAG  >  NC_000913/13985‑14281

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