Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,435,829 A→T 8.0% C48S (TGC→AGC)  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,8290AT 8.0% 425.4 / 23.8 138C48S (TGC→AGC) PP_1256alpha‑ketoglutarate semialdehyde dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (61/66);  new base T (5/6);  total (66/72)
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.69e-01

Only 100 of 141 total aligned reads displayed.

CGCCTTGGGCAGGCCGGTTTCGATCACGGCGCGGTCGATCAGGGCATCGCCGAGGGCCTCGATCTGTGTGGCGATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACGGATGGCGTCGCGGCTGCCGGTTACCGGCCTCTGGCCGATCAGCAGGTTGCCTGTGAGGGGCATGTAGCCT  >  NC_002947/1435680‑1435977
                                                                                                                                                     |                                                                                                                                                    
cgcCTTGGGCAGGCCGGTTTCGATCACGGCGCGGTCGATCAGGGCATCGCCGAGGGCCTCGATCTGTGTGGCGATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCa                                                                                                                                                      <  2:2944489/150‑1 (MQ=255)
   cTTGGGCAGGCCGGTTTCGATCACGGCGCGGTCGATCAGGGCATCGCCGAGGGCCTCGATCTGTGTGGCGATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGc                                                                                                                                                   >  2:1929419/1‑150 (MQ=255)
        gCAGGCCGGTTTCGATCACGGCGCGGTCGATCAGGGCATCGCCGAGGGCCTCGATCTGTGTGGCGATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCTGCGCGCTGGCCTggg                                                                                                                                              <  2:1691771/150‑1 (MQ=255)
              cGGTTTCGATCACGGCGCGGTCGATCAGGGCATCGCCGAGGGCCTCGATCTGTGTGGCGATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACgt                                                                                                                                        <  1:1929419/150‑1 (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:2375175/150‑1 (MQ=255)
                                ggTCGATCAGGGCATCGCCGACGGCCTCGATCTGTGTGGCGATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGa                                                                                                                      >  1:91024/1‑150 (MQ=255)
                                 gTCGATCAGGGCATCGCCGAGGGCCTCGATCTGTGTGGCGCTGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAg                                                                                                                     <  2:2491245/150‑1 (MQ=255)
                                    gATCAGGGCATCGCCGAGGGCCTCGATCTGTGTGGCGATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTa                                                                                                                  >  1:207458/1‑150 (MQ=255)
                                     aTCAGGGCATCGCCGAGGGCCTCGATCTGTGTGGCGATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCCCGTGTTCGCCTGTGCCGCCGAGGTAg                                                                                                                 <  1:1187602/150‑1 (MQ=255)
                                        aGGGCATCGCCGAGGGCCTCGATCTGTGTGGCGATGGATTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCTGCCGAGGTAGGcc                                                                                                              <  1:2334856/150‑1 (MQ=255)
                                         gggCATCGCCGAGGGCCTCGATCTGTGTGGCGATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCtg                                                                                                                               >  2:2696553/1‑132 (MQ=255)
                                         gggCATCGCCGAGGGCCTCGATCTGTGTGGCGATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCtg                                                                                                                               <  1:2696553/132‑1 (MQ=255)
                                           gCATCGCCGAGGGCCTCGATCTGTGTGGCGATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGc                                                                                                           >  1:208190/1‑150 (MQ=255)
                                             aTCGCCGAGGGCCTCGATCTGTGTGGCGATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTc                                                                                                         >  1:2604719/1‑150 (MQ=255)
                                                 ccGAGGGCCTCGATCTGTGTGGCGATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAgg                                                                                                                <  2:579617/139‑1 (MQ=255)
                                                 ccGAGGGCCTCGATCTGTGTGGCGATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAgg                                                                                                                >  1:580136/1‑139 (MQ=255)
                                                 ccGAGGGCCTCGATCTGTGTGGCGATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAgg                                                                                                                >  1:579617/1‑139 (MQ=255)
                                                 ccGAGGGCCTCGATCTGTGTGGCGATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAAc                                                                                                     <  2:2604719/150‑1 (MQ=255)
                                                   gAGGGCCTCGATCTGTGTGGCGATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCACGTTCCGCCGAGGTAGGCCGGCTCCAACGt                                                                                                   <  1:82428/150‑1 (MQ=255)
                                                     gggCCTCGATCTGTGTGGCGATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCt                                                                                                 <  2:1170834/150‑1 (MQ=255)
                                                     gggCCTCGATCTGTGTGGCGATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCt                                                                                                 >  1:1622556/1‑150 (MQ=255)
                                                          tCGATCTGTGTGGCGATGGCTTCAAGGACCTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCg                                                                                            <  1:1451434/150‑1 (MQ=255)
                                                            gATCTGTGTGGCGATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCggt                                                                                          <  2:2799957/150‑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)
                                                                          tGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATgg                                                                            <  1:125903/150‑1 (MQ=255)
                                                                           ggCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGc                                                                           <  2:2066495/150‑1 (MQ=255)
                                                                             cTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGtt                                                                                      <  2:2993692/137‑1 (MQ=255)
                                                                                 aaGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACGGAt                                                                     >  1:2947305/1‑150 (MQ=255)
                                                                                 aaGGAACTCTGCTAGTTGGTCGAGCGATGTTTCGCGGTAGGCATAGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACGGAt                                                                     <  2:752365/150‑1 (MQ=255)
                                                                                   ggAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGTCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACGGATgg                                                                   >  2:2913947/1‑150 (MQ=255)
                                                                                         ctGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACGGATGGCGTcgc                                                             >  2:1704132/1‑150 (MQ=255)
                                                                                          tGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGCGCTCGCCTGTGCCGCCGCGGTCGGAAGCATCCGACGGCTGACCGGGGGTTGGGTCGATGGCACGGATGGCGTcgcg                                                            >  1:2983225/1‑150 (MQ=255)
                                                                                             tCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCTGCGCGCTGGCCTGGGcg                                                                                                                                            >  2:1406892/1‑66 (MQ=37)
                                                                                               gttgttCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACGGATGGCGTCGCGGCTGc                                                       >  1:455682/1‑150 (MQ=255)
                                                                                                  gttCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACGGATGGCGTCGCGGCTGCCgg                                                    <  1:347987/150‑1 (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                                                                                                                 >  1:2141107/3‑85 (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:121093/1‑150 (MQ=255)
                                                                                                                         gCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACGGATGGCGTCGCGGCTGCCGGTTACCGGCCTCTGGCCGATcagc                             <  2:935189/150‑1 (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)
                                                                                                                                 tGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCg                                                                                            <  2:961045/79‑1 (MQ=255)
                                                                                                                                 tGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCg                                                                                            >  1:961045/1‑79 (MQ=255)
                                                                                                                                 tGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACGGATGGCGTCGCGGCTGCCGGTTACCGGCCTCTGGCCGATCAGCAGGTTGcc                     <  1:1857861/150‑1 (MQ=255)
                                                                                                                                  gccgccCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGt                                                                                  >  2:1961593/1‑88 (MQ=255)
                                                                                                                                  gccgccCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGt                                                                                  <  1:1961593/88‑1 (MQ=255)
                                                                                                                                  gccgccCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACGGATGGCGTCGCGGCTGCCGGTTACCGGCCTCTGGCCGATCAGCAGGTTGCCt                    <  1:618647/150‑1 (MQ=255)
                                                                                                                                    cgccCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACgg                                                                       >  2:424629/1‑97 (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                                                                                          <  2:2121229/74‑1 (MQ=38)
                                                                                                                                        cAGGCCTGCGCGCTGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACGGATGGCGTCGCGGCTGCCGGTTACCGGCct                                          >  1:1204604/1‑122 (MQ=255)
                                                                                                                                        cAGGCCTGCGCGCTGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACGGATGGCGTCGCGGCTGCCGGTTACCGGCct                                          <  2:1204604/122‑1 (MQ=255)
                                                                                                                                        cAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACGGATGGCGTCGCGGCTGCCGGTTACCGGCCTCTGGCCGATCAGCAGGTTGCCtgtg                 <  2:167422/147‑1 (MQ=255)
                                                                                                                                        cAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACGGATGGCGTCGCGGCTGCCGGTTACCGGCCTCTGGCCGATCAGCAGGTTGCCtgtg                 >  1:167422/1‑147 (MQ=255)
                                                                                                                                          ggCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACGGATGGCGTCGCGGCTGCCGGTTACCGGCCTCTGGCCGATCAGCAGGTTGCCTGTGAgggg            <  2:3030360/150‑1 (MQ=255)
                                                                                                                                           gCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACGGATGGCGTCGCGGCTGCCGGTTACCGGCCTCTGGCCGATCAGCAGGTTGCCTGTGAGGGGc           >  2:2814253/1‑150 (MQ=255)
                                                                                                                                             cAGCGCGCAGGCCTGGGCCACGTGTTGGTCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACGGATGGCGTCGCGGCTGCCGGTTACCGGCCTCTGGCCGATCAGCAGGTTGCCTGTGAGGGGCAt         <  1:95/150‑1 (MQ=255)
                                                                                                                                             cAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACGGATGGCGTCGCGGCTGCCGGTTACCGGCCTCTGGCCGATCAGCAGGTTGCCTGTGAGGGGCAt         >  2:2017978/1‑150 (MQ=255)
                                                                                                                                              aGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACGGATGGCGTCGCGGCTGCCGGTTACCGGCCTCTGGCCGATCAGCAGGTTGCCTGTGAGGGGCATg        >  1:932914/1‑150 (MQ=255)
                                                                                                                                              aGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACGGATGGCGTCGCGGCTGCCGGTTACCGGCCTCTGGCCGATCAGCAGGTTGCCTGTGAGGGGCATg        >  2:1292190/1‑150 (MQ=255)
                                                                                                                                               gcgcgcAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACGGATGGCGTCGCGGCTGCCGGTTACCGGCCTCTGGCCGATCAGCAGGTTGCCTGTGAGGGGCATGt       <  2:289907/150‑1 (MQ=255)
                                                                                                                                                 gcgcTGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTggg                                                                                   <  2:1468536/72‑1 (MQ=255)
                                                                                                                                                 gcgcTGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTggg                                                                                   >  1:1468536/1‑72 (MQ=255)
                                                                                                                                                  cgcAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACGGATGGCGTCGCGGCTGCCGGTTATCGACATCTGTGCGATGAGCAGTTTTCCTGTTGGGGTCATGTAGc    <  2:1274323/150‑1 (MQ=255)
                                                                                                                                                  cgcAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTAGGTCGATGGCACGGATGGCGTCGCGGCTGCCGGTTACCGGCCTCTGGCCGATCAGCAGGTTGCCTGTGAGGGGCATGTAGc    <  1:1556921/150‑1 (MQ=255)
                                                                                                                                                    cAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACGGATGGCGTCGCGGCTGCCGGTTACCGGCCTCTGGCCGATCAGCAGGTTGCCTGTGAGGGGCATGTAGcct  <  1:1975485/150‑1 (MQ=255)
                                                                                                                                                     aGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACGGATGGCGTCGCGGCTGCCGGTTACCGGCCTCTGGCCGATCAGCAGGTTg                       >  1:2486133/1‑128 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CGCCTTGGGCAGGCCGGTTTCGATCACGGCGCGGTCGATCAGGGCATCGCCGAGGGCCTCGATCTGTGTGGCGATGGCTTCAAGGAACTCTGCTCGTTGTTCGAGCGATGTTTCGCGGTAGGCATCGAATGCCGCCCAGGCCAGCGCGCAGGCCTGGGCCACGTGTTCGCCTGTGCCGCCGAGGTAGGCCGGCTCCAACGTCTGACCGGTGGTTGGGTCGATGGCACGGATGGCGTCGCGGCTGCCGGTTACCGGCCTCTGGCCGATCAGCAGGTTGCCTGTGAGGGGCATGTAGCCT  >  NC_002947/1435680‑1435977

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