Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,917,404 A→T 7.1% F25I (TTC→ATC)  ccmA ← cytochrome c biogenesis ABC transporter ATP‑binding protein CcmA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,917,4040AT 7.1% 302.6 / 11.2 99F25I (TTC→ATC) ccmAcytochrome c biogenesis ABC transporter ATP‑binding protein CcmA
Reads supporting (aligned to +/- strand):  ref base A (42/49);  new base T (3/4);  total (46/53)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01

Only 100 of 102 total aligned reads displayed.

GCATGACGCTGTTCTGCCAGGGGCTTGCCACCCAGCAGAATCTGCCCGGCCGTCGGCTGCATCAGCCCGGCCAGCAGGCGCAACAGGCTGGTCTTGCCGCTGCCGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCC‑‑CAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCACGGTTTCGCCTTATCTTCAGTAGGGCTTTGAAGCAGCCAACACCGCTCCAGAACCCTTGTGTCTCAAGTCGCCCCG  >  NC_002947/4917256‑4917551
                                                                                                                                                      |                                                                                                                                                   
gCATGACGCTGTTCTGCCAGGGGCTTGCCACCCAGCAGAATCTGCCCGGCCGTCGGCTGCATCAGCCCGGCCAGCAGGCGCAACAGGCTGGTCTTGCCGCTGCCGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCC‑‑CAGCTCGaaa                                                                                                                                                    <  2:1422264/150‑1 (MQ=255)
     aCGCTGTTCTGCCAGGGGCTTGCCACCCAGCAGAATCTGCCCGGCCGTCGGCTGCATCAGCCCGGCCAGCAGGCGCAACAGGCTGGTCTTGCCGCTGCCGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCC‑‑CAGCTCGAAATCGAg                                                                                                                                               <  1:1128028/150‑1 (MQ=255)
       gCTGTTCTGCCAGGGGCTTGCCACCCAGCAGAATCTGCCCGGCCGTCGGCTGCATCAGCCCGGCCAGCAGGCGCAACAGGCTGGTCTTGCCGCTGCCGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCC‑‑CAGCTCGAAATCGAGCt                                                                                                                                             >  2:1682331/1‑150 (MQ=255)
         tGTTCTGCGAGGGGCTTGCCACCCAGCAGAAGCTGCCCGGCCGTCGGCTGCATCAGCCCGGCCAGCAGGCGCAACAGGCTGGTCTTGCCGCTGCCGTTGGGGNCACTGANCTGCAGCATGTCNCCGGCGCC‑‑CAGCTCGAAATCGAGCTGt                                                                                                                                           <  1:1068597/150‑1 (MQ=255)
         tGTTCTGCCAGGGGCTTGCCACCCAGCAGAATCTGCCCGGCCGTCGGCTGCATCAGCCCGGCCAGCAGGCGCAACAGGCTGGTCTTGCCGCTGCCGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCC‑‑CAGCTCGAAATCGAGCTGt                                                                                                                                           <  1:756634/150‑1 (MQ=255)
           ttCTGCCGGGGGCGTGCCACCCAGCGGAATCTGCCCGGCCGTGGGCTTCATCAGCCCGGCGAGCAGGTGCAACAGGCTGGTCTTGCTGCTGCCGTTGGGGCCACTGATCTGCAGCATGGCGCCGGCGCC‑‑CAGCTCGATGTCGAGCttggc                                                                                                                                         <  2:255187/150‑5 (MQ=255)
           ttCTGCCAGGGGCTTGCCACCCAGCAGAATCTGCCCGGCCGTCGGCTGCATCAGCCCGGCCAGCAGGCGCAACAGGCTGGTCTTGCCGCTGCCGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCC‑‑CAGCTCGAAATCGAGCTGTTc                                                                                                                                         <  1:1673732/150‑1 (MQ=255)
           ttCTGCCAGGGGCTTGCCACCCAGCAGAATCTGCCCGGCCGTCGGCTGCATCAGCCCGGCCAGCAGGCGCAACAGGCTGGTCTTGCCGCTGCCGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCC‑‑CAGCTCGAAATCGAGCTGTTc                                                                                                                                         >  1:150006/1‑150 (MQ=255)
             cTGCCAGGGGCTTGCCACCCAGCAGAATCTGCCCGGCCGTCGGCTGCATCAGCCCGGCCAGCAGGCGCAACAGGCTGGTCTTGCCGCTGCCGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCC‑‑CAGCTCGAAATCGAGCTGTTCGa                                                                                                                                       >  1:825398/1‑150 (MQ=255)
               gCCAGGGGCTTGCCACCCAGCAGAATCTGCCCGGCCGTCGGCTGCATCAGCCCGGCCAGCAGGCGCAACAGGCTGGTCTTGCCGCTGCCGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCC‑‑CAGCTCGaaa                                                                                                                                                    <  2:376437/135‑1 (MQ=255)
               gCCAGGGGCTTGCCACCCAGCAGAATCTGCCCGGCCGTCGGCTGCATCAGCCCGGCCAGCAGGCGCAACAGGCTGGTCTTGCCGCTGCCGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCC‑‑CAGCTCGaaa                                                                                                                                                    >  1:376437/1‑135 (MQ=255)
                   ggggCTTGCCACCCAGCAGAATCTGCCCGGCCGTCGGCTGCATCAGCCCGGCCAGCAGGCGCAACAGGCTGGTCTTGCCGCTGCCGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCC‑‑CAGCTCGAAATCGAGCTGTTCGAAcagca                                                                                                                                 <  2:150006/150‑1 (MQ=255)
                    gggCTTGCCCCCCAGCAGAATCTGCCCGGCCGTCGGCTGCATCAGCCCGGCCAGCAGGCGCAACAGGCTGGTCTTGCCGCTGCCGTTGGGGCCACTGATTTGCAGCATGTCGCCGGCACC‑‑CAGCTCAAAGGCAAGCTGTTCGGAcaggag                                                                                                                                >  2:771575/1‑150 (MQ=255)
                        ttGCCACCCAGCAGAATCTGCCCGGCCGTCGGCTGCATCAGCCCGGCCAGCAGGCGCAACAGGCTGGTCTTGCCGCTGCCGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCC‑‑CAGCTCGAAATCGAGCTGTTCGAACAGCAGgcgc                                                                                                                            <  2:1723395/150‑1 (MQ=255)
                             aCCCAGCAGAATCTGCCCGGCCGTCGGCTGCATCAGCCCGGCCAGCAGGCGCAACAGGCTGGTCTTGCCGCTGCCGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCC‑‑CAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTc                                                                                                                       >  2:541524/1‑150 (MQ=255)
                             aCCCAGCAGAATCTGCCCGGCCGTCGGCTGCATCAGCCCGGCCAGCAGGCGCAACAGGCTGGTCTTGCCGCTGCCGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCC‑‑CAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTc                                                                                                                       <  1:1211321/150‑1 (MQ=255)
                             aCCCAGCAGAATCTGCCCGGCCGTCGGCTGCATCAGCCCGGCCAGCAGGCGCAACAGGCTGGTCTTGCCGCTGCCGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCC‑‑CAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCACTc                                                                                                                       <  1:1934386/150‑1 (MQ=255)
                                  gcagAATCTGCCCGGCCGTCGGCTGCATCAGCCCGGCCAGCAGGCGCAACAGGCTGGTCTTGCCGCTGCCGTTAGGGCCACTGATCTGCAGCATGTCGCCGGCGCC‑‑CAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCt                                                                                                                  >  2:232857/1‑150 (MQ=255)
                                    agAATCTGCCCGGCCGTCGGCTGCATCAGCCCGGCCAGCAGGCGCAACAGGCTGGTCTTGCCGCTGCCGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCC‑‑CAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCg                                                                                                                >  1:1028577/1‑150 (MQ=255)
                                    agAATCTGCCCGGCCGTCGGCTGCATCAGCCCGGCCAGCAGGCGCAACAGGCTGGTCTTGCCGCTGCCGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCC‑‑CAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCg                                                                                                                >  1:1028576/1‑150 (MQ=255)
                                    agAATCTGCCCGGCCGTCGGCTGCATCAGCCCGGCCAGCAGGCGCAACAGGCTGGTCTTGCCGCTGCCGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCC‑‑CAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCg                                                                                                                <  2:869253/150‑1 (MQ=255)
                                      aaTCTGCCCGGCCGTCGGCTGCATCAGCCCGGCCAGCAGGCGCAACAGGCTGGTCTTGCCGCTGCCGTTGGGGCCACTGATCTGCAGCATGTCGCCCGGCTG‑‑CATTGCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGca                                                                                                              <  1:1320647/150‑1 (MQ=255)
                                        tCTGCCCGGCCGTCGGCTGCATCAGCCCGGCCAGCAGGCGCAACAGGCTGGTCTTGCCGCTGCCGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCC‑‑CAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGcagg                                                                                                            >  1:1281248/1‑150 (MQ=255)
                                             ccGGCCGTCGGCTGCATCAGCCCGGCCAGCAGGCGCAACAGGCTGGTCTTGCCGCTGCCGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCC‑‑CAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGcaggccagg                                                                                                       >  2:1427669/1‑150 (MQ=255)
                                              cGGCCGTCGGCTGCATCAGCCCGGCCAGCAGGCGCAACAGGCTGGTCTTGCCGCTGCCGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCC‑‑CAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGcaggccaggc                                                                                                      >  2:880689/1‑150 (MQ=255)
                                              cGGCCGTCGGCTGCATCAGCCCGGCCAGCAGGCGCAACAGGCTGGTCTTGCCGCTGCCGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCC‑‑CAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGcaggccaggc                                                                                                      <  1:1421206/150‑1 (MQ=255)
                                                  cGTCGGCTGCATCAGCCCGGCCAGCAGGCGCAACAGGCTGGTCTTGCCGCTGCCGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCC‑‑CAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCcgc                                                                                                  <  1:1093984/150‑1 (MQ=255)
                                                   gTCGGCTGCATCAGCCCGGCCAGCAGGCGCAACAGGCTGGTCTTGCCGCTGCCGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCC‑‑CAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCcgcg                                                                                                 <  2:1099967/150‑1 (MQ=255)
                                                       gCTGCATCAGCCCGGCCAGCAGGCGCAACAGGCTGGTCTTGCCGCTGCCGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCC‑‑CAGCTCGAAATCGAGCTGTTCGAAca                                                                                                                                    <  2:268326/111‑1 (MQ=255)
                                                       gCTGCATCAGCCCGGCCAGCAGGCGCAACAGGCTGGTCTTGCCGCTGCCGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCC‑‑CAGCTCGAAATCGAGCTGTTCGAAca                                                                                                                                    >  1:268326/1‑111 (MQ=255)
                                                         tGCATCAGCCCGGCCAGCAGGCGCAACAGGCTGGTCTTGCCGCTGCCGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCC‑‑CAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTgg                                                                                           <  2:825398/150‑1 (MQ=255)
                                                         tGCATCAGCCCGGCCAGCAGGCGCAACAGGCTGGTCTTGCCGCTGCCGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCC‑‑CAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTgg                                                                                           <  1:541524/150‑1 (MQ=255)
                                                         tGCATCAGCCCGGCCAGCAGGCGCAACAGGCTGGTCTTGCCGCTGCCGTTGGGGCCACTGATCTGCAGCATATGGCCGGCGCC‑‑CAGCTCGACATAGAGCTGTTCTAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTgg                                                                                           >  1:1400250/1‑150 (MQ=255)
                                                             tCAGCCCGGCCAGCAGGCGCAACAGGCTGGTCTTGCCGCTGCCGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCC‑‑CAGCTCGAAATCGAGc                                                                                                                                              <  2:475558/95‑1 (MQ=255)
                                                             tCAGCCCGGCCAGCAGGCGCAACAGGCTGGTCTTGCCGCTGCCGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCC‑‑CAGCTCGAAATCGAGc                                                                                                                                              >  1:475558/1‑95 (MQ=255)
                                                             tCAGCCCGGCCAGCAGGCGCAACAGGCTGGTCTTGCCGCTGCCGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCC‑‑CAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGt                                                                                       >  1:1158569/1‑150 (MQ=255)
                                                                gCCCGGCCAGCAGGCGCAACAGGCTGGTCTTGCCGCTGCCGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCC‑‑CAGCTCGATTTCGAGCTgg                                                                                                                                           <  1:938626/95‑2 (MQ=255)
                                                                gCCCGGCCAGCAGGCGCAACAGGCTGGTCTTGCCGCTGCCGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCC‑‑CAGCTCGATTTCGAGCTgg                                                                                                                                           >  2:938626/1‑94 (MQ=255)
                                                                   cGGCCAGCAGGCGCAACAGGCTGGTCTTGCCGCTGCCGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCC‑‑CAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAggg                                                                                 >  2:462540/1‑150 (MQ=255)
                                                                   cGGCCAGCAGGCGCAACAGGCTGGTCTTGCCGCTCCCGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCC‑‑CAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCACGCGGCTTGGAGGTGAAggg                                                                                 >  2:384530/1‑150 (MQ=255)
                                                                      ccAGCAGGCGCAACAGGCTGGTCTTGCCGCTGCCGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCC‑‑CAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCa                                                                              <  1:880689/150‑1 (MQ=255)
                                                                       cagcagGCGCAACAGGCTGGTCTTGCCGCTGCCGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCC‑‑CAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCAc                                                                             <  1:1682331/150‑1 (MQ=255)
                                                                        agcagGCGCAACAGGCTGGTCTTGCCGCTGCCGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCC‑‑CAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCCGGCCCGCGGCTTGGAGGTGAAgg                                                                                  >  2:451459/1‑144 (MQ=255)
                                                                        agcagGCGCAACAGGCTGGTCTTGCCGCTGCCGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCC‑‑CAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAgg                                                                                  <  1:451459/144‑1 (MQ=255)
                                                                          cagGCGCAACAGGCTGGTCTTGCCGCTGCCGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCC‑‑CAGCTCGATTTCGAGCTgggc                                                                                                                                         <  2:1919105/87‑4 (MQ=39)
                                                                          cagGCGCAACAGGCTGGTCTTGCCGCTGCCGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCC‑‑CAGCTCGATTTCGAGCTgggc                                                                                                                                         >  1:1919105/1‑84 (MQ=39)
                                                                          cagGCGCAACAGGCTGGTCTTGACGCTGCCGTTTTGGCCATGGATCTGCAGCATGTCGCCGGCGCC‑‑CAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCACGGt                                                                          <  1:821049/150‑1 (MQ=255)
                                                                             gcgcAACAGGCTGGTCTTGCCGCTGCCGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCC‑‑CAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCGGCGGCTTGGGGGTGAAGGGTCACGGTTTc                                                                       >  1:1093732/1‑150 (MQ=255)
                                                                             gcgcAACAGGCTGGTCTTGCCGCTGCCGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCC‑‑CAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCACGGTTTc                                                                       <  2:1860926/150‑1 (MQ=255)
                                                                             gcgcAACAGGCTGGTCTTGCCGCTGCCGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCC‑‑CAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGGGAAGGGTCACGGTTTc                                                                       >  1:922307/1‑150 (MQ=255)
                                                                             ccgcAACAGTCCGTGCTTGCCGCTGCCGTTGGGGCCACTGATCTGCAGCAAGTTCCTGTGACTGGAGTTCAGAAGTGGGTCTTCGAACTTTACGCGCCTGTCGCGCTTGCTGGCCTCGCTCGCTGCTTGGAGGTGAAGGGTCACGGTTTc                                                                       <  1:925613/149‑1 (MQ=255)
                                                                               gcAACAGGCTGGTCTTGCCGCTGCCGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCC‑‑CAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAggg                                                                                 <  2:831765/138‑1 (MQ=255)
                                                                               gcAACAGGCTGGTCTTGCCGCTGCCGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCC‑‑CAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCACGGTTTCGc                                                                     <  1:1427669/150‑1 (MQ=255)
                                                                                  aCAGGCTGGTCTTGCCGCTGCCGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCC‑‑CAGCTCGAAATCGAGCTGTTCGaa                                                                                                                                      >  1:1563059/1‑82 (MQ=255)
                                                                                  aCAGGCTGGTCTTGCCGCTGCCGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCC‑‑CAGCTCGAAATCGAGCTGTTCGaa                                                                                                                                      <  2:1563059/82‑1 (MQ=255)
                                                                                   cAGGCTGGTCTTGCCGCTGCCGTTGGCGCGACTGATCTGCAGCATGTCGCCGGGGCC‑‑CAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCACGTTTTCGCCTTa                                                                 <  1:1858554/150‑1 (MQ=255)
                                                                                     ggCTGGTCTTGCCGCTGCCGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCC‑‑CAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCACGGtt                                                                         <  1:284830/140‑1 (MQ=255)
                                                                                     ggCTGGTCTTGCCGCTGCCGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCC‑‑CAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCACGGtt                                                                         >  2:284830/1‑140 (MQ=255)
                                                                                           tCTTGCCGCTGCCGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCC‑‑CAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCACGGttt                                                                        <  1:233776/135‑1 (MQ=255)
                                                                                           tCTTGCCGCTGCCGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCC‑‑CAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCACGGttt                                                                        >  2:233776/1‑135 (MQ=255)
                                                                                              tGCCGCTGCCGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCC‑‑CAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCACGGTTTCGCCTTATCTTCAGTAgg                                                      <  2:1158569/150‑1 (MQ=255)
                                                                                               gCCGCTGCCGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCC‑‑CAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCACGGTTTCGCCTTATCTTCAGTAggg                                                     >  1:799694/1‑150 (MQ=255)
                                                                                                ccGCTGCCGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCC‑‑CAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCACGGTTTCGCCTTATCTTCAGTAGGGc                                                    >  1:1497574/1‑150 (MQ=255)
                                                                                                  gCTGCCGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCC‑‑CAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAg                                                                                         <  2:1767099/111‑1 (MQ=255)
                                                                                                  gCTGCCGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCC‑‑CAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAg                                                                                         >  1:1767099/1‑111 (MQ=255)
                                                                                                   cTGCCGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCC‑‑CAGATCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCACGGTTTCGCCTTATCTTCAGTAGGGCttt                                                 <  2:1400250/150‑1 (MQ=255)
                                                                                                      ccGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCC‑‑CAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCACGGTTTCGCCTTATCTTCAGTAGGGCTTTGaa                                              <  1:1562642/150‑1 (MQ=255)
                                                                                                        gTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCC‑‑CAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCACGGTTTCGCCTTATCTTCAGTAGGGCTTTGAagc                                            <  2:1028577/150‑1 (MQ=255)
                                                                                                        gTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCC‑‑CAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCACGGTTTCGCCTTATCTTCAGTAGGGCTTTGAagc                                            <  2:1028576/150‑1 (MQ=255)
                                                                                                         ttGGGGCCACTGATCTGCAGCATGTCGCCGGCGCC‑‑CAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCACGGTTTCGCCTTATCTTCAGTAGGGCTTTGAagca                                           <  2:1497574/150‑1 (MQ=255)
                                                                                                         ttGGGGCCACTGATCTGCAGCATGTCGCCGGCGCC‑‑CAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCACGGTTTCGCCTTATCTTCAGTAGGGCTTTGAagca                                           <  2:799694/150‑1 (MQ=255)
                                                                                                           ggggCCACTGATCTGCAGCATGTCGCCGGCGCC‑‑CAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCACGGTTTCGCCTTATCTTCAGTAGGGCTTTGAagcagc                                         >  1:357742/1‑150 (MQ=255)
                                                                                                           ggggCCACTGATCTGCAGCATGTCGCCGGCGCC‑‑CAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCACGGTTTCGCCTTATCTTCAGTAGGGCTTTGAagcagc                                         <  1:1835987/150‑1 (MQ=255)
                                                                                                           ggggCCACTGATCTGCAGCATGTCGCCGGCGCC‑‑CAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCACGGTTTCGCCTTATCTTCAGTAGGGCTTTGAagcagc                                         <  1:1448903/150‑1 (MQ=255)
                                                                                                           ggggCCACTGATCTGCAGCATGTCGCCGGCGCC‑‑CAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCACGGTTTCGCCTTATCTTCAGTAGGGCTTTGAagcagc                                         <  1:290586/150‑1 (MQ=255)
                                                                                                           ggggCCACTGATCTGCAGCATGTCGCCGGCGCC‑‑CAGCTCGAAATCGAGCGGAGCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCACGGTTTCGCCTTATCTTCAGTAGTGCTTTGAagcagc                                         >  2:1528840/1‑150 (MQ=255)
                                                                                                            gggCCACTGATCTGCAGCATGTCGCCGGCGCC‑‑CAGCTCGAAATCGAGCTGTTCGaa                                                                                                                                      >  1:376169/1‑56 (MQ=255)
                                                                                                            gggCCACTGATCTGCAGCATGTCGCCGGCGCC‑‑CAGCTCGAAATCGAGCTGTTCGaa                                                                                                                                      <  2:376169/56‑1 (MQ=255)
                                                                                                            gggCCACTGATCTGCAGCATGTCGCCGGCGCC‑‑CAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCACGGTTTCGCCTTATCTTCAGTAGGGCATTGAAGCAGcc                                        >  1:298661/1‑150 (MQ=255)
                                                                                                              gCCACTGATCTGCAGCATGTCGCCGGCGCC‑‑CCGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGCGGTCAAGGGTCACGGTTTCGCCTTATCTTCAGTAGGGCTTTGAAGCAGCCaa                                      <  2:273153/150‑1 (MQ=255)
                                                                                                                    gATCTGCAGCATGTCGCCGGCGCC‑‑CAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAg                                                                                   >  2:1763247/1‑99 (MQ=255)
                                                                                                                    gATCTGCAGCATGTCGCCGGCGCC‑‑CAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAg                                                                                   <  1:1763247/99‑1 (MQ=255)
                                                                                                                     aTCTGCAGCATGTCGCCGGCGCC‑‑CAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCACGGTTTCGCCTTATCTTCAGTa                                                        <  2:43622/125‑1 (MQ=255)
                                                                                                                     aTCTGCAGCATGTCGCCGGCGCC‑‑CAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCACGGTTTCGCCTTATCTTCAGTa                                                        >  1:43622/1‑125 (MQ=255)
                                                                                                                         gcagcaTGTCGCCGGCGCC‑‑CAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCACGGTTTCGCCTTATCTTCAGTAGGGCTTTGAAGCAGCCAACACCGCTCCAg                           >  2:1600933/1‑150 (MQ=255)
                                                                                                                          cagcaTGTCGCCGGCGCC‑‑CAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCACGGTTTCGCCTTATCTTCAGTAGGGCTTTGAAGCAGCCAACACCGCTCCAGa                          >  2:1569240/1‑150 (MQ=255)
                                                                                                                                  cGCCGGCGCC‑‑CAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCACGGTTTCGCCTTATCTTCAGTAGGGCTTTGAATCAGCCAACACCGCTCCAGAACCCTtgt                  >  2:317432/1‑150 (MQ=255)
                                                                                                                                  cGCCGGCGCC‑‑CAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCACGGTTTCGCCTTATCTTCAGTAGGGCTTTGAAGCAGCCAACACCGCTCCAGAACCCTtgt                  <  1:1600933/150‑1 (MQ=255)
                                                                                                                                           gaacAGCTCGATTTCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCACGGTTTCGCCTTATCTTCAGTAGGGCTTTGAag                                             <  2:320455/113‑1 (MQ=255)
                                                                                                                                           gaacAGCTCGATTTCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCACGGTTTCGCCTTATCTTCAGTAGGGCTTTGAag                                             >  1:320455/4‑116 (MQ=255)
                                                                                                                                           c‑‑cAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCACGGTTTCGCCTTATCTTCAGTAGGGCTTTGAagca                                           >  2:1095239/1‑116 (MQ=255)
                                                                                                                                           c‑‑cAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCACGGTTTCGCCTTATCTTCAGTAGGGCTTTGAagca                                           <  1:1095239/116‑1 (MQ=255)
                                                                                                                                           c‑‑cAGCTCGAAATCGAGCTGTTCGAACAGCAGGCACCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCACGGTTTCGCCTTATCTTCAGTAGGGCTTTGAagca                                           >  2:1029129/1‑116 (MQ=255)
                                                                                                                                           c‑‑cAGCTCGAAATCGAGCTGTTCGAACAGCAGGCACCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCACGGTTTCGCCTTATCTTCAGTAGGGCTTTGAagca                                           <  1:1029129/116‑1 (MQ=255)
                                                                                                                                              cAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCACGGTTTCGCCTTATCTTCAGTAGGGCTTTGAAGCAGCCAACACCGCTCCAGAACCCTTGTGTCTCAAGTc        >  1:1665687/1‑150 (MQ=255)
                                                                                                                                                gCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCACGGTTTCGCCTTATCTTCAGTAGGGCTTTGAAGCAGCCAACACCGCTCCAGAACCCTTGTGTCTCAAGTCGc      <  1:384530/150‑1 (MQ=255)
                                                                                                                                                 cTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCACGGTTTCGCCTTATCTTCAGTAGGGCTTTGAAGCAGCCAACACCGCTCCAGAACCCTTGTGTCTCAAGTCGcc     >  2:593661/1‑150 (MQ=255)
                                                                                                                                                 cTCGAAATCGAGCTGTTCGAACAGCAGGCGCAACGCGCTCCCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCCCAGTGTCGCCTTATCTTCAGTAGGGCTTTGAAGCAGCCAACACCGCTCCAGAACCCTTGTGTCTCAAGTCGcc     >  1:646864/1‑150 (MQ=255)
                                                                                                                                                    gAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAgg                                                                                  >  1:987361/1‑70 (MQ=255)
                                                                                                                                                    gAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCACGGTTTCGCCTTATCTTCAGTAGGGCTTTGAAGCAGCCAACACCGCTCCAGAACCCTTGTGTCTCAAGTCGCCCcg  >  1:1548285/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                   
GCATGACGCTGTTCTGCCAGGGGCTTGCCACCCAGCAGAATCTGCCCGGCCGTCGGCTGCATCAGCCCGGCCAGCAGGCGCAACAGGCTGGTCTTGCCGCTGCCGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCC‑‑CAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCACGGTTTCGCCTTATCTTCAGTAGGGCTTTGAAGCAGCCAACACCGCTCCAGAACCCTTGTGTCTCAAGTCGCCCCG  >  NC_002947/4917256‑4917551

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