Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,953,446 G→A 6.5% R46R (CGC→CGT) ‡ mnmC ← methyltransferase/FAD‑dependent oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,953,4460GA 6.5% 565.1 / 19.5 184R46R (CGC→CGT) ‡mnmCmethyltransferase/FAD‑dependent oxidoreductase
Reads supporting (aligned to +/- strand):  ref base G (84/88);  new base A (3/9);  total (87/97)
Fisher's exact test for biased strand distribution p-value = 1.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.40e-01

Only 100 of 186 total aligned reads displayed.

TCGACGCTGACGAAGTGCAGGCGCGCGTCGCTGTGTGCGTGCTGGTCGAACAGCTGCCAGGCGCAGAAAAAATTCATGCCGGTGCCAAAACCGGTTTCGCCGATCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGTTGACTGCGAAATAGACGTCGTCATAGTGCCGCGAATGGGGGCGGCCCTGGTCGTCCCAGTCGATCTGGGCGTGCTGGAGGAGGGTGGGCATGGTGGGCTCGG  >  NC_002947/1953297‑1953594
                                                                                                                                                     |                                                                                                                                                    
tcGACGCTGACGAAGTGCAGGCGCGCGTCGCTGTGTGCGTGCTGGTCGAACAGCTGCCAGGCGCAGAAAAAATTCATGCCGGTGCCAAAACCGGTTTCGCCGATCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAgg                                                                                                                                                      >  2:2427280/1‑150 (MQ=255)
     gCTGACGAAGTGCAGGCGCGCGTCGCTGTGTGCGTGCTGGTCGAACAGCTGCCAGGCGCAGAAAAAATTCATGCCGGTGCCAAAACCGGTTTCGCCGATCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGt                                                                                                                                                 >  2:582968/1‑150 (MQ=255)
     gCTGACGAAGTGCAGGCGCGCGTCGCTGTGTGCGTGCTGGTCGAACAGCTGCCAGGCGCAGAAAAAATTCATGCCGGTGCCAAAACCGGTTTCGCCGATCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGt                                                                                                                                                 >  1:2408600/1‑150 (MQ=255)
       tGACGAAGTGCAGGCGCGCGTCGCTGTGTGCGTGCTGGTCGAACAGCTGCCAGGCGCAGAAAAAATTCATGCCGGTGCCAAAACCGGTTTCGCCGATCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCt                                                                                                                                               >  1:2902046/1‑150 (MQ=255)
        gACGAAGTGCAGGCGCGCGTCGCTGTGTGCGTGCTGGTCGAACAGCTGCCAGGCGCAGAAAAAATTCATGCCGGTGCCAAAACCGGTTTCGCCGATCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTg                                                                                                                                              >  1:2973289/1‑150 (MQ=255)
        gACGAAGTGCAGGCGCGCGTCGCTGTGTGCGTGCTGGTCGAACAGCTGCCAGGCGCAGAAAAAATTCATGCCGGTGCCAAAACCGGTTTCGCCGATCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTg                                                                                                                                              >  1:4756560/1‑150 (MQ=255)
         aCGAAGTGCAGGCGCGCGTCGCTGTGTGCGTGCTGGTCGAACAGCTGCCAGGCGCAGAAAAAATTCATGCCGGTGCCAAAACCGGTTTCGCCGATCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTgg                                                                                                                                             >  1:1856747/1‑150 (MQ=255)
         aCGAAGTGCAGGCGCGCGTCGCTGTGTGCGTGCTGGTCGAACAGCTGCCAGGCGCAGAAAAAATTCATGCCGGTGCCAAAACCGGTTTCGCCGATCACCGTGCAGGCGTGCGGCGCCAGGTCGGCAAAACGCTCGGCCAGGCGGGTCTgg                                                                                                                                             >  1:1856871/1‑150 (MQ=255)
             aGTGCAGGCGCGCGTCGCTGTGTGCGTGCTGGTCGAACAGCTGCCAGGCGCAGAAAAAATTCATGCCGGTGCCAAAACCGGTTTCGCCGATCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGa                                                                                                                                         <  1:2095341/150‑1 (MQ=255)
             aGTGCAGGCGCGCGTCGCTGTGTGCGTGCTGGTCGAACAGCTGCCAGGCGCAGAAAAAATTCATGCCGGTGCCAAAACCGGTTTCGCCGATCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGa                                                                                                                                         <  1:1311873/150‑1 (MQ=255)
             aGTGCAGGCGCGCGTCGCTGTGTGCGTGCTGGTCGAACAGCTGCCAGGCGCAGAAAAAATTCATGCCGGTGCCAAAACCGGTTTCGCCGATCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGACCCGCCTGGCCGa                                                                                                                                         <  1:4354648/150‑1 (MQ=255)
             aGTGCAGGCGCGCGTCGCTGTGTGCGTGCTGGTCGAACAGCTGCCAGGCGCAGAAAAAATTCATGCCGGTGCCAAAACCGGTTTCGCCGATCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGACCCGCCTGGCCGa                                                                                                                                         <  1:731245/150‑1 (MQ=255)
              gTGCAGGCGCGCGTCGCTGTGTGCGTGCTGGTCGAACAGCTGCCAGGCGCAGAAAAAATTCATGCCGGTGCCAAAACCGGTTTCGCCGATCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGACCCGCCTGGCCGAg                                                                                                                                        <  1:2626432/150‑1 (MQ=255)
                  aGGCGCGCGTCGCTGTGTGCGTGCTGGTCGAACAGCTGCCAGGCGCAGAAAAAATTCATGCCGGTGCCAAAACCGGTTTCGCCGATCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGACCCGCCTGGCCGa                                                                                                                                         <  1:225249/145‑1 (MQ=255)
                  aGGCGCGCGTCGCTGTGTGCGTGCTGGTCGAACAGCTGCCAGGCGCAGAAAAAATTCATGCCGGTGCCAAAACCGGTTTCGCCGATCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGACCCGCCTGGCCGa                                                                                                                                         >  2:225249/1‑145 (MQ=255)
                        gcgTCGCTGTGTGCGTGCTGGTCGAACAGCTGCCAGGCGCAGAAAAAATTCATGCCGGTGCCAAAACCGGTTTCGCCGATCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGc                                                                                                                              <  1:5335857/150‑1 (MQ=255)
                          gTCGCTGTGTGCGTGCTGGTCGAACAGCTGCCAGGCGCAGAAAAAATTCATGCCGGTGCCAAAACCGGTTTCGCCGATCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCtt                                                                                                                            <  1:2921371/150‑1 (MQ=255)
                           tCGCTGTGTGCGTGCTGGTCGAACAGCTGCCAGGCGCAGAAAAAATTCATGCCGGTGCCAAAACCGGTTTCGCCGATCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTg                                                                                                                           <  1:2512555/150‑1 (MQ=255)
                             gCTGTGTGCGTGCTGGTCGAACAGCTGCCAGGCGCAGAAAAAATTCATGCCGGTGCCAAAACCGGTTTCGCCGATCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGt                                                                                                                         <  1:1033106/150‑1 (MQ=255)
                              cTGTGTGCGTGCTGGTCGAACAGCTGCCAGGCGCAGAAAAAATTCATGCCGGTGCCAAAACCGGTTTCGCCGATCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGtt                                                                                                                        >  2:1925419/1‑150 (MQ=255)
                               tgtgtgCGTGCTGGTCGAACAGCTGCCAGGCGCAGAAAAAATTCATGCCGGTGCCAAAACCGGTTTCGCCGATCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGttt                                                                                                                       <  1:3818231/150‑1 (MQ=255)
                                  gtgcgtgcTGGTCGAACAGCTGCCAGGCGCAGAAAAAATTCATGCCGGTGCCAAAACCGGTTTCGCCGATCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTttctt                                                                                                                    >  2:3949518/1‑150 (MQ=255)
                                  gtgcgtgcTGGTCGAACAGCTGCCAGGCGCAGAAAAAATTCATGCCGGTGCCAAAACCGGTTTCGCCGATCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTttctt                                                                                                                    >  1:2197672/1‑150 (MQ=255)
                                     cgtgcTGGTCGAACAGCTGCCAGGCGCAGAAAAAATTCATGCCGGTGCCAAAACCGGTTTCGCCGATCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTc                                                                                                                                                >  2:5112399/1‑119 (MQ=255)
                                      gtgcTGGTCGAACAGCTGCCAGGCGCAGAAAAAATTCATGCCGGTGCCAAAACCGGTTTCGCCGATCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGAt                                                                                                                <  1:4057258/150‑1 (MQ=255)
                                       tgcTGGTCGAACAGCTGCCAGGCGCAGAAAAAATTCATGCCGGTGCCAAAACCGGTTTCGCCGATCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATg                                                                                                               <  1:2743029/150‑1 (MQ=255)
                                            gTCGAACAGCTGCCAGGCGCAGAAAAAATTCATGCCGGTGCCAAAACCGGTTTCGCCGATCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAATCGCTCGGCCAGGCGGGTCTGGCCGAGGACCACGTGCTTGGTTTCTTCGATGCCTTc                                                                                                          >  2:1758795/1‑150 (MQ=255)
                                            gTCGAACAGCTGCCAGGCGCAGAAAAAATTCATGCCGGTGCCAAAACCGGTTTCGCCGATCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGACCACGTGCTTGGTTTCTTCGATGCCTTc                                                                                                          >  1:3151697/1‑150 (MQ=255)
                                            gTCGAACAGCTGCCAGGCGCAGAAAAAATTCATGCCGGTGCCAAAACCGGTTTCGCCGATCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTc                                                                                                          >  2:1412880/1‑150 (MQ=255)
                                            gTCGAACAGCTGCCAGGCGCAGAAAAAATTCATGCCGGTGCCAAAACCGGTTTCGCCGATCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTc                                                                                                          >  2:4661505/1‑150 (MQ=255)
                                             tCGAACAGCTGCCAGGCGCAGAAAAAATTCATGCCGGTGCCAAAACCGGTTTCGCCGATCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCg                                                                                                         >  1:945543/1‑150 (MQ=255)
                                              cgcaCAGCTGCCACGCGCAGAAAAAATTCATGCCGGTGCCAAAACCGGTTTCGCCGATCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGt                                                                                                        <  2:4653151/147‑1 (MQ=255)
                                               gAACAGCTGCCAGGCGCAGAAAAAATTCATGCCGGTGCCAAAACCGGTTTCGCCGATCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGACCACGTGCTTGGTTTCTTCGATGCCTTCGtt                                                                                                       >  1:1469411/1‑150 (MQ=255)
                                               gAACAGCTGCCAGGCGCAGAAAAAATTCATGCCGGTGCCAAAACCGGTTTCGCCGATCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGtt                                                                                                       >  1:1756122/1‑150 (MQ=255)
                                               gAACAGCTGCCAGGCGCAGAAAAAATTCATGCCGGTGCCAAAACCGGTTTCGCCGATCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGtt                                                                                                       <  2:3652100/150‑1 (MQ=255)
                                                 aCAGCTGCCAGGCGCAGAAAAAATTCATGCCGGTGCCAAAACCGGTTTCGCCGATCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGTTGa                                                                                                     <  1:2427280/150‑1 (MQ=255)
                                                  cAGCTGCCAGGCGCAGAAAAAATTCATGCCGGTGCCAAAACCGGTTTCGCCGATCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGACCACGTGCTTGGTTTCTTCGATGCCTTCGTTGAc                                                                                                    >  1:628137/1‑150 (MQ=255)
                                                     cTGCCAGGCGCAGAAAAAATTCATGCCGGTGCCAAAACCGGTTTCGCCGATCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGACCCGCCTGGCCGAg                                                                                                                                        >  1:392796/1‑111 (MQ=255)
                                                     cTGCCAGGCGCAGAAAAAATTCATGCCGGTGCCAAAACCGGTTTCGCCGATCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGACCCGCCTGGCCGAg                                                                                                                                        <  2:392796/111‑1 (MQ=255)
                                                      tGCCAGGCGCAGAAAAAATTCATGCCGGTGCCAAAACCGGTTTCGCCGATCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGTTGACTGCg                                                                                                >  2:170450/1‑150 (MQ=255)
                                                        ccAGGCGCAGAAAAAATTCATGCCGGTGCCAAAACCGGTTTCGCCGATCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGTTGACTGCGaa                                                                                              <  2:3728315/150‑1 (MQ=255)
                                                        ccAGGCGCAGAAAAAATTCATGCCGGTGCCAAAACCGGTTTCGCCGATCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGTTGACTGCGaa                                                                                              <  2:2197672/150‑1 (MQ=255)
                                                        ccAGGCGCAGAAAAAATTCATGCCGGTGCCAAAACCGGTTTCGCCGATCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGTTGACTGCGaa                                                                                              >  2:3273595/1‑150 (MQ=255)
                                                          aGGCGCAGAAAAAATTCATGCCGGTGCCAAAACCGGTTTCGCCGATCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGTTGACTGCGAAAt                                                                                            >  2:4592357/1‑150 (MQ=255)
                                                            gcgcAGAAAAAATTCATGCCGGTGCCAAAACCGGTTTCGCCGATCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGACCACGTGCTTGGTTTCTTCGATGCCTTCGTTGACTGCGAAATAg                                                                                          >  2:3165051/1‑150 (MQ=255)
                                                            gcgcAGAAAAAATTCATGCCGGTGCCAAAACCGGTTTCGCCGATCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGTTGACTGCGAAATAg                                                                                          >  1:4029916/1‑150 (MQ=255)
                                                              gcAGAAAAAATTCATGCCGGTGCCAAAACCGGTTTCGCCGATCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGTTGACTGCGAAATAGAc                                                                                        <  2:431672/150‑1 (MQ=255)
                                                                aGAAAAAATTCATGCCGGTGCCAAAACCGGTTTCGCCGATCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGTTGACTGCGAAATAGAcgt                                                                                      >  1:3658003/1‑150 (MQ=255)
                                                                 gAAAAAATTCATGCCGGTGCCAAAACCGGTTTCGCCGATCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGTTGACTGCGAAATAGAcgtc                                                                                     <  1:5246334/150‑1 (MQ=255)
                                                                     aaaTTCATGCAGGTGCCAAAACCGGTTTCGCCGATCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGTTGACTGCGAAATAGACGTCGTCa                                                                                 >  2:3680482/1‑150 (MQ=255)
                                                                           aTGCCGGTGCCAAAACCGGTTTCGCCGATCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGTTGACTGCGAAATAGACGTCGTCATAGTGc                                                                           >  2:2753996/1‑150 (MQ=255)
                                                                                 gTGCCAAAACCGGTTTCGCCGATCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGACCACGTGCTTGGTTTCTTCGATGCCTTCGTTGACTGCGAAATAGACGTCGTCATAGTGCCGCGaa                                                                     >  1:4973481/1‑150 (MQ=255)
                                                                                  tGCCAAAACCGGTTTCGCCGATCACCGCGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGTTGACTGCGAAATAGACGTCGTCATAGTGCCGCGAAt                                                                    <  1:3959579/150‑1 (MQ=255)
                                                                                   gCCAAAACCGGTTTCGCCGATCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAac                                                                                                                                    <  2:2476631/85‑1 (MQ=255)
                                                                                     cAAAACCGGTTTCGCCGATCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGTTGACTGCGAAATAGACGTCGTCATAGTGCCGCGAATggg                                                                 >  2:3858799/1‑150 (MQ=255)
                                                                                       aaaCCGGTTTCGCCGATCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGTTGACTGCGAAATAGACGTCGTCATAGTGCCGCGAATggggg                                                               <  1:699712/150‑1 (MQ=255)
                                                                                        aaCCGGTTTCGCCGATCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAacac                                                                                                                                  >  1:3349139/1‑82 (MQ=255)
                                                                                               ttCGCCGATCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGTTGACTGCGAAATAGACGTCGTCATAGTGCCGCGAATGGGGGCGGCCCTg                                                       <  1:3583484/150‑1 (MQ=255)
                                                                                                 cGCCGATCACCGTGCAGGCGTGCGGCGCAAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGTTGACTGCGAAATAGACGTCGTCATAGTGCCGCGAATGGGGGCGGCCCTGgt                                                     <  2:631981/150‑1 (MQ=255)
                                                                                                  gCCGATCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGACCCCGTGCTTGGGTTCTTCGAGGCCTTCGTTGACTGCGAAATAGACGTCGTCATAGTGCCGCGAATGGGGGCGGCCCTGgtc                                                    <  2:2777595/150‑1 (MQ=255)
                                                                                                  gCCGATCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGTTGa                                                                                                     <  1:4255091/101‑1 (MQ=255)
                                                                                                  gCCGATCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGTTGa                                                                                                     >  2:4255091/1‑101 (MQ=255)
                                                                                                  gCCGATCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGTTGACTGCGAAATAGACGTCGTCATAGTGCCGCGAATGGGGGCGGCCCTGgtc                                                    >  2:3392504/1‑150 (MQ=255)
                                                                                                  gCCGATCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGTTGACTGCGAAATAGACGTCGTCATAGTGCCGCGAATGGGGGCGGCCCTGgtc                                                    >  1:2391259/1‑150 (MQ=255)
                                                                                                  gCCGATCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGTTGACTGCGAAATAGACGTCGTCATAGTGCCGCGAATGGGGGCGGCCCTGgtc                                                    >  2:2839281/1‑150 (MQ=255)
                                                                                                   ccGATCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGa                                                                                                                 <  1:3014844/88‑1 (MQ=255)
                                                                                                   ccGATCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGa                                                                                                                 >  2:3014844/1‑88 (MQ=255)
                                                                                                     gATCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGTTGACTGCGAAATAGACGTCGTCATAGTGCCGCGAATGGGGGCGGCCCTGgtcgtc                                                 <  2:1856747/150‑1 (MQ=255)
                                                                                                      aTCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGTTGACTGCGAAATAGACGTCGTCATAGTGCCGCGAATGGGGGCGGCCCTGGTCGTcc                                                <  2:2008833/150‑1 (MQ=255)
                                                                                                      aTCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGTTGACTGCGAAATAGACGTCGTCATAGTGCCGCGAATGGGGGCGGCCCTGGTCGTcc                                                <  2:4746796/150‑1 (MQ=255)
                                                                                                        cACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGTTGACTGCGAAATAGACGTCGTCATAGTGCCGCGAATGGGGGCGGCCCTGGTCGTCCCa                                              <  2:1044218/150‑1 (MQ=255)
                                                                                                         aCCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGTTGACTGCGAAATAGACGTCGTCATAGTGCCGCGAATGGGGGCGGCCCTGGTCGTCCCAg                                             <  1:516094/150‑1 (MQ=255)
                                                                                                           cGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGTTGACTGCGAAATAGACGTCGTCATAGTGCCGCGAATGGGGGCGGCCCTGGTCGTCCCAGTc                                           >  2:605917/1‑150 (MQ=255)
                                                                                                                  gCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGTTGACTGCGAAATAGACGTCGTCATAGTGCCGCGAATGGGGGCGGCCCTGGTCGTCCCGGTCGATCTgg                                    >  1:1398526/1‑150 (MQ=255)
                                                                                                                  gCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGTTGACTGCGAAATAGACGTCGTCATAGTGCCGCGAATGGGGGCGGCCCTGGTCGTCCCAGTCGATCTgg                                    <  2:748852/150‑1 (MQ=255)
                                                                                                                       cggcgCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGTTGACTGCGAAATAGACGTCGTCATAGTGCCGCGAATGGGGGCGGCCCTGGTCGTCCCAGTCGATCTGGGCGTg                               <  2:5215869/150‑1 (MQ=255)
                                                                                                                        ggcgCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGTTGACTGCGAAATAGACGTCGTCATAGTGCCGCGAATGGGGGCGGCCCTGGTCGTCCCAGTCGATCTGGGCGTGc                              <  1:1592516/150‑1 (MQ=255)
                                                                                                                           gcCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGTTGACTGCGAAATAGACGTCGTCATAGTGCCGCGAATGGGGGCGGCCCTGGTCGTCCCAGTCGATCTGGGCGTGCTgg                           <  2:3799808/150‑1 (MQ=255)
                                                                                                                            ccAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGTTGACTGCGAAATAGACGTCGTCATAGTGCCGCGAATGGGGGCGGCCCTGGTCGTCCCAGTCGATCTGGGCGTGCTgga                          <  1:4999347/150‑1 (MQ=255)
                                                                                                                               ggTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGTTGACTGCGAAATAGACGTCGTCATAGTGCCGCGAATGGGGGCGGCCCTGGTCGTCCCAGTCGATCTGGGCGTGCTggagga                       <  2:256423/150‑1 (MQ=255)
                                                                                                                                gTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGTTGACTGCGAAATAGACGTCGTCATAGTGCCGCGAATGGGGGCGGCCCTGGTCGTCCCAGTCGATCt                                      <  2:499879/134‑1 (MQ=255)
                                                                                                                                gTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGTTGACTGCGAAATAGACGTCGTCATAGTGCCGCGAATGGGGGCGGCCCTGGTCGTCCCAGTCGATCTGGGCGTGCTGgaggag                      >  2:4204130/1‑150 (MQ=255)
                                                                                                                                 tttgCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTTCTTGGTTTCTTCGATGCCTTCGTTGACTGAGAAATAGACGTCGTCATAGTGCCGCGAATGGTGGCGGCCCTGGTCGTCCCAGTCGATCTGGGCGTGCTGGaggagg                     <  1:3786386/147‑1 (MQ=255)
                                                                                                                                      aaaaCGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGTTGACTGCGAAATAGACGTCGTCATAGTGCCGCGAATGGGGGCGGCCCTGGTCGTCCCAGTCGATCTGGGCGTGCTGGAGGAGGGTggg                >  2:2800982/1‑150 (MQ=255)
                                                                                                                                      aaaaCGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGTTGACTGCGAAATAGACGTCGTCATAGTGCCGCGAATGGGGGCGGCCCTGGTCGTCCCAGTCGATCTGGGCGTGCTGGAGGAGGGTggg                >  1:1991447/1‑150 (MQ=255)
                                                                                                                                      aaaaCGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGTTGACTGCGAAATAGACGTCGTCATAGTGCCGCGAATGGGGGCGGCCCTGGTCGTCCCAGTCGATCTGGGCGTGCTGGAGGAGGGTggg                >  2:2854045/1‑150 (MQ=255)
                                                                                                                                      aaaaCGATCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGTTGACCGCGAAATAGACGTCGTCATAGTGCCGCGAATGGGGGCGGCCCTGGTCGTCCCAGTCGATCTGGGCGTGCTGGAGGAGGGTggg                >  2:478357/1‑150 (MQ=255)
                                                                                                                                       aaaCGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGTTGACTGCGAAATAGACGTCGTCATAGTGCCGCGAATGGGGGCGGCCCTGGTCGTCCCAGTCGATCTGGGCGTGCTGGaggagg                     <  1:2459151/144‑1 (MQ=255)
                                                                                                                                       aaaCGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGTTGACTGCGAAATAGACGTCGTCATAGTGCCGCGAATGGGGGCGGCCCTGGTCGTCCCAGTCGATCTGGGCGTGCTGGaggagg                     >  2:2459151/1‑144 (MQ=255)
                                                                                                                                        aaCGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGTTGACTGCGAAATAGACGTCGTCATAGTGCCGCGAATGGGGGCGGCCCTGGTCGTCCCAGTCGATCTGGGCGTGCTGGAGGAGGGTGGGCa              >  1:460897/1‑150 (MQ=255)
                                                                                                                                        aaCGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGTTGACTGCGAAATAGACGTCGTCATAGTGCCGCGAATGGGGGCGGCCCTGGTCGTCCCAGTCGATCTGGGCGTGCTGGAGGAGGGTGGGCa              >  2:2599673/1‑150 (MQ=255)
                                                                                                                                          cGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGTTGACTGCGAAATAGACGTCGTCATAGTGCCGCGAATGGGGGCGGCCCTGGTCGTCCCAGTCGATCTGGGCGTGCTGGAGGAGGGTGGGCAtg            >  1:3923853/1‑150 (MQ=255)
                                                                                                                                           gCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGTTGACTGCGAAATAGACGTCGTCATAGTGCCGCGAATGGGGGCGGCCCTGGTCGTCCCAGTCGATCTGGGCGTGCTGGAGGAGGGTGGGCAtgg           >  1:4586981/1‑150 (MQ=255)
                                                                                                                                           gCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGTTGACTGCGAAATAGACGTCGTCATAGTGCCGCGAATGGGGGCGGCCCTGGTCGTCCCAGTCGATCTGGGCGTGCTGGAGGAGGGTGGGCAtgg           >  2:1340541/1‑150 (MQ=255)
                                                                                                                                            cTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGTTGACTGCGAAATAGACGTCGTCATAGTGCCGCGAATGGGGGCGGCCCTGGTCGTCCCAGTCGATCTGGGCGTGCTGGAGGAGGGTGGGCAtggt          <  1:2805138/150‑1 (MQ=255)
                                                                                                                                             tCGGCCAGACCCGCCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGTTGACTGCGAAATAGACGTCGTCATAGTGCCGCGAATGGGggc                                                              <  2:1598761/97‑1 (MQ=255)
                                                                                                                                              cGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGTTGACTGCGAAATAGACGTCGTCATAGTGCCGCGAATGGGGGCGGCCCTGGTCGTCCCAGTCGATCTGGGCGTGCTGGAGGAGGGTGGGCAtggtgg        <  1:4407943/150‑1 (MQ=255)
                                                                                                                                               ggCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGTTGACTGCGAAATAGACGTCGTCATAGTGCCGCGAATGGGGGCGGCCCTGGTCGTCCCAGTCGATCTGGGCGTGCTGGAGGAGGGTGGGCATGGTggg       >  1:1165220/1‑150 (MQ=255)
                                                                                                                                                  cAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGTTGACTGCGAAATAGACGTCGTCATAGTGCCGCGAATGGGGGCGGCCCTGgtcgt                                                  >  1:4750846/1‑104 (MQ=255)
                                                                                                                                                   aGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGTTGACTGCGAAATAGACGTCGTCATAGTGCCGCGAATGGGGGCGGCCCTGGTCGTCCCAGTCGATCTGGGCGTGCTGGAGGAGGGTGGGCATGGTGGGCTCg   >  1:321893/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TCGACGCTGACGAAGTGCAGGCGCGCGTCGCTGTGTGCGTGCTGGTCGAACAGCTGCCAGGCGCAGAAAAAATTCATGCCGGTGCCAAAACCGGTTTCGCCGATCACCGTGCAGGCGTGCGGCGCCAGGTTGGCAAAACGCTCGGCCAGGCGGGTCTGGCCGAGGAACACGTGCTTGGTTTCTTCGATGCCTTCGTTGACTGCGAAATAGACGTCGTCATAGTGCCGCGAATGGGGGCGGCCCTGGTCGTCCCAGTCGATCTGGGCGTGCTGGAGGAGGGTGGGCATGGTGGGCTCGG  >  NC_002947/1953297‑1953594

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: