Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,917,404 A→T 5.8% 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 5.8% 226.5 / 2.1 76F25I (TTC→ATC) ccmAcytochrome c biogenesis ABC transporter ATP‑binding protein CcmA
Reads supporting (aligned to +/- strand):  ref base A (32/39);  new base T (2/3);  total (34/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCATGACGCTGTTCTGCCAGGGGCTTGCCACCCAGCAGAATCTGCCCGGCCGTCGGCTGCATCAGCCCGGCCAGCAGGCGCAACAGGCTGGTCTTGCCGCTGCCGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCCCAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCACGGTTTCGCCTTATCTTCAGTAGGGCTTTGAAGCAGCCAACACCGCTCCAGAACCCTTGTGTCTCAAGTCGCCCCGC  >  NC_002947/4917256‑4917552
                                                                                                                                                     |                                                                                                                                                    
gCATGACGCTGTTCTGCCAGGGGCTTGCCACCCAGCAGAATCTGCCCGGCCGTCGGCTGCATCAGCCCGGCCAGCAGGCGCAACAGGCTGGTCTTGCCGCTGCCGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCCCAGCTCGaaa                                                                                                                                                     <  1:31179/150‑1 (MQ=255)
    gACGCTGTTCTGCCAGGGGCTTGCCACCCAGCAGAATCTGCCCGGCCGTCGGCTGCATCAGCCCGGCCAGCAGGCGCAACAGGCTGGTCTTGCCGCTGCCGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCCCAGCTCGAAATCGa                                                                                                                                                 >  1:1771461/1‑150 (MQ=255)
    gACGCTGTTCTGCCAGGCAATCTCCACCCAGCAGAATCTGCCCGGCCGTCGGCTGCATAAGCCCGGCCAGAAGGCTCAACAGGCTGGTCTTGCCGCTGCCGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCCCAGCTCGAAATCGa                                                                                                                                                 <  1:1038159/150‑1 (MQ=255)
     aCGCTGTTCTGCCAGGGGCTTGCCACCCAGCAGAATCTGCCCGGCCGTCGGCTGCATCAGCCCGGCCAGCAGGCGCAACAGGCTGGTCTTGCCGCTGCCGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCCCAGCTCGAAATCGAg                                                                                                                                                <  2:1673845/150‑1 (MQ=255)
      cGCTGTTCTGCCAGGGGCTTGCCACCCAGCAGAATCTGCCCGGCCGTCGGCTGCATCAGCCCGGCCAGCAGGCGCAACAGGCTGGTCTTGCCGCTGCCGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCCCAGCTCGAAATCGAGc                                                                                                                                               <  1:806812/150‑1 (MQ=255)
       cgtGTTCCGCCTGGGTCGTACCACCCAGCAGAATCTGCCCGGCCGTCGGCTGCATCAGCCCGGCCAGCAGGCGCAACAGGCTGGTCTTGCTGCTGCCGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCCCAGCTCGAAATCGAGCt                                                                                                                                              <  2:870776/148‑1 (MQ=255)
         tGTTCTGCCAGGGGCTGGCCACCCAGCAGAATCTGCCCGGCCGTCGGCTGCATCAGCCCGGCCAGCAGGCGCAACAGGCTGGTCTTGCCGCTGCCGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCCCAGCTCGATTTCGAGCTgg                                                                                                                                            <  2:458308/150‑2 (MQ=255)
           ttCTGCCAGGGGCTTGCCACCCAGCAGAATCTGCCCGGCCGTCGGCTGCATCAGCCCGGCCAGCAGGCGCAACAGGCTGGTCTTGCCGCTGCCGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCCCAGCTCGAAATCGAGCTNTTc                                                                                                                                          <  2:1202620/150‑1 (MQ=255)
              tGCCAGGGGCTTGCCACCCAGCAGAATCTGCCCGGCCGTCGGCTGCATCCGCCCGGCCAGCAGGCGCAACAGGCTGGTCTTGCCGCTGCCGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCCCAGCTCGAAATCGAGCTGTTCGaa                                                                                                                                       <  1:1366080/150‑1 (MQ=255)
                     ggCTTGCCACCCAGCAGAATCTGCCCGGCCGTCGGCTGCATCAGCCCGGCCAGCAGGCGCAACAGGCTGGTCTTGCCGCTGCCGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCCCAGCTCGAAATCGAGCTGTTCGAACAGCAgg                                                                                                                                <  1:1109091/150‑1 (MQ=255)
                      gCTTGCCACCCAGCAGAATCTGCCCGGCCGTCGGCTGCATCAGCCCGGCCAGCAGGCGCAACAGGCTGGTCTTGCCGCTGCCGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCCCAGCTCGAAATCGAGCTGTTCGAACAGCAGgc                                                                                                                               <  2:1737211/150‑1 (MQ=255)
                        ttGCCACCCAGCAGAATCTGCCCGGCCGTCGGCTGCATCAGCCCGGCCAGCAGGCGCAACAGGCTGGTCTTGCCGCTGCCGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCCCAGCTCGAAATCGAGCTGTTCTAACAGCAGgcgc                                                                                                                             >  2:516819/1‑150 (MQ=255)
                               ccAGCAGAATCTGCCCGGCCGTCGGCTGCATCAGCCCGGCCAGCAGGCGCAACAGGCTGGTCTTGCCGCTGCCGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCCCAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTcgc                                                                                                                      <  1:426204/150‑1 (MQ=255)
                               ccAGCAGAATCTGCCCGGCCGTCGGCTGCATCAGCCCGGCCAGCAGGCGCAACAGGCTGGTCTTGCCGCTGCCGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCCCAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTcgc                                                                                                                      <  2:1560537/150‑1 (MQ=255)
                                 agcagAATCTGCCCGGCCGTCGGCTGCATCAGCCCGGCCAGCAGGCGCAACAGGCTGGTCTTGCCGCTGCCGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCCCAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCgcgc                                                                                                                    <  2:1127752/150‑1 (MQ=255)
                                  gcagAATCTGCCCGGCCGTCGGCTGCATCAGCCCGGCCAGCAGGCGCAACAGGCTGGTCTTGCCGCTGCCGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCCCAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCt                                                                                                                   >  2:1314973/1‑150 (MQ=255)
                                        tCTGCCCGGCCGTCGGCTGCATCAGCCCGGCCAGCAGGCGCAACAGGCTGGTCTTGCCGCTGCCGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCCCAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGcagg                                                                                                             >  2:1526308/1‑150 (MQ=255)
                                            cccGGCCGTCGGCTGCATCAGCCCGGCCAGCAGGCGCAACAGGCTGGTCTTGCCGCTGCCGTTGGGGCCACTCATCTGCAGCATGTCGCCGGCGCCCAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGcaggccag                                                                                                         <  1:751493/150‑1 (MQ=255)
                                              cGGCCGTCGGCTGCATCAGCCCGGCCAGCAGGCGCAACAGGCTGGTCTTGCCGCTGCCGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCCCAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGcaggccaggc                                                                                                       <  2:379656/150‑1 (MQ=255)
                                                  cGTCGGCTGCATCAGCCCGGCCAGCAGGCGCAACAGGCTGGTCTTGCCGCTGCCGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCCCAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCcgc                                                                                                   >  2:192754/1‑150 (MQ=255)
                                                    tCGGCTGCATCAGCCCGGCCAGCAGGCGCAACAGGCTTGTCTTGCCGCTGCCGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCCCAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCgg                                                                                                 >  2:917661/1‑150 (MQ=255)
                                                     cGGCTGCATCAGCCCGGCCAGCAGGCGCAACAGGCTGGTCTTGCCGCTGCCGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCCCAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGc                                                                                                <  2:297602/150‑1 (MQ=255)
                                                       gCTGCATCAGCCCGGCCAGCAGGCGCAACAGGCTGGTCTTGCCGCTGCCGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCCCAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCtt                                                                                              >  2:1807889/1‑150 (MQ=255)
                                                            aTCAGCCCGGCCAGCAGGCGCAACAGGCTGGTCTTGCCGCTGCCGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCCCAGCTCGAAATCGAGCTGTTCAAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAgg                                                                                         >  1:1551645/1‑150 (MQ=255)
                                                                gCCCGGCCAGCAGGCGCAACAGGCTGGTCTTGCCGCTGCCGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCCCAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGcag                                                                                                              >  2:587957/1‑125 (MQ=255)
                                                                gCCCGGCCAGCAGGCGCAACAGGCTGGTCTTGCCGCTGCCGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCCCAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGcag                                                                                                              <  1:587957/125‑1 (MQ=255)
                                                                gCCCGGCCAGCAGGCGCAACAGGCTGGTCTTGCCGCTGCCGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCCCAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGaa                                                                                     >  2:554033/1‑150 (MQ=255)
                                                                 cccGGCCAGCAGGCGCAACAGGCTGGTCTTGCCGCTGCCGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCCCAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAg                                                                                    <  1:554033/150‑1 (MQ=255)
                                                                  ccGGCCAGCAGGCGCAACAGGCTGGTCTTGCCGCTGCCGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCCCAGCTCGAAATCGAGCTGTTCGCGCAGCAGGCGCCAGTCGCCCTCCCAGGCCAGGTCCGCGCCCTGGAGGTGAAgg                                                                                   >  2:116636/1‑150 (MQ=255)
                                                                  ccGGCCAGCAGGCGCAACAGGCTGGTCTTGCCGCTGCCGACGCGGTCCACGCGGCTTACCATGTCGCCGGCGCCCAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAgg                                                                                   <  2:50052/150‑1 (MQ=255)
                                                                      ccAGCAGGCGCAACAGGCTGGTCTTGCCGCTGCCGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCCCAGCTCGAAATCGCGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCa                                                                               <  2:138030/150‑1 (MQ=255)
                                                                      ccAGCAGGCGCAACAGGCTGGTCTTGCCGCTGCCGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCCCAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCa                                                                               <  2:1079080/150‑1 (MQ=255)
                                                                       cagcagGCGCAACAGGCTGGTCTTGCCGCTGCCGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCCCAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCAc                                                                              <  2:1786193/150‑1 (MQ=255)
                                                                        agcagGCGCAACAGGCTGGTCTTGCCGCTGCCGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCCCAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCACg                                                                             >  1:10534/1‑150 (MQ=255)
                                                                          cagGCGCAACAGGCTGGTCTTGCCGCTGCCGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCCCAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCACGGt                                                                           >  1:1636002/1‑150 (MQ=255)
                                                                            ggCGCAACAGGCTGGTCTTGCCGCTGCCGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCCCAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCACGGttt                                                                         >  1:1028684/1‑150 (MQ=255)
                                                                               gcAACAGGCTTGTCTTGCCGCTGCCGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCCCAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCACGGTTTCGc                                                                      <  1:917661/150‑1 (MQ=255)
                                                                               gcAACAGGCTGGTCTTGCCGCTGCCGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCCCAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCACGGTTTCGc                                                                      <  1:848310/150‑1 (MQ=255)
                                                                               gcAACAGGCTGGTCTTGCCGCTGCCGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCCCAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCACGGTTTCGc                                                                      <  2:1771461/150‑1 (MQ=255)
                                                                                 aaCAGGCTGGTCTTGCCGCTGCCGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCCCAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCACGGTTTCGCCt                                                                    >  2:1109668/1‑150 (MQ=255)
                                                                                   cAGGCTGGTCTTGCCGCTGCCGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCCCAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGTTGAAGGGTCACGGTTTCGCCTTa                                                                  <  1:258386/150‑1 (MQ=255)
                                                                                      gCTGGTCTTGCCGCTGCCGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCCCAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCACGGTTTCGCCTTATCt                                                               >  1:958393/1‑150 (MQ=255)
                                                                                        tcgTCTTGCCGCTGCCGTTGGCGCCACTGACCTGCAGACCAAACCCGTCCCACCATTTGTTACCAAGCTGTTCGAACAGCAGACGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCACGGTTTCGCATTATCTTc                                                             >  2:544105/3‑150 (MQ=255)
                                                                                         ggTCTTGCCGCTGCCGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCCCAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCACGGTTTCGCCTTATCTTCa                                                            <  1:192754/150‑1 (MQ=255)
                                                                                          ttCTTGCCGCTGCCGTTGGGGCACCTGATCTGCCGCATCCCGCCGGCGCCCAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCACGGTTTCGCCTTATCTTCAg                                                           <  2:958393/149‑1 (MQ=255)
                                                                                          gTCTTGCCGCTGCCGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCCCAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCACGGTTTCGCCTTATCTTCAg                                                           >  1:234150/1‑150 (MQ=255)
                                                                                            cTTGCCGCTGCCGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCCCAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCACGGTTTCGCCTTATCTTCAGTa                                                         <  1:1526308/150‑1 (MQ=255)
                                                                                            cTTGCCGCTGCCGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCCCAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCACGGTTTCGCCTTATCTTCAGTa                                                         <  2:1694427/150‑1 (MQ=255)
                                                                                             ttGCCGCTGCCGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCCCAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCACGGTTTCGCGTTATCTTCAGTAg                                                        >  1:388795/1‑150 (MQ=255)
                                                                                                     gCCGTTGGGGCCACTGATGTGCAGCATGTCGCCGGCGCCCAGCTCGAAATCGCGCTGTTCGAACAGCAGGCGCGAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCACGGTTTCGCCTTATCTTCAGTAGGGCTTTGa                                                <  2:358984/150‑1 (MQ=255)
                                                                                                      ccGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCCCAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCACGGTTTCGCCTTATCTTCAGTAGGGCTTTGaa                                               <  2:1279758/150‑1 (MQ=255)
                                                                                                      ccGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCCCAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCACGGTTTCGCCTTATCTTCAGTAGGGCTTTGaa                                               >  1:1757612/1‑150 (MQ=255)
                                                                                                        gTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCCCAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCACGGTTTCGCCTTATCTTCAGTAGGGCTTTGAagc                                             >  2:1122252/1‑150 (MQ=255)
                                                                                                         ttGGGGCCACTGATCTGCAGCATGTCGCCGGCGCCCAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCACGGATCTACTTTATCTTCAGTAGGGCTTTGAagca                                            <  1:527739/150‑1 (MQ=255)
                                                                                                         ttGGGGCCACTGATCTCCAGCATGTCGCCGGCGCCCAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCACGGTTTCGCCTTATCTTCAGTAGGGCTTTGAagca                                            <  2:969358/150‑1 (MQ=255)
                                                                                                           ggggCCACTGATCTGCAGCATGTCGCCGGCGCCCAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCACGGTTTCGCCTTATCTTCAGTAGGGCTTTGAagcagc                                          >  2:2145/1‑150 (MQ=255)
                                                                                                               gCCACTGATCTGCAGCATGTCGCCGGCGCCCAGCTCGAAATCGAGCTGTTCTAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCACGGTTTCGCCTTATCTTCAGTAg                                                        <  1:1900047/133‑1 (MQ=255)
                                                                                                               gCCACTGATCTGCAGCATGTCGCCGGCGCCCAGCTCGAAATCGAGCTGTTCTAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCACGGTTTCGCCTTATCTTCAGTAg                                                        >  2:1900047/1‑133 (MQ=255)
                                                                                                                  aCTGATCTGCAGCATGTCGCCGGCGCCCAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCACGGTTTCGCCTTATCTTCAGTAGGGCTTTGAAGCAGCCAacac                                    >  2:1333723/1‑150 (MQ=255)
                                                                                                                  aCTGATCTGCAGCATGTCGCCGGCGCCCAGCAGGTTTTCGTGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCACGGTTTCGCCTTATCTTCAGTAGGGCTTTGAAGCAGCCAacac                                    <  2:151014/150‑1 (MQ=255)
                                                                                                                       tCTGCCGCATGTCGCCTGCGCCCAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCACGGTTTCGCCTTATCTTCAGTAGGGCTTTGAAGCAGCCAACACCGCTc                               <  2:527129/150‑1 (MQ=255)
                                                                                                                         tGCAGCATGTCGCCGGCGCCCAGCTCGAAATCGAGCTGTTCGCAAAGAAGAACGTAGTCGCGATCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCACGGTTTCGCCTTATCTTCAGTAGGGCTTTGAAGCAGCCAACACCGCTCCa                             >  1:1470509/1‑150 (MQ=255)
                                                                                                                         tGCAGCATGTCGCCGGCGCCCAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGa                                                                                           <  1:1387784/88‑1 (MQ=255)
                                                                                                                         tGCAGCATGTCGCCGGCGCCCAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGa                                                                                           >  2:1387784/1‑88 (MQ=255)
                                                                                                                                  tCGCCGGCGCCCAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCACGGTTTCGCCTTATCTTCAGTAGGGCTTTGAAGCAGCCAACACCGCTCCAGAACCCTtg                    >  2:1758805/1‑150 (MQ=255)
                                                                                                                                   cGCCGGCGCCCAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCACGGTTTCGCCTTATCTTCAGTAGGGCTTTGAAGCAGCCAACACCGCTCCAGAACCCTtgt                   <  2:1028684/150‑1 (MQ=255)
                                                                                                                                   cGCCGGCGCCCAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCACGGTTTCGCCTTATCTTCAGTAGGGCTTTGAAGCAGCCAACACCGCTCCAGAACCCTtgt                   <  1:1122252/150‑1 (MQ=255)
                                                                                                                                      cGGCGCCCAGCTCGCAATCCAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCACGGTTTCGCCTTATCTTCAGTAGGGCTTTGAAGCAGCCAACACCGCTCCAGAACCCTTGTGtc                <  2:1883814/150‑1 (MQ=255)
                                                                                                                                       ggCGCCCAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCACGGTTTCGCCTTATCTTCAGTAGGGCTTTGAAGCAGCCAACACCGCTCCAGAACCCTTGTGtct               >  1:1900807/1‑150 (MQ=255)
                                                                                                                                            ccAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCACGGTTTCGCCTTATCTTCAGTAGGGCTTTGAAGCAGCCAACACCGCTCCAGAACCCTTGTGTCTCAAGt          >  2:633899/1‑150 (MQ=255)
                                                                                                                                            acAGCTCGATTTCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCACGGTTTCGCCTTATCTTCAGTAGGGCTTTGaa                                               >  2:1151864/2‑113 (MQ=255)
                                                                                                                                            acAGCTCGATTTCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCACGGTTTCGCCTTATCTTCAGTAGGGCTTTGaa                                               <  1:1151864/112‑1 (MQ=255)
                                                                                                                                             cAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCACGGTTTCGCCTTATCTTCAGTAGGGCTTTGAAGCAGCCAACACCGCTCCAGAACCCTTGTGTCTCAAGTc         >  2:319418/1‑150 (MQ=255)
                                                                                                                                                cTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCACGGTTTCGCCTTATCTTCAGTAGGGCTTTGAAGCAGCCAACACCGCTCCAGAACCCTTGTGTCTCAAGTCGcc      >  2:1309247/1‑150 (MQ=255)
                                                                                                                                                 tCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCACGGTTTCGCCTTATCTTCAGTAGGGCTTTGAAGCAGCCAACACCGCTCCAGAACCCTTGTGTCTCAAGTCGccc     >  1:726227/1‑150 (MQ=255)
                                                                                                                                                 tCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCACGGTTTCGCCTTATCTTCAGTAGGGCTTTGAAGCAGCCAACACCGCTCCAGAACCCTTGTGTCTCAAGTCGccc     >  1:1544804/1‑150 (MQ=255)
                                                                                                                                                  cGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCACGGTTTCGCCTTATCTTCAGTAGGGCTTTGAAGCAGCCAACACCGCTCCAGAACCCTTGTGTCTCAAGTCGcccc    <  1:1314973/150‑1 (MQ=255)
                                                                                                                                                    aaaTCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCACGGTTTCGCCTTATCTTCAGTAGGGCTTTGAAGCAGCCAACACCGCTCCAGAACCCTTGTGTCTCAAGTCGCCCcgc  <  2:1757612/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GCATGACGCTGTTCTGCCAGGGGCTTGCCACCCAGCAGAATCTGCCCGGCCGTCGGCTGCATCAGCCCGGCCAGCAGGCGCAACAGGCTGGTCTTGCCGCTGCCGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCCCAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCACGGTTTCGCCTTATCTTCAGTAGGGCTTTGAAGCAGCCAACACCGCTCCAGAACCCTTGTGTCTCAAGTCGCCCCGC  >  NC_002947/4917256‑4917552

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: