Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 651,426 T→A 5.1% K414M (AAG→ATG)  dsbD‑I ← thiol/disulfide interchange protein 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947651,4260TA 5.1% 268.0 / 3.8 90K414M (AAG→ATG) dsbD‑Ithiol/disulfide interchange protein 1
Reads supporting (aligned to +/- strand):  ref base T (34/51);  new base A (3/2);  total (37/53)
Fisher's exact test for biased strand distribution p-value = 3.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGCCGGCAGGCTGCCGCTGGCGGCCAGTGACGCCCCCGGCAGTGGGCGCAGTGGGTCGCCCTGGCCAGTGAGGGCGCCGTACCAGCAAGCCAGCGCGTAGACCAGCAAGGCCAGACCGACCAGTTGCGCCAGGCGTTGGGCCGTCGTCTTGACGACGAATTCCAGCGCACCGAGGAACAGTGCCACACCCGCAGCGAGGAAACCGACCAGCAGCAAGGTCGCAGGGCCCGGCAGTACACGGCTGAGCAAGCCGATGGCCAGCCCCAGCAGCAGTACGCCAATGGCGTTTTTCACC  >  NC_002947/651277‑651572
                                                                                                                                                     |                                                                                                                                                  
ctgcCGGCAGGCTGCCGCTGGCGGCCAGTGACGCCCCCGGCAGTGGGCGCAGTGGGTCGCCCTGGCCAGTGAGGGCGCCGTACCAGCAAGCCAGCGCGTAGACCAGCAAGGCCAGACCGACCAGTTGCGCCAGGCGTTGGGCCGTCGTCt                                                                                                                                                    <  1:159253/150‑1 (MQ=255)
ctgcCGGCAGGCTGCCGCTGGCGGCCAGTGACGCCCCCGGCAGTGGGCGCAGTGGGTCGCCCTGGCCAGTGAGGGCGCCGCCTGATCAAGCCAGCGCGTAGACCAGCAAGGCCAGACCGACCAGTTGCGCCAGGCGTTGGGCCGTCGTCt                                                                                                                                                    >  1:118603/1‑150 (MQ=255)
 tgcCGGCAGGCTGCCGCTGGCGGCCAGTGACGCCCCCGGCAGTGGGCGCAGTGGGTCGCCCTGGCCAGTGAGGGCGCCGTACCAGCAAGCCAGCGCGTAGACCAGCAAGGCCAGACCGACCAGTTGCGCCAGGCGTTGGGCCGTCGTCtt                                                                                                                                                   >  1:1446687/1‑150 (MQ=255)
   ccGGCAGGCTGCCGCTGGCGGCCAGTGACGCCCCCGGCAGTGGGCGCAGTGGGTCGCCCTGGCCAGTGAGGGCGCCGTACCAGCAAGCCAGCGCGTAGACCAGCAAGGCCAGACCGACCAGTTGCGCCAGGCGTTGGGCCGTCGTCTTga                                                                                                                                                 <  2:1924273/150‑1 (MQ=255)
         ggCTGCCGCTGGCGGCCAGTGACGCCCCCGGCAGTGGGCGCAGTGGGTCGCCCTGGCCAGTGAGGGCGCCGTACCAGCAAGCCAGCGCGTAGACCAGCAAGGCCAGACCGACCAGTTGCGCCAGGCGTTGGGCCGTCGTCAAGacgacgg                                                                                                                                           <  2:1337906/150‑2 (MQ=255)
          gCTGCCGCTGGCGGCCAGTGACGCCCCCGGCAGTGGGCGCAGTGGGTCGCCCTGGCCAGTGAGGGCGCCGTACCAGCAAGCCAGCGCGTAGACCAGCAAGGCCAGACCGACCAGTTGCGCCAGGCGTTGGGCCGTCGTCTTGACGACGaa                                                                                                                                          >  2:1545332/1‑150 (MQ=255)
           cTGCCGCTGGCGGCCAGTGACGCCCCCGGCAGTGGGCGCAGTGGGTCGCCCTGGCCAGTGAGGGCGCCGTACCAGCAAGCCAGCGCGTAGACCAGCAAGGCCAGACCGACCAGTTGCGCCAGGCGTTGGGCCGTCGTCAAGacgacggc                                                                                                                                          >  2:1184075/1‑147 (MQ=255)
           cTGCCGCTGGCGGCCAGTGACGCCCCCGGCAGTGGGCGCAGTGGGTCGCCCTGGCCAGTGAGGGCGCCGTACCAGCAAGCCAGCGCGTAGACCAGCAAGGCCAGACCGACCAGTTGCGCCAGGCGTTGGGCCGTCGTCAAGacgacggc                                                                                                                                          <  1:1184075/149‑3 (MQ=255)
           cTGCCGCTGGCGGCCAGTGACCCCCCCGGCAGTGGGCGCAGTGGGTCGCCCTGGCCAGTGAGGGCGCCGTACCAGCAAGCCAGCGCGTAGACCAGCAAGGACAGACCGAACAGTTGCGCCAGGGGTTGGGCCGCAGTCTTGACGAAGAAt                                                                                                                                         >  2:665180/1‑150 (MQ=255)
               cGCTGGCGGCCAGTGAGGCCCCCGGTAGTGGGCGCAGTGGGTCGCCCTGGCCAGTGAGGGCGCCGTACCAGCAAGCCTGCGCGTAGACCAGCAAGGCCAGACCGACCAGTTGCGCCAGGCGTTGGGCCGTCGTCTTGACGACGAATTCCa                                                                                                                                     <  2:396219/150‑1 (MQ=255)
               cGCTGGCGGCCAGTGACGCCCCCGGCAGTGGGCGCAGTGGGTCGCCCTGGCCAGTGAGGGCGCCGTACCAGCAAGCCAGCGCGTAGACCAGCAAGGCCAGACCGACCAGTTGCGCCAGGCGTTGGGCCGTCGTCTTGACGACGAATTCCa                                                                                                                                     <  1:740950/150‑1 (MQ=255)
               cGCTGGCGGCCAGTGACGCCCCCGGCAGTGGGCGCAGTGGGTCGCCCTGGCCAGTGAGGGCGCCGTACCAGCAAGCCAGCGCGTAGACCAGCAAGGCCAGACCGACCAGTTGCGCCAGGCGTTGGGCCGTCGTCTTGACGACGAATTCCa                                                                                                                                     >  1:1636242/1‑150 (MQ=255)
                gCTGGCGGCCAGTGACGCCCCCGGCAGTGGGCGCAGTGGGTCGCCCTGGCCAGTGAGGGCGCCGTACCAGCAAGCCAGCGCGTAGACCAGCAAGGCCAGACCGACCAGTTGCGCCAGGCGTTGGGCCGTCGTCTTGACGACGAATTCCAg                                                                                                                                    <  1:1709482/150‑1 (MQ=255)
                 cTGGCGGCCAGTGACGCCCCCGGCAGTGGGCGCAGTGGGTCGCCCTGGCCAGTGAGGGCGCCGTACCAGCAAGCCAGCGCGTAGACCAGAAAAGCCAGACCGACCAGTTGCGCCAGGCGTTGGGCCGTCGTCTTGACGACGAATTCCAgc                                                                                                                                   <  2:1606277/150‑1 (MQ=255)
                      ggcCAGTGACGCCCCCGGCAGTGGGCGCAGTGGGTCGCCCTGGCCAGTGAGGGCGCCGTACCAGCAAGCCAGCGCGTAGACCAGCAAGGCCAGACCGACCAGTTGCGCCAGGCGTTGGGCCGTCGTCTTGACGACGAATTCCAGCGCAcc                                                                                                                              >  1:1499243/1‑150 (MQ=255)
                       gcCAGTGACGCCCCCGGCAGTGGGCGCAGTGGGTCGCCCTGGCCAGTGAGGGCGCCGTACCAGCAAGCCAGCGCGTAGACCAGCAAGGCCAGACCGACCAGTTGCGCCAGGCGTTGGGCCGTCGTCTTGACGACGAATTCCAGCGCCCCg                                                                                                                             >  2:1112/1‑150 (MQ=255)
                                gCCCCCGGCAGTGGGCGCAGTGGGTCGCCCTGGCCAGTGAGGGCGCCGTACCAGCAAGCCAGCGCGTAGACCAGCAAGGCCAGACCGACCAGTTGCGCCAGGCGTTGGGCCGTCGTCTTGACGACGAATTCCAGCGCACCGAGGAACAGt                                                                                                                    >  2:1331927/1‑150 (MQ=255)
                                  ccccGGCAGTGGGCGCAGTGGGTGGCCCTGTCCAGTGAGGGCGCCGTACCAGCAAGCCAGCGCGTAGACCAGCAAGGCCAGACCGACCAGTTGCGCCAGGCGTTGGGCCGTCGTCTTGACGACGGATTCCAGCGCACCGAGGAACGGTGc                                                                                                                  >  2:948572/1‑150 (MQ=255)
                                   cccGGCAGTGGGCGCAGTGGGTCGCCCTGGCCAGTGAGGGCGCCGTACCAGCAAGCCAGCGCGTAGACCAGCAAGGCCAGACCGACCAGTTGCGCCAGGCGTTGGGCCGTCGTCTTGACGACGAATTCCAGCGCACCGAGGAACAGTGcc                                                                                                                 <  1:1763142/150‑1 (MQ=255)
                                    cgggCAGTGGGCGCAGTGGGTCGCCCTGGCCAGTGAGGGCGCCGTACCAGCAAGCCAGCGCGTAGACCAGCAAGGCCAGACCGACCAGTTGCGCCAGGCGTTGGGCCGTCGTCTTGACGACGAATTCCAGCGCACCGAGGAACAGTGCca                                                                                                                <  2:636368/148‑1 (MQ=255)
                                      ggCAGTGGGCGCAGTGGGTCGCCCTGGCCAGGGAGGGCGCCGTCGCTGCCGATCTGGGTGTAGACCAGCAAGGCCAGACCGACCAGTTGCGCCAGGCGTTGGGCCGTCGTCTTGACGACGAATTCCAGCGCACCGAGGAACAGTGCcaca                                                                                                              <  2:529016/150‑1 (MQ=255)
                                             ggCGCAGTGGGTCGCCCTGGCCAGTGAGGGCGCCGTACCAGCAAGCCAGCGCGTAGACCAGCAAGGCCAGACCGACCAGTTGCGCCAGGCGTTGGGCCGTCGTCTTGACGACGAATTCCAGCGCACCGAGGAACAGTGCCACACCCGCAg                                                                                                       <  2:181225/150‑1 (MQ=255)
                                              gcgcAGTGGGTCGCCCTGGCCAGTGAGGGCGCCGTACCAGCAAGCCAGCGCGTAGACCAGCAAGGCCAGACCGACCAGTTGCGCCAGGCGTTGGGCCGTCGTCTTGACGACGAATTCCAGCGCACCGAGGAACAGTGCCACACCCGCAGc                                                                                                      <  2:94496/150‑1 (MQ=255)
                                               ggcAGTGGGTCGCCCTGGCCAGTGAGGGCGCCGTACCAGCAAGCCAGCGCGTATACCAGCAAGGCCAGACCGACCAGTTGAGCCAGGCGTTGGGCCGTCGTCTTGACGACGAATTCCAGCGCACCGAGGAACAGTGCCACACCCGCAGCg                                                                                                     >  1:1471011/2‑150 (MQ=255)
                                                 cAGTGGGTCGCCCTGGCCAGTGAGGGCGCCGTACCAGCAAGCCAGCGCGTAGACCATCAAGGCCAGACCGACCAGTTGCGCCAGCCGTTTGGCCGTCGTCTTGACGACGAATTCCAGCGCACCGAGGAACAGTGCCACACCCGCAGCGAg                                                                                                   <  1:927193/150‑1 (MQ=255)
                                                 cAGTGGGTCGCCCTGGCCAGTGAGGGCGCCGTACCAGCAAGCCAGCGCGTAGACCAGCAAGGCCAGACCGACCAGTTGCGCCAGGCGTTGGGCCGTCGTTTTGACGACGAATTCCAGCGCACCGAGGAACAGTGCCACACCCGCAGCGAg                                                                                                   <  2:265945/150‑1 (MQ=255)
                                                 cAGTGGGTCGCCCTGGCCAGTGAGGGCGCCGTACCAGCAAGCCAGCGCGTAGACCAGCAAGGCCAGACCGACCAGTTGCGCCAGGCGTTGGGCCGTCGTCTTGACGACGAATTCCAGCGCACCGAGGAACAGTGCCACACCCGCAGCGAg                                                                                                   <  1:717905/150‑1 (MQ=255)
                                                 cAGTGGGCACCCCTGGCCAGTGAGGGCGCCGTACCAGCAAGCCAGCGCGTAGACCAGCAAGGCCAGACCGACCAGTTGCGCCAGGCGTTGGGCCGTCGTCTTGACGACGAATTCCAGCGCACCGAGGAACAGTGCCACACCCGCAGCGAg                                                                                                   <  1:1020274/150‑1 (MQ=255)
                                                      ggTCGCCCTGGCCAGTGAGGGCGCCGTACCAGCAAGCCAGCGCGTAGACCAGCAAGGCCAGACCGACCAGTTGCGCCAGGCGTTGGGCCGTCGTCTTGACGACGAATTCCAGCGCACCGAGGAACAGTGCCACACCCGCAGCGa                                                                                                    <  1:898813/144‑1 (MQ=255)
                                                      ggTCGCCCTGGCCAGTGAGGGCGCCGTACCAGCAAGCCAGCGCGTAGACCAGCAAGGCCAGACCGACCAGTTGCGCCAGGCGTTGGGCCGTCGTCTTGACGACGAATTCCAGCGCACCGAGGAACAGTGCCACACCCGCAGCGa                                                                                                    >  2:898813/1‑144 (MQ=255)
                                                        tCGCCCTGGCCAGTGAGGGCGCCGTACCAGCAAGCCAGCGCGTAGACCAGCAAGGCCAGACCGACCAGTTGCGCCAGGCGTTGGGCCGTCGTCTTGACGACGAAt                                                                                                                                         <  1:266823/105‑1 (MQ=255)
                                                        tCGCCCTGGCCAGTGAGGGCGCCGTACCAGCAAGCCAGCGCGTAGACCAGCAAGGCCAGACCGACCAGTTGCGCCAGGCGTTGGGCCGTCGTCTTGACGACGAAt                                                                                                                                         >  2:266823/1‑105 (MQ=255)
                                                          gCCCTGGCCAGTGAGGGCGCCGTACCAGCAAGCCAGCGCGTAGACCAGCAAGGCCAGACCGACCAGTTGCGCCAGGCGTTGGGCCGTCGTCTTGACGACGAATTCCAGCGCACCGAGGAACAGTGCCACACCCGCAGCGAGGAAACCGAc                                                                                          <  2:1636242/150‑1 (MQ=255)
                                                          gCCCTGGCCAGTGAGGGCGCCGTACCAGCAAGCCAGCGCGTAGACCAGCAAGGCCAGACCGACCAGTTGCGCCAGGCGTTGGGCCGTCGTCTTGACGACGAATTCCAGCGCACCGAGGAACAGTGCCACACCCGCAGCGAGGAAACCGAc                                                                                          <  2:48028/150‑1 (MQ=255)
                                                           cccTGGCCAGTGAGGGCGCCGTACCAGCAAGCCAGCGCGTAGACCAGCAAGGCCAGACCGACCAGTTGCGCCAGGCGTTGGGCCGTCGTCTTGACGACGAATTCCAGCGCACCGAGGAACAGTGCCACACCCGCAGCGAg                                                                                                   >  2:1871110/1‑140 (MQ=255)
                                                           cccTGGCCAGTGAGGGCGCCGTACCAGCAAGCCAGCGCGTAGACCAGCAAGGCCAGACCGACCAGTTGCGCCAGGCGTTGGGCCGTCGTCTTGACGACGAATTCCAGCGCACCGAGGAACAGTGCCACACCCGCAGCGAg                                                                                                   <  1:1871110/140‑1 (MQ=255)
                                                            ccTGGCCAGTGAGGGCGCCGTACCAGCAAGCCAGCGCGTAGACCAGCAAGGCCAGACCGACCAGTTGCGCCAGGCGTTGGGCCGTCGTCTTga                                                                                                                                                 >  1:718627/1‑93 (MQ=255)
                                                            ccTGGCCAGTGAGGGCGCCGTACCAGCAAGCCAGCGCGTAGACCAGCAAGGCCAGACCGACCAGTTGCGCCAGGCGTTGGGCCGTCGTCTTga                                                                                                                                                 <  2:718627/93‑1 (MQ=255)
                                                               ggCCAGTGAGGGCGCCGTACCAGCAAGCCAGCGCGTAGACCAGCAAGGCCAGACCGACCAGTTGCGCCAGGCGTTGGGCCGTCGTCTTGACGACGAATTCCAGCGCACCGAGGAACAGTGCCACACCCGCAGCGAGGAAACCGACcagca                                                                                     <  2:1499243/150‑1 (MQ=255)
                                                                    gcggcggCGCCGCGCCAGCAAGCCAGCGCGTAGACCAGCAAGGCCAGACCGACCAGTTGCGCCAGGCGTTGGGCCGTCGTCTTGACGACGAATTCCAGCGCACCGAGGAACAGTGCCACACCCGCATCGAGGAAACCGACCAGCAGCAAg                                                                                <  2:942216/145‑1 (MQ=255)
                                                                    gTGAGGGCGCCGTACCAGCAAGCCAGCGCGTAGACCAGCAAGGCCAGACCGACCAGTTGCGCCAGGCGTTGGGCCGTCGTCTTGACGACGAATTCCAGCGCACCGAGGAACAGTGCCACAcc                                                                                                            <  1:1692371/122‑1 (MQ=255)
                                                                    gTGAGGGCGCCGTACCAGCAAGCCAGCGCGTAGACCAGCAAGGCCAGACCGACCAGTTGCGCCAGGCGTTGGGCCGTCGTCTTGACGACGAATTCCAGCGCACCGAGGAACAGTGCCACAcc                                                                                                            >  2:1692371/1‑122 (MQ=255)
                                                                    gTGAGGGCGCCGTACCAGCAAGCCAGCGCGTAGACCAGCAAGGCCAGACCGACCAGTTGCGCCAGGCGTTGGGCCGTCGTCTTGACGACGAATTCCAGCGCACCGAGGAACAGTGCCACACCCGCAGCGAGGAAACCGACCAGCAGCAAg                                                                                <  2:1752762/150‑1 (MQ=255)
                                                                    gTGAGGGCGCCGTACCAGCAAGCCAGCGCGTAGACCAGCAAGGCCAGACCGACCAGTTGCGCCAGGCGTTGGGCCGTCATCTTGACGACGAATTCCGGTGCGCGCATATACATTGCCACGCCCGCAGCGAGGAAACCGACCAGCAGCAAg                                                                                <  1:799443/150‑1 (MQ=255)
                                                                         ggCGCCGTACCAGCAAGCCAGCGCGTAGACCAGCAAGGCCAGACCGACCAGTTGCGCCAGGCGTTGGGCCGTCGTCTTGACGACGAATTCCAGCGCa                                                                                                                                >  2:1543452/1‑97 (MQ=255)
                                                                         ggCGCCGTACCAGCAAGCCAGCGCGTAGACCAGCAAGGCCAGACCGACCAGTTGCGCCAGGCGTTGGGCCGTCGTCTTGACGACGAATTCCAGCGCa                                                                                                                                <  1:1543452/97‑1 (MQ=255)
                                                                         ggCGCCGTACCAGCAAGCCAGCGCGTAGACCAGCAAGGCCAGACCGACCAGTTGCGCCAGGCGTTGGGCCGTCGTCTTGACGACGAATTCCAGCGCACCGAGGAACAGTGCCACACCCGCAGCGAGGAAACCGACCAGCAGCAAGGTCGc                                                                           >  1:360726/1‑150 (MQ=255)
                                                                            gcCGTACCAGCAAGCCAGCGCGTAGACCAGCAAGGCCAGACCGACCAGTTGCGCCAGGCGTTGGGCCGTCGTCTTGACGACGAATTCCAGCGCACCGAGGAACAGTGCCACACCCGCAGCGAGGAAACCGACCAGCAGCAAGGTCGCAgg                                                                        >  2:1089661/1‑150 (MQ=255)
                                                                             ccGTACCAGCAAGCCAGCGCGTAGACCAGCAAGGCCAGACCGACCAGTTGCGCCAGGCGTTGGGCCGTCGTCTTGACGACGAATTCCAGCGCACCGAGGAACAGTGCCACACCCGCAGCGAGGAAACCGACCAGCAGCAAGGTCGCAggg                                                                       <  1:1545332/150‑1 (MQ=255)
                                                                                 aCCAGCAAGCCAGCGCGTAGACCAGCAAGGCCAGACCGACCAGTTGCGCCAGGCGTTGGGCCGTCGTCTTGACGACGAATTCCAGCGCACCGAGGAACAGTGCCACACCCGCAGCGAGGAAACCGACCAGCAGCAAGGTCGCAGGGCCCg                                                                   <  2:1446687/150‑1 (MQ=255)
                                                                                  ccAGCAAGCCAGCGCGTAGACCAGCAAGGCCAGACCGACCAGTTGCGCCAGGCGTTGGGCCGTCGTCTTGACGACGAATTCCAGCGCACCGAGGAACAGTGCCACACCCGCAGCGAGGAAACCGACCAGCAGCAAGGTCGCAGGGCCCgg                                                                  <  2:607866/150‑1 (MQ=255)
                                                                                      cAAGCCAGCGCGTAGACCAGCAAGGCCAGACCGACCAGTTGCGCCAGGCGTTGGGCCGTCGTCTTGACGACGAATTCCAGCGCACCGAGGAACAGTGCCACACCCGCAGCGAGGAAACCGACCAGCAGCAAGGTCGCAGGGCCCGGCAGt                                                              >  1:1233684/1‑150 (MQ=255)
                                                                                       aaGCCAGCGCGTAGACCAGCAAGGCCAGACCGACCAGTTGCGCCAGGCGTTGGGCCGTCGTCTTGACGACGAATTCCAGCGCACCGAGGAACAGTGCCACACCCGCAGCGAGGAAACCGACCAGCAGCAAGGTCGCAGGGCCCGGCAGTa                                                             <  2:656221/150‑1 (MQ=255)
                                                                                          ccAGCGCGTAGACCAGCAAGGCCAGACCGACCAGTTGCGCCAGGCGTTGGGCCGTCGTCTTGACGACGAATTCCAGCGCACCGAGGAACAGTGCCACACCCGCAGCGAGGAAACCGACCAGCAGCAAGGTCGCAGGGCCCGGCAGTacac                                                          >  2:1172312/1‑150 (MQ=255)
                                                                                          ccAGCCCGTAGAGCACAAAGGCCAGACCGACCAGTTGCGCCACGAGATGGGCCGACGTCATGACGACGAATTCCAGAACACCGAGACCCAGTGCCACACCCCCAGCGAGGAAACCGACCAGCAGCAAGGTAGCAGGGCCCGGCAGTacac                                                          >  2:786141/1‑150 (MQ=255)
                                                                                            aGCGCGTAGACCAGCAAGGCCAGACCGACCAGTTGCGCCAGGCGTTGGGCCGTCGTCTTGACGACGAATTCCAGCGCACCGAGGAACAGTGCCACACCCGCAGCGCGGAAACCGACCAGCAGCAAGGTCGCCGGGCCCGGCCGTACACgg                                                        >  2:466984/1‑150 (MQ=255)
                                                                                             gcgcGTAGACCAGCAAGGCCAGACCGACCAGTTGCGCCAGGCGTTGGGCCGTCGTCTTGATGACGAATTCCAGCGCACCGAGGAACAGTGCCACACCCGCAGCGAGGAAACCGACCAGCAGCAAGGTCGCAGGGCCCGGCAGTACACGGc                                                       >  1:1908555/1‑150 (MQ=255)
                                                                                                ggTAGACCAGCAAGGCCAGACCGACCAGTCGCGCCAGGCGTTGGGCCGTCGTCTTGACGACGAATTCCAGCGCACCGAGGAACAGTGCCACACCCGCAGCGAGGAAACCGACCAGCAGCAAGGTCGCAGGGCCCGGCAGTACACGGCTGa                                                    <  2:1235761/149‑1 (MQ=255)
                                                                                                cgTAGACCAGCAAGGCCAGCCCGACCAGTTGCGCCAGGCGTTGGGCCCGCGTCTTGACGACGAATTCCAGCGTACCGAGGAACAGTGCCACACCCGCAGCGATGAAACCGACCAGCAGCAAGGTCGCAGGGCCCGGCAGTACACGGCTGa                                                    <  1:665180/150‑1 (MQ=255)
                                                                                                 gTAGACCAGCAAGGCCAGACCGACCAGTTGCGCGAGGCGTTGGGCCGTCGTCTTGACGACGAATTCCAGCGCACCGAGGAACAGTGCCACCGCCGCAGCGAGGAAACCGACCAGCAGCAAGGTCGCAGGGCCCGGCAGTACACGGCTGAg                                                   <  2:166870/150‑1 (MQ=255)
                                                                                                 gTAGACCAGCAAGGCCAGACCGACCAGTTGCGCCAGGCGTTGGGCCGTCGTCTTGACGACGAATTCCAGCGCACCGAGGAACTGTGCCACACCCGCAGCGAGGAAACCGACCAGCAGCAAGGTCGCAGGGCCCGGCAGTACACGGCTGAg                                                   <  2:514819/150‑1 (MQ=255)
                                                                                                      ccAGCAAGGCCAGACCGGCCAGTTGCGCCAGGCGTTGGGCCGTCGTCTTGACGACGAATTCCAGCGCACCGAGGAACAGTGCCACACCCGCAGCGAGGAAACCGACCAGCAGCAAGGTCGCAGGGCCCGGCAGTACACGGCTGAGCAAGc                                              <  2:1560879/150‑1 (MQ=255)
                                                                                                      ccAGCAAGGCCAGACCGACCAGTTGCGCCAGGCGTTGGGCCGTCGTCTTGACGACGAATTCCAGCGCACCGAGGAACAGTGCCACACCCGCAGCGAGGAAACCGACCAGCAGCAAGGTCGCAGGGCCCGGCAGTACACGGCTGAGCAAGc                                              <  1:1089661/150‑1 (MQ=255)
                                                                                                       cAGCAAGGCCAGACCGACCAGTTGCGCCAGGCGTTGGGCCGTCGTCTTGACGACGAATTCCAGCGCACCGAGGAACAGTGCCACACCCGCAGCGAGGAAACCGACCAGCAGCAAGGTCGCAGGGCCCGGCAGTACACGGCTGAGCAAGcc                                             >  1:842262/1‑150 (MQ=255)
                                                                                                         gCAAGGCCAGACCGACCAGTTGCGCCAGGCGTTGGGCCGTCGTCTTGACGACGAATTCCAGCGCACCGAGGAACAGTGcc                                                                                                                 >  1:38337/1‑80 (MQ=255)
                                                                                                         gCAAGGCCAGACCGACCAGTTGCGCCAGGCGTTGGGCCGTCGTCTTGACGACGAATTCCAGCGCACCGAGGAACAGTGcc                                                                                                                 <  2:38337/80‑1 (MQ=255)
                                                                                                         gCAAGGCCAGACCGACCAGTTGCGCCAGGCGTTGGGCCGTCGTCTTGACGACGAATTCCAGCGCACCGAGGAACAGTGCCACACCCGCAGCGAGGAAACCGACCAGCAGCAAGGTCGCAGGGCCCGGCAGTACACGGCTGAGCAAGCCGa                                           >  1:1424366/1‑150 (MQ=255)
                                                                                                          cAAGGCCAGACCGACCAGTTGCGCCAGGCGTTGGGCCGTCGTCTTGACGACGAATTCCAGCGCACCGAGGAACAGTGCCACACCCGCAGCGAGGAAACCGACCAGCAGCAAGGTCGCAGGGCCCGGCAGTACACGGCTGAGCAAGCCGAt                                          >  1:352496/1‑150 (MQ=255)
                                                                                                          cAAGGCCAGACCGACCAGTTGCGCCAGGCGTTGGGCCGTCGTCTTGACGACGAATTCCAGCGCACCGAGGAACAGTGCCACACCCGCAGCGAGGAAACCGACCAGCAGCAAGGTCGCAGGGCCCGGCAGTACACGGCTGAGCAAGCCGAt                                          <  2:1429553/150‑1 (MQ=255)
                                                                                                               ccAGACCGACCAGTTGCGCCAGGCGTTGGGCCGTCGTCTTGACGACGAATTCCAGCGCACCGAGGAACAGTGCCACACCCGCAGCGAGGAAACCGACCAGCAGCAAGGTCGCAGGGCCCGGCAGTACACGGCTGAGCAAGCCGATGGCCa                                     >  1:1021234/1‑150 (MQ=255)
                                                                                                                 aGACCGACCAGTTGCGCCAGGCGTTGGGCCGTCGTCTTGACGACGAATTCCAGCGCACCGAGGAACAGTGCCACACCCGCAGCGAGGAAACCGACCAGCAGCAAGGTCGCAGGGCCCGGCAGTACACGGCTGAGCAAGCCGATGGCCAGc                                   >  1:1787046/1‑150 (MQ=255)
                                                                                                                    ccgaccAGTTGCGCCAGGCGTTGGGCCGTCGTCTTGACGACGAATTCCAGCGCACCGAGGAACAGTGCCACACCCGCAGCGACGAAACCGACCAGCAGCAAGGTCGCAGGGCCCGGCAGTACACGGCTGAGCAAGCCGATGGCCAGcccc                                <  2:101322/150‑1 (MQ=255)
                                                                                                                     cgaccAGTTGCGCCAGGCGTTGGGCCGTCGTCTTGACGACGAATTCCAGCGCACCGAGGAACAGTGCCACACCCGCAGCGAGGAAACCGACCAGCAGCAAGGTCGCAGGGCCCGGCAGTACACGGCTGAGCAAGCCGATGGCCAGCCCca                               <  1:410455/150‑1 (MQ=255)
                                                                                                                        ccAGTTGCGCCAGGCGTTGGGCCGTCGTCTTGACGACGAATTCCAGCGCACCGAGGAACAGTGCCACACCCGCAGCGAGGAAACCGACCAGCAGCAAGGTCGCAGGGCCCgg                                                                  >  1:1133435/1‑112 (MQ=255)
                                                                                                                        ccAGTTGCGCCAGGCGTTGGGCCGTCGTCTTGACGACGAATTCCAGCGCACCGAGGAACAGTGCCACACCCGCAGCGAGGAAACCGACCAGCAGCAAGGTCGCAGGGCCCgg                                                                  <  2:1133435/112‑1 (MQ=255)
                                                                                                                         cAGTTGCGCCAGGCGTTGGGCCGTCGTCTTGACGACGAATTCCAGCGCACCGAGGAACAGTGCCACACCCGCAGCGAGGAAACCGACCAGCAGCAAGGTCGCAGGGCCCGGCAGTACACGGCTGAGCAAGCCGATGGCCAGCCCcagcag                           <  2:1471011/150‑1 (MQ=255)
                                                                                                                           tttGCGCCAGGCGCTGGGCCGTCGTCTTGACGACGAATTCCAGCGCACCGAGGAACAGTGCCACACCCGCAGCGAGGAACCCGACCAGCAGCAAGGTCGCAGGGCCCGGCAGTACACGGCTGAGCAAGCCGATGGCCAGCCCcagcagca                         <  1:856762/149‑1 (MQ=255)
                                                                                                                           gTTGCGCCAGGCGTTGGGCCGTCGTCTTGACGACGAATTCCAGCGCACCGAGGAACAGTGCCACACCCGCAGCGAGGAAACCGACCAGCAGCAAGGTCGCAGGGCCCGGCAGTACACGGCTGAGCAAGCCGATGGCCAGCCCcagcagca                         >  2:1211199/1‑150 (MQ=255)
                                                                                                                            ttGCGCCAGGCGTTGGGCCGTCGTCTTGACGACGAATTCCAGCGCACCGAGGAACAGTGCCACACCCGCAGCGAGGAAACCGACCAGCAGCAAGGTCGCAGGGCCCGGCAGTACACGGCTGAGCAAGCCGATGGCCAGCCCcagcagcag                        <  2:1021234/150‑1 (MQ=255)
                                                                                                                                 ccAGGCGTTGGGCCGTCGTCTTGACGACGAATTCCAGCGCACCGAGGAACAGTGCCACACCCGCAGCGAGGAAACTGACCAGCAGCAAGGTCGCAGGGCCCGGCAGTACACGGCTGAGCAAGCCGATGGCCAGCCCCAGCAGCAGTACGc                   >  1:602123/1‑150 (MQ=255)
                                                                                                                                 ccAGGCGTTGGGCCGTCGTCTTGACGACGAATTCCAGCGCACCGAGGAACAGTGCCACACCCGCAGCGAGGAAACCGACCAGCAGCAAGGTCGCAGGGCCCGGCAGTACACGGCTGAGCAAGCCGATGGCCa                                     >  1:1211544/1‑132 (MQ=255)
                                                                                                                                 ccAGGCGTTGGGCCGTCGTCTTGACGACGAATTCCAGCGCACCGAGGAACAGTGCCACACCCGCAGCGAGGAAACCGACCAGCAGCAAGGTCGCAGGGCCCGGCAGTACACGGCTGAGCAAGCCGATGGCCa                                     <  2:1211544/132‑1 (MQ=255)
                                                                                                                                 ccAGGCGTTGGGCCGTCGTCTTGACGACGAATTCCAGCGCACCGAGGAACAGTGCCACACCCGCAGCGAGGAAACCGACCAGCAGCAAGGTCGCAGGGCCCGGCAGTACACGGCTGAGCAAGCCGATGGCCAGCCCCAGCAGCAGTACGc                   >  2:979239/1‑150 (MQ=255)
                                                                                                                                 ccAGGCGTTGGGCCGTCGTCTTGACGACGAATTCCAGCGCACCGAGGAACAGTGCCACACCCGCAGCGAGGAAACCGACCAGCAGCAAGGTCGCAGGGCCCGGCAGTACACGGCTGAGCAAGCCGATGGCCAGCCCCAGCAGCAGTACGc                   >  1:1502562/1‑150 (MQ=255)
                                                                                                                                 ccAGGCGTTGGGCCGTCGTCTTGACGACGAATTCCAGCGCACCGAGGAACAGTGCCACACCCGCAGCGAGGAAACCGACCAGCAGCAAGGTCGCAGGGCCCGGCAGTACACGGCTGAGCAAGCCGATGGCCAGCCCCAGCAGCAGTACGc                   >  2:1017357/1‑150 (MQ=255)
                                                                                                                                  cAGGCGTTGGGCCGTCGTCTTGACGACGAATTCCAGCGCACCGAGGAACAGTGCCACACCCGCAGCGAGGAAACCGACCAGCAGCAAGGTCGCAGGGCCCGGCAGTACACGGCTGAGCAAGCCGATGGCCAGCCCCAGCAGCAGTACGcc                  <  2:352496/150‑1 (MQ=255)
                                                                                                                                      cGTTGGGCCGTCGTCTTGACGACGAATTCCAGCGCACCGAGGAACAGTGCCACACCCGCAGCGAGGAAACCGACCAGCAGCAAGGTCGCAGGGCCCGGCAGTACACGGCTGAGCAAGCCGATGGCCAGCCCCAGCAGCAGTACGCCAATg              <  1:1211199/150‑1 (MQ=255)
                                                                                                                                       gTTGGGCCGTCGTCTTGACGACGAATTCCAGCGCACCGAGGAACAGTGCCACACCCGCAGCGAGGAAACCGACCAGCAGCaa                                                                                 <  2:1375405/82‑1 (MQ=255)
                                                                                                                                       gTTGGGCCGTCGTCTTGACGACGAATTCCAGCGCACCGAGGAACAGTGCCACACCCGCAGCGAGGAAACCGAACAGCATCaa                                                                                 >  1:1375405/1‑82 (MQ=255)
                                                                                                                                            ttcgtcgtCAAGACGACGAATTCCAGCGCACCGAGGAACAGTGCCACACCCGCAGCGAGGAAACCGACCAGCAGCAAGGTCGCAGGGCCCGGCAGTACACTGATGAACAAACCGATGGCCAGCCCCAGCAGCATTACGCCAATGGCGttt        >  2:1534706/3‑150 (MQ=255)
                                                                                                                                                tcgtcTTGGCGACGAATGCCAGCGCACCGGGGAACAGTGCCACACCCGCAGCGAGGAAACCGACCAGCAGCAAGGTCGCAGGGCCCGGCAGTACACGGCTGAGCAAGCCGATGGCCAGCCCCAGCAGCAGTACGCCAATGGCGTTTTTCa    <  2:172452/150‑1 (MQ=255)
                                                                                                                                                 cgtcTTGACGACGAATGCCAGCGCACCGAGGAACAGTGCCACACCCGCAGCGAGGAAACCGACCAGCAGCAAGGTCGCAGGGCCCGGCAGTACACGGCTGAGCAAGCCGATGGCCAGCCCCAGCAGCAGTACGCCAATGGCGTTTTTCAc   <  1:775125/150‑1 (MQ=255)
                                                                                                                                                  gtcTTGACGACGAATTCCAGCGCACCGAGGAACAGTGCCACACCCGCAGCGAGGAAACCCACCAGCAGCAAGGTCGCAGGGCCCGGCAGTACAGGGCTGAGCAAGCCGATGGCCAGCCCCAGCAGCAGTACGCCAATGGCGTTTTTCAcc  <  1:786141/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
CTGCCGGCAGGCTGCCGCTGGCGGCCAGTGACGCCCCCGGCAGTGGGCGCAGTGGGTCGCCCTGGCCAGTGAGGGCGCCGTACCAGCAAGCCAGCGCGTAGACCAGCAAGGCCAGACCGACCAGTTGCGCCAGGCGTTGGGCCGTCGTCTTGACGACGAATTCCAGCGCACCGAGGAACAGTGCCACACCCGCAGCGAGGAAACCGACCAGCAGCAAGGTCGCAGGGCCCGGCAGTACACGGCTGAGCAAGCCGATGGCCAGCCCCAGCAGCAGTACGCCAATGGCGTTTTTCACC  >  NC_002947/651277‑651572

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: