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.

GCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAAC‑GGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCAT  >  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)
                                     tttCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGacggaac‑g                                                                                                               <  1:241787/150‑1 (MQ=255)
                                        cAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAAC‑GGATc                                                                                                           >  2:277743/1‑151 (MQ=255)
                                         aGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAAC‑GGATcc                                                                                                          >  3:492876/1‑151 (MQ=255)
                                          gCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAAC‑GGATCCt                                                                                                         <  1:1754615/151‑1 (MQ=255)
                                           cGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAAC‑GGATCCTg                                                                                                        <  1:904612/151‑1 (MQ=255)
                                            gAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTTCAATATAAAAACAGTACCGCCCGACGGGAC‑GGATCCTgg                                                                                                       >  4:365509/1‑151 (MQ=255)
                                            gAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAAC‑GGATCCTgg                                                                                                       >  4:1589634/1‑151 (MQ=255)
                                                tGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAAC‑GGATCCCGGCttt                                                                                                   >  2:95007/1‑151 (MQ=255)
                                                 gATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAAC‑GGATCCTGGCtttt                                                                                                  >  2:49812/1‑151 (MQ=255)
                                                 gATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAAC‑GGATCCTGGCtttt                                                                                                  <  4:1316442/151‑1 (MQ=255)
                                                    aCAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAAC‑GGATCCTGGCTTttgtt                                                                                               <  1:771916/151‑1 (MQ=255)
                                                     cAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAAC‑GGATCCTGGCTTttgttg                                                                                              <  3:352881/151‑1 (MQ=255)
                                                      aaCTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAAC‑GGATCCTGGCTTttgtt                                                                                               >  4:1791045/1‑149 (MQ=255)
                                                       aCTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAAC‑GGATCCTGGCTTTTGTTGaa                                                                                            <  1:1595432/151‑1 (MQ=255)
                                                       aCTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAAAAGTACCGCCCCACGGAAA‑GGATCCCGGCTTtttttg                                                                                              >  1:1666237/1‑149 (MQ=255)
                                                           tGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTTGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAAC‑GGATCCTGGCTTTTGTTGAActct                                                                                        >  1:447485/1‑151 (MQ=255)
                                                            gCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAAC‑GGATCCTGGCTTTTGTTGAAct                                                                                          <  4:1486990/148‑1 (MQ=255)
                                                             cTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAAC‑GGATCCTGGCTTTTGTTGAACTCtt                                                                                       >  1:144063/1‑150 (MQ=255)
                                                                ccGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAAC‑GGATCCTGGCTTTTGTTGAACTCTTc                                                                                      <  1:27780/148‑1 (MQ=255)
                                                                  gggAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAAC‑GGATCCTGGCTTTTGTTGAACTCTTCGCGGa                                                                                 <  1:198397/151‑1 (MQ=255)
                                                                    gAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAAC‑GGATCCTGGCTTTTGTTGAACTCTTCGCGGaa                                                                                >  4:1322247/1‑150 (MQ=255)
                                                                     aaaTAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAAC‑GGATCCTGGCTTTTGTTGAACTCTTCGCGGAAca                                                                              <  1:1534393/151‑2 (MQ=255)
                                                                     aaaTAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAAC‑GGATCCTGGCTTTTGTTGAACTCTTCGCGGAACt                                                                              <  4:298899/151‑1 (MQ=255)
                                                                      aaTAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAAC‑GGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTg                                                                             >  2:150360/1‑151 (MQ=255)
                                                                      aaTAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAAC‑GGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTg                                                                             <  1:361382/151‑1 (MQ=255)
                                                                      aaTAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAAC‑GGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTg                                                                             <  3:329956/151‑1 (MQ=255)
                                                                       aTAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAAC‑GGATCCTGGCTTTTGTTGAACTCTTCGCGGaa                                                                                >  4:1277676/1‑147 (MQ=255)
                                                                       aTAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAAC‑GGATCCTGGCTTTTGTTGAACTCTTCGCGGaa                                                                                >  3:42481/1‑147 (MQ=255)
                                                                         aGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAAC‑GGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGAt                                                                           <  1:380404/150‑1 (MQ=255)
                                                                          gggTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAAC‑GGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGAt                                                                           >  2:317996/1‑149 (MQ=255)
                                                                          gggTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAAC‑GGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGAt                                                                           >  2:468862/1‑149 (MQ=255)
                                                                          gggTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCCACGGAAC‑GGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGAt                                                                           >  1:94887/1‑149 (MQ=255)
                                                                                ttcgTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAAC‑GGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAAcc                                                                   >  4:100041/4‑151 (MQ=255)
                                                                                 ttGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAAC‑GGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAAcc                                                                   <  4:862096/150‑1 (MQ=255)
                                                                                 ttGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAAC‑GGATCCTGGCTTTTGTTGAAATCTTCGCGGAACTGATAGTAAAcc                                                                   >  4:1250358/1‑150 (MQ=255)
                                                                                  tGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAAC‑GGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCt                                                                  <  3:20724/150‑1 (MQ=255)
                                                                                     aGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAAC‑GGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGgc                                                              <  2:31079/151‑1 (MQ=255)
                                                                                      gCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAAC‑GGATCCTGGCCTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGgcg                                                             >  4:1065184/1‑151 (MQ=255)
                                                                                         gCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAAC‑GGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCaa                                                          <  1:608140/151‑1 (MQ=255)
                                                                                          ccGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAAC‑GGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAAt                                                         <  1:596183/151‑1 (MQ=255)
                                                                                           cGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAAAAGTACCGCCCGACGGGAC‑GGATCCTGGCCGTTGTTTAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAAtt                                                        >  4:1674023/1‑151 (MQ=255)
                                                                                            gAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAAC‑GGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAAttt                                                       <  4:277913/151‑1 (MQ=255)
                                                                                             aaCGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAAC‑GGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTg                                                      <  1:1895441/151‑1 (MQ=255)
                                                                                             aaCGGCACGTTAAAATCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTTAACTATTAAAACAAAACCACCCCACGGGAC‑GGACTCTGGCCTTTTTTGAACTCCTCTCGGGACTGATAGTAAACCGTTGCGCCAttt                                                       >  3:404958/1‑150 (MQ=255)
                                                                                                ggCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAAC‑GGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGttt                                                   <  1:1484831/151‑1 (MQ=255)
                                                                                                 gCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCCAATATAAAAACAGTACCGCCCGACGGAAC‑GGATCCTGGCTTTTGTTGAACTCTTCGCGGGACTGATAGTAAACCTCGGCGCAATTTGttt                                                   >  1:642164/1‑150 (MQ=255)
                                                                                                 gCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAAC‑GGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTc                                                  <  4:144109/151‑1 (MQ=255)
                                                                                                    cGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCACGACGGAAC‑GGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTct                                                 <  1:1586340/149‑2 (MQ=255)
                                                                                                      ttAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAAC‑CGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGa                                              <  4:1859296/150‑1 (MQ=255)
                                                                                                       tAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGCTCCAATATAAAAAAATTACCGCCCTAACGAACGGGATCTTGCCTTTTTTTG‑ACTCTTCGCGGAACTGAGAGTAAAACCCGGCGCAATTTTTTTCGGGaa                                             >  4:407303/1‑150 (MQ=255)
                                                                                                        aaaCTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAAC‑GGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGaaca                                           <  4:1166897/151‑1 (MQ=255)
                                                                                                        aaaCTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGGACCGCCCGACGGGAC‑GGATCCTTGGTTTTGTTGAACTCTTCGCGGCACTGATAGTAAACCTCGGCGCAATTTGTTTCgggga                                             >  4:1622152/1‑147 (MQ=255)
                                                                                                          aCTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAAC‑GGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGaacaa                                          >  1:1763878/1‑150 (MQ=255)
                                                                                                            tCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAAC‑GGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGaacaaa                                         >  3:495121/1‑149 (MQ=255)
                                                                                                                cgcgCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAAC‑GGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTc                                   <  4:1786425/151‑1 (MQ=255)
                                                                                                                   gcAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAAC‑GGATCCTGGCCTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTcgc                                 >  4:1661516/1‑150 (MQ=255)
                                                                                                                      gATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAAC‑GGGTCCTGGCTTTTGTTTAACTCTTCGCGGAACTGATAGTAAAACTCGGCGCAATTTTTTTCGGGGACAAACAGCCcgcggg                              >  1:1029045/1‑148 (MQ=255)
                                                                                                                      gATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAAC‑GGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGGTAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTg                              >  2:429110/1‑150 (MQ=255)
                                                                                                                          gTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAAC‑GGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGcc                         <  4:1037867/151‑1 (MQ=255)
                                                                                                                          gTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTAACGCCCGACGGAAC‑GGATCCTGGCTTTTGTTGAACTCCTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCgg                           >  4:75297/1‑149 (MQ=255)
                                                                                                                           tAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAATACCGCCCGACGGAAC‑GGATCCCGGCTTTTGTTGAAACCTTCGCCGAACTGATAGTAAAACTCCGCGCAATTTGGTTCGGGGACAAACAGCTCGCGTGCGGcc                         >  1:1810832/1‑150 (MQ=255)
                                                                                                                            aaCGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAAC‑GGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCt                        <  4:653192/150‑1 (MQ=255)
                                                                                                                            aaCGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAAC‑GGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTg                       <  1:1548303/151‑1 (MQ=255)
                                                                                                                              cGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAAC‑GGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTgta                     <  4:1526864/151‑1 (MQ=255)
                                                                                                                                 cacaGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAAC‑GGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTgtacgt                  <  1:667820/151‑1 (MQ=255)
                                                                                                                                 cacaGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAAC‑GGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTgtacgt                  <  3:226271/151‑1 (MQ=255)
                                                                                                                                 cacaGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAAC‑GGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTgtacgt                  <  1:1347038/151‑1 (MQ=255)
                                                                                                                                  acaGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAAC‑GGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTgtacgta                 <  1:640162/151‑1 (MQ=255)
                                                                                                                                        cGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAAC‑GGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACtcatc           <  2:231790/151‑1 (MQ=255)
                                                                                                                                         gTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAAC‑GGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACtcatca          <  4:1770102/151‑1 (MQ=255)
                                                                                                                                          ttGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAAC‑GGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACtcatc           >  2:395053/1‑149 (MQ=255)
                                                                                                                                           tGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAAC‑GGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAg         >  4:624059/1‑150 (MQ=255)
                                                                                                                                           tGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAAC‑GGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAg         <  4:1100844/150‑1 (MQ=255)
                                                                                                                                              aaCCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAAC‑GGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTAc      <  1:1189665/150‑1 (MQ=255)
                                                                                                                                                ccGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAAC‑GGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGc    >  1:649575/1‑150 (MQ=255)
                                                                                                                                                 cGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAAC‑GGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCAt  <  1:1005974/151‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAAC‑GGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTACGCAT  >  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: ‑