Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,731,930 A→T 6.8% L661Q (CTG→CAG)  PP_1526 ← beta‑(1‑3)‑glucosyl transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,731,9300AT 6.8% 220.5 / 6.5 73L661Q (CTG→CAG) PP_1526beta‑(1‑3)‑glucosyl transferase
Reads supporting (aligned to +/- strand):  ref base A (31/37);  new base T (2/3);  total (33/40)
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.11e-01

CCCGATGCGGCACGATGATCATCGCCGCCGACCACAACAGCGCGCCGATGGTGAAGAAGATGTTCATGCCATCGGCGATCCACGGTAGCCAGCCGGCCAGGAAGTGGTAGCGCTGGCCACGGGTCAGCTGGCTGCCTTTGCCGCGCAGCAGCGCGCCGGCGTGGTGCTTGATGATCTGGATGGCGCCGTAGGCCCAGCGGAAGCGTTGCTTCTTGAAGTCGATGAAGGTGTCGGGCATCAGGCCCTTGCCGTAGCTGTTGTGGGCGTAGGCGGCGGACAGGCCTTTCTC  >  NC_002947/1731781‑1732069
                                                                                                                                                     |                                                                                                                                           
cccGATGCGGCACGATGATCATCGCCGCCGACCACAACAGCGCGCCGATGGTGAAGAAGATGTTCATGCCATCGGCGATCCACGGTAGCCAGCCGGCCAGGAAGTGGTAGCGCTGGCCACGGGTCAGCTGGCTGCCTTTGCCGCgcagca                                                                                                                                             >  1:207321/1‑150 (MQ=255)
  cGATGCGGCACGATGATCATCGCCGCCGACCACAACAGCGCGCCGATGGTGAAGAAGATGTTCATGCCATCGGCGATCCACGGTAGCCAGCCGGCCAGGAAGTGGTAGCGCTGGCCACGGGTCAGCTGGCTGCCTTTGCCGCGCagcagc                                                                                                                                           <  1:494825/150‑1 (MQ=255)
   gATGCGGCACGATGATCATCGCCGCCGACCACAACAGCGCGCCGATGGTGAAGAAGATGTTCATGCCATCGGCGATCCACGGTAGCCAGCCGGCCAGGAAGTGGTAGCGCTGGCCACGGGTCAGCTGGCTGCCTTTGCCGCGCAGCAgcg                                                                                                                                          >  2:1044396/1‑150 (MQ=255)
         gCACGATGATCATCGCCGCCGACCACAACCGCGCGCCGATGGTGAAGAAGATGTTCATGCCATCGGCGATCCACGGTAGCCAGCCGGCCAGGAAGTGGTAGCGCTGGCCACGGGTCAGCTGGCTGCCTTTGCCGCGCAGCAGCGCGCCgc                                                                                                                                    >  2:384409/1‑149 (MQ=255)
         gCACGATGATCATCGCCGCCGACCACAACAGCGCGCCGATGGTGAAGAAGATGTTCATGCCATCGGCGATCCACGGTAGCCAGCCGGCCAGGAAGTGGTAGCGCTGGCCACGGGTCAGCTGGCTGCCTTTGCCGCGCAGCAGCGCGCCgg                                                                                                                                    >  2:449204/1‑150 (MQ=255)
         gCACGATGATCATCGCCGCCGACCACAACAGCGCGCCGATGGTGAAGAAGATGTTCATGCCATCGGCGATCCACGGTAGCCAGCCGGCCAGGAAGTGGTAGCGCTGGCCACGGGTCAGCTGGCTGCCTTTGCCGCGCAGCAGCGCGCCgg                                                                                                                                    >  1:1276370/1‑150 (MQ=255)
         gCACGATGATCATCGCCGCCGACCACAACAGCGCGCCGATGGTGAAGAAGATGTTCATGCCATCGGCGATCCACGGTAGCCAGCCGGCCAGGAAGTGGTAGCGCTGGCCACGGGTCAGCTGGCTGCCTTTGCCGCGCAGCAGCGCGCCgg                                                                                                                                    >  1:504127/1‑150 (MQ=255)
          cACGATGATCATCGCCGCCGACCACAACAGCGCGCCGATGGTGAAGAAGATGTTCATGCCATCGGCGATCCACGGTAGCCAGCCGGCCAGGAAGTGGTAGCGCTGGCCACGGGTCAGCTGGCTGCCTTTGCCGCGCTGCTgcgcggcaaa                                                                                                                                   <  2:287838/150‑6 (MQ=255)
           aCGATGATCATCGCCGCCGACCACAACAGCGCGCCGATGGTGAAGAAGATGTTCATGCCATCGGCGATCCACGGTAGCCAGCCGGCCAGGAAGTGGTAGCGCTGGCCACGGGTCAGCTGGCTGCCTTTGCCGCGCAGCAGCGCGCCGGCg                                                                                                                                  <  2:828364/150‑1 (MQ=255)
            cGATGATCATCGCCGCCGACCACAACAGCGCGCCGATGGTGAAGAAGATGTTCATGCCATCGGCGATCCACGGTAGCCAGCCGGCCAGGAAGTGGTAGCGCTGGCCACGGGTCAGCTGGCTGCCTTTGCCGCGCAGCAGCGCGCCGGCgt                                                                                                                                 <  2:1011090/150‑1 (MQ=255)
             gatgatCATCGCCGCCGACCACAACAGCGCGCCGATGGTGAAGAAGATGTTCATGCCATCGGCGATCCACGGTAGCCAGCCGGCCAGGAAGTGGTAGCGCTGGCCAAGGGTCGGCGTGCGGCCTTTGCCGCGCAGCAGCGCGCCGGCgtg                                                                                                                                >  2:473949/1‑150 (MQ=255)
               tgatCATCGCCGCCGACCACAACAGCGCGCCGATGGTGAAGAAGATGTTCATGCCATCGGCGATCCACGGTAGCCAGCCGGCCAGGAAGTGGTAGCGCTGGCCACGGGTCAGCTGGCTGCCTTTGCCGCGCAGCAGCGCGCCGGCgtggt                                                                                                                              <  2:414766/150‑1 (MQ=255)
                  tcatcGCCGCCGACCACAACAGCGCGCCGATGGTGAAGAAGATGTTCATGCCATCGGCGATCCACGGTAGCCAGCCGGCCAGGAAGTGGTAGCGCTGGCCACGGGTCAGCTGGCTGCCTTTGCCGCGCAGCAGCGCGCCGGCGTGGTGCt                                                                                                                           <  2:1276370/150‑1 (MQ=255)
                     tcGCCGCCGACCACAACAGCGCGCCGATGGTGAAGAAGATGTTCATGCCATCGGCGATCCACGGTAGCCAGCCGGCCAGGAAGTGGTAGCGCTGGCCACGGGTCAGCTGGCTGCCTTTGCCGCgcagca                                                                                                                                             <  2:952173/129‑1 (MQ=255)
                     tcGCCGCCGACCACAACAGCGCGCCGATGGTGAAGAAGATGTTCATGCCATCGGCGATCCACGGTAGCCAGCCGGCCAGGAAGTGGTAGCGCTGGCCACGGGTCAGCTGGCTGCCTTTGCCGCgcagca                                                                                                                                             >  1:952173/1‑129 (MQ=255)
                     tcGCCGCCGACCACAACAGCGCGCCGATGGTGAAGAAGATGTTCATGCCATCGGCGATCCACGGTAGCCAGCCGGCCAGGAAGTGGTAGCGCTGGCCACGGGTCAGCTGGCTGCCTTTGCCGCGCAGCAGCGCGCCGGCGTGGTGCTtga                                                                                                                        <  2:207321/150‑1 (MQ=255)
                      cgccgcCGACCACAACAGCGCGCCGATGGTGAAGAAGATGTTCATGCCATCGGCGATCCACGGTAGCCAGCCGGCCAGGAAGTGGTAGCGCTGGCCACGGGTCAGCTGGCTGCCTTTGCCGCGCAGCAGCGCGCCGGCGTGGTGCTtgat                                                                                                                       <  1:230351/150‑1 (MQ=255)
                             gACCACAACAGCGCGCCGATGGTGAAGAAGAGGTTCATGCCATCGGCGATCCACGGTAGCCAGCCGGCCAGGAAGTGGTAGCGCTGGCCACGGGTCAGCTGGCTGCCTTTGCCGCGCAGCAGCGCGCCGGCGTGGTGCTTGATGATCTgg                                                                                                                >  2:876797/1‑150 (MQ=255)
                               ccACAACAGCGCGCCGATGGTGAAGAAGATGTTCATGCCATCGGCGATCCACGGTAGCCAGCCGGCCAGGAAGTGGTAGCGCTGGCCACGGGTCAGCTGGCTGCCTTTGCCGCGCAGCAGCGCGCCGGCGTGGTGCTTGATGATCTGGAt                                                                                                              <  1:1044396/150‑1 (MQ=255)
                                 acaacaGCGCGCCGATGGTGAAGAAGATGTTCATGCCATCGGCGATCCACGGTAGCCAGCCGGCCAGGAAGTGGTAGCGCTGGCCACGGGTCAGCTGGCTGCCTTTGCCGCGCAGCAGCGCGCCGGCGTGGTGCTTGATGATCTGGATgg                                                                                                            >  2:207163/1‑150 (MQ=255)
                                         gcgcCGATGGTGAAGAAGATGTTCATGCCATCGGCGATCCACGGTAGCCAGCCGGCCAGGAAGTGGTAGCGCTGGCCACGGGTCAGCTGGCTGCCTTTGCCGCGCAGCAGCGCGCCGGCGTGGTGCTTGATGATCTGGATGGCGCCGTAg                                                                                                    >  2:1053696/1‑150 (MQ=255)
                                              gATGGTGAAGAAGATGTTCATGCCATCGGCGATCCACGGTAGCCAGCCGGCCAGGAAGTGGTAGCGCTGGCCACGGGTCAGCTGGCTGCCTTTGCCGCGCAGCAGCGCGCCGGCGTGGTGCTTGATGATCTGGATGGCGCCGTAGGCCCa                                                                                               <  1:328202/150‑1 (MQ=255)
                                                tGGTGAAGAAGATGTTCATGCCATCGGCGATCCACGGTAGCCAGCCGGCCAGGAAGTGGTAGCGATGGCCACGGCTGAGCTGGCTGCCTTTGCCGCGCAGCAGCGCGCCGGCGTGGTGCTTGATGATCTGGATGGCGCCGTAGGCCCAGc                                                                                             <  1:1242133/150‑1 (MQ=255)
                                                  gTGAAGAAGATGTTCATGCCATCGGCGATCCACGGTAGCCAGCCGGCCAGGAAGTGGTAGCGCTGGCCACGGGTCAGCTGGCTGCCTTTGCCGCGCAGCAGCGCGCCGGCGTGGTGCTTGATGATCTGGATGGCGCCGTAGGCCCAGCgg                                                                                           >  2:1403119/1‑150 (MQ=255)
                                                   tGAAGAAGATGTTCATGCCATCGGCGATCCACGGTAGCCAGCCGGCCAGGAAGTGGTAGCGCTGGCGCCGGGTCACCTGGCTGCCTTTGCCGCGCAGCAGCGCGCCGGCGCGGTGCTTGATGATCTGGATCACGCCATAGGCCCAGCggc                                                                                          >  2:121512/1‑149 (MQ=255)
                                                      agaagaTGTTCATGCCATCGGCGATCCACGGTAGCCAGCCGGCCAGGAAGTGGTAGCGCTGGCCACGGGTCAGCTGGCTGCCTTTGCCGCGCAGCAGCGCGCCGGCGTGGTGCTTGATGATCTGGATGGCGCCGTAGGCCCAGCGGAAGc                                                                                       <  2:494280/150‑1 (MQ=255)
                                                       gaagaTGTTCATGCCATCGGCGATCCACGGTAGCCAGCCGGCCAGGAAGTGGTAGCGCTGGCCACGGGTCAGCTGGCTGCCTTTGCCGCGCAGCAGCGCGCCGGCGTGGTGCTTGATGATCTGGATGGCGCCGTAGGCCCAGCGGAAGCg                                                                                      >  2:192708/1‑150 (MQ=255)
                                                        aagaTGTTCATGCCATCGGCGATCCACGGTAGCCAGCCGGCCAGGAAGTGGTAGCGCTGGCCACGGGTCAGCTGGCTGCCTTTGCCGCGCAGCAGCGCGCCGGCGTGGTGCTTGATGATCTGGATGGCGCCGTAGGCCCAGCGGAAGCGt                                                                                     >  1:148878/1‑150 (MQ=255)
                                                        aagaTGTTCATGCCATCGGCGATCCACGGTAGCCAGCCGGCCAGGAAGTGGTAGCGCTGGCCACGGGTCAGCTGGCTGCCTTTGCCGCGCAGCAGCGCGCCGGCGTGGTGCTTGATGATCTGGATGGCGCCGTAGGCCCAGCGGAAGCGt                                                                                     >  1:1068375/1‑150 (MQ=255)
                                                         agaTGTTCATGCCATCGGCGATCCACGGTAGCCAGCCGGCCAGGAAGTGGTAGCGCTGGCCACGGGTCAGCTGGCTGCCTTTGCCGCGCAGCAGCGCGCCGGCGTGGTGCTTGATGATCTGGATGGCGCCGTAGGCCCAGCGGAAGCGtt                                                                                    <  1:769907/150‑1 (MQ=255)
                                                         agaTGTTCATGCCATCGGCGATCCACGGTAGCCAGCCGGCCAGGAAGTGGTAGCGCTGGCCACGGGTCAGCTGGCTGCCTTTGCCGCGCAGCAGCGCGCCGGCGTGGTGCTTGATGATCTGGATGGCGCCGTAGGCCCAGCGGAAGCGtt                                                                                    <  2:1068375/150‑1 (MQ=255)
                                                              ttCATGCCATCGGCGATCCACGGTAGCCAGCCGGCCAGGAAGTGGTAGCGCTGGCCACGGGTCAGCTGGCTGCCTTTGCCGCGCAGCAGCGCGCCGGCGTGGTGCTTGATGATCTGGATGGCGCCGTAGGCCCAGCGGAAGCGTTGcttc                                                                               >  1:1338072/1‑150 (MQ=255)
                                                                cATGCCATCGGCGATCCACGGTAGCCAGCCGGCCAGGAAGTGGTAGCGCTGGCCACGGGTCAGCTGGCTGCCTTTGCCGCGCAGCAGCGCGCCGGCGTGGTGCTTGATGATCTGGATGGCGCCGTAGGCCCAGCGGAAGCGTTGcttctt                                                                             <  2:498336/150‑1 (MQ=255)
                                                                   gCCATCGGCGATCCACGGTAGCCAGCCGGCCAGGAAGTGGTAGCGCTGGCCACGGGTCAGCTGGCTGCCTTTGCCGCGCAGCAGCGCGCCGGCGTGGTGCTTGATGATCTGGCAGGCGCCGTAGGCCCAGCGGAAGCGTTGCTTCTTGaa                                                                          <  1:473949/150‑1 (MQ=255)
                                                                     cATCGGCGATCCACGGTAGCCAGCCGGCCAGGAAGTGGTAGCGCTGGCCACGGGTCAGCTGGCTGCCTTTGCCGCGCAGCAGCGCGCCGGCGTGGTGCTTGATGATCTGGATGGCGCCGTAGGCCCAGCGGAAGCGTTGCTTCTTGAAGt                                                                        >  1:139251/1‑150 (MQ=255)
                                                                      aTCGGCGATCCACGGTAGCCAGCCGGCCAGGAAGTGGTAGCGCTGGCCACGGGTCAGCTGGCTGCCTTTGCAGTGCAGCAGCGCGCCGCCGTGGTGCTTGATGATCTGGATGGCGCCGTAGGCCCAGCGGAAGCGTTGCTTCTTGAAGTc                                                                       >  1:806661/1‑150 (MQ=255)
                                                                       tCGGCGATCTACGGTAGCCAGCCGGCCAGGAAGTGGTAGCGCTGGCCACGGGTCAGCTGGCTGCCTTTGCCGCGCAGCAGCGCGCCGGCGTGGTGCTTGATGATCTGGATg                                                                                                             <  2:951385/111‑1 (MQ=255)
                                                                       tCGGCGATCCACGGTAGCCAGCCGGCCAGGAAGTGGTAGCGCTGGCCACGGGTCAGCTGGCTGCCTTTGCCGCGCAGCAGCGCGCCGGCGTGGTGCTTGATGATCTGGATg                                                                                                             >  1:951385/1‑111 (MQ=255)
                                                                        cGGCGATCCACGGTAGCCAGCCGGCCAGGAAGTGGTAGCGCTGGCCACGGGTCAGCTGGCTGCCTTTGCCGCGCAGCAGCGCGCCGGCGTGGTGCTTGATGATCTGGATGGCGCCGTAGGCCCAGCGGAAGCGTTGCTTCTTGAAGTCGa                                                                     <  1:1062784/150‑1 (MQ=255)
                                                                          gCGATCCACGGTAGCCAGCCGGCCAGGAAGTGGTAGCGCTGGCCACGGGTCAGCTGGCTGCCTTTGCCGCGCAGCAGCGCGCCGGCGTGGTGCTTGATGATCTGGATGGCGCCGTAGGCCCAGCGGAAGCGTTGCTTCTTGAAGTCGATg                                                                   >  2:418575/1‑150 (MQ=255)
                                                                            gATCCACGGTAGCCAGCCGGCCAGGAAGTGGTAGCACTGGCCACGGGTCAGTTGTCATCGTTTGTCTCGCATCAGCGCGCCGGCGTGGTGCTTGATGATCTGGATGGCGACGTAGGCCCAGAGGAAGCGTTGCTTCTTGAAGTCGATGaa                                                                 <  2:1379126/150‑1 (MQ=255)
                                                                             atctaaggTCGCCCGCCTGCCAGGAAGTGGTAGCGCTGGCCACGGGTCAGCTCGAAGCCCTTGCCCCGCATCAGCGCGCCGGCGTGGTGCTTGATGCTCTAGATGGCACCGTTTGCCTAGTGTAAGCGTTGCTTCTTGAAGTCGATGAAg                                                                <  1:678599/144‑1 (MQ=255)
                                                                             aTCCACGGTAGCCAGCCGGCCAGGAAGTGGTAGCGCTGGCCACGGGTCAGCTGGCTGCCTTTGCCGCGCAGCAGCGCGCCGGCGTGGTGCTTGATGATCTGGATGGCGCCGTAGGCCCAGCGGAAGCGTTGCTTCTTGAAGTCGATGAAg                                                                <  2:421297/150‑1 (MQ=255)
                                                                             aTCCACGGTAGCCAGCCGGCCAGGAAGTGGTAGCGCTGGCCACGGGTCAGCTGGCTGCCTTTGCCGCGCAGCAGCGCGCCGGCGTGGTGCTTGATGATCTGGATGGCGCCGTAGGCCCAGCGGAAGCGTTGCTTCTTGAAGTCGATGAAg                                                                <  2:139251/150‑1 (MQ=255)
                                                                                cACGGTAGCCAGCCGGCCAGGAAGTGGTAGCGCTGGCCACGGGTCAGCTGGCTGCCTTTGCCGCGCAGCAGCGCGCCGGCGTGGTGCTtg                                                                                                                         >  2:133296/1‑90 (MQ=255)
                                                                                cACGGTAGCCAGCCGGCCAGGAAGTGGTAGCGCTGGCCACGGGTCAGCTGGCTGCCTTTGCCGCGCAGCAGCGCGCCGGCGTGGTGCTtg                                                                                                                         <  1:133296/90‑1 (MQ=255)
                                                                                 aCGGTAGCCAGCCGGCCAGGAAGTGGTAGCGCTGGCCACGGGTCAGCTGGCTGCCTTTGCCGCGCAGCAGCGCGCCGGCGTGGTGCTTGATGATCTGGATGGCGCCGTAGGCCCAGCGGAAGCGTTGCTTCTTGAAGTCGATGAAGgtgt                                                            >  2:291791/1‑150 (MQ=255)
                                                                                  cGGTAGCCAGCCGGCCAGGGAGTGGTAGCGCTGGCCACGGGTCAGCGGGCTGCCTTTGCCGCGCAGCAGCGCGCCGGCGTGGTTCTTGATGATCTGGATGGCGCCGTAGGCCCAGCGGAAGCGTTGCTTCTTGTAGTCGATGACGGCGTc                                                           >  2:678599/1‑150 (MQ=255)
                                                                                  cGGTAGCCAGCCGGCCAGGAAGTGGTAGCGCTGGCCACGGGTCAGCTGGCTGCCTTTGCCGCGCAGCAGCGCGCCGGCGTGGTGCTTGATGATCTGGATGGGGCCGTAGGCCCAGCGGAAGCGTTGCTTCTTGAAGTCGAGGAAGGTGTc                                                           <  2:763871/150‑1 (MQ=255)
                                                                                    gTAGCCAGCCGGCCAGGAAGTGGTAGCGCTGGCCACGGGTCAGCTGGCTGCCTTTGCCGCGCAGCAGCGCGCCGGCGTGGTGCTTGATGATCTGGATGGCGCCGTAGGCCCAGCGGAAGCGTTGCTTCTTGAAGTCGATGAAGGTGTCgg                                                         <  2:18031/150‑1 (MQ=255)
                                                                                         cagccGGCCAGGAAGTGGTAGCGCTGGCCACGGGTCAGCTGGCTGCCTTTGCCGCGCAGCAGCGCGCCGGCGTGGTGCTTGATGATCTGGATGGCGCCGTAGGCCCAGCGGAAGCGTTGCTTCTTGAAGTCGATGAAGGTGTCGGGCATc                                                    >  2:608873/1‑150 (MQ=255)
                                                                                                  aGGAAGTGGTAGCGCTGGCCACGGGTCAGCTGGCTGCCTTTGCCGCGCAGCAGCGCGCCGGCGTGGTGCTTGATGATCTGGATGGCGCCGTAGGCCCAGCGGAAGCGTTGCTTCTTGAAGTCGATGAAGGTGTCGGGCATCAGGCCCTTg                                           >  1:1025781/1‑150 (MQ=255)
                                                                                                    gAAGTGGTAGCGCTGGCCACGGGTCAGCTGGCTGCCTTTGCCGCGCAGCAGCGCGCCGGCGTGGTGCTTGATGATCTGGATGGCGCCGTAGGCCCAGCGGAAGCGTTGCTTCTTGAAGTCGa                                                                     <  2:1051654/122‑1 (MQ=255)
                                                                                                    gAAGTGGTAGCGCTGGCCACGGGTCAGCTGGCTGCCTTTGCCGCGCAGCAGCGCGCCGGCGTGGTGCTTGATGATCTGGATGGCGCCGTAGGCCCAGCGGAAGCGTTGCTTCTTGAAGTCGa                                                                     >  1:1051654/1‑122 (MQ=255)
                                                                                                    gAAGTGGTAGCGCTGGCAACGGGTCAGCTGGCTGCCTTTGCCGCGCAGCAGCGCGCCGGTGTGGTGCTTGATGATCTGGATGGCGCCGTAGGCCCAGCGGAAGCGTTGCTTCTTGAAGTCGATGAAGGTGTCGGGCATCAGGCCCTTGcc                                         <  1:127128/150‑1 (MQ=255)
                                                                                                          gTAGCGCTGGCCACGGGTCAGCTGGCTGCCTTTGCCGCGCAGCAGCGCGCCGGCGTGGTGCTTGATGATCTGGATGGCGCCGTAGGCCCAGCGGAAGCGTTGCTTCTTGAAGTCGATGAAGGTGTCGGGCATCAGGCCCTTGCCGTAGCt                                   <  1:876797/150‑1 (MQ=255)
                                                                                                           tAGCGCTGGCCACGGGTCAGCTGGCTGCCTTTGCCGCGCAGCAGCGCGCCGGCGTGGTGCTTGATGATCTGGATGGCGCCGTAGGCCCAGCGGAAGCGTTGCTTCTTGAAGTCGATGAAGGTGTc                                                           >  2:735978/1‑125 (MQ=255)
                                                                                                           tAGCGCTGGCCACGGGTCAGCTGGCTGCCTTTGCCGCGCAGCAGCGCGCCGGCGTGGTGCTTGATGATCTGGATGGCGCCGTAGGCCCAGCGGAAGCGTTGCTTCTTGAAGTCGATGAAGGTGTc                                                           <  1:735978/125‑1 (MQ=255)
                                                                                                             gcgcTGGCCACGGGTCAGCTGGCTGCCTTTGCCGCGCAGCAGCGCGCCGGCGTGGTGCTTGATGATCTGGATGGCGCCGTAGGCCCAGCGGAAGCGTTGCTTCTTGAAGTCGATGAAGGTGTCGGGCATCAGGCCCTTGCCGTAGCtgtt                                <  1:211100/150‑1 (MQ=255)
                                                                                                             gcgcGGGCCACGGGTCAGCTGGCTGCCTTTGCCGCGCAGCAGCGCGCCGGCGTGGTGCTTGATGATCTGGATGGCGCCGTAGGCCCAGCGGAAGCGTTGCTTCTTGAAGTCGATGAAGGTGTCGGGCATCAGGCCCTTGCCGTAGCtgtt                                <  2:535103/150‑1 (MQ=255)
                                                                                                              cgcTGGCCACGGGTCAGCTGGCTGCCTTTGCCGCGCAGCAGCGCGCCGGCGTGGTGCTTGATGATCTGGATGGCGCCGTAGGCCCAGCGGAAGCGTTGCTTCTTGAAGTCGATGAAGGTGTCGGGCATCAGGCCCTTGCCGTAGCtgttg                               <  1:121512/150‑1 (MQ=255)
                                                                                                                 tGGCCACGGGTCAGCTGGCTGCCTTTGCCGCGCAGCAGCGCGCCGGCGTGGTGCTTGATGATCTGGATGGCGCCGTAGGCCCAGCGGAAGCGTTGCTTCTTGAAGTCGATGAAGGTGTCGGGCATCAGGCCCTTGCCGTAGCTGTTGTgg                            <  1:608873/150‑1 (MQ=255)
                                                                                                                   gCCACGGGTCAGCTGGCTGCCTTTGCCGCGCAGCAGCGCGCCGGCGTGGTGCTTGATGATCTGGATGGCGCCGTAGGCCCAGCGGAAGCGTTGCTTCTTGAAGTCGATGAAGGTGTCGGGCATCAGGCCCTTGCCGTAGCTGTTGTGGGc                          <  1:192708/150‑1 (MQ=255)
                                                                                                                    ccACGGGTCAGCTGGCTGCCTTTGCCGCGCAGCAGCGCGCCGGCGTGGTGCTTGATGATCTGGATGCTCATCTAGGCCCAGCGGAAGCGTTGCTTCTTGAAGTCGATGAAGGTGTCGGGCATCAGGCCCTTGCCGTAGCTGTTGTGGGCg                         <  1:449204/150‑1 (MQ=255)
                                                                                                                        gggTCAGCTGGCTGCCTTTGCCGCGCAGCAGCGCGCCGGCGTGGTGCTTGATGATCTGGATGGCGCCGTAGGCCCAGCGGAAGCGTTGCTTCTTGAAGTCGATGAAGGTGTCGGGCATCAGGCCCTTGCCGTAGCTGTTGTGGGCGTAgg                     >  2:818283/1‑150 (MQ=255)
                                                                                                                             aGCTGGCTGCCTTTGCCGCGCAGCAGCGCGCCGGCGTGGTGCTTGATGATCTGGATGGCGCCGTAGGCCCAGCGGAAGCGTTGCTTCTTGAAGTCGATGAAGGTGTCGGGCATCAGGCCCTTGCCGTAGCTGTTGTGGGCGTAGgcggcg                <  1:207163/150‑1 (MQ=255)
                                                                                                                                tggctgCCTTTGCCGCGCAGCAGCGCGCCGGCGTGGTGCTTGATGATCTGGATGGCGCCGTAGGCCCAGCGGAAGCGTTGCTTCTTGAAGTCGATGAAGGTGTCGGGCATCAGGCCCTTGCCGTAGCTGTTGTGGGCGTAGGCGGCGGAc             >  2:834848/1‑150 (MQ=255)
                                                                                                                                 ggctgCCTTTGCCGCGCAGCAGCGCGCCGGCGTGGTGCTTGATGATCTGGATGGCGCCGTAGGCCCAGCGGAAGCGTTGCTTCTTGAAGTCGATGAAGGTGTCGGGCATCAGGCCCTTGCCGTAGCTGTTGTGGGCGTAGGCGGCGGACa            <  1:989340/150‑1 (MQ=255)
                                                                                                                                  gctgCCTTTGCCGCGCAGCAGCGCGCCGGCGTGGTGCTTGATGATCTGGATGGCGCCGTAGGCCCAGCGGAAGCGTTGCTTCTTGAAGTCGATGAAGGTGTCGGGCATCAGGCCCTTGCCGTAGCTGTTGTGGGCGTAGGCGGCGGACAg           >  1:1086627/1‑150 (MQ=255)
                                                                                                                                    tgCCTTTGCCGCGCAGCAGCGCGCCGGCGTGGTGCTTGATGATCTGGATGGCGCCGTAGGCCCAGCGGAAGCGTTGCTTCTTGAAGTCGATGAAGGTGTCGGGCATCAGGCCCTTGCCGTAGCTGTTGTGGGCGTAGGCGGCGGACAGGc         >  1:1008750/1‑150 (MQ=255)
                                                                                                                                       cTTTGCCGCGCAGCAGCGCGCCGGCGTGGTGCTTGATGATCTGGATGGCGCCGTAGGCCCAGCGGAAGCGTTGCTTCTTTAAGTCGATGAAGGTGTCGGGCATCAGGCCCTTGCCGTAGCTGTTGTGGGCGTAGGCGGCGGACAGGCCtt      >  1:208978/1‑150 (MQ=255)
                                                                                                                                        gccggcgcgCTGCTGCGCGCCGGCGTGGTGCTTGATGATCTGGATGGCGCCGTAGGCCCAGCGGAAGCGTTGCTTCTTGAAGTCGATGAAGGTGTCGGGCATCAGGCCCTTGCCGTAGCtgtt                                <  2:413843/118‑1 (MQ=255)
                                                                                                                                        gccggcgcgCTGCTGCGCGCCGGCGTGGTGCTTGATGATCTGGATGGCGCCGTAGGCCCAGCGGAAGCGTTGCTTCTTGAAGTCGATGAAGGTGTCGGGCATCAGGCCCTTGCCGTAGCtgtt                                >  1:413843/6‑123 (MQ=255)
                                                                                                                                          tGCCGCGCAGCAGCGCGCCGGCGTGGTGCTTGATGATCTGGATGGCGCCGTAGGCCCAGCGGAAGCGTTGCTTCTTGAAGTCGATGAAGGTGTCGGGCATCAGGCCCTTGCCGTAGCTGTTGTGGGCGTAGGCGGCGGACAGGCCTTtct   <  2:504127/150‑1 (MQ=255)
                                                                                                                                          cggcgcgCTGCTGCGCGCCGGCGTGGTGCTTGATGATCTGGATGGCGCCGTAGGCCCAGCGGAAGCGTTGCTTCTTGAAGTCGATGAAGGTGTCGGGCATCAGGCCCTTGCCGTAGCtgttgt                              <  2:1251505/120‑1 (MQ=255)
                                                                                                                                          cggcgcgCTGCTGCGCGCCGGCGTGGTGCTTGATGATCTGGATGGCGCCGTAGGCCCAGCGGAAGCGTTGCTTCTTGAAGTCGATGAAGGTGTCGGGCATCAGGCCCTTGCCGTAGCtgttgt                              >  1:1251505/4‑123 (MQ=255)
                                                                                                                                           gCCGCGCAGCAGCGCGCCGGCGTGGTGCTTGATGATCTGGATGGCGCCGTAGGCCCAGCGGAAGCGTTGCTTCTTGAAGTCGATGAAGGTGTCGGGCATCAGGCCCTTGCCGTAGCTGTTGTGGGCGTAGGCGGCGGACAGGCCTTtctc  <  1:291791/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                           
CCCGATGCGGCACGATGATCATCGCCGCCGACCACAACAGCGCGCCGATGGTGAAGAAGATGTTCATGCCATCGGCGATCCACGGTAGCCAGCCGGCCAGGAAGTGGTAGCGCTGGCCACGGGTCAGCTGGCTGCCTTTGCCGCGCAGCAGCGCGCCGGCGTGGTGCTTGATGATCTGGATGGCGCCGTAGGCCCAGCGGAAGCGTTGCTTCTTGAAGTCGATGAAGGTGTCGGGCATCAGGCCCTTGCCGTAGCTGTTGTGGGCGTAGGCGGCGGACAGGCCTTTCTC  >  NC_002947/1731781‑1732069

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