Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,435,822 A→T 6.7% L50Q (CTG→CAG)  PP_1256 ← alpha‑ketoglutarate semialdehyde dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,435,8220AT 6.7% 417.6 / 15.8 133L50Q (CTG→CAG) PP_1256alpha‑ketoglutarate semialdehyde dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (62/62);  new base T (4/5);  total (66/67)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 136 total aligned reads displayed.

GAATGCGCGCCTTGGGCAGGCCGGTTTCGATCACGGCGCGGTCGATCAGGGCATCGCCGAGGGCCTCGATCTGTGTGGCGATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACGGATGGCGTCGCGGCTGCCGGTTACCGGCCTCTGGCCGATCAGCAGGTTGCCTGTGAGGGGCATG  >  NC_002947/1435673‑1435971
                                                                                                                                                     |                                                                                                                                                     
gAATGCGCGCCTTGGGCAGGCCGGTTTCGATCACGGCGCGGTCGATCAGGGCATCGCCGAGGGCCTCGATCTGTGTGGCGATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCAAGGCCa                                                                                                                                                       >  2:413268/1‑150 (MQ=255)
   tGCGCGCCTTGGGCAGGCCGGTTTCGATCACGGCGCGGTCGATCAGGGCATCGCCGAGGGCCTCGATCTGTGTGGCGATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAgcg                                                                                                                                                    >  2:182860/1‑150 (MQ=255)
      gcgcCTTGGGCAGGCCGGTTTCGATCACGGCGCGGTCGATCAGGGCATCGCCGAGGGCCTCGATCTGTGTGGCGATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAgcgcgc                                                                                                                                                 >  1:246117/1‑150 (MQ=255)
       cgcCTTGGGCAGGCCGGTTTCGATCACGGCGCGGTCGATCAGGGCATCGCCGAGGGCCTCGATCTGTGTGGCGATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCa                                                                                                                                                <  2:2944489/150‑1 (MQ=255)
              ggCAGGCCGGTTTCGATCACGGCGCGGTCGATCAGGGCATCGCCGAGGGCCTCGATCTGTGTGGCGATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTgg                                                                                                                                         <  2:1936264/150‑1 (MQ=255)
               gCAGGCCGGTTTCGATCACGGCGCGGTCGATCAGGGCATCGCCGAGGGCCTCGATCTGTGTGGCGATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCTGCGCGCTGGCCTggg                                                                                                                                        <  2:1691771/150‑1 (MQ=255)
                   gctggTTTCGATCACGGCGCGGTCGATCGGGGCATCGCCGAGGGCCGCGATCTGTGTGGCGATGGCTTCAAGGAACTCTGCTCGTTGGTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCAc                                                                                                                                    <  1:453370/147‑1 (MQ=255)
                     cGGTTTCGATCACGGCGCGGTCGATCAGGGCATCGCCGAGGGCCTCGATCTGTGTGGCGATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACgt                                                                                                                                  <  1:1929419/150‑1 (MQ=255)
                     cGGTTTCGATCACGGCGCGGTCGATCAGGGCATCGCCGAGGGCCTCGATCTGTGTGGCGATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCCAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACgt                                                                                                                                  <  1:2485399/150‑1 (MQ=255)
                        tttCGATCACGGCGCGGTCGATCAGGGCATCGCCGAGGGCCTCGATCTGTGTGGCGATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTg                                                                                                                                          >  1:944816/1‑139 (MQ=255)
                        tttCGATCACGGCGCGGTCGATCAGGGCATCGCCGAGGGCCTCGATCTGTGTGGCGATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGtt                                                                                                                               <  1:182860/150‑1 (MQ=255)
                           cGATCACGGCGCGGTCGATCAGGGCATCGCCGAGGGCCTCGATCTGTGTGGCGATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGc                                                                                                                            >  2:95/1‑150 (MQ=255)
                           cGATCACGGCGCGGTCGATCAGGGCATCGCCGAGGGCCTCGATCTGTGTGGCGATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGc                                                                                                                            <  2:2429143/150‑1 (MQ=255)
                            gATCACGGCGCGGTCGATCAGGGCATCGCCGAGGGCCTCGATCTGTGTGGCGATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGcc                                                                                                                           >  2:1095948/1‑150 (MQ=255)
                                 cGGCGCGGTCGATCAGGGCATCGCCGAGGGCCTCGATCTGTGTGGCGATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTgc                                                                                                                      <  1:1861178/150‑1 (MQ=255)
                                       ggTCGATCAGGGCATCGCCGAGGGCCTCGATCTGTGTGGCGATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGa                                                                                                                <  2:2602550/150‑1 (MQ=255)
                                       ggTCGATCAGGGCATCGCCGAGGGCCTCGATCTGTGTGGCGATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGa                                                                                                                <  2:2375175/150‑1 (MQ=255)
                                       ggTCGATCAGGGCATCGCCGACGGCCTCGATCTGTGTGGCGATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGa                                                                                                                >  1:91024/1‑150 (MQ=255)
                                           gATCAGGGCATCGCCGAGGGCCTCGATCTGTGTGGCGATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTa                                                                                                            >  1:207458/1‑150 (MQ=255)
                                            aTCAGGGCATCGCCGAGGGCCTCGATCTGTGTGGCGATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCCCGTGTTCGCCTGTGCCGCCGAGGTAg                                                                                                           <  1:1187602/150‑1 (MQ=255)
                                                gggCATCGCCGAGGGCCTCGATCTGTGTGGCGATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCtg                                                                                                                         <  1:2696553/132‑1 (MQ=255)
                                                gggCATCGCCGAGGGCCTCGATCTGTGTGGCGATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCtg                                                                                                                         >  2:2696553/1‑132 (MQ=255)
                                                gggCATCGCCGAGGGCCTCGATCTGTGTGGCGATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCg                                                                                                       >  2:2984268/1‑150 (MQ=255)
                                                 ggCATCGCCGAGGGCCTCGATCTGTGTGGCGATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCgg                                                                                                      >  2:1451434/1‑150 (MQ=255)
                                                  gCATCGCCGAGGGCCTCGATCTGTGTGGCGATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGc                                                                                                     >  1:208190/1‑150 (MQ=255)
                                                   cATCGCCGAGGGCCTCGATCTGTGTGGCGATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCt                                                                                                    <  2:448577/150‑1 (MQ=255)
                                                        ccGAGGGCCTCGATCTGTGTGGCGATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAgg                                                                                                          <  2:579617/139‑1 (MQ=255)
                                                        ccGAGGGCCTCGATCTGTGTGGCGATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAgg                                                                                                          >  1:579617/1‑139 (MQ=255)
                                                            gggCCTCGATCTGTGTGGCGATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCt                                                                                           >  1:1622556/1‑150 (MQ=255)
                                                            gggCCTCGATCTGTGTGGCGATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCt                                                                                           <  2:1170834/150‑1 (MQ=255)
                                                                cgcGATCTGTGTGGCGATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGAGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGAcc                                                                                       <  2:2806661/148‑1 (MQ=255)
                                                                           tgGCGATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGt                                                                            <  1:676124/150‑1 (MQ=255)
                                                                           tgGCGATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCAGAGGAAGGCCGGCTCCAACGTGTGCCCGGTGGTTGGGt                                                                            >  1:2789048/1‑150 (MQ=255)
                                                                             gCGATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCt                                                                                                                          <  2:572632/102‑1 (MQ=255)
                                                                             gCGATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCt                                                                                                                          >  1:572632/1‑102 (MQ=255)
                                                                              cGATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGa                                                                         >  2:2765339/1‑150 (MQ=255)
                                                                               gATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGAt                                                                        <  1:617185/150‑1 (MQ=255)
                                                                                aTGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATg                                                                       >  2:2046624/1‑150 (MQ=255)
                                                                                 tGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATgg                                                                      <  1:125903/150‑1 (MQ=255)
                                                                                  ggCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGc                                                                     <  2:2066495/150‑1 (MQ=255)
                                                                                    cTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGtt                                                                                >  1:2993692/1‑137 (MQ=255)
                                                                                     ttCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAgg                                                                                                          <  1:231813/110‑1 (MQ=255)
                                                                                     ttCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACg                                                                  >  1:697901/1‑150 (MQ=255)
                                                                                      tCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACgg                                                                 <  1:1558541/150‑1 (MQ=255)
                                                                                        aaGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACGGAt                                                               >  1:2947305/1‑150 (MQ=255)
                                                                                        aaGGAACTCTGCTAGTTGGTCGAGCGATGTTTCGCGGTAGGCATAGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACGGAt                                                               <  2:752365/150‑1 (MQ=255)
                                                                                          ggAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGTCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACGGATgg                                                             >  2:2913947/1‑150 (MQ=255)
                                                                                             aCTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGAt                                                                        >  2:1209626/1‑136 (MQ=255)
                                                                                                ctGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACGGATGGCGTcgc                                                       >  2:1704132/1‑150 (MQ=255)
                                                                                                 tGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGCGCTCGCCTGTGCCGCCGCGGTCGGAAGCATCCGACGGCTGACCGGGGGTTGGGTCGATGGCACGGATGGCGTcgcg                                                      >  1:2983225/1‑150 (MQ=255)
                                                                                                    tCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCTGCGCGCTGGCCTGGGcg                                                                                                                                      <  1:1406892/67‑2 (MQ=37)
                                                                                                    tCGTTGTTCGAGAGAAGTTTCGCGGTAGGCATCGCATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCCCGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCCACGTCTGACCGGAGGCTGGGTCGATGGCACGAATGGCGTCGCGGcc                                                   >  2:1401305/1‑149 (MQ=255)
                                                                                                          ttCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACGGATGGCGTCGCGGCTGCCGGt                                             >  1:2582265/1‑150 (MQ=255)
                                                                                                          ttCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACGGATGGCGTCGCGGCTGCCGGt                                             >  1:289907/1‑150 (MQ=255)
                                                                                                           tCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACGGATGGCGTCGCGGCTGCCGGtt                                            >  2:1474248/1‑150 (MQ=255)
                                                                                                            cGAGCGATGTTTCGCGGTTGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACGGATGGCGTCGCGGCTGCCGGTTa                                           <  2:1446856/150‑1 (MQ=255)
                                                                                                            cGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGAt                                                                        >  2:605539/1‑121 (MQ=255)
                                                                                                            cGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGAt                                                                        <  1:605539/121‑1 (MQ=255)
                                                                                                             gggCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAg                                                                                                           <  2:2141107/83‑1 (MQ=255)
                                                                                                             gAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACGGATGGCGTAGCGGCTGCCGGTTAc                                          >  2:2138653/1‑150 (MQ=255)
                                                                                                              aGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCGAACGTCTGACCGGTGGTGTGGTCGATGGCGCATAGTGCGTCGCGGCTGCCGGTTAcc                                         >  1:2307272/1‑150 (MQ=255)
                                                                                                              aGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACGGATGGCGTCGCGGCTGCCGGTTAcc                                         >  2:774760/1‑150 (MQ=255)
                                                                                                               gCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACGGATGGCGTCGCGGCTGCCGGTTACCg                                        <  1:1095948/150‑1 (MQ=255)
                                                                                                               gCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGAGGTTGGGTCGATGGCACGGATGGCGTCGCGGCTGCCGGTTACCg                                        >  1:2324195/1‑150 (MQ=255)
                                                                                                                cGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACGGATGGCGTCGCGGCTGCCGGTTACCgg                                       <  2:2582265/150‑1 (MQ=255)
                                                                                                                 gATGTTTCGCGGTAGGCATCGAATGCCGACCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACGGATGGCGTCGCGGCTGCCGGTTACCGGc                                      >  1:1274323/1‑150 (MQ=255)
                                                                                                                      ttCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACGGATGGCGTCGCGGCTGCCGGTTACCGGCCTCTg                                 <  1:2649935/150‑1 (MQ=255)
                                                                                                                       tCGCGGTAGGCATCGAATGCCGCCCCGGCCAGCGCGCAGGACTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTCGGTCGCTCGCAAGCATGGCGTCGCTGCTGCTGGTTCCCGGGctcttg                                <  2:2491313/150‑3 (MQ=255)
                                                                                                                        cgcgGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACGGATGGCGTCGCGGCTGCCGGTTACCGGCCTCTGGc                               <  1:413268/150‑1 (MQ=255)
                                                                                                                              aGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACGGATGGCGTCGCGGCTGCCGGTTACCGGCCTCTGGCCGATca                         >  2:1975485/1‑150 (MQ=255)
                                                                                                                               ggCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACGGATGGCGTCGCGGCTGCCGGTTACCGGCCTCTGGCCGATcag                        >  1:749571/1‑150 (MQ=255)
                                                                                                                                gCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACGGATGGCGTCGCGGCTGCCGGTTACCGGCCTCTGGCCGATcagc                       <  2:935189/150‑1 (MQ=255)
                                                                                                                                gCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACGGATGGCGTCGCGGCTGCCGGTTACCGGCCTCTGGCCGATcagc                       >  2:121093/1‑150 (MQ=255)
                                                                                                                                 cATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACGGATGGCGTCGCGGCTGCCGGTTACCGGCCTCTGGCCGATcagca                      >  2:1337857/1‑150 (MQ=255)
                                                                                                                                  aTCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACGGATGGCGTCGCGGCTGCCGGTTa                                           <  1:928913/128‑1 (MQ=255)
                                                                                                                                  aTCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACGGATGGCGTCGCGGCTGCCGGTTa                                           >  2:928913/1‑128 (MQ=255)
                                                                                                                                  aTCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACGGATGGCGTCGCGGCTGCCGGTTACCGGCCTCTGGCCGATcagcag                     >  2:1857861/1‑150 (MQ=255)
                                                                                                                                    cGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACGGATGGCGTCGCGGCTGCCGGTTACCGGCCTCTGGCCGATCAGCAGGt                   <  2:207458/150‑1 (MQ=255)
                                                                                                                                      aaTGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACGGATGGCGTCGCGGCTGCCGGTTACCGGCCTCTGGCCGATCAGCAGGTTg                 >  1:576092/1‑150 (MQ=255)
                                                                                                                                       aTGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACGGATGGCGTCGCGGCTGCCGGTTACCGGCCTCTGGCCGa                            >  2:853051/1‑138 (MQ=255)
                                                                                                                                       aTGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACGGATGGCGTCGCGGCTGCCGGTTACCGGCCTCTGGCCGa                            <  1:853051/138‑1 (MQ=255)
                                                                                                                                        tGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCg                                                                                      >  1:961045/1‑79 (MQ=255)
                                                                                                                                        tGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCg                                                                                      <  2:961045/79‑1 (MQ=255)
                                                                                                                                        tGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACGGATGGCGTCGCGGCTGCCGGTTACCGGCCTCTGGCCGATCAGCAGGTTGcc               <  1:1857861/150‑1 (MQ=255)
                                                                                                                                         gccgccCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGt                                                                            <  1:1961593/88‑1 (MQ=255)
                                                                                                                                         gccgccCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGt                                                                            >  2:1961593/1‑88 (MQ=255)
                                                                                                                                           cgccCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACgg                                                                 <  1:424629/97‑1 (MQ=255)
                                                                                                                                           cgccCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACGGATGGCGTCGCGGCTGCCGGTTACCGGCCTCTGGCCGATCAGCAGGTTGCCtgt            >  2:279438/1‑150 (MQ=255)
                                                                                                                                             cccAGGCCTGCGCGCTGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACGGATGGCGTCGCGGCTGCCg                                               <  1:2887884/113‑1 (MQ=255)
                                                                                                                                             cccAGGCCTGCGCGCTGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGGGGTTGGGTCGATGGCACGGATGGCGTCGCGGCTGCCg                                               >  2:2887884/1‑113 (MQ=255)
                                                                                                                                             cccAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACGGATGGCGTCGCGGCTGCCGGTTACCGGCCTCTGGCCGATCAGCAGGTTGCCTGTGa          >  2:1243253/1‑150 (MQ=255)
                                                                                                                                               cAGGCCTGCGCGCTGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCggt                                                                                    >  1:2121229/1‑74 (MQ=38)
                                                                                                                                               cAGGCCTGCGCGCTGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCggt                                                                                    <  2:2121229/74‑1 (MQ=38)
                                                                                                                                               cAGGCCTGCGCGCTGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACGGATGGCGTCGCGGCTGCCGGTTACCGGCct                                    >  1:1204604/1‑122 (MQ=255)
                                                                                                                                               cAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACGGATGGCGTCGCGGCTGCCGGTTACCGGCCTCTGGCCGATCAGCAGGTTGCCtgtg           >  1:167422/1‑147 (MQ=255)
                                                                                                                                               cAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACGGATGGCGTCGCGGCTGCCGGTTACCGGCCTCTGGCCGATCAGCAGGTTGCCtgtg           <  2:167422/147‑1 (MQ=255)
                                                                                                                                                 ggCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACGGATGGCGTCGCGGCTGCCGGTTACCGGCCTCTGGCCGATCAGCAGGTTGCCTGTGAgggg      <  2:3030360/150‑1 (MQ=255)
                                                                                                                                                  gCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACGGATGGCGTCGCGGCTGCCGGTTACCGGCCTCTGGCCGATCAGCAGGTTGCCTGTGAGGGGc     >  2:2814253/1‑150 (MQ=255)
                                                                                                                                                    cAGCGCGCAGGCCTGGGCCACGTGTTGGTCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACGGATGGCGTCGCGGCTGCCGGTTACCGGCCTCTGGCCGATCAGCAGGTTGCCTGTGAGGGGCAt   <  1:95/150‑1 (MQ=255)
                                                                                                                                                     aGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACGGATGGCGTCGCGGCTGCCGGTTACCGGCCTCTGGCCGATCAGCAGGTTGCCTGTGAGGGGCATg  >  2:1292190/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GAATGCGCGCCTTGGGCAGGCCGGTTTCGATCACGGCGCGGTCGATCAGGGCATCGCCGAGGGCCTCGATCTGTGTGGCGATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACGGATGGCGTCGCGGCTGCCGGTTACCGGCCTCTGGCCGATCAGCAGGTTGCCTGTGAGGGGCATG  >  NC_002947/1435673‑1435971

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑