Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,173,926 G→C 5.7% R78P (CGG→CCG)  PP_2784 → short‑chain dehydrogenase/reductase family oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,173,9260GC 5.7% 441.6 / 11.6 133R78P (CGG→CCG) PP_2784short‑chain dehydrogenase/reductase family oxidoreductase
Reads supporting (aligned to +/- strand):  ref base G (77/48);  new base C (5/3);  total (82/51)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 133 total aligned reads displayed.

GCACCTTTTGATCTTGCTTGATCGGGATCTGGCAACCTTGCAGGGGTGGGTCACCGAGGGCGAATTTGCCGCTCGTATCGAAACCCATCAAGCGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGCGACGAGAACGAACCCGAAGACCTCGACAACTGGCACAAGGTCATTTCGATCAACCTCAACGGCACCTTTTACGTCACCTCGTTGTGCCTGCCACTGATGGCGGAT  >  NC_002947/3173780‑3174074
                                                                                                                                                  |                                                                                                                                                    
gCACCTTTTGATCTTGCTTGATCGGGATCTGGCAACCTTGCAGGGGTGGGTCACCGAGGGCGAATTTGCCGCTCGTATCGAAACCCATCAAGCGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGt                                                                                                                                                   <  1:3480080/150‑1 (MQ=255)
  aCCTTTTGATCTTGCTTGATCGGGATCTGGCAACCTTGCAGGGGTGGGTCACCGAGGGCGAATTTGCCGCTCGTATCGAAACCCATCAAGCGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCg                                                                                                                                                 >  2:2738211/1‑150 (MQ=255)
   ccTTTTGATCTTGCTTGATCGGGATCTGGCAACCTTGCAGGGGTGGGTCACCGAGGGCGAATTTGCCGCTCGTATCGAAACCCATCAAGCGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCgg                                                                                                                                                >  1:1213984/1‑150 (MQ=255)
    cTTTTGATCTTGCTTGATCGGGATCTGGCAACCTTGCAGGGGTGGGTCACCGAGGGCGAATTTGCCGCTCGTATCGAAACCCATCAAGCGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGc                                                                                                                                               <  1:983328/150‑1 (MQ=255)
    cTTTTGATCTTGCTTGATCGGGATCTGGCAACCTTGCAGGGGTGGGTCACCGAGGGCGAATTTGCCGCTCGTATCGAAACCCATCAAGCGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGc                                                                                                                                               >  1:1002699/1‑150 (MQ=255)
    cTTTTGATCTTGCTTGATCGGGATCTGGCAACCTTGCAGGGGTGGGTCACCGAGGGCGAATTTGCCGCTCGTATCGAAACCCATCAAGCGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGc                                                                                                                                               >  1:3002570/1‑150 (MQ=255)
    cTTTTGATCTTGCTTGATCGGGATCTGGCAACCTTGCAGGGGTGGGTCACCGAGGGCGAATTTGCCGCTCGTATCGAAACCCATCAAGCGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGc                                                                                                                                               >  1:2703025/1‑150 (MQ=255)
    cTTTTGATCTTGCTTGATCGGGATCTGGCAACCTTGCAGGGGTGGGTCACCGAGGGCGAATTTGCCGCTCGTATCGAAACCCATCAAGCGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGc                                                                                                                                               <  1:2297757/150‑1 (MQ=255)
    cTTTTGATCTTGCTTGATCGGGATCTGGCAACCTTGCAGGGGTGGGTCACCGAGGGCGAATTTGCCGCTCGTATCGAAACCCATCAAGCGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCCGGTCGGc                                                                                                                                               <  2:3215222/150‑1 (MQ=255)
    cTTTTGATCTTGCTTGATCGGGATCTGGCAACCTTGCAGGGGTGGGTCACCGAGGGCGAATTTGCCGCTCGTATCGAAACCCATCAAGCGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCCGGTCGGc                                                                                                                                               <  2:3214900/150‑1 (MQ=255)
     ttttGATCTTGCTTGATCGGGATCTGGCAACCTTGCAGGGGTGGGTCACCGGGGGCGAATTTGCCGCTCGTATCGAAACCCATCAAGCGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCt                                                                                                                                              >  2:467891/1‑150 (MQ=255)
     ttttGATCTTGCTTGATCGGGATCTGGCAACCTTGCAGGGGTGGGTCACCGAGGGCGAATTTGCCGCTCGTATCGAAACCCATCAAGCGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCt                                                                                                                                              >  1:5079053/1‑150 (MQ=255)
     ttttGATCTTGCTTGATCGGGATCTGGCAACCTTGCAGGGGTGGGTCACCGAGGGCGAATTTGCCGCTCGTATCGAAACCCATCAAGCGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCt                                                                                                                                              >  1:1560241/1‑150 (MQ=255)
     ttttGATCTTGCTTGATCGGGATCTGGCAACCTTGCAGGGGTGGGTCACCGAGGGCGAATTTGCCGCTCGTATCGAAACCCATCAAGCGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCt                                                                                                                                              >  1:1560140/1‑150 (MQ=255)
       ttGATCTTGCTTGATCGGGATCTGGCAACCTTGCAGGGGTGGGTCACCGAGGGCGAATTTGCCGCTCGTATCGAAACCCATCAAGCGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCttt                                                                                                                                            >  2:4160588/1‑150 (MQ=255)
        tGATCTTGCTTGATCGGGATCTGGCAACCTTGCAGGGGTGGGTCACCGAGGGCGAATTTGCCGCTCGTATCGAAACCCATCAAGCGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTg                                                                                                                                           <  2:1368724/150‑1 (MQ=255)
        tGATCTTGCTTGATCGGGATCTGGCAACCTTGCAGGGGTGGGTCACCGAGGGCGAATTTGCCGCTCGTATCGAAACCCATCAAGCGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTg                                                                                                                                           <  2:918579/150‑1 (MQ=255)
        tGATCTTGCTTGATCGGGATCTGGCAACCTTGCAGGGGTGGGTCACCGAGGGCGAATTTGCCGCTCGTATCGAAACCCATCAAGCGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTATg                                                                                                                                           >  1:2284654/1‑150 (MQ=255)
         gATCTTGCTTGATCGGGATCTGGCAACCTTGCAGGGGTGGGTCACCGAGGGCGAATTTGCCGCTCGTATCGAAACCCATCAAGCGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGt                                                                                                                                          <  1:2753037/150‑1 (MQ=255)
          aTCTTGCTTGATCGGGATCTGGCAACCTTGCAGGGGTGGGTCACCGAGGGCGAATTTGCCGCTCGTATCGAAACCCATCAAGCGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTc                                                                                                                                         <  1:371552/150‑1 (MQ=255)
             ttgcttgATCGGGATCTGGCAACCTTGCAGGGGTGGGTCACCGAGGGCGAATTTGCCGCTCGTATCGAAACCCATCAAGCGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGAt                                                                                                                                      >  2:2419356/1‑150 (MQ=255)
               gcttgATCGGGATCTGGCAACCTTGCAGGGGTGGGTCACCGAGGGCGAATTTGCCGCTCGTATCGAAACCCATCGAGCGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGAtgt                                                                                                                                    >  2:3535550/1‑150 (MQ=255)
                 ttgATCGGGATCTGGCAACCTTGCAGGGGTGGGTCACCGAGGGCGAATTTGCCGCTCGTATCGAAACCCATCAAGCGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGc                                                                                                                                  <  1:4694488/150‑1 (MQ=255)
                  tgATCGGGATCTGGCAACCTTGCAGGGGTGGGTCACCGAGGGCGAATTTGCCGCTCGTATCGAAACCCATCAAGCGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCt                                                                                                                                 >  1:4910675/1‑150 (MQ=255)
                    aTCGGGATCTGGCAACCTTGCAGGGGTGGGTCACCGAGGGCGAATTTGCCGCTCGTATCGAAACCCATCAAGCGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGATCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTgg                                                                                                                               >  2:4924905/1‑150 (MQ=255)
                      cGGGATCTGGCAACCTTGCAGGGGTGGGTCACCGAGGGCGAATTTGCCGCTCGTATCGAAACCCATCAAGCGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTc                                                                                                                             <  2:3498943/150‑1 (MQ=255)
                         gATCTGGCAACCTTGCAGGGGTGGGTCACCGAGGGCGAATTTGCCGCTCGTATCGAAACCCATCAAGCGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAAc                                                                                                                          >  1:1693807/1‑150 (MQ=255)
                             tGGCAACCTTGCAGGGGTGGGTCACCGAGGGCGAATTTGCCGCTCGTATCGAAACCCATCAAGCGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAgcg                                                                                                                      >  2:1784732/1‑150 (MQ=255)
                                 aaCCTTGCAGGGGTGGGTCACCGAGGGCGAATTTGCCGCTCGTATCGAAACCCATCAAGCGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCgg                                                                                                                  <  2:1454653/150‑1 (MQ=255)
                                  aCCTTGCAGGGGTGGGTCACCGAGGGCGAATTTGCCGCTCGTATCGAAACCCATCAAGCGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCggg                                                                                                                 <  2:2602444/150‑1 (MQ=255)
                                  aCCTTGCAGGGGTGGGTCACCGAGGGCGAATTTGCCGCTCGTATCGAAACCCATCAAGCGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCggg                                                                                                                 <  1:3943340/150‑1 (MQ=255)
                                       gCAGGGGTGGGTCACCGAGGGCGAATTTGCCGCTCGTATCGAAACCCATCAAGCGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTg                                                                                                            >  2:4071511/1‑150 (MQ=255)
                                                 gtCACCGAGGGCGAATTTGCCGCTCGTATCGAAACCCATCAAGCGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGCGACGAgaac                                                                                                  >  2:1380696/1‑150 (MQ=255)
                                                     ccGAGGGCGAATTTGCCGCTCGTATCGAAACCCATCAAGCGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGCGACGAgaacgaac                                                                                              <  2:2175638/150‑1 (MQ=255)
                                                     ccGAGGGCGAATTTGCCGCTCGTATCGAAACCCATCAAGCGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGCGACGAgaacgaac                                                                                              >  1:1878835/1‑150 (MQ=255)
                                                      cGAGGGCGAATTTGCCGCTCGTATCGAAACCCATCAAGCGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGCGACGAGAACGAAcc                                                                                             >  1:832919/1‑150 (MQ=255)
                                                          ggCGAATTTGCCGCTCGTATCGAAACCCATCAAGCGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGCGACGAGAACGAACCCGaa                                                                                         <  1:1489463/150‑1 (MQ=255)
                                                             gAATTTGCCGCTCGTATCGAAACCCATCAAGCGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGCGACGAGAACGAACCCGAAGAc                                                                                      <  2:2225065/150‑1 (MQ=255)
                                                             gAATTTGCCGCTCGTATCGAAACCCATCAAGCGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGCGACGAGAACGAACCCGAAGAc                                                                                      >  1:376868/1‑150 (MQ=255)
                                                                tttGCCGCTCGTATCGAAACCCATCAAGCGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCt                                                                                                                                 >  1:390492/1‑104 (MQ=255)
                                                                tttGCCGCTCGTATCGAAACCCATCAAGCGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCt                                                                                                                                 <  2:390492/104‑1 (MQ=255)
                                                                tttGCCGCTCGTATCGAAACCCATCAAGCGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGCGACGAGAACGAACCCGAAGACCTc                                                                                   <  2:4313229/150‑1 (MQ=255)
                                                                        tCGTATCGAAACCCATCAAGCGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGCGACGAGAACGAACCCGAAGACCTCGACAACTg                                                                           >  1:4949105/1‑150 (MQ=255)
                                                                          gTATCGAAACCCATCAAGCGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGCGACGAGAACGAACCCGAAGACCTCGACAACTGGc                                                                         <  2:1360685/150‑1 (MQ=255)
                                                                           tATCGAAACCCATCAAGCGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGCGACGAGAACGAACCCGAAGACCTCGACAACTGGca                                                                        >  1:3075915/1‑150 (MQ=255)
                                                                              cGAAACCCATCAAGCGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGCGACGAGAACGAACCCGAAGACCTCGACAACTGGCACaa                                                                     <  2:1663617/150‑1 (MQ=255)
                                                                               gAAACCCATCAAGCGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGCGACGAGAACGAACCCGAAGACCTCGACAACTGGCACAAg                                                                    >  1:4719921/1‑150 (MQ=255)
                                                                                  aCCAATCAAGCGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGCGACGAGAACGAACCCGAAGACCTCGACAACTGGCACAAGGTc                                                                 >  2:3747618/1‑150 (MQ=255)
                                                                                      aTCAAGCGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGCGACGAGAACGAACCCGAAGACCTCGACACCTGGCACAAGGTCAttt                                                             >  2:3851604/1‑150 (MQ=255)
                                                                                      aTCAAGCGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGCGACGAGAACGAACCCGAAGACCTCGACAACTGGCACAAGGTCAttt                                                             >  2:4148422/1‑150 (MQ=255)
                                                                                      aTCAAGCGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGCGACGAGAACGAACCCGAAGACCTCGACAACTGGCACAAGGTCAttt                                                             >  1:4512624/1‑150 (MQ=255)
                                                                                      aTCAAGCGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGCGACGAGAACGAACCCGAAGACCTCGACAACTGGCACAAGGTCAttt                                                             >  1:3920095/1‑150 (MQ=255)
                                                                                          aGCGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGCGACGAGAACGAACCCGAAGACCTCGACAACTGGCACAAGGTCATTTCGAt                                                         >  1:3906217/1‑150 (MQ=255)
                                                                                           gCGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGCGACGAGAACGAACCCGAAGACCTCGACAACTGGCACAAGGTCATTTCGATc                                                        >  2:5026138/1‑150 (MQ=255)
                                                                                            cGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGCGACGAGAACGAACCCGAAGACCTCGACAACTGGCACAAGGTCATTTCGATCa                                                       <  1:1597143/150‑1 (MQ=255)
                                                                                            cGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGCGACGAGAACGAACCCGAAGACCTCGACAACTGGCACAAGGTCATTTCGATCa                                                       <  1:1597204/150‑1 (MQ=255)
                                                                                             gAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGCGACGAGAACGAACCCGAAGACCTCGACAACTGGCACAAGGTCATTTCGATCaa                                                      >  1:2004841/1‑150 (MQ=255)
                                                                                             gAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGCGACGAGAACGAACCCGAAGACCTCGACAACTGGCACAAGGTCATTTCGATCaa                                                      >  1:2004540/1‑150 (MQ=255)
                                                                                              aaCATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGCGACGAGAACGAACCCGAAGACCTCGACAACTGGCACAAGGTCATTTCGATCAAc                                                     <  1:4820614/150‑1 (MQ=255)
                                                                                              aaCATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGCGACGAGAACGAACCCGAAGACCTCGACAACTGGCACAAGGTCATTTCGATCAAc                                                     <  1:467891/150‑1 (MQ=255)
                                                                                               aCATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGc                                                                                                                                  >  2:541132/1‑72 (MQ=255)
                                                                                                  tCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGCGACGAGAACGAACCCGAAGACCTCGACAACTGGCACAAGGTCATTTCGATCAACCTCa                                                 >  1:5257/1‑150 (MQ=255)
                                                                                                  tCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGCGACGAGAACGAACCCGAAGACCTCGACAACTGGCACAAGGTCATTTCGATCAACCTCa                                                 >  1:979710/1‑150 (MQ=255)
                                                                                                   cGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGCGACGAGAACGAACCCGAAGACCTCGACAACTGGCACAAGGTCATTTCGATCAACCTCaa                                                <  1:1435484/150‑1 (MQ=255)
                                                                                                    gCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGCGACGAGAACGAACCCGAAGACCTCGACAACTGGCACAAGGTCATTTCGATCAACCTCAAc                                               >  1:1048632/1‑150 (MQ=255)
                                                                                                    gCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGCGACGAGAACGAACCCGAAGACCTCGACAACTGGCACAAGGTCATGTCGATCAACCTCAAc                                               >  1:3902180/1‑150 (MQ=255)
                                                                                                       gATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGCGACGAGAACGAACCCGAAGACCTCGACAACTGGCACAAGGTCATTTCGAt                                                         <  2:4083889/137‑1 (MQ=255)
                                                                                                       gATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGCGACGAGAACGAACCCGAAGACCTCGACAACTGGCACAAGGTCATTTCGAt                                                         >  1:4083889/1‑137 (MQ=255)
                                                                                                       gATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGCGACGAGAACGAACCCGAAGACCTCGACAACTGGCACAAGGTCATTTCGATCAACCTCAACGGc                                            >  1:4532894/1‑150 (MQ=255)
                                                                                                           tcGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGCGACGAGAACGAACCCGAAGACCTCGACAACTGGCACAAGGTCATTTCGATCAACCTCAACGGCACCt                                        <  1:548183/150‑1 (MQ=255)
                                                                                                            cGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGCGACGAGAACGAACCCGAAGACCTCGACAACTGGCACAAGGTCATTTCGATCAACCTCAACGGCACCtt                                       <  1:1369613/150‑1 (MQ=255)
                                                                                                            cGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGCGACGAGAACGAACCCGAAGACCTCGACAACTGGCACAAGGTCATTTCGATCAACCTCAACGGCACCtt                                       <  1:1187283/150‑1 (MQ=255)
                                                                                                             gCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGCGACGAGAACGAACCCGAAGACCTCGACAACTGGCACAAGGTCATTTCGATCAACCTCAACGGCACCttt                                      <  2:4830868/150‑1 (MQ=255)
                                                                                                             gCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGCGACGAGAACGAACCCGAAGACCTCGACAACTGGCACAAGGTCATTTCGATCAACCTCAACGGCACCttt                                      >  2:4169565/1‑150 (MQ=255)
                                                                                                             gCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGCGACGAGAACGAACCCGAAGACCTCGACAACTGGCACAAGGTCATTTCGATCAACCTCAACGGCACCttt                                      >  2:1588009/1‑150 (MQ=255)
                                                                                                             gCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGCGACGAGAACGAACCCGAAGACCTCGACAACTGGCACAAGGTCATTTCGATCAACCTCAACGGCACCttt                                      <  1:2537039/150‑1 (MQ=255)
                                                                                                                  ccTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGCGACGAGAACGAACCCGAAGACCTCGACAACTGGCACAAGGTCATTTCGATCAACCTCAACGGCACCTTTTACGt                                 >  1:667371/1‑150 (MQ=255)
                                                                                                                       aaCTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGCGACGAGAACGAACCCGAAGACCTCGACAACTGGCACAAGGTCATTTCGATCAACCTCAACGGCACCTTTTACGTCACCt                            >  2:482994/1‑150 (MQ=255)
                                                                                                                       aaCTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGCGACGAGAACGAACCCGAAGACCTCGACAACTGCCACAAGGTCATTTCGATCAACCTCAACGGCACCTTTTACGTCACCt                            >  1:4958934/1‑150 (MQ=255)
                                                                                                                         cTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGCGACGAGAACGAACCCGAAGACCTCGACAACTGGCACAAGGTCATTTCGATCAACCTCAACGGCACCTTTTACGTCACCTCg                          >  1:1865866/1‑150 (MQ=255)
                                                                                                                          ttCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGCGACgaga                                                                                                    >  2:2232742/1‑75 (MQ=255)
                                                                                                                          ttCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGCGACGAGAACGAACCCGAAGACCTCGACAACTGGCACAAGGTCATTTCGATCAACCTCAACGGCACCTTTTACGTCACCTCGt                         >  1:2862450/1‑150 (MQ=255)
                                                                                                                           tCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGCGACGAGAACGAACCCGAAGACCTCGACAACTGGCACAAGGTCATTTCGATCAACCTCAACGGCACCTTTTACGTCACCTCGtt                        <  1:3143954/150‑1 (MQ=255)
                                                                                                                            cTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGCGACGAGAACGAACCCGAAGACCTCGACAACTGGCACAAGGTCATTTCGATCAACCTCAACGGCACCTTTTACGTCACCTCGTtg                       <  1:4875673/150‑1 (MQ=255)
                                                                                                                             tGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGCGACGAGAACGAACCCGAAGACCTCGACAACTGGCACAAGGTCATTTCGATCAACCTCAACGGCACCTTTTACGTCACCTCGTtgt                      >  1:1325403/1‑150 (MQ=255)
                                                                                                                             tGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGCGACGAGAACGAACCCGAAGACCTCGACAACTGGCACAAGGTCATTTCGATCAACCTCAACGGCACCTTTTACGTCACCTCGTtgt                      <  2:1238320/150‑1 (MQ=255)
                                                                                                                                 cAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGCGACGAGAACGAACCCGAAGACCTCGACAACTGGCACAAGGTCATTTCGATCAACCTCAACGGCACCTTTTACGTCACCTCGTTGtgcct                  >  1:1837229/1‑150 (MQ=255)
                                                                                                                                 cAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGCGACGAGAACGAACCCGAAGACCTCGACAACTGGCACAAGGTCATTTCGATCAACCTCAACGGCACCTTTTACGTCACCTCGTTGtgcct                  >  1:1837066/1‑150 (MQ=255)
                                                                                                                                  aaGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGCGACGAGAACGAACCCGAAGACCTCGACAACTGGCACAAGGTCATTTCGATCAACCTCAACGGCACCTTTTACGTCACCTCGTTGtgcctg                 >  1:432012/1‑150 (MQ=255)
                                                                                                                                  aaGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGCGACGAGAACGAACCCGAAGACCTCGACAACTGGCACAAGGTCATTTCGATCAACCTCAACGGCACCTTTTACGACACCTCGTTGtgcctg                 >  2:3257853/1‑150 (MQ=255)
                                                                                                                                    gggCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGCGACGAGAACGAACCCGAAGACCTCGACAACTGGCACAAGGTCATTTCGATCAACCTCAACGGCACCTTTTACGTCACCTCGTTGtgcctgcc               >  1:980287/1‑150 (MQ=255)
                                                                                                                                      gCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGCGACGAGAACGAACCCGAAGACCTCGACAACTGGCACAAGGTCATTTCGATCAACCTCAACGGCACCTTTTACGTCACCTCGTTGTGCCTGCCAc             <  2:92003/150‑1 (MQ=255)
                                                                                                                                      acaaagCCGACCCGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGCGACGAGAACGAACCCGAAGACCTCGACAACTGGCACAAGGTCATTTCGATCAACCTCAACGGCACCTTTTACGTCACCTCGTTGTGCCTGCCAc             >  2:4445804/6‑150 (MQ=255)
                                                                                                                                      acaaagCCGACCCGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGCGACGAGAACGAACCCGAAGACCTCGACAACTGGCACAAGGTCATTTCGATCAACCTCAACGGCACCTTTTACGTCACCTCGTTGTGCCTGCCAc             >  1:4934287/6‑150 (MQ=255)
                                                                                                                                       ccTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGCGACGAGAACGAACCCGAAGACCTCGACAACTGGCACAAGGTCATTTCGATCAACCTCAACGGCACCTTTTACGTCACCTCGTTGTGCCTGCCACt            <  1:513337/150‑1 (MQ=255)
                                                                                                                                        aaagCCGACCCGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGCGACGAGAACGAACCCGAAGACCTCGACAACTGGCACAAGGTCATTTCGATCAACCTCAACGGCACCTTTTACGTCACCTCGTTGTGCCTGCCACTg           >  2:4346126/4‑150 (MQ=255)
                                                                                                                                         tGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGCGACGAGAACGAACCCGAAGACCTCGACAACTGGCACAAGGTCATTTCGATCAACCTCAACGGCACCTTTTACGTCACCTCGTTGTGCCTGCCACTGa          <  2:3775148/150‑1 (MQ=255)
                                                                                                                                         tGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGCGACGAGAACGAACCCGAAGACCTCGACAACTGGCACAAGGTCATTTCGATCAACCTCAACGGCACCTTTTACGTCACCTCGTTGTGCCTGCCACTGa          >  2:1368436/1‑150 (MQ=255)
                                                                                                                                             cGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGCGACGAGAACGAACCCGAAGACCTCGACAACTGGCACAAGGTCATTTCGATCAACCTCAACGGCACCTTTTACGTCACCTCGTTGTGCCTGCCACTGATGGc      >  2:5027605/1‑150 (MQ=255)
                                                                                                                                                 cGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGCGACGAGAACGAACCCGAAGACCTCGACAACTGGCACAAGGTCATTTCGATCAACCTCAACGGCACCTTTTACGTCACCTCGTTGTGCCTGCCACTGATGGCGGAt  >  2:3600920/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
GCACCTTTTGATCTTGCTTGATCGGGATCTGGCAACCTTGCAGGGGTGGGTCACCGAGGGCGAATTTGCCGCTCGTATCGAAACCCATCAAGCGAACATCGCCGATCTCGCCAGCCTGCAACTTCTGTTCAAGGGCCTGGCCGACCGGGTCGGCTTTGTCGATGTGCTGGTCAACAGCGCCGGGGTCTGCGACGAGAACGAACCCGAAGACCTCGACAACTGGCACAAGGTCATTTCGATCAACCTCAACGGCACCTTTTACGTCACCTCGTTGTGCCTGCCACTGATGGCGGAT  >  NC_002947/3173780‑3174074

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