Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,515,565 A→G Y117H (TAC→CAC)  dam ← DNA adenine methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,515,5650AG97.8% 656.8 / 0.5 182Y117H (TAC→CAC) damDNA adenine methyltransferase
Reads supporting (aligned to +/- strand):  ref base A (3/1);  new base G (92/85);  total (96/86)
Fisher's exact test for biased strand distribution p-value = 6.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.93e-01

Only 100 of 182 total aligned reads displayed.

GCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCAT  >  NC_000913/3515417‑3515712
                                                                                                                                                    |                                                                                                                                                    
gCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGc                                                                                                                                                    >  4:213350/1‑151 (MQ=255)
    aaGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGtt                                                                                                                                                >  1:1294026/1‑151 (MQ=255)
    aaGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGt                                                                                                                                                 <  4:623276/150‑1 (MQ=255)
      gACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCgg                                                                                                                                                  >  4:549145/1‑147 (MQ=255)
      gACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCa                                                                                                                                              <  3:57676/151‑1 (MQ=255)
        cTCACAATAGAAAAAGGCATTCTGCGCTTTTTCACCGAAGCGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCa                                                                                                                                              <  1:1689359/149‑1 (MQ=255)
         tCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCaa                                                                                                                                             >  1:1388795/1‑149 (MQ=255)
          cacaATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAata                                                                                                                                          <  1:1755072/151‑1 (MQ=255)
             aaTAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAAtata                                                                                                                                        >  1:1522841/1‑150 (MQ=255)
                 gAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATaaaaa                                                                                                                                    >  4:1738961/1‑150 (MQ=255)
                 gAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATaaaaa                                                                                                                                    >  2:161990/1‑150 (MQ=255)
                  aaaaaGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACa                                                                                                                                  <  4:671064/151‑1 (MQ=255)
                   aaaaGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATaaaaa                                                                                                                                    >  4:583375/1‑148 (MQ=255)
                    aaaGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGt                                                                                                                                >  1:485271/1‑151 (MQ=255)
                      aGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTAc                                                                                                                              <  4:493811/151‑1 (MQ=255)
                      aGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAATa                                                                                                                               >  2:258676/1‑150 (MQ=255)
                          aTTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGc                                                                                                                           <  4:1069900/150‑1 (MQ=255)
                               gcgcTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGacg                                                                                                                     <  1:486189/151‑1 (MQ=255)
                               gcgcTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGa                                                                                                                       <  1:1271773/149‑1 (MQ=255)
                                cgcTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGacgg                                                                                                                    >  1:55513/1‑151 (MQ=255)
                                cgcTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCAACGCGCAGATTGTAACGACACAGGCCGTTGTAACCGGAGAGGATCGAAAATTAAAACAATACCGaccaatgg                                                                                                                    >  1:904233/1‑146 (MQ=255)
                                cgcTTTTTCAGCGAAGTGATACAACGCTGCTTCCAGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGacgg                                                                                                                    <  3:125808/151‑1 (MQ=255)
                                  cTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGacgga                                                                                                                   >  2:67316/1‑150 (MQ=255)
                                    ttttCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGacggaac                                                                                                                 <  1:866669/150‑1 (MQ=255)
                                     tttCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGacggaacg                                                                                                               <  1:241787/150‑1 (MQ=255)
                                        cAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATc                                                                                                           >  2:277743/1‑151 (MQ=255)
                                         aGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATcc                                                                                                          >  3:492876/1‑151 (MQ=255)
                                          gCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCt                                                                                                         <  1:1754615/151‑1 (MQ=255)
                                           cGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTg                                                                                                        <  1:904612/151‑1 (MQ=255)
                                            gAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTTCAATATAAAAACAGTACCGCCCGACGGGACGGATCCTgg                                                                                                       >  4:365509/1‑151 (MQ=255)
                                            gAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTgg                                                                                                       >  4:1589634/1‑151 (MQ=255)
                                                tGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCCGGCttt                                                                                                   >  2:95007/1‑151 (MQ=255)
                                                 gATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCtttt                                                                                                  >  2:49812/1‑151 (MQ=255)
                                                 gATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCtttt                                                                                                  <  4:1316442/151‑1 (MQ=255)
                                                    aCAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTttgtt                                                                                               <  1:771916/151‑1 (MQ=255)
                                                     cAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTttgttg                                                                                              <  3:352881/151‑1 (MQ=255)
                                                      aaCTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTttgtt                                                                                               >  4:1791045/1‑149 (MQ=255)
                                                       aCTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGaa                                                                                            <  1:1595432/151‑1 (MQ=255)
                                                       aCTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAAAAGTACCGCCCCACGGAAAGGATCCCGGCTTtttttg                                                                                              >  1:1666237/1‑149 (MQ=255)
                                                           tGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTTGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAActct                                                                                        >  1:447485/1‑151 (MQ=255)
                                                            gCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAAct                                                                                          <  4:1486990/148‑1 (MQ=255)
                                                             cTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCtt                                                                                       >  1:144063/1‑150 (MQ=255)
                                                                ccGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTc                                                                                      <  1:27780/148‑1 (MQ=255)
                                                                  gggAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGa                                                                                 <  1:198397/151‑1 (MQ=255)
                                                                    gAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGaa                                                                                >  4:1322247/1‑150 (MQ=255)
                                                                     aaaTAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAAca                                                                              <  1:1534393/151‑2 (MQ=255)
                                                                     aaaTAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACt                                                                              <  4:298899/151‑1 (MQ=255)
                                                                      aaTAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTg                                                                             >  2:150360/1‑151 (MQ=255)
                                                                      aaTAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTg                                                                             <  1:361382/151‑1 (MQ=255)
                                                                      aaTAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTg                                                                             <  3:329956/151‑1 (MQ=255)
                                                                       aTAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGaa                                                                                >  4:1277676/1‑147 (MQ=255)
                                                                       aTAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGaa                                                                                >  3:42481/1‑147 (MQ=255)
                                                                         aGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGAt                                                                           <  1:380404/150‑1 (MQ=255)
                                                                          gggTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGAt                                                                           >  2:317996/1‑149 (MQ=255)
                                                                          gggTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGAt                                                                           >  2:468862/1‑149 (MQ=255)
                                                                          gggTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCCACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGAt                                                                           >  1:94887/1‑149 (MQ=255)
                                                                                ttcgTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAAcc                                                                   >  4:100041/4‑151 (MQ=255)
                                                                                 ttGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAAcc                                                                   <  4:862096/150‑1 (MQ=255)
                                                                                 ttGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAAATCTTCGCGGAACTGATAGTAAAcc                                                                   >  4:1250358/1‑150 (MQ=255)
                                                                                  tGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCt                                                                  <  3:20724/150‑1 (MQ=255)
                                                                                     aGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGgc                                                              <  2:31079/151‑1 (MQ=255)
                                                                                      gCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCCTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGgcg                                                             >  4:1065184/1‑151 (MQ=255)
                                                                                         gCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCaa                                                          <  1:608140/151‑1 (MQ=255)
                                                                                          ccGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAAt                                                         <  1:596183/151‑1 (MQ=255)
                                                                                           cGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAAAAGTACCGCCCGACGGGACGGATCCTGGCCGTTGTTTAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAAtt                                                        >  4:1674023/1‑151 (MQ=255)
                                                                                            gAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAAttt                                                       <  4:277913/151‑1 (MQ=255)
                                                                                             aaCGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTg                                                      <  1:1895441/151‑1 (MQ=255)
                                                                                             aaCGGCACGTTAAAATCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTTAACTATTAAAACAAAACCACCCCACGGGACGGACTCTGGCCTTTTTTGAACTCCTCTCGGGACTGATAGTAAACCGTTGCGCCAttt                                                       >  3:404958/1‑150 (MQ=255)
                                                                                                ggCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGttt                                                   <  1:1484831/151‑1 (MQ=255)
                                                                                                 gCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCCAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGGACTGATAGTAAACCTCGGCGCAATTTGttt                                                   >  1:642164/1‑150 (MQ=255)
                                                                                                 gCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTc                                                  <  4:144109/151‑1 (MQ=255)
                                                                                                    cGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCACGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTct                                                 <  1:1586340/149‑2 (MQ=255)
                                                                                                      ttAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACCGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGa                                              <  4:1859296/150‑1 (MQ=255)
                                                                                                       tAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGCTCCAATATAAAAAAATTACCGCCCTAACGAACGGGATCTTGCCTTTTTTTGACTCTTCGCGGAACTGAGAGTAAAACCCGGCGCAATTTTTTTCGGGaa                                             >  4:407303/1‑150 (MQ=255)
                                                                                                        aaaCTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGaaca                                           <  4:1166897/151‑1 (MQ=255)
                                                                                                        aaaCTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGGACCGCCCGACGGGACGGATCCTTGGTTTTGTTGAACTCTTCGCGGCACTGATAGTAAACCTCGGCGCAATTTGTTTCgggga                                             >  4:1622152/1‑147 (MQ=255)
                                                                                                          aCTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGaacaa                                          >  1:1763878/1‑150 (MQ=255)
                                                                                                            tCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGaacaaa                                         >  3:495121/1‑149 (MQ=255)
                                                                                                                cgcgCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTc                                   <  4:1786425/151‑1 (MQ=255)
                                                                                                                   gcAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCCTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTcgc                                 >  4:1661516/1‑150 (MQ=255)
                                                                                                                      gATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGGTCCTGGCTTTTGTTTAACTCTTCGCGGAACTGATAGTAAAACTCGGCGCAATTTTTTTCGGGGACAAACAGCCcgcggg                              >  1:1029045/1‑148 (MQ=255)
                                                                                                                      gATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGGTAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTg                              >  2:429110/1‑150 (MQ=255)
                                                                                                                          gTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGcc                         <  4:1037867/151‑1 (MQ=255)
                                                                                                                          gTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTAACGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCCTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCgg                           >  4:75297/1‑149 (MQ=255)
                                                                                                                           tAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAATACCGCCCGACGGAACGGATCCCGGCTTTTGTTGAAACCTTCGCCGAACTGATAGTAAAACTCCGCGCAATTTGGTTCGGGGACAAACAGCTCGCGTGCGGcc                         >  1:1810832/1‑150 (MQ=255)
                                                                                                                            aaCGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCt                        <  4:653192/150‑1 (MQ=255)
                                                                                                                            aaCGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTg                       <  1:1548303/151‑1 (MQ=255)
                                                                                                                              cGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTgta                     <  4:1526864/151‑1 (MQ=255)
                                                                                                                                 cacaGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTgtacgt                  <  1:667820/151‑1 (MQ=255)
                                                                                                                                 cacaGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTgtacgt                  <  3:226271/151‑1 (MQ=255)
                                                                                                                                 cacaGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTgtacgt                  <  1:1347038/151‑1 (MQ=255)
                                                                                                                                  acaGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTgtacgta                 <  1:640162/151‑1 (MQ=255)
                                                                                                                                        cGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACtcatc           <  2:231790/151‑1 (MQ=255)
                                                                                                                                         gTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACtcatca          <  4:1770102/151‑1 (MQ=255)
                                                                                                                                          ttGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACtcatc           >  2:395053/1‑149 (MQ=255)
                                                                                                                                           tGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAg         >  4:624059/1‑150 (MQ=255)
                                                                                                                                           tGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAg         <  4:1100844/150‑1 (MQ=255)
                                                                                                                                              aaCCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTAc      <  1:1189665/150‑1 (MQ=255)
                                                                                                                                                ccGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGc    >  1:649575/1‑150 (MQ=255)
                                                                                                                                                 cGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCAt  <  1:1005974/151‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCAT  >  NC_000913/3515417‑3515712

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