Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,304,875 C→T 5.6% intergenic (+482/+233) eco → / ← mqo ecotin, a serine protease inhibitor/malate dehydrogenase, FAD/NAD(P)‑binding domain

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,304,8750CT 5.6% 242.7 / 5.5 93intergenic (+482/+233)eco/mqoecotin, a serine protease inhibitor/malate dehydrogenase, FAD/NAD(P)‑binding domain
Reads supporting (aligned to +/- strand):  ref base C (58/27);  new base T (3/2);  total (64/29)
Fisher's exact test for biased strand distribution p-value = 6.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.01e-01

GGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACG  >  NC_000913/2304726‑2305023
                                                                                                                                                     |                                                                                                                                                    
ggCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTc                                                                                                                                                      >  1:3853608/1‑150 (MQ=40)
ggCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCTCTGTGAAGTGCTCCACCCCGTAGGTc                                                                                                                                                      <  1:2550916/150‑1 (MQ=25)
 gCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTCGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCg                                                                                                                                                     >  2:636183/1‑150 (MQ=255)
  ccTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCgg                                                                                                                                                    >  1:1282881/1‑150 (MQ=255)
   cTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGa                                                                                                                                                   <  1:2608363/150‑1 (MQ=255)
   cTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGa                                                                                                                                                   >  1:189803/1‑150 (MQ=255)
    tACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGAt                                                                                                                                                  >  1:3830560/1‑150 (MQ=255)
    tACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGGTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGTAGTGCTCCACCCCGTAGGGGGGAt                                                                                                                                                  >  1:501638/1‑150 (MQ=255)
       acCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAg                                                                                                                                               >  2:4155133/1‑150 (MQ=255)
        ccGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAgg                                                                                                                                              >  2:2203497/1‑150 (MQ=255)
         cGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGc                                                                                                                                             >  2:80390/1‑150 (MQ=255)
          gCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCg                                                                                                                                            <  2:2073485/150‑1 (MQ=255)
                aaGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTAcg                                                                                                                                      <  2:3118918/150‑1 (MQ=255)
                aaGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCATACACAGCTGAGCCGTGCGCCACCCCGTAGGTCGGATCAGGCGGATAcg                                                                                                                                      >  2:1181993/1‑150 (MQ=14)
                 aGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTAcgc                                                                                                                                     >  1:2269703/1‑150 (MQ=255)
                  gTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCCAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTAcgcc                                                                                                                                    <  1:2495317/150‑1 (MQ=17)
                     cTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGACTACACCCCTGTGACGTGCGCCACCCCGTAGGTCGGATAAGGCGGAGAcgcagaa                                                                                                                                 >  1:2049516/1‑146 (MQ=14)
                      tCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCaa                                                                                                                                >  1:2792646/1‑149 (MQ=255)
                      tCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCAt                                                                                                                                >  2:3932726/1‑150 (MQ=255)
                      tCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCCCTGTGAAGTGCACCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCaa                                                                                                                                >  1:2928990/1‑149 (MQ=255)
                      tCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTATTATCAGGCCTACCCAGCTGGGAAGTGCTACACCCCGTAGGTCGGATAAGGCGGTTACGCCGCAt                                                                                                                                >  1:1046280/1‑150 (MQ=14)
                      tCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCg                                                                                                                                            >  2:300713/1‑138 (MQ=17)
                      tCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCg                                                                                                                                            <  1:300713/138‑1 (MQ=17)
                      ggcACCCCGTAGGTATGATAAGGCGGTGACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCAt                                                                                                                                <  1:4045605/148‑1 (MQ=255)
                       ccACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGATGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATc                                                                                                                               >  2:2680066/1‑150 (MQ=255)
                        cACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGCAGTGCACCACCCCGTTGGGCGGATAAGGCGGTTACGCCGCAAcc                                                                                                                              >  1:2845765/1‑150 (MQ=14)
                        cACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATcc                                                                                                                              >  2:1487279/1‑150 (MQ=255)
                        cACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATcc                                                                                                                              >  1:3270541/1‑150 (MQ=255)
                        cACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACACCGTAGTTCGGTTAAGGCGGTTACGCCGCATcc                                                                                                                              >  1:2307385/1‑150 (MQ=255)
                        cACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCGCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCAGcc                                                                                                                              >  1:757648/1‑150 (MQ=255)
                        cACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCGACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATcc                                                                                                                              >  1:2004348/1‑150 (MQ=255)
                        cACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCAAACACCGCTGTGAAGTGCTCCCCCCCGTAGGTCGGATAAGGCGCTTACGCCGCAtac                                                                                                                              >  2:875270/1‑148 (MQ=255)
                         aCCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCTCCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCg                                                                                                                             >  2:4035195/1‑150 (MQ=255)
                         aCCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCg                                                                                                                             >  1:1255715/1‑150 (MQ=255)
                          ccccGTAGGTTGGATAATGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGa                                                                                                                            <  2:874671/150‑1 (MQ=255)
                          ccccGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGa                                                                                                                            <  2:1897333/150‑1 (MQ=255)
                          ccccGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGa                                                                                                                            <  1:3567965/150‑1 (MQ=255)
                          ccccGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGa                                                                                                                            <  1:2644119/150‑1 (MQ=255)
                           cccGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGAc                                                                                                                           >  2:1319503/1‑150 (MQ=255)
                                   ttGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTGTTATAGGGCCTACACGGCTGTGAGGTGCTCCACCCCGTAGGACTGATAAGTCGGTTACGCCGCATCCGACATCTAACg                                                                                                                   >  2:2850564/1‑150 (MQ=14)
                                   ttGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCCCCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCTAACg                                                                                                                   >  1:682097/1‑150 (MQ=255)
                                   ttGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACg                                                                                                                   >  1:1735323/1‑150 (MQ=39)
                                   ttGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACg                                                                                                                   >  2:764135/1‑150 (MQ=39)
                                    tGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATc                                                                                                                               <  1:4134539/137‑1 (MQ=25)
                                    tGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATc                                                                                                                               >  2:4134539/1‑137 (MQ=25)
                                    tGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGc                                                                                                                  <  2:284935/150‑1 (MQ=38)
                                       aTAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTCCACCGCTGAGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCa                                                                                                               >  2:2441836/1‑150 (MQ=25)
                                       aTAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCa                                                                                                               >  2:2441796/1‑150 (MQ=35)
                                       aTAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTAGGATAATGCGGTTACGCCGCATCCGACATCTAACGCCCa                                                                                                               >  2:2441912/1‑150 (MQ=25)
                                          aGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGc                                                                                                            >  1:2186785/1‑150 (MQ=34)
                                            gCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCg                                                                                                          <  1:3166935/150‑1 (MQ=34)
                                            gCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGGTACGCCGCATCCGACATCTAACGCCCAAGCCg                                                                                                          >  1:1809941/1‑150 (MQ=21)
                                             cGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCATACACCGCTGTGACGTGCTCCACCCCGTTGGTCGGATCAGGCGGTTACGCCGCATCCGACATCTAACTCCccggccgg                                                                                                         >  2:3952546/1‑150 (MQ=11)
                                               gTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGa                                                                                                                                                   <  1:2369328/106‑1 (MQ=35)
                                               gTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGa                                                                                                                                                   >  2:2369328/1‑106 (MQ=35)
                                               gTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCCACATCTAACGCCCAAGCCGGtt                                                                                                       >  1:3472904/1‑150 (MQ=21)
                                               gTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTAGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCGAGCCGGtt                                                                                                       >  1:1644627/1‑150 (MQ=21)
                                                ttACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTg                                                                                                      <  1:1487279/150‑1 (MQ=33)
                                                 tACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGc                                                                                                     <  1:2167371/150‑1 (MQ=32)
                                                     ccgcTTGCGACCTCCTCCGTCCGAGCCGGTTGCGTGAGGCGACGCTTGATCGTCTTATCAGGCATTGACCGCGTTGAAGTGGACCTGCCCCTAGGTCGGATAAGGGGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTgt                                                                                                 <  2:1473672/150‑2 (MQ=2)
                                                          tCCGCCCTCCTACGTCCGTGCCGGTTGCCTGATGCGTCGCGGGCGGGTCTTATCAGGCCTACACCGCTGGGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGa                                                                                            <  1:2474414/150‑1 (MQ=11)
                                                          tCCGACATCCACCGTCCGCGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGa                                                                                            <  1:2474399/150‑1 (MQ=17)
                                                            cGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAATGCGGTTACGCCACAACCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACg                                                                                          >  1:1175626/1‑150 (MQ=12)
                                                            cGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGGCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGGTGCCTGATGCGACg                                                                                          >  1:1369877/1‑150 (MQ=17)
                                                            cGACATCCAACGTCCGAGCCGGTTGCCAGATGCGACGCGGGCGCGTCTTATCAGGCCTACACCGCGCGGAACTGCGCCCCCCCGTAGGTCGGCTCAGGCCGGTACGCCACCACCGACCTCGATAGCCCAAGGCGGGTGGCTGATGCGACg                                                                                          >  1:1175572/1‑150 (MQ=2)
                                                             gACATCCAACGTCCGAGCCGGTTTCCTGATGCGACGCTGGCGCGTCTTTTCAGGCCTACACCGCTGTGAAGTTCTCCACCCCGTAGGTCGGATAATGCGGTTACGCCGCATCCGACATCTAAAGCCCAAGACGGTTTCCTGATGCGACGc                                                                                         >  1:2073828/1‑150 (MQ=11)
                                                             gACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGc                                                                                         >  1:3913404/1‑150 (MQ=32)
                                                             gACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGc                                                                                         >  2:2237827/1‑150 (MQ=32)
                                                              aCATCCAAAGTCCGAGCCGGTTGCCTGATGCGACTATGGATGGTCTTATGTGTGCTAAAACGCTTTGAATTGCTCCACCCCGTAGGTCTGAAAAGGCGGTTACGCCGCATCCGACATCAAACGCCCTAGCCGGTTGCCTGATGCGACGCt                                                                                        <  1:2638777/150‑1 (MQ=11)
                                                                 tCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGGTAAGGCGGTTAAGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGgc                                                                                     >  1:257260/1‑150 (MQ=17)
                                                                 tCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGgc                                                                                     <  1:809978/150‑1 (MQ=32)
                                                                  ccAACGTCCGAGCCGGTTTCCTGATGCGACGCTGGCGCGTCTTATCAGGCCGACCCCGCTGTGTAGTGCTCCACCCCGGAGGTCGGATAAGGCGGGTACGCCGCATCCGACAGCTAACGCCCCAGCCGGTTGCCTGATGCGACGCTGgcg                                                                                    >  1:44674/1‑150 (MQ=11)
                                                                  ccAACGTCCGAGCCGGTTGCCAGATGCGACGCTGGCGCGGCTTATCAGGCCTACACCGCTGTGTAGTGCACCACCCCGGAGGTCGGATAAGGCGGTTACGCCGCCGCCGACCTCTAACGCCCCAGCCGGTTGCCTGATGCGACGCTGgcg                                                                                    >  1:44675/1‑150 (MQ=11)
                                                                   cAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCTTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGgcgc                                                                                   <  1:3355451/150‑1 (MQ=17)
                                                                   cAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGgcgc                                                                                   <  2:2274972/150‑1 (MQ=32)
                                                                   cAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGgcgc                                                                                   <  1:3355544/150‑1 (MQ=32)
                                                                     aCGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACGCCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGt                                                                                 >  2:4106997/1‑150 (MQ=17)
                                                                     aCGTCCGAGCCGGTTGCCTGACGCGCCGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGt                                                                                 <  1:3758044/150‑1 (MQ=17)
                                                                      cGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTc                                                                                <  1:2022388/150‑1 (MQ=32)
                                                                      cGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTc                                                                                <  1:1319503/150‑1 (MQ=32)
                                                                        tCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGAGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCtt                                                                              >  2:3015533/1‑150 (MQ=11)
                                                                                                                                          aCCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCa            >  1:3946225/1‑150 (MQ=11)
                                                                                                                                           ccccGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCaa           >  1:72062/1‑150 (MQ=31)
                                                                                                                                           ccccGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCaa           >  1:2887058/1‑150 (MQ=31)
                                                                                                                                            cccGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATCAGGCGGATACGCCGCATCCAAc          >  2:3671159/1‑150 (MQ=17)
                                                                                                                                            cccGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAAc          >  1:3788719/1‑150 (MQ=32)
                                                                                                                                            cccGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGCCCTACACCCCTCAGAAGTGCGCCACCCCGTATGGTGGATAAGGCGGAGACGCCGCATCCAAc          >  2:3373431/1‑150 (MQ=11)
                                                                                                                                             ccGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACa         >  2:3489890/1‑150 (MQ=32)
                                                                                                                                             ccGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACa         >  2:3489855/1‑150 (MQ=32)
                                                                                                                                                 aGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCCAGCCGGTTGCCTGAAGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTa     <  2:2465556/150‑1 (MQ=18)
                                                                                                                                                   gTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGATGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGAGGTTACGCCGCATCCAACATCTAAc   >  2:577848/1‑150 (MQ=21)
                                                                                                                                                    tCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACg  >  1:2552469/1‑150 (MQ=34)
                                                                                                                                                    tCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACg  >  1:2552317/1‑150 (MQ=34)
                                                                                                                                                    tCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACg  <  1:1104392/150‑1 (MQ=34)
                                                                                                                                                    tCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCGTATCAGGCCTCCCCCGCTGGGGCGTGCTCCACCCCGTAGGTGGGAGAAGGCGGTGACGCCGCCTccaccagccaccg  >  1:3318865/1‑143 (MQ=11)
                                                                                                                                                     cGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACg  >  1:3062879/1‑149 (MQ=34)
                                                                                                                                                     cGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACg  <  2:3062879/149‑1 (MQ=34)
                                                                                                                                                     cGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACg  <  2:3062874/149‑1 (MQ=34)
                                                                                                                                                     |                                                                                                                                                    
GGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCGACATCCAACGTCCGAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCTAACGCCCAAGCCGGTTGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTACACCGCTGTGAAGTGCTCCACCCCGTAGGTTGGATAAGGCGGTTACGCCGCATCCAACATCTAACG  >  NC_000913/2304726‑2305023

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