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,5650AG98.9% 338.7 / ‑3.9 93Y117H (TAC→CAC) damDNA adenine methyltransferase
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (41/51);  total (41/52)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.99e-01

GCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTA  >  NC_000913/3515417‑3515707
                                                                                                                                                    |                                                                                                                                              
gCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGc                                                                                                                                              >  2:874975/1‑151 (MQ=255)
  gTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCgg                                                                                                                                            <  1:938816/151‑1 (MQ=255)
   tAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGt                                                                                                                                           >  2:906463/1‑151 (MQ=255)
     aGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGtt                                                                                                                                          >  1:481536/1‑150 (MQ=255)
       aCTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTc                                                                                                                                         >  2:771115/1‑149 (MQ=255)
        cTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCaa                                                                                                                                       >  2:1077311/1‑150 (MQ=255)
         tCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTc                                                                                                                                         >  2:213671/1‑147 (MQ=255)
         tCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCaa                                                                                                                                       >  1:133393/1‑149 (MQ=255)
         tCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTAAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCaaa                                                                                                                                      >  2:224616/1‑150 (MQ=255)
          cacaATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCaaa                                                                                                                                      >  2:413351/1‑149 (MQ=255)
          cacaATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCaaa                                                                                                                                      <  2:1295353/149‑1 (MQ=255)
          cacaATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAat                                                                                                                                     >  2:189865/1‑150 (MQ=255)
           acaATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAatat                                                                                                                                   <  1:1074982/151‑1 (MQ=255)
             aaTAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAAtata                                                                                                                                  >  2:62713/1‑150 (MQ=255)
                aGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATaaaaa                                                                                                                              <  1:1090856/151‑1 (MQ=255)
                aGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATaaaa                                                                                                                               <  1:1120823/150‑1 (MQ=255)
                 gAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAAc                                                                                                                             <  1:305631/151‑1 (MQ=255)
                  aaaaaGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATaaaaa                                                                                                                              >  2:870526/1‑149 (MQ=255)
                  aaaaaGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAAc                                                                                                                             <  2:902423/150‑1 (MQ=255)
                  aaaaaGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACa                                                                                                                            <  2:1160584/151‑1 (MQ=255)
                  aaaaaGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACa                                                                                                                            <  1:983399/151‑1 (MQ=255)
                   aaaaGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACa                                                                                                                            >  1:156330/1‑150 (MQ=255)
                   aaaaGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAg                                                                                                                           <  2:781486/151‑1 (MQ=255)
                     aaGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTa                                                                                                                         <  1:521741/151‑1 (MQ=255)
                     aaGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTa                                                                                                                         <  1:101102/151‑1 (MQ=255)
                      aGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTa                                                                                                                         >  2:825884/1‑150 (MQ=255)
                      aGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTa                                                                                                                         >  2:517609/1‑150 (MQ=255)
                        gCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCg                                                                                                                      <  2:540226/151‑1 (MQ=255)
                        gCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCg                                                                                                                      <  2:876505/151‑1 (MQ=255)
                         cATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTTAAATATAAAAAAAGTaaccc                                                                                                                     >  1:1243616/1‑147 (MQ=255)
                         cATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGc                                                                                                                     <  1:70127/151‑1 (MQ=255)
                          aTTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGcc                                                                                                                    <  1:701286/151‑1 (MQ=255)
                             cTGCGCTTTTTCAGCGAAGTGATAAAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGCTACAATATAAAAAAAGTAcccccc                                                                                                                   >  1:701475/1‑149 (MQ=255)
                               gcgcTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGac                                                                                                                <  1:93534/150‑1 (MQ=255)
                               gcgcTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGac                                                                                                                <  1:1262315/150‑1 (MQ=255)
                               gcgcTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGa                                                                                                                 <  2:749446/149‑1 (MQ=255)
                               gcgcTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGa                                                                                                                 <  1:1251650/149‑1 (MQ=255)
                                    ttttCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGacggaac                                                                                                           <  1:832883/150‑1 (MQ=255)
                                       tCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGacggaacgga                                                                                                        <  2:1077494/150‑1 (MQ=255)
                                        cAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGAt                                                                                                       <  2:137765/150‑1 (MQ=255)
                                        cAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAAAAGTACCGCCCGACGGGACGGGTc                                                                                                      >  1:390622/1‑151 (MQ=255)
                                           cGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCCAATATAAAAAAAGGACCGCCCCACGGAACGGAtgc                                                                                                     >  2:914571/1‑147 (MQ=255)
                                            gAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGGTCCAATATAAAAACAGTAACCCCCCacgggacgggcccccg                                                                                                  >  1:692827/1‑144 (MQ=255)
                                               gTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGAGCCCGGCtt                                                                                               >  1:748940/1‑151 (MQ=255)
                                                 gATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCtttt                                                                                             >  2:543400/1‑151 (MQ=255)
                                                    aCAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTttgtt                                                                                          <  2:1085344/151‑1 (MQ=255)
                                                         tctGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTTTTGAAc                                                                                      >  2:119932/1‑150 (MQ=255)
                                                          ctGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCCAATATAAAAACAGTACCCCCCGACGGAAAGGGTCCTGGCTTtttttgga                                                                                       >  1:937330/1‑146 (MQ=255)
                                                            gCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCtt                                                                                  <  1:1078977/151‑1 (MQ=255)
                                                            gCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCtt                                                                                  <  2:776583/151‑1 (MQ=255)
                                                               tCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTcg                                                                                <  2:788328/150‑1 (MQ=255)
                                                                  gggAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGa                                                                            <  1:330988/151‑1 (MQ=255)
                                                                   ggAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTTAAATATAAAAACAGTACCCCCCGCCGGAACGGATCCTCGCTTTTGTTGGACTCTTCGCgg                                                                             >  2:86543/1‑149 (MQ=255)
                                                                   ggAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGaa                                                                           <  1:859339/151‑1 (MQ=255)
                                                                    gAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAAAATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCCCGGaa                                                                           >  1:126472/1‑150 (MQ=255)
                                                                    gAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGCGGCGGTTCCAAAATAAAAAACATACCGCCCGACGGAACGGACACCGGCTTATGGTGAAAACTTTCCGGaa                                                                           >  2:892236/1‑150 (MQ=255)
                                                                      aaTAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTg                                                                        <  2:66374/151‑1 (MQ=255)
                                                                        tAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCCAATATAAAAACAATACCGCCCCACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGGACTGa                                                                       >  1:216609/1‑150 (MQ=255)
                                                                         aGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGAt                                                                      <  2:1236196/150‑1 (MQ=255)
                                                                          gggTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAg                                                                    <  2:49899/151‑1 (MQ=255)
                                                                            gTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTACGCGGAACTGATAGt                                                                   <  1:311675/150‑1 (MQ=255)
                                                                              tttttGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTaa                                                                 <  2:138130/150‑1 (MQ=255)
                                                                               ttttGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTaaa                                                                >  2:398082/1‑150 (MQ=255)
                                                                                  tGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAAcc                                                              >  2:795363/1‑149 (MQ=255)
                                                                                   gTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTc                                                            <  2:1135272/150‑1 (MQ=255)
                                                                                     aGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCgg                                                          <  2:363955/150‑1 (MQ=255)
                                                                                     aGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCgg                                                          >  2:1291811/1‑150 (MQ=255)
                                                                                           cGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAAtt                                                   <  1:1252658/151‑1 (MQ=255)
                                                                                            gagcggctggtgaaaCTCACCTCGCAGATTGTAACGACAAAGGACGTTGTAACGGTGGCGGTTCAAATATAAAAACAGAACCCCCCGACGGAACGGATGCTGGCTTTGGTTGATCTCTTCGCGGAACTGATAGTAAACCTGGGCGCAAttt                                                  <  1:396247/139‑1 (MQ=255)
                                                                                             aaCGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTg                                                 <  2:944788/151‑1 (MQ=255)
                                                                                              aCGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCCAATATAAAAACAGTACCGCCCGACGGAACGGATCCCGGGTTTTTTTGAAATCTTCGCGGCACTGCTAGTAAACCTCGGCGCAACTTGt                                                >  1:362637/1‑151 (MQ=255)
                                                                                                   aCGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCGTTTGTTGAAATCTTCGCGGGACTGATAGTAAACCTCGGCGCAATTTGTTTCgg                                           >  2:605331/1‑151 (MQ=255)
                                                                                                     gTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGa                                         <  2:770262/151‑1 (MQ=255)
                                                                                                     gTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCATCGCGGAAATGATAGTAAACCTCGGCGCAATTTGTTTCGGGa                                         <  1:927761/151‑1 (MQ=255)
                                                                                                      ttAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGGACGGAACCTGGCTTTTGTTGAAATCTTCGCGGAACTGATAGTAAAACTCGGCGCAATTTGTTTTggg                                          >  1:556690/1‑149 (MQ=255)
                                                                                                        aaaCTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGaaca                                      <  2:1125633/151‑1 (MQ=255)
                                                                                                          aCTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGaacaaa                                    <  1:275739/151‑1 (MQ=255)
                                                                                                          aCTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGaacaaa                                    <  1:1223691/151‑1 (MQ=255)
                                                                                                           cTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGAACCTGGGTTTTGTTGAAATCTTCGCGGAACTGATAGTAAACCTCGGCGCCAATTGTTTTGGGaacaaa                                    >  1:875766/1‑150 (MQ=255)
                                                                                                             cACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAg                                 <  1:62714/151‑1 (MQ=255)
                                                                                                                cgcgCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTc                              <  2:701509/151‑1 (MQ=255)
                                                                                                                cgcgCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTc                              <  2:50214/151‑1 (MQ=255)
                                                                                                                cgcgCAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCCACGGAAAGGATCCTGGCCCTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAAAAGCTc                              >  2:1159372/1‑151 (MQ=255)
                                                                                                                  cgcAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTcgc                            <  2:886866/151‑1 (MQ=255)
                                                                                                                    cAGATTGTAACGACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGt                          <  2:512080/151‑1 (MQ=255)
                                                                                                                            aaCGACACAGGCCGTTGTAACCGTGGCGGGTCAAATATAAAAAAAGTACCGCCCCACGGAAAGGGACCTGGCTTTTGTTGGACTCTTCGCGGAACTGATAGTAAAACTCCGCGCAATTTGTTTCGGGAACCAACAGCTCGCGTGCGGcc                    >  1:1161280/1‑149 (MQ=255)
                                                                                                                               gACACAGGCCGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTgtac               <  1:221389/151‑1 (MQ=255)
                                                                                                                                       ccGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACtcat       <  1:363308/151‑1 (MQ=255)
                                                                                                                                        cGTTGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACtcat       <  1:914530/150‑1 (MQ=255)
                                                                                                                                        cGTTGTAACCGTGGCGGTTCAAATATAAAAAAAGTACCGCCCGACGGAAAGGATCCTGGCTTTTGTTTAAATCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACtca        >  1:654833/1‑149 (MQ=255)
                                                                                                                                          ttGTAACCGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCCTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACtcatc      >  2:1074040/1‑149 (MQ=255)
                                                                                                                                           tGTAACCGTGGCGGTTCAAATATAAAAAAAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAAATCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCCCGCGTGCGGCCTGTACGTACtcat       >  1:1033384/1‑147 (MQ=255)
                                                                                                                                                ccGTGGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTTGGGAACAAACAGCTCGCGTGCGGCCTGTACGTActcctccgga  >  1:609381/1‑143 (MQ=255)
                                                                                                                                                    |                                                                                                                                              
GCGTAAGACTCACAATAGAAAAAGGCATTCTGCGCTTTTTCAGCGAAGTGATACAACTCTGCTTCCGGGAAATAGGGTTTTTTGTAGCGGCCGAACGGCACGTTAAACTCACCGCGCAGATTGTAACGACACAGGCCGTTGTAACCGTAGCGGTTCAAATATAAAAACAGTACCGCCCGACGGAACGGATCCTGGCTTTTGTTGAACTCTTCGCGGAACTGATAGTAAACCTCGGCGCAATTTGTTTCGGGAACAAACAGCTCGCGTGCGGCCTGTACGTACTCATCAGTA  >  NC_000913/3515417‑3515707

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