Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,960,926 G→C 6.8% G501G (GGC→GGG) ‡ recC ← exodeoxyribonuclease V subunit RecC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,960,9260GC 6.8% 442.0 / 17.1 133G501G (GGC→GGG) ‡recCexodeoxyribonuclease V subunit RecC
Reads supporting (aligned to +/- strand):  ref base G (56/67);  new base C (5/4);  total (61/72)
Fisher's exact test for biased strand distribution p-value = 7.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.95e-01

Only 100 of 134 total aligned reads displayed.

ATAAGGTAGAACCGATTGCCACTCGCCCTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGACACAAAACGACTGTCAGGCAGTGATAACAGGCTGATAAACGCTTC  >  NC_000913/2960779‑2961075
                                                                                                                                                   |                                                                                                                                                     
  aaGGTAGAACCGATTGCCACTCGCCCTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATTTATGCCCCAACGAATGCCgg                                                                                                                                                   <  1:2410717/150‑1 (MQ=255)
    ggTAGAACCGATTGCCACTCGCCCTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGAt                                                                                                                                                 <  1:497046/150‑1 (MQ=255)
       aGAACCGATTGCCACTGGCCCTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCg                                                                                                                                              <  2:1525633/150‑1 (MQ=255)
         aaCCGATTGCCACTCGCCCTGCGCGCACTCCATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTAGAACTCTGGCGTCGACTTAAAAATAAAAAAAATTGTCTTCATCTATGCCCCAACGAATGCCGGATTCCtt                                                                                                                                            >  2:2345708/1‑150 (MQ=255)
          aCCGATTGCCACTCGCCCTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTg                                                                                                                                           >  1:833743/1‑150 (MQ=255)
          aCCGATTGCCACTCGCCCTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTGGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTg                                                                                                                                           <  2:148491/150‑1 (MQ=255)
           ccGATTGCCACTCGCCCTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATCCGGCATTCGTTgg                                                                                                                                          <  2:1110619/150‑2 (MQ=255)
             gATTGCCACTCGCCCTGCGCGCGCTCCATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGGTGTCCGGTGGCGGGGAGGTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATCCGGCATTCGTTgggg                                                                                                                                        <  2:1625081/150‑4 (MQ=255)
                 gCCACTCGCCCTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACt                                                                                                                                    <  1:595421/150‑1 (MQ=255)
                 gCCACTCGCCCTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACt                                                                                                                                    <  2:1518781/150‑1 (MQ=255)
                 gCCACTCGCCCTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACCTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACt                                                                                                                                    >  2:862834/1‑150 (MQ=255)
                 gCCACTCGCCCTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGACTTACAGCTAAATAAAGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACt                                                                                                                                    >  2:1328967/1‑150 (MQ=255)
                   cACTCGCCCTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATTGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTgg                                                                                                                                  <  2:204074/150‑1 (MQ=255)
                   cACTCGCCCTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTgg                                                                                                                                  >  1:1157493/1‑150 (MQ=255)
                   cACTCGCCCTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTgg                                                                                                                                  <  2:2226387/150‑1 (MQ=255)
                   cACTAGCCCTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTgg                                                                                                                                  <  2:2235438/150‑1 (MQ=255)
                    aCTCGCCCTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGc                                                                                                                                 >  1:811548/1‑150 (MQ=255)
                    aCTCGCCCTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGc                                                                                                                                 >  1:1693795/1‑150 (MQ=255)
                    aCTCGCCCTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGc                                                                                                                                 >  2:524905/1‑150 (MQ=255)
                      tCGCCCTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTACAGCTAGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGCCCCACTGGCGt                                                                                                                               >  1:741736/1‑150 (MQ=255)
                       cGCCCTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTa                                                                                                                              <  1:298282/150‑1 (MQ=255)
                        gCCCTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTaa                                                                                                                             >  1:839155/1‑150 (MQ=255)
                         cccTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTaaa                                                                                                                            <  1:2444846/150‑1 (MQ=255)
                         cccTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTaaa                                                                                                                            >  1:1026743/1‑150 (MQ=255)
                          ccTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAAt                                                                                                                           <  2:1614394/150‑1 (MQ=255)
                          ccTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGGGTGTTGTCCGGTGGCGGGGAGTGTCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAAt                                                                                                                           <  1:946139/150‑1 (MQ=255)
                               gcgcTCTCCATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTg                                                                                                                                   <  1:2008754/137‑1 (MQ=255)
                               gcgcTCTCCATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGc                                                                                                                      >  2:971861/1‑150 (MQ=255)
                               gcgcTCTCCATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGc                                                                                                                      >  1:1605155/1‑150 (MQ=255)
                                             gTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGt                                                                                                        >  1:260404/1‑150 (MQ=255)
                                             gTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGt                                                                                                        <  2:2098830/150‑1 (MQ=255)
                                              tAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTg                                                                                                       <  1:1988329/150‑1 (MQ=255)
                                               aGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGa                                                                                                      <  1:1067737/150‑1 (MQ=255)
                                               aGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGa                                                                                                      <  2:1683631/150‑1 (MQ=255)
                                                 cccAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATg                                                                                                    <  2:610445/150‑1 (MQ=255)
                                                 cccAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATg                                                                                                    <  1:1394543/150‑1 (MQ=255)
                                                   caacaaCATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGAAAATAACGCAGCCCTTCTTCGGTGATGTc                                                                                                  <  2:1449031/150‑1 (MQ=255)
                                                    aacaacATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCa                                                                                                 <  2:1327483/150‑1 (MQ=255)
                                                     acaacaTACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCaa                                                                                                >  1:1666091/1‑150 (MQ=255)
                                                           tACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACcgcg                                                                                          <  1:543292/150‑1 (MQ=255)
                                                             cgcgTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGcgcc                                                                                        <  1:1484601/150‑1 (MQ=255)
                                                                 tCAGGCCAAATCGCCAGGTGTGTAGTACGCTGTCGTCGACTTAAAAATAAAAAAAGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCa                                                                                    >  2:1899780/1‑150 (MQ=255)
                                                                  cAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAg                                                                                   >  1:1389834/1‑150 (MQ=255)
                                                                  cAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAg                                                                                   <  2:2071614/150‑1 (MQ=255)
                                                                     gCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCAc                                                                                >  1:431114/1‑150 (MQ=255)
                                                                     gCCAAATCGCCAGGTGTGTTGTACTCTGTCGTCGACTTAAAACTAGCAAAAATTGTCGTCATCTAAGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCAc                                                                                >  2:1372434/1‑150 (MQ=255)
                                                                      ccAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCAcc                                                                               <  2:35308/150‑1 (MQ=255)
                                                                      ccAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCAcc                                                                               >  1:2086434/1‑150 (MQ=255)
                                                                      ccAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCAcc                                                                               <  2:667341/150‑1 (MQ=255)
                                                                       cAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCg                                                                              <  1:31058/150‑1 (MQ=255)
                                                                         aaTCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGc                                                                            <  2:839155/150‑1 (MQ=255)
                                                                         aaTCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCATCTTCGGAGATGTCAAACCGCGCCGCCAGCACCCGc                                                                            >  2:187440/1‑150 (MQ=255)
                                                                             gCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACAt                                                                        >  1:555308/1‑150 (MQ=255)
                                                                                   tgtgtTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCa                                                                  >  1:653996/1‑150 (MQ=255)
                                                                                      gttgtCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACg                                                               >  2:1427188/1‑150 (MQ=255)
                                                                                      gttgtCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACg                                                               >  1:72489/1‑150 (MQ=255)
                                                                                      gttgtCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACg                                                               >  1:930670/1‑150 (MQ=255)
                                                                                          tCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTTGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAg                                                           <  1:816960/150‑1 (MQ=255)
                                                                                          tCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAg                                                           <  1:119687/150‑1 (MQ=255)
                                                                                            cGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGcc                                                         >  2:1181899/1‑149 (MQ=255)
                                                                                                ggCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATc                                                     >  1:828633/1‑150 (MQ=255)
                                                                                                    gggAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCctct                                                 <  1:1429791/150‑1 (MQ=255)
                                                                                                     ggAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCTGCACATCCAGCAACGCCAGCACATCCTCTg                                                <  1:862834/150‑1 (MQ=255)
                                                                                                     ggAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTg                                                <  2:942921/150‑1 (MQ=255)
                                                                                                      gAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGa                                               <  2:409083/150‑1 (MQ=255)
                                                                                                      gAGTTCCAGCTCGCGAACGTTGTCGTAATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGa                                               <  2:202744/150‑1 (MQ=255)
                                                                                                       aGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGac                                              <  2:1157493/150‑1 (MQ=255)
                                                                                                       aGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGac                                              <  2:833743/150‑1 (MQ=255)
                                                                                                       aGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGac                                              <  1:2242025/150‑1 (MQ=255)
                                                                                                       aGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGac                                              <  1:1728092/150‑1 (MQ=255)
                                                                                                        gTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAg                                                           >  1:2394322/1‑136 (MQ=255)
                                                                                                        gTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAg                                                           <  2:2394322/136‑1 (MQ=255)
                                                                                                           ccAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGACACaa                                          >  1:2263794/1‑150 (MQ=255)
                                                                                                            cAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGACACaaa                                         >  1:1022585/1‑150 (MQ=255)
                                                                                                            cAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGACACaaa                                         >  1:1217271/1‑150 (MQ=255)
                                                                                                              gctggggAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAAGGCAGCCCTTCTTCGGGGATGTCAAACCGCGCCGCCAGCTCCGGCACATCCAGCAACGCCAGCACATCGTCTGACACAAAAc                                       <  2:100258/144‑1 (MQ=255)
                                                                                                               cTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGACACAAAACg                                      >  2:2273044/1‑150 (MQ=255)
                                                                                                                  gcgAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGACACAAAAc                                       >  2:411546/1‑146 (MQ=255)
                                                                                                                  gcgAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGACACAAAAc                                       <  1:411546/146‑1 (MQ=255)
                                                                                                                   cgAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGACACAAAACGACTg                                  >  2:663347/1‑150 (MQ=255)
                                                                                                                    gAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGACACAAAACGACTGt                                 >  2:2251863/1‑150 (MQ=255)
                                                                                                                            tcgtcATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCCGCAACGCCAGCACATCCTCTGCCACAAAACGACTGTCAGGCAGt                         >  1:2399056/1‑150 (MQ=255)
                                                                                                                             cgtcATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGACACAAAACGACTGTCAGGCAGTg                        >  2:530892/1‑150 (MQ=255)
                                                                                                                             cgtcATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGACACAAAACGACTGTCAGGCAGTg                        <  2:1125844/150‑1 (MQ=255)
                                                                                                                              gtcATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACCTCCAGCAACGCCAGCACATCCTCTGACACAAAACGACTGTCAGGCAGTGa                       >  1:1942739/1‑150 (MQ=255)
                                                                                                                              gtcATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGACACAAAACGACTGTCAGGCAGTGa                       <  2:1418643/150‑1 (MQ=255)
                                                                                                                               gcATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGACACAAAACGACTGTCAGGCAGTGAt                      <  2:2071607/149‑1 (MQ=255)
                                                                                                                                cATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGACACAAAACGACTGTCAGGCAGTGATa                     <  1:1328967/150‑1 (MQ=255)
                                                                                                                                cATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGACACAAAACGACTGTCAGGCAGTGATa                     <  1:2400782/150‑1 (MQ=255)
                                                                                                                                  tCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGACACAAAACGACTGTCAGGCAGTGATAAc                   >  1:2319095/1‑150 (MQ=255)
                                                                                                                                      tGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGACACAAAACGACTGTCAGGCAGTGATAACAGGc               <  2:72489/150‑1 (MQ=255)
                                                                                                                                       gCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGACACAAAACGACTGTCAGGCAGTGATAACAGGCt              >  1:894733/1‑150 (MQ=255)
                                                                                                                                        ccccAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGACACAAAACGACTGTCAGGCAGTGATAACAGGCTg             >  1:323887/1‑150 (MQ=255)
                                                                                                                                        ccccAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGACACAAAACGACTGTCAGGCAGTGATAACAGGCTg             >  1:680060/1‑150 (MQ=255)
                                                                                                                                         gtcAACGAATCCGGCATTCGTTGCCCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCCCATCCAGCAACGCCAGCACATCCTCTGACACAAAACGACTGTCAGGCAGGGATAACAGGATGa            >  2:2016910/3‑150 (MQ=255)
                                                                                                                                         gtcAACGAATCCGGCATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGACACAAAACGACTGTCAGGCAGTGATAACAGGCTGa            >  2:1167845/3‑150 (MQ=255)
                                                                                                                                          ccAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGACACAAAACGACTGTCAGGCAGTGATAACAGGCTGAt           >  1:1897434/1‑150 (MQ=255)
                                                                                                                                           cAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGACACAAAACGACTGTCAGCCTGTTATCACTGCCTGaca          <  1:864825/150‑3 (MQ=255)
                                                                                                                                               gAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGACACAAAACGACTGTCAGGCAGTGATAACAGGCTGATAAACg      <  2:1605155/150‑1 (MQ=255)
                                                                                                                                                   gCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGACACAAAACGACTGTCAGGCAGTGATAACAGGCTGATAAACGCTTc  <  1:1427188/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
ATAAGGTAGAACCGATTGCCACTCGCCCTGCGCGCTCTCCATCGCGTAGCCCAACAACATACGCGTCAGGCCAAATCGCCAGGTGTGTTGTCCGGTGGCGGGGAGTTCCAGCTCGCGAACGTTGTCGTCATCTATGCCCCAACGAATGCCGGATTCGTTGACCCACTGGCGTAAATAACGCAGCCCTTCTTCGGTGATGTCAAACCGCGCCGCCAGCACCGGCACATCCAGCAACGCCAGCACATCCTCTGACACAAAACGACTGTCAGGCAGTGATAACAGGCTGATAAACGCTTC  >  NC_000913/2960779‑2961075

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