Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,342,193 A→G T95A (ACG→GCG)  melA → alpha‑galactosidase, NAD(P)‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,342,1930AG93.7% 264.8 / 1.2 95T95A (ACG→GCG) melAalpha‑galactosidase, NAD(P)‑binding
Reads supporting (aligned to +/- strand):  ref base A (4/2);  new base G (42/47);  total (46/49)
Fisher's exact test for biased strand distribution p-value = 4.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.87e-01

TGGAAGAGTCGCATATTGTGGTGCGTAAGCTGATGGATTCAGCAGGGGCCAGCGGCAAAATCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCACGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGTATTATGCGCGCGCTACGTACCATTCCGCATCTGTGGCAAATTTGCGAGGACATGACGGAAGTCTGCC  >  NC_000913/4342047‑4342334
                                                                                                                                                  |                                                                                                                                             
tGGAAGAGTCGCATATTGTGGTGCGTCAGCTGATGGATTCAGCAGGGGCCAGCGGCAAAATCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGATATTAACCTTGCAgggg                                                                                                                                           >  2:818279/1‑150 (MQ=255)
 ggAAGAGTCGCATATTGTGGTGCGTAAGCTGATGGATTCAGCAGGGGCCAGCGGCAAAATCCCCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCGCGGTg                                                                                                                                          <  1:937915/151‑1 (MQ=255)
 ggAAGAGTCGCATATTGTGGTGCGTAAGCTGATGGATTCAGCAGGGGCCAGCGGCAAAATCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCGCgg                                                                                                                                            >  2:1341021/1‑149 (MQ=255)
 ggAAGAGTCGCATATTGTGGTGCGTAAGCTGATGGATTCAGCAGGGGCCAGCGGCAAAATCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCGCGGt                                                                                                                                           <  1:682153/150‑1 (MQ=255)
 ggAAGAGTCGCATATTGTGGTGCGTAAGCTGATGGATTCAGCAGGGGCCAGCGGCAAAATCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCGCGGTg                                                                                                                                          <  2:1341404/151‑1 (MQ=255)
  gAAGAGTCGCATATTGTGGTGCGTAAGCTGATGGATTCAGCAGGGGCCAGCGGCAAAATCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCGCGGTGa                                                                                                                                         <  2:1455384/151‑1 (MQ=255)
   aaGAGTCGCATATTGTGGTGCGTAAGCTGATGGATTCAGCAGGGGCCAGCGGCAAAATCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACATTACGCgg                                                                                                                                            >  1:770612/1‑147 (MQ=255)
     gagTCGCATATTGTGGTGCGTAAGCTGATGGATTCAGCAGGGGCCAGCGGCAAAATCCCCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCGCGGTGACt                                                                                                                                       <  1:50915/150‑1 (MQ=255)
     gagTCGCATATTGTGGTGCGTAAGCTGATGGATTCAGCAGGGGCCAGCGGCAAAATCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCGCGGTGACTg                                                                                                                                      <  1:208328/151‑1 (MQ=255)
     gagTCGCATATTGTGGTGCGTAAGCTGATGGATTCAGCAGGGGCCAGCGGCAAAATCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCGCGGTGACTg                                                                                                                                      <  2:1136216/151‑1 (MQ=255)
     gagTCGCATATTGTGGTGCGTAAGCTGATGGATTCAGCAGGGGCCAGCGGCAAAATCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAAACTTTCGCGGGGACt                                                                                                                                       >  2:5670/1‑150 (MQ=255)
      agTCGCATATTGTGGTGCGTAAGCTGATGGATTCAGCAGGGGCCAGCGGCAAAATCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCACGGTGACTg                                                                                                                                      >  2:716282/1‑150 (MQ=255)
        tCGCATATTGTGGTGCGTAAGCTGATGGATTCAGCAGGGGCCAGCGGCAAAATCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATTAACCTTTCGCGGTggc                                                                                                                                        >  1:743059/1‑144 (MQ=255)
         cGCATATTGTGGTGCGTAAGCTGATGGATTCAGCAGGGGCCAGCGGCAAAATCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCGCGGTGACTGAttt                                                                                                                                  <  1:734077/151‑1 (MQ=255)
         cGCATATTGTGGTGCGTAAGCTGATGGATTCAGCAGGGGCCAGCGGCAAAATCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCGCGGTGACTGAtt                                                                                                                                   <  1:249715/150‑1 (MQ=255)
            atatTGTGGTGCGTAAGCTGATGGATTCAGCAGGGGCCAGCGGCAAAATCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTTCGCGGGGGCTGATTTCg                                                                                                                                >  2:1714262/1‑150 (MQ=255)
            atatTGTGGTGCGTAAGCTGATGGATTCAGCAGGGGCCAGCGGCAAAATCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCGCGGGGACTGAttt                                                                                                                                  >  2:329115/1‑148 (MQ=255)
                   ggtgCGTAAGCTGATGGATTCAGCAGGGGCCAGCGGCAAAATCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCGCGGTGACTGATTTCGAGGTCTg                                                                                                                         <  1:776017/150‑1 (MQ=255)
                   ggtgCGTAAGCTGATGGATTCAGCAGGGGCCAGCGGCAAAATCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCGCGGTGACTGATTTCGAGGTCTg                                                                                                                         <  2:56881/150‑1 (MQ=255)
                    gtgCGTAAGCTGATGGATTCAGCAGGGGCCAGCGGCAAAATCCCCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCGCGGTGACTGATTTCGAGGTCTGTa                                                                                                                       <  1:220088/151‑1 (MQ=255)
                     tgCGTAAGCTGATGGATTCAGCAGGGGCCAGCGGCAAAATCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCCTGCGCGGTGGCTGATTTCGAGGTCTGt                                                                                                                        >  2:1320708/1‑149 (MQ=255)
                      gCGTAAGCTGATGGATTCAGCGGGGGCCAGCGGCAAAATCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCGCGGTGACTGATTTCGAGGTCTGTAAg                                                                                                                     <  1:1209662/151‑1 (MQ=255)
                         tAAGCTGATGGATTCAGCAGGGGCCAGCGGCAAAATCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTTTGAAACTTTCGCGGTGACTGGTTTTGGGGTCTGTAAACgg                                                                                                                  >  2:1034579/1‑151 (MQ=255)
                          aaGCTGATGGATTCAGCAGGGGCCAGCGGCAAAATCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCGCGGTGACTGATTTCGAGGTCTGTAAGCgg                                                                                                                  <  1:539466/150‑1 (MQ=255)
                          aaGCTGATGGATTCAGCAGGGGCCAGCGGCAAAATCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCGCGGGGACTGATTTTGAGGGCTGTAAGCgg                                                                                                                  >  2:239408/1‑150 (MQ=255)
                            gCTGATGGATTCAGCAGGGGCCAGCGGCAAAATCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCGCGGTGACTGATTTCGAGGTCTGTAAGCGGCAt                                                                                                               <  1:377954/151‑1 (MQ=255)
                             cTGATGGATTCAGCAGGGGCCAGCGGCAAAATCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCGCGGTGACTGATTTCGAGGTCTGTAAGCGGCAt                                                                                                               >  1:1522282/1‑150 (MQ=255)
                              tGATGGATTCAGCAGGGGCCAGCGGCAAAATCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCGCGGTGACTGATTTCGAGGTCTGTAAGCGGCATgg                                                                                                             >  2:1223836/1‑151 (MQ=255)
                               gATGGATTCAGCAGGGGCCAGCGGCAAAATCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCGCGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGt                                                                                                            <  1:1658297/151‑1 (MQ=255)
                               gATGGATTCAGCAGGGGCCAGCGGCAAAATCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCGCGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGt                                                                                                            <  2:1671165/151‑1 (MQ=255)
                                aTGGATTCAGCAGGGGCCAGCGGCAAAATCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCGCGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGt                                                                                                            <  1:82698/150‑1 (MQ=255)
                                aTGGATTCAGCAGGGGCCAGCGGCAAAATCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCGCGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTc                                                                                                           <  1:1261977/151‑1 (MQ=255)
                                aTGGATTCAGCAGGGGCCAGCGGCAAAATCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCGCGGGGACTGATTTCGAGGTCTGTAAGCGGCACgg                                                                                                             >  2:1508953/1‑149 (MQ=255)
                                 tGGATTCAGCAGGGGCCAGCGGCAAAATCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTTTTAAACTTGCGCGGGGGCTGGTTTTCGGGGCTGTAAGCgg                                                                                                                  >  1:225910/1‑143 (MQ=255)
                                 tGGATTCAGCAGGGGCCAGCGGCAAAATCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCGCGGTGACTGATTTCGAGGTCTGTAAGCGGCATgg                                                                                                             <  2:425975/148‑1 (MQ=255)
                                  ggATTCAGCAGGGGCCAGCGGCAAAATCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCGCGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTccg                                                                                                         <  1:368660/151‑3 (MQ=255)
                                    aTTCAGCAGGGGCCAGCGGCAAAATCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCGCGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGa                                                                                                       <  1:417588/151‑1 (MQ=255)
                                       cagaagTGGCGAGCGGCAAAACCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCTGATTTTGTCGTTGTGGCATTTCAGATTGGCGGTTATGAGCCTTGCGCGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGaa                                                                                                      <  2:1217068/149‑1 (MQ=255)
                                        agcagGGGCCAGCGGCAAAATCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCGCGGTGACTGATTTCGTGGTCTGTAAGCGGCATGGGCTGGAACa                                                                                                    >  1:465613/1‑150 (MQ=255)
                                        agcagGGGCCAGCGGCAAAACCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCGCGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACaa                                                                                                   <  2:1059599/151‑1 (MQ=255)
                                            ggggCCAGCGGCAAAATCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCGCGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCa                                                                                               <  1:860924/151‑1 (MQ=255)
                                              ggCCAGCGGCAAAATCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTTCGCGGTGACTGATTTCGAGGTCTGTAAACGGCATGGGTTGGAACAAACCAtt                                                                                             >  2:1707546/1‑151 (MQ=255)
                                                  aGCGGCAAAATCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCGCGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGcc                                                                                          <  2:1076672/150‑1 (MQ=255)
                                                  aGCGGCAAAATCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCGCGGGGACTGATTTTGAGGGCTTTAAGCGGCATTGTCTGGGACAAAACCATGCCg                                                                                         >  1:1401455/1‑151 (MQ=255)
                                                  aGCGGCAAAATCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGCTATGAAACCTGCACGGTGGCTGGTTTCGGGGTCTGTAAGCGGCATGGGCTGGAAAAAAACATTTcc                                                                                          >  1:229007/1‑150 (MQ=255)
                                                   gCGGCAAAATCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCGCGGTGACTGATTTTGAGGGCTGTAAGCGGCATGGGCTGGGACAAAACAATGCCg                                                                                         >  1:156566/1‑150 (MQ=255)
                                                   gCGGCAAAATCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCGCGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCg                                                                                         <  1:994985/150‑1 (MQ=255)
                                                    cGGCAAAATCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATTAACCCTGCACCGGGACTGATTTTGAGGGCTGTAAGCGGCATGGGCTGGAACAAACCATTTcc                                                                                          >  1:1479862/1‑148 (MQ=255)
                                                     ggCAAAATCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGGAACATGCGCGGGGTCTGGTTTCGAGGTCTGTAAGCGGCATGGTGTGGGACAAAACATTGcc                                                                                          >  1:598432/1‑147 (MQ=255)
                                                        aaaaTCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCGCGGTGACTGATTTTGAGGTCTGTTAGCGGCATGGTCTGGGACAAACCATTGCCCGTACg                                                                                    >  2:1415406/1‑150 (MQ=255)
                                                        aaaaTCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAAACCTGCGCGGTGACTGATTTTGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTTCCCAt                                                                                       >  2:812333/1‑147 (MQ=255)
                                                          aaTCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAAACTTGCGCGGTGGCTGATTTTGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCAATGCCGATACggt                                                                                  >  2:15727/1‑148 (MQ=255)
                                                                  gCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCGCGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCgg                                                                          <  1:1496668/150‑1 (MQ=255)
                                                                    cacaCCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTTTGAAACTTGCGCGGGGGCTGATTTCGAGGTCTGTAAGAGGCATGGTCTTGAACAAAACCTTTCCGAGACGTTGGGGCCGGGCg                                                                       >  2:1232958/1‑151 (MQ=255)
                                                                      cacCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCGCGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGt                                                                     <  1:1488833/151‑1 (MQ=255)
                                                                           aaCAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAAACCTGCGCGGGGGCTGGTTTTGGGGGCTGTAAACGGCATGGTCTGGAACAAAACATTTCCGATAAGTTGGGGGCGGGCGGgatt                                                                  >  1:1363021/1‑149 (MQ=255)
                                                                            aCAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCGCGGGGACTGATTTTGAGGTCTGTAAGCGGCCTGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGgattat                                                                >  2:1368462/1‑150 (MQ=255)
                                                                            aCAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCGCGGGGACTGATTTTGAGGTCTGTAAGCGGCATGGGCTGGGAAAAACCATTGCCGATACGTTTGGGCCGGGCGGgatt                                                                  >  2:867234/1‑148 (MQ=255)
                                                                                aaagaaGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATTAAACTTGCGCGGGGACTGATTTCGAGGGCTGTAAGCGGCATGGGCTGGGACAAAACATTGCCGATACGTTGGGGGCGGGGGGGATTATgcgc                                                            >  2:1462217/1‑150 (MQ=255)
                                                                                   gaagCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCGCGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGTATTATgcgcgcgc                                                        <  2:379129/151‑1 (MQ=255)
                                                                                       ccTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCGCGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGTATTATGCGCGCGCTa                                                      <  2:598486/149‑1 (MQ=255)
                                                                                       ccTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCGCGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGTATTATGCGCGCGCTACg                                                    <  1:5667/151‑1 (MQ=255)
                                                                                       ccTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCGCGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGTATTATGCGCGCGCTACg                                                    <  2:397551/151‑1 (MQ=255)
                                                                                         ttAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGTGCGGGGAATGATTTCGGGGGCTGTAAGCGGGATTGGCTGGGACAAACCATTGCCGGTACGTTTGGGCCGGGCGGGATTATTGGCGCGCTAc                                                     >  2:1568690/1‑148 (MQ=255)
                                                                                         ttAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCGCGGTGGCTGATTTTGAGGGCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGGCGGGCGGGATTATGCGCGCGCTACGt                                                   >  2:201459/1‑150 (MQ=255)
                                                                                          tAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCGCGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGTATTATGCGCGCGCTACGt                                                   <  1:627932/149‑1 (MQ=255)
                                                                                          tAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCGCGGGGGCTGGTTTTTAGGGCTGTAAGCGGCATGGTCTGGGAAAAACCATTTCCGGTACGTTTGGGGCGGGCGGGATTTTTCGCGCCCTACGTa                                                  >  1:1667735/1‑150 (MQ=255)
                                                                                           agagGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCGCGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGTATTATGCGCGCGCTACGTAc                                                 <  1:794565/150‑1 (MQ=255)
                                                                                            gagGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAAACTTGCGCGGGGGCTGATTTCGGGGGCTGTAAGCGGCATGGTCTGGAACAAAACATTGCCGATACGTTGGGGCCGGGCGGGATTATGCGCGCGCCACGTAcc                                                >  1:1059151/1‑150 (MQ=255)
                                                                                              ggATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATTAACCTTTCGCCGTGGCTGGTTTCGTGGTGTGTAAGGGGAATGGTCTGGAACAAAACATTGCCGACACGTTGGGGCCCGGGGGGATTTTTCGGGCGCTACGCACCaat                                             >  1:1357342/1‑149 (MQ=255)
                                                                                                aTGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCCTGCGCGGTGGCTGATTTCGAGGTCTGTAAGCGGCATGGGCTGGGACAAACCATTGCCGGTACGTTGGGGCCGGGCGGGATTATGCGCGCGCTACGTACCCTTcc                                           >  1:115712/1‑151 (MQ=255)
                                                                                                   ccGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCGCGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGTATTATGCGCGCGCTACGTACCATTCCGCa                                        <  2:233222/151‑1 (MQ=255)
                                                                                                    cGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTTTTAACCCCCCGCGGGGACCTGTTTTTAGGGCTGTAAGCGGCAAGGGCTGGAACAAAACATTGCCGCAACCCTGGGGGCGGGCGGGATTATgcgcgcgcgcgctacca                                               >  1:1593924/1‑134 (MQ=255)
                                                                                                     gATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTTCGCGGGGGCTGATTTTGAGGGCTGTTAGCGGCATGGGCTGGGACAAAACATTGCCGATACGGTTGGGGCGGGGGGGATTATTCGGGCGGCACGTACCATTCCGcct                                       >  1:1549611/1‑148 (MQ=255)
                                                                                                     gATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCGCGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGTATTATGCGCGCGCTACGTACCATTCCGCATc                                      <  2:1417787/151‑1 (MQ=255)
                                                                                                     gATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCCTTCGCGGGGACTTATTTCGAGGTCTGTAAGCGGCATGGGCTGGAACAAACCAATGCCGGTAAGTTTGGGCCGGGGGGGATTATGCGCGCGCTACGTACCATTCCGc                                         >  2:369952/1‑148 (MQ=255)
                                                                                                             cGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCGCGGGGGCTGATTTCGAGGTCTGTAAACGGCATGGTCTGGAACAAACCATTTCCGATACGTTGGGGGCGGGCGGGATTATGCGCGCGCTACGTACCAATTCGCATCt                                     >  2:1076997/1‑144 (MQ=255)
                                                                                                              gtggtgGCATTTCAGATTGGCGGTTATGAACCTTGCGCGGGGGCTGATTTCGTGGTCTGTAAGCGGCATGGTCTGGGAAAAAACATTGCCGGTACGTTGGGGCCGGGCGGGATTATGCGCGCGCTACCTACCATTCCGCATCTGTgg                                 >  1:1055263/1‑147 (MQ=255)
                                                                                                              gtggtgGCATTTCAGATTGGCGGTTATGAACCTTGCACGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGTATTATGCGCGCGCTACGTACCATTCCGCATCTGTGGCaaa                             <  2:58709/151‑1 (MQ=255)
                                                                                                                ggtggCATTTCAGATTGGCGGTTATGAACCTTGCGCGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGTATTATGCGCGCGCTACGTACCATTCCGCATCTGTGGCAAAtt                           <  1:239386/151‑1 (MQ=255)
                                                                                                                   ggCATTTCAGATTGGCGGTTATGAACCTTGCGCGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGTATTATGCGCGCGCTACGTACCATTCCGCATCTGTGGCAAATTTg                         <  2:156266/150‑1 (MQ=255)
                                                                                                                   ggCATTTCAGATTGGCGGTTATGAACCTTGCACGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGTATTATGCGCGCGCTACGTACCATTCCGCATCTGTGGCAAATTTGc                        <  2:1201859/151‑1 (MQ=255)
                                                                                                                          cAGATTGGCGGTTATGAACCTTGCGCGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGTATTATGCGCGCGCTACGTACCATTCCGCATCTGTGGCAAATTTGCGAGGAc                  <  2:983195/150‑1 (MQ=255)
                                                                                                                           aGATTGGCGGTTATGAACCTTGCGCGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGTATTATGCGCGCGCTACGTACCATTCCGCATCTGTGGCAAATTTGCGAGGACAt                <  1:369924/151‑1 (MQ=255)
                                                                                                                            gATTGGCGGTTATGAACCTTGCGCGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGTATTATGCGCGCGCTACGTACCATTCCGCATCTGTGGCAAATTTGCGAGGACATg               <  1:1265037/151‑1 (MQ=255)
                                                                                                                            gATTGGCGGTTATGAACCTTGCGCGGGGACTGATTTCGAGGTCTGTAAACGGCCTGGTCTGGGACAAACCATTGCCGATACGTTGGGGCCGGGCGGGATTATGCGCGCGCTACGTACCATTCCTCCTCTGTGGCAAATTTTTGAgggcctg               >  1:1340625/1‑146 (MQ=255)
                                                                                                                             aTTGGCGGTTATGAACCTTGCGCGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGTATTATGCGCGCGCTACGTACCATTCCGCATCTGTGGCAAATTTGCGAGGACATGa              <  2:734047/151‑1 (MQ=255)
                                                                                                                                ggCGGTTATGAACCTTGCGCGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGTATTATGCGCGCGCTACGTACCATTCCGCATCTGTGGCAAATTTGCGAGGACATGACg            <  1:773309/150‑1 (MQ=255)
                                                                                                                                 gCGGATATGAACCTTGCGCGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGTATTATGCGCGCGCTACGTACCATTCCGCATCTGTGGCAAATTTGCGAGGACATGACgg           >  1:461800/1‑150 (MQ=255)
                                                                                                                                   ggTTATGAACCTTGCGCGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGTATTATGCGCGCGCTACGTACCATTCCGCATCTGTGGCAAATTTGCGAGGACATGACGGAAg        <  1:1462118/151‑1 (MQ=255)
                                                                                                                                   ggTTATGAACCTTGCGCGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGTATTATGCGCGCGCTACGTACCATTCCGCATCTGTGGCAAATTTGCGAGGACATGACGGAAg        <  2:743118/151‑1 (MQ=255)
                                                                                                                                      tATGAACCTTGCGCGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGTATTATGCGCGCGCTACGTACCATTCCGCATCTGTGGCAAATTTGCGAGGACATGACGGAAGTCt     >  1:197964/1‑151 (MQ=255)
                                                                                                                                         gAACCTTGCGCGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGTATTATGCGCGCGCTACGTACCATTCCGCATCTGTGGCAAATTTGCGAGGACATGACGGAAGTCTGcc  >  1:1309993/1‑151 (MQ=255)
                                                                                                                                         gAACCTTGCGCGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGTATTATGCGCGCGCTACGTACCATTCCGCATCTGTGGCAAATTTGCGAGGACATGACGGAAGTCTGcc  >  1:46698/1‑151 (MQ=255)
                                                                                                                                         gAACCTTGCGCGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGTATTATGCGCGCGCTACGTACCATTCCGCATCTGTGGCAAATTTGCGAGGACATGACGGAAGTCTGc   <  2:1619298/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                             
TGGAAGAGTCGCATATTGTGGTGCGTAAGCTGATGGATTCAGCAGGGGCCAGCGGCAAAATCACCTGCCACACCCAACAGAAAGAAGCCTTAGAGGATGCCGATTTTGTCGTGGTGGCATTTCAGATTGGCGGTTATGAACCTTGCACGGTGACTGATTTCGAGGTCTGTAAGCGGCATGGTCTGGAACAAACCATTGCCGATACGTTGGGGCCGGGCGGTATTATGCGCGCGCTACGTACCATTCCGCATCTGTGGCAAATTTGCGAGGACATGACGGAAGTCTGCC  >  NC_000913/4342047‑4342334

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