Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 14,119 G→A 7.7% intergenic (+40/‑49) 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,1190GA 7.7% 474.0 / 18.7 156intergenic (+40/‑49)dnaK/dnaJchaperone Hsp70, with co‑chaperone DnaJ/chaperone Hsp40, DnaK co‑chaperone
Reads supporting (aligned to +/- strand):  ref base G (69/74);  new base A (3/9);  total (72/84)
Fisher's exact test for biased strand distribution p-value = 1.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.61e-01

Only 100 of 162 total aligned reads displayed.

AGCAGCAACATGCCCAGCAGCAGACTGCCGGTGCTGATGCTTCTGCAAACAACGCGAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCGTGAAATCAGAAAGGCCTACAAACGCCTGGCCATGAAAT  >  NC_000913/13970‑14261
                                                                                                                                                     |                                                                                                                                              
 gcagcaACATGCCCAGCAGCAGACTGCCGGTGCTGATGCTTCTGCAAACAACGCGAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAggg                                                                                                                                               <  2:2102596/150‑1 (MQ=255)
 gcagcaACATGCCCAGCAGCAGACTGCCGGTGCTGATGCTTCTGCAAACAACGCGAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAggg                                                                                                                                               >  2:1567945/1‑150 (MQ=255)
  cagcaACATGCCCAGCAGCAGACTGCCGGTGCTGATGCTTCTGCAAACAACGCGAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAgggg                                                                                                                                              <  2:1771553/150‑1 (MQ=255)
      aaCATGCCCAGCAGCAGACTGCCGGTGCTGATGCTTCTGCAAACAACGCGAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAAtt                                                                                                                                          <  1:790154/150‑1 (MQ=255)
       aCATGCCCAGCAGCAGACTGCCGGTGCTGATGCTTCTGCAAACAACGCGAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAAttt                                                                                                                                         <  1:469951/150‑1 (MQ=255)
         aTGCCCAGCAGCAGACTGCCGGTGCTGATGCTTCTGCAAACAACGCGAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTcc                                                                                                                                       <  2:1252524/150‑1 (MQ=255)
         aTGCCCAGCAGCAGACTGCCGGTGCTGATGCTTCTGCAAACAACGCGAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTcc                                                                                                                                       <  2:1031868/150‑1 (MQ=255)
          tGCCCAGCAGCAGACTGCCGGTGCTGATGCTTCTGCAAACAACGCGAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCct                                                                                                                                      >  2:713566/1‑150 (MQ=255)
            cccAGCAGCAGACTGCCGGTGCTGATGCTTCTGCAAACAACGCGAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCctct                                                                                                                                    >  2:1430283/1‑150 (MQ=255)
             ccAGCAGCAGACTGCCGGTGCTGATGCTTCTGCAAACAACGCGAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGGAGAGGAAATTCccctt                                                                                                                                   <  2:284552/150‑2 (MQ=255)
                 cagcagACTGCCGGTGCTGATGCTTCTGCAAACAACGCGAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGGAGAGGAAATTCCCCTTCGcc                                                                                                                               <  2:151939/150‑1 (MQ=255)
                 cagcagACTGCCGGTGCTGATGCTTCTGCAAACAACGCGAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGGAGAGGAAATTCCCCTTCGcc                                                                                                                               <  2:2086167/150‑1 (MQ=255)
                 cagcagACTGCCGGTGCTGATGCTTCTGCAAACAACGCGAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGcc                                                                                                                               <  2:1866768/150‑1 (MQ=255)
                 cagcagACTGCCGGTGCTGATGCTTCTGCAAACAACGCGAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGcc                                                                                                                               <  1:2031945/150‑1 (MQ=255)
                 aagcagACTGGCTGTGCTGATGCTTCTGCAAACAACGAGAAATATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGCCACGGGCGGAGAGGAAATTCCCCTTCGcc                                                                                                                               <  2:142073/149‑1 (MQ=255)
                  agcagACTGCCGGTGCTGATGCTTCTGCAAACAACGCGAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGccc                                                                                                                              <  2:486030/150‑1 (MQ=255)
                   gcagACTGCCGGTGCTGATGCTTCTGCAAACAACGCGAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGGAGAGGAAATTCCCCTTCGCCCg                                                                                                                             <  2:256899/150‑1 (MQ=255)
                   gcagACTGCCGGTGCTGATGCTTCTGCAAACAACGCGAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGGAGAGGAAATTCCCCTTCGCCCg                                                                                                                             <  1:720423/150‑1 (MQ=255)
                    cagACTGCCGGTGCTGATGCTTCTGCAAACAACGCGAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGCCACGGCCGAAGGGGAATTACCTCACCGCCCGt                                                                                                                            >  1:849843/1‑150 (MQ=255)
                      gACTGCCGGTGCTGATGCTTCTGCAAACAACGCGAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGc                                                                                                                          <  2:1891697/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)
                          gCCGGTGCTGATGCTTCTGCAAACAACGCGAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTc                                                                                                                      >  2:345828/1‑150 (MQ=255)
                           ccGGTGCTGATGCTTCTGCAAACAACGCGAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCa                                                                                                                     >  2:1035533/1‑150 (MQ=255)
                               ggCTGATGCTTCTGCAAACAACGCGAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTa                                                                                                                 <  1:490820/149‑1 (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)
                                  tGATGCTTCTGCAAACAACGCGAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAggg                                                                                                              <  2:245217/150‑1 (MQ=255)
                                   gATGCTTCTGCAAACAACGCGAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAgggg                                                                                                             >  1:1006644/1‑150 (MQ=255)
                                     tGCTTCTGCAAACAACGCGAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACAGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCa                                                                                                           >  2:1524600/1‑150 (MQ=255)
                                         tCTGCAAACAACGCGACAGATGACGATGTTGGCGACACTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCACTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATGTCCTCTCCGCCCGTGCATTCATCTAGGGGCAAttt                                                                                                       <  2:383160/150‑1 (MQ=255)
                                          cTGCAAACAACGCGAAAGATGACGATGTTGTCGAAGCTGAATTTCAAGAAGTCAAAGACAAAAAATAAACGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCACCGCCCGTCCATTCATAAAGGCGCAATTTa                                                                                                      >  2:106359/1‑150 (MQ=255)
                                           tGCAAACAACGCGAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTaa                                                                                                     >  1:255340/1‑150 (MQ=255)
                                            gCAAACAACGCGAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTaaa                                                                                                    >  1:1700215/1‑150 (MQ=255)
                                            gCAAACAACGCGAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCGATAAGCGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATATAGGGGCAATTTaaa                                                                                                    >  2:425216/1‑150 (MQ=255)
                                             cAAACAACGCGAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAGGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGAGCGAATGGGAATATCCTCTCCGCCCGTGCATTCATATAGGGGCAATTTaaaa                                                                                                   >  1:1116685/1‑150 (MQ=255)
                                              aaaCAACGCGAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTaaaaa                                                                                                  <  1:1855089/150‑1 (MQ=255)
                                               aacaacGCGAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTaaaaaa                                                                                                 >  1:655516/1‑150 (MQ=255)
                                                acaacGCGAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGTGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTaaaaatg                                                                                                >  1:255474/1‑148 (MQ=255)
                                                    cgcgAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAGAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATgg                                                                                            <  2:1772216/150‑1 (MQ=255)
                                                    cgcgAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCTCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATgg                                                                                            <  1:1998309/150‑1 (MQ=255)
                                                      cgAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCTTTCATCTAGGGGCAATTTAAAAAAGATGGca                                                                                          >  1:415493/1‑149 (MQ=255)
                                                             tGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAg                                                                                   >  1:1482045/1‑150 (MQ=255)
                                                                    gttgtCGACTCTGAATTTGAATAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGAttatta                                                                            <  1:853715/150‑1 (MQ=255)
                                                                     ttgtCGACGCTGATTTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATTCTGACACGGTCGAATTGGAATTTCCTCTCAGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGGTAAGCAAGATTATTAc                                                                           <  2:1941945/150‑1 (MQ=255)
                                                                     ttgtCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTAc                                                                           >  2:1741267/1‑150 (MQ=255)
                                                                        tCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACgag                                                                        <  1:515595/150‑1 (MQ=255)
                                                                          gACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGTCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGAt                                                                      <  1:1095570/150‑1 (MQ=255)
                                                                           aCGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACTGGCGAAGGGGAATTTCCTCTCCGCACGTGCATTCATCTAGGGGCAATTTAAAAAAGATTGCTAATCAAGATTATTACGAGAtt                                                                     >  2:72731/1‑150 (MQ=255)
                                                                           aCGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGAtt                                                                     >  2:515861/1‑150 (MQ=255)
                                                                           aCGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGAtt                                                                     <  1:1567945/150‑1 (MQ=255)
                                                                              cTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTa                                                                  >  1:1404959/1‑150 (MQ=255)
                                                                                gAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAgg                                                                >  1:571040/1‑150 (MQ=255)
                                                                                gAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAgg                                                                <  1:1430283/150‑1 (MQ=255)
                                                                                gAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAgg                                                                <  1:1374262/150‑1 (MQ=255)
                                                                                gAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAgg                                                                >  2:862335/1‑150 (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)
                                                                                                aaGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCg                                                >  2:2162815/1‑150 (MQ=255)
                                                                                                    cAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAag                                            <  2:238973/150‑1 (MQ=255)
                                                                                                    cAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAag                                            >  1:1293736/1‑150 (MQ=255)
                                                                                                    cAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTACAAAACAGCGGaac                                            >  1:1759314/1‑149 (MQ=255)
                                                                                                        aaaTAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTACGCGTTTCCAAAACAGCGGAAGAGCg                                        <  2:117195/150‑1 (MQ=255)
                                                                                                         aataatCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCGt                                       <  1:862335/150‑1 (MQ=255)
                                                                                                          ataatCGCCCTATAAACGTGTAATTATACTGACACGGGCGAAGGGGAATTTGCTCTCCGCCCGTGCATTCATCTATGGGCAATTTATAAATTATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCGTg                                      <  2:55836/150‑1 (MQ=255)
                                                                                                            aatCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCGTGaa                                    <  2:655516/150‑1 (MQ=255)
                                                                                                             atCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCGTGaaa                                   <  2:1498463/150‑1 (MQ=255)
                                                                                                              tCGCCCTATAAACGGGTTATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCACGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGATATTTTAGGCGTTTCCAAAACAGCGGAACAGCATGAAAt                                  >  1:1769365/1‑150 (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)
                                                                                                                        aaCGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCGTGAAATCAGAAAGGcc                        >  2:775975/1‑150 (MQ=255)
                                                                                                                           gggTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGCTTATTACGAGATTTTAGGCGTTTCCAAACCAGCGGAAGAGCGTGAATTCAGAAAGGCCTAc                     >  1:1310450/1‑150 (MQ=255)
                                                                                                                           gggTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCGTGAAATCAGAAAGGCCTAc                     <  2:828365/150‑1 (MQ=255)
                                                                                                                            ggTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATATTAGGCGTTTCCAAAACAGCGGAAGAGCGTGAAATCAGAAAGGCCTACa                    <  1:1550794/150‑1 (MQ=255)
                                                                                                                              tAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCGTGAAATCAGAAAGGCCTACaaa                  <  2:1006644/150‑1 (MQ=255)
                                                                                                                              tAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCGTGAAATCAGAAAGGCCTACaaa                  <  1:1781796/150‑1 (MQ=255)
                                                                                                                                 ttATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCGTGAAATCAGAAAGGCCTACAAACGc               >  2:474890/1‑150 (MQ=255)
                                                                                                                                    taCTTACACGTGCGAAGGGGAATTACCTCTCCGCCGGTGCATTCATATAGGGGCAATATAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCGTGAAATCAGAAAGGCCTACAAACGCCTg            <  2:1456584/150‑1 (MQ=255)
                                                                                                                                     aCTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCGTGAAATCAGAAAGGCCTACAAACGCCTgg           <  1:1028234/150‑1 (MQ=255)
                                                                                                                                        tgcacGGGCGGAGAGGAAATTCCCCTTCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCGTGAAATCAGAAAGGCCTACAAACGCCTGGCCa        >  1:900920/3‑150 (MQ=255)
                                                                                                                                        gACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCGTGAAATCAGAAAGGCCTACAAACGCCTGGCCa        >  2:273068/1‑150 (MQ=255)
                                                                                                                                          cacGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCGTGAAATCAGAAAGGCCTACAAACGCCTGGCCATg      <  2:594084/150‑1 (MQ=255)
                                                                                                                                          cacGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCCAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGGGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCGTGAAATCAGAAAGGCCTACAAACGCCTGGCCATg      >  1:1682618/1‑150 (MQ=255)
                                                                                                                                           acGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCGTGAAATCAGAAAGGCCTACAAACGCCTGGCCATGa     <  1:2150934/150‑1 (MQ=255)
                                                                                                                                              ggCGGAGAGGAAATTCCCCTTCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCGTGAAATCAGAAAGGCCTACAAACGCCTGGCCATg      <  2:367724/146‑1 (MQ=255)
                                                                                                                                              ggCGGAGAGGAAATTCCCCTTCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCGTGAAATCAGAAAGGCCTACAAACGCCTGGCCATg      >  1:367724/1‑146 (MQ=255)
                                                                                                                                              ggCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCGTGAAATCAGAAAGGCCTACAAACGCCTGGCCATGAAAt  >  1:1651906/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                              
AGCAGCAACATGCCCAGCAGCAGACTGCCGGTGCTGATGCTTCTGCAAACAACGCGAAAGATGACGATGTTGTCGACGCTGAATTTGAAGAAGTCAAAGACAAAAAATAATCGCCCTATAAACGGGTAATTATACTGACACGGGCGAAGGGGAATTTCCTCTCCGCCCGTGCATTCATCTAGGGGCAATTTAAAAAAGATGGCTAAGCAAGATTATTACGAGATTTTAGGCGTTTCCAAAACAGCGGAAGAGCGTGAAATCAGAAAGGCCTACAAACGCCTGGCCATGAAAT  >  NC_000913/13970‑14261

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: