Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,978,415 G→A 5.2% R754C (CGC→TGC)  cobN ← aerobic cobaltochelatase subunit CobN

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,978,4150GA 5.2% 412.0 / 7.1 136R754C (CGC→TGC) cobNaerobic cobaltochelatase subunit CobN
Reads supporting (aligned to +/- strand):  ref base G (71/57);  new base A (3/4);  total (74/62)
Fisher's exact test for biased strand distribution p-value = 7.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.71e-01

Only 100 of 140 total aligned reads displayed.

ACCGGTTGCCACTCACTCAGGTCGGGTAACTGTGCCGTGCCCTCCAATGCCTGCCCGATCACCTGCAGTGCCAGCAGCTCAAGGCGCTCGCGGGTGTCGCCACAGGTGCGCCAAGGTTCATTGCTCATGGCCTGCAGCTGTTCAGGGCGTATGCCCTGCCACGGCTGCCCAAGGTCGCAGTCGAGCGGGTCCCAGCCCGGCACCAGTGCTTTGGCCAAAGCCCGCAGCAGGCTGGCATTACCCCCACGACCGTCGCCACGCTCGACCCGCAACAACGCCAGCAGCGTATCCAGGC  >  NC_002947/3978267‑3978561
                                                                                                                                                    |                                                                                                                                                   
ccctgTTGCCACTCACTCAGGTCGGGTAACTGTGCCGTGCCCTCCAATGCCTGCCCGATCACCTGCAGTGCCAGCAGCTCAAGGCGCTCGCGGGTGTCGCCACAGGTGCGCCAAGGTTCATTGCTCATGGCCTGCAGCTGTTCAGGGCGt                                                                                                                                                    >  1:1264895/5‑150 (MQ=255)
aCCGGTTGCCACTCACTCAGGTCGGGTAACTGTGCCGTGCCCTCCAATGCCTGCCCGATCACCTGCAGTGCCAGCAGCTCAAGGCGCTCGCGGGTGTCGCCACAGGTGCGCCAAGGTTCATTGCTCATGGCCTGCAGCTGTTCAGGGCGt                                                                                                                                                    >  1:1768759/1‑150 (MQ=255)
   ggTTGCCACTCACTCAGGTCGGGTAACTGTGCCGTGCCCTCCAATGCCTGCCCGATCACCTGCAGTGCCAGCAGCTCAAGGCGCTCGCGGGTGTCGCCACAGGTGCGCCAAGGTTCATTGCTCATGGCCTGCAGCTGTTCAGGGCGTATg                                                                                                                                                 >  1:2346880/1‑150 (MQ=255)
    gTTGCCACTCACTCAGGTCGGGTAACTGTGCCGTGCCCTCCAATGCCTGCCCGATCACCTGCAGTGCCAGCAGCTCAAGGCGCTCGCGGGTGTCGCCACAGGTGCGCCAAGGTTCATTGCTCATGGCCTGCAGCTGTTCAGGGCGTATGc                                                                                                                                                >  2:2147397/1‑150 (MQ=255)
     ttGCCACTCACTCAGGTCGGGTAACTGTGCCGTGCCCTCCAATGCCTGCCCGATCACCTGCAGTGCCAGCAGCTCAAGGCGCTCGCGGGTGTCGCCACAGGTGCGCCAAGGTTCATTGCTCATGGCCTGCAGCTGTTCAGGGCATACGcc                                                                                                                                               <  2:2016186/150‑1 (MQ=255)
       gCCACTCACTCAGGTCGGGTAACTGTGCCGTGCCCTCCAATGCCTGCCCGATCACCTGCAGTGCCAGCAGCTCAAGGCGCTCGCGGGTGTCGCCACAGGTGCGCCAAGGTTCATTGCTCATGGCCTGCAGCTGTTCAGGGCATACGcc                                                                                                                                               <  2:1699281/148‑1 (MQ=255)
       gCCACTCACTCAGGTCGGGTAACTGTGCCGTGCCCTCCAATGCCTGCCCGATCACCTGCAGTGCCAGCAGCTCAAGGCGCTCGCGGGTGTCGCCACAGGTGCGCCAAGGTTCATTGCTCATGGCCTGCAGCTGTTCAGGGCATACGcc                                                                                                                                               >  1:1699281/1‑148 (MQ=255)
          actcactcAGGTCGGGTAACTGTGCCGTGCCCTCCAATGCCTGCCCGATCACCTGCAGTGCCAGCAGCTCAAGGCGCTCGCGGGTGTCGCCACAGGTGCGCCAAGGTTCATTGCTCATGGCCTGCAGCTGTTCAGGGCATACGccc                                                                                                                                              >  1:426182/1‑146 (MQ=255)
                tcaGGTCGGGTAACTGTGCCGTGCCCTCCAATGCCTGCCCGATCACCTGCAGTGCCAGCAGCTCAAGGCGCTCGCGGGTGTCGCCACAGGTGCGCCAAGGTTCATTGCTCATGGCCTGCAGCTGTTCAGGGCGTATGCCCTGCCACgg                                                                                                                                      >  2:1175728/1‑148 (MQ=255)
                tcaGGTCGGGTAACTGTGCCGTGCCCTCCAATGCCTGCCCGATCACCTGCAGTGCCAGCAGCTCAAGGCGCTCGCGGGTGTCGCCACAGGTGCGCCAAGGTTCATTGCTCATGGCCTGCAGCTGTTCAGGGCGTATGCCCTGCCACgg                                                                                                                                      <  1:1175728/148‑1 (MQ=255)
                tcaGGTCGGGTAACTGTGCCGTGCCCTCCAATGCCTGCCCGATCACCTGCAGTGCCAGCAGCTCAAGGCGCTCGCGGGTGTCGCCACAGGTGCGCCAAGGTTCATTGCTCATGGCCTGCAGCTGTTCAGGGCGTATGCCCTGCCACGGCt                                                                                                                                    >  1:2203254/1‑150 (MQ=255)
                tcaGGTCGGGTAACTGTGCCGTGCCCTCCAATGCCTGCCCGATCACCTGCAGTGCCAGCAGCTCAAGGCGCTCGCGGGTGTCGCCACAGGTGCGCCAAGGTTCATTGCTCATGGCCTGCAGCTGTTCAGGGCGTATGCCCTGCCACGGCt                                                                                                                                    >  1:1042480/1‑150 (MQ=255)
                 caGGTCGGGTAACTGTGCCGTGCCCTCCAATGCCTGCCCGATCACCTGCAGTGCCAGCAGCTCAAGGCGCTCGCGGGTGTCGCCACAGGTGCGCCAAGGTTCATTGCTCATGGCCTGCAGCTGTTCAGGGCGTATGCCCTGCCACGGCTg                                                                                                                                   <  1:1806775/150‑1 (MQ=255)
                  gggTCGGGTAACTGTGCCGTGCCCTCCAATGCCTTCCCGATCACCTGCAGTGCCAGCAGCTCAAGGCGCTCGCGGGTGTCGCCACAGGTGCGCCAAGGTTCATTGCTCATGGCCTGCAGCTGTTCAGGGCGTATGCCCTGCCACGGCTGc                                                                                                                                  <  2:737310/149‑1 (MQ=255)
                  aGGTCGGGTAACTGTGCCGTGCCCTCCAATGCCTGCCCGATCACCTGCAGTGCCAGCAGCTCAAGGCGCTCGCGGGTGTCGCCACAGGTGCGCCAAGGTTCATTGCTCATGGCCTGCAGCTGTTCAGGGCGTATGCCCTGCCACGGCTGc                                                                                                                                  <  1:360795/150‑1 (MQ=255)
                      cGGGTAACTGTGCCGTGCCCTCCAATGCCTGCCCGATCACCTGCAGTGCCAGCAGCTCAAGGCGCTCGCGGGTGTCGCCACAGGTGCGCCAAGGTTCATTGCTCATGGCCTGCAGCTGTTCAGGGCGTATGCCCTGCCACGGCTGCCCaa                                                                                                                              <  2:419840/150‑1 (MQ=255)
                        ggTAACTGTGCCGTGCCCTCCAATGCCTGCCCGATCACCTGCAGTGCCAGCAGCTCAAGGCGCTCGCGGGTGTCGCCACAGGTGCGCCAAGGTTCATTGCTCATGGCCTGCAGCTGTTCAGGGCGTATGCCCTGCCACGGCTGCCCAAgg                                                                                                                            <  1:1328073/150‑1 (MQ=255)
                          caacgctgccgtcccCTCCACTGCCTGCCCGATCACCTGCAGTGCCAGCAGCTCAAGGCGCTCGCGGGTGTCGCCACAGGTGCGCCAAGGTTCATTGCTCATGGCCTGCAGCTGTTCAGGGCGTATGCCCTGCCACGGCTGCCCAAGGTc                                                                                                                          <  2:10616/144‑1 (MQ=255)
                              tgtgCCGTGCCCTCCAATGCCTGCCCGATCACCTGCAGTGCCAGCAGCTCAAGGCGCTCGCGGGTGTCGCCACAGGTGCGCCAAGGTTCATTGCTCATGGCCTGCAGCTGTTCAGGGCGTATGCCCTGCCACGGCTGCCCAAGGTc                                                                                                                          <  2:2298335/146‑1 (MQ=255)
                              tgtgCCGTGCCCTCCAATGCCTGCCCGATCACCTGCAGTGCCAGCAGCTCAAGGCGCTCGCGGGTGTCGCCACAGGTGCGCCAAGGTTCATTGCTCATGGCCTGCAGCTGTTCAGGGCGTATGCCATGCGACGTCTGTCCAAGGTCGCAg                                                                                                                      >  2:1365104/1‑150 (MQ=255)
                                tgccgtgccCTCCAATGCCTGCCCGATCACCTGCAGTGCCAGCAGCTCAAGGCGCTCGCGGGTGTCGCCACAGGTGCGCCAAGGTTCATTGCTCATGGCCTGCAGCTGTTCAGGGCGTATGCCCTGCCACGGCTGCCCAAGGTCGCAGTc                                                                                                                    <  2:450326/150‑1 (MQ=255)
                                  ccgtgccCTCCAATGCCTGCCCGATCACCTGCAGTGCCAGCAGCTCAAGGCGCTCGCGGGTGTCGCCACAGGTGCGCCAAGGTTCATTGCTCATGGCCTGCAGCTGTTCAGGGCGTATGCCCTGCCACGGCTGCCCAAGGTCGCAGTCGa                                                                                                                  <  1:2147397/150‑1 (MQ=255)
                                  ccgtgccCTCCAATGCCTGCCCGATCACCTGCAGTGCCAGCAGCTCAAGGCGCTCGCGGGTGTCGCCACAGGTGCGCCAAGGTTCATTGCTCATGGCCTGCAGCTGTTCAGGGCGTATGCCCTGCCACGGCTGCCCAAGGTCGCAGTCGa                                                                                                                  <  1:1564894/150‑1 (MQ=255)
                                   cgtgccCTCCAATGCCTGCCCGATCACCTGCAGTGCCAGCAGCTCAAGGCGCTCGCGGGTGTCGCCACAGGTGCGCCAAGGTTCATTGCTCATGGCCTGCAGCTGTTCAGGGCg                                                                                                                                                     <  2:2066327/114‑1 (MQ=255)
                                   cgtgccCTCCAATGCCTGCCCGATCACCTGCAGTGCCAGCAGCTCAAGGCGCTCGCGGGTGTCGCCACAGGTGCGCCAAGGTTCATTGCTCATGGCCTGCAGCTGTTCAGGGCg                                                                                                                                                     >  1:2066327/1‑114 (MQ=255)
                                   cgtgccCTCCAATGCCTGCCCGATCACCTGCAGTGCCAGCAGCTCAAGGCGCTCGCGGGTGTCGCCACAGGTGCGCCAAGGTTCATTGCTCATGGCCTGCAGCTGTTCAGGGCGTATGCCCTGCCACGGCTGCCCAAGGTCGCAGTCGAg                                                                                                                 >  1:1130073/1‑150 (MQ=255)
                                       cccGCCAATGCCTGCCCGATCACCTGCAGTGCCAGCAGCTCAAGGCGCTCGCGGGTGTCGCCACAGGTGCGCCAAGGTTCATTGCTCATGGCCTGCAGCTGTTCAGGGCGTATGCCCTGCCACGGCTGCCCAAGGTCGCAGTCGAGCggg                                                                                                             <  2:1150124/150‑1 (MQ=255)
                                        ccTCCAATGCCTGCCCGATCACCTGCAGTGCCAGCAGCTCAAGGCGCTCGCGGGTGTCGCCACAGGTGCGCCAAGGTTCTAGCCTCATGGCCTGCAGCTGTTCAGGGCGTATGCCCTGCCACGGCTGCCCAAGGTCGCAGTCGAGCGGGt                                                                                                            >  1:1972008/1‑150 (MQ=255)
                                        ccTCCAATGCCTGCCCGATCACCTGCAGTGCCAGCAGCTCAAGGCGCTCGCGGGTGTCGACACAGGTGCGCCAAGGTTCATTGCTCATGGCCTGCAGCTGTTCAGGGCGTATGCCCTGCCACGGCTGCCCAAGGTCGCAGTCGAGCGGGt                                                                                                            >  1:1611875/1‑150 (MQ=255)
                                         cTCCAATGCCTGCCCGATCACCTGCAGTGCCAGCAGCTCAAGGCGCGCGCGGGTGTCGCCACAGGTGCGCCAAGGTTCATTGCTCATGGCCTGCAGCTGTTCAGGGCGTATGCCCTGCCACGGCTGCCCAAGGTCGCAGTCGAGCgg                                                                                                              <  1:143335/147‑1 (MQ=255)
                                          tCCAATGCCTGCCCGATCACCTGCAGTGCCAGCAGCTCAAGGCGCTCGCGGGTGTCGCCACAGGTGCGCCAAGGTTCATTGCTCATGGCCTGCAGCTGTTCAGGGCGTATGCCCTGCCACGGCt                                                                                                                                    >  1:1651574/1‑124 (MQ=255)
                                          tCCAATGCCTGCCCGATCACCTGCAGTGCCAGCAGCTCAAGGCGCTCGCGGGTGTCGCCACAGGTGCGCCAAGGTTCATTGCTCATGGCCTGCAGCTGTTCAGGGCGTATGCCCTGCCACGGCt                                                                                                                                    <  2:1651574/124‑1 (MQ=255)
                                          tCCAATGCCTGCCCGATCACCTGCAGTGCCAGCAGCTCAAGGCGCTCGCGGGTGTCGCCACAGGTGCGCCAAGGTTCATTGCTCATGGCCTGCAGCTGTTCAGGGCGTATGCCCTGCCACGGCTGCCCAAGGTCGCAGTCGAGCGGGTcc                                                                                                          >  1:1329671/1‑150 (MQ=255)
                                           ccAACGCCTGCCCGATCCGCTGCAGCGCCAGCAGCTCAAGGCGCTCGCGGGTGTCGCCACAGGTGCGCCAAGGTTCATTGCTCATGGCCTGCAGCTGTCCAGGGCGGATGCCCTGCCACGGCTGCCCAAGGTCGCAGGCGAGCGGGTccc                                                                                                         <  2:737468/150‑1 (MQ=255)
                                            cAATGCCTGCCCGATCACCTGCAGTGCCAGCAGCTCAAGGCGCTCGCGGGTGTCGCCACAGGTGCGCCAAGGTTCATTGCTCATGGCCTGCAGCTGTTCAGGGCGTATGCCCTGCCACGGCTGCCCAAGGTCGCAGTCGAGCGGGTCCCa                                                                                                        >  1:612456/1‑150 (MQ=255)
                                                gcctgccCGATCACCTGCAGTGCCAGCAGCTCAAGGCGCTCGCGGGTGTCGCCACAGGTGCGCCAAGGTTCATTGCTCATGGCCTGCAGCTGTTCAGGGCGTATGCCCTGCCACGGCTGCCCAAGGTCGCAGTCGAGCGGGTCCCAGccc                                                                                                    <  2:316712/150‑1 (MQ=255)
                                                gcctgccCGATCACCTGCAGTGCCAGCAGCTCAAGGCGCTCGCGGGTGTCGCCACAGGTGCGCCAAGGTTCATTGCTCATGGCCTGCAGCTGTTCAGGGCGTATGCCCTGCCACGGCTGCCCAAGGTCGCAGTCGAGCGGGTCCCAGccc                                                                                                    <  2:1329671/150‑1 (MQ=255)
                                                  ctgccCGATCACCTGCAGTGCCAGCAGCTCAAGGCGCTCGCGGGTGTCGCCACAGGTGCGCCAAGGTTCATTGCTCACGGCCTGCAGCTGTTCAGGGCGTATGCCCTGCCACGGCTGCCAAAGGTCGCAGTCGAGCGGGTCCCAGCCCgg                                                                                                  <  1:2286967/150‑1 (MQ=255)
                                                   tgccCGATCACCTGCAGTGCCAGCAGCTCAAGGCGCTCGCGGGTGTCGCCACAGGTGCGCCAAGGTTCATTGCTCATGGCCTGCAGCTGTTCAGGGCGTATGCCCTGCCACGGCTGCCCAAGGTCGCAGTCGAGCGGGTCCCAGCCCGGc                                                                                                 >  2:265497/1‑150 (MQ=255)
                                                   tgccCGATCACCTGCAGTGCCAGCAGCTCAAGGCGCTCGCGGGTGTCGCCACAGGTGCGCCAAGGTTCATTGCTCATGGCCTGCAGCTGTTCAGGGCGTATGCCCTGCCACGGCTGCCCAAGGTCGCAGTCGAGCGGGTCCCAGCCCGGc                                                                                                 >  2:1876509/1‑150 (MQ=255)
                                                     cccGATCACCTGCAGTGCCAGCAGCTCAAGGCGCTCGCGGGTGTCGCCACAGGTGCGCCAAGGTTCATTGCTCATGGCCTGCAGCTGTTCAGGGCGTATGCCCTGCCACGGCTGCCCAAGGTCGCAGTCGAGCGGGTCCCAGCCCGGCAc                                                                                               >  2:1422332/1‑150 (MQ=255)
                                                      ccGATCACCTGCAGTGCCAGCAGCTCAAGGCGCTCGCGGGTGTCGCCACAGGTGCGCCAAGGTTCATTGCTCATGGCCTGCAGCTGTTCAGGGCGTATGCCCTGCCACGGCTGCCCAAGGTCGCAGTCGAGCGGGTCCCAGCCCGGCAcc                                                                                              <  2:2321378/150‑1 (MQ=255)
                                                       cGATCACCTGCAGTGCAAGCAGCTCAAGGCGCTCGCGGGTGTCGCCACAGGTGCGCCAAGGTTCATTGCTCATGGCCTGCAGCTGTTCAGGGCGTATGCCCTGCCACGGCTGCCCAAGGTCGCAGTCGAGCGGGTCCCAGCCCGGCACCa                                                                                             <  2:596430/150‑1 (MQ=255)
                                                        agtCACCTGCAGTGCCAGCAGCTCAAGGCGCTCGCGGGTGTCGCCACAGGTGCGCCAAGGTTCATTGCTCATGGCCTGCAGCTGTTCAGGGCGTATGCCCTGCCACGGCTGCCCAAGGTCGCAGTCGAGCGGGTCCCAGCCCGGCACCAg                                                                                            >  2:171117/3‑150 (MQ=255)
                                                              cTGCAGTGCCAGCAGCTCAGCGCGCTCGCCCGTGTCGCCACAGGTGCGCCAAGGTGCATTGCTCATGGCCTGCAGCTGTTCAGGGCCTATGCCCTGCCACGGCTGCCCAAGGTCGCAGTCGAGCGGGTCCCAGCCCGGCACCAGTGCttt                                                                                      <  2:390724/150‑1 (MQ=255)
                                                                 cAGTGCCAGCAGCTCAAGGCGCTCGCGGGTGTCGCCACAGGTGCGCCAAGGTTCATTGCTCATGGCCTGCAGCTGTTCAGGGCGTATGCCCTGCCACGGCTGCCCAAGGTCGCAGTCGAGCGGGTCCCAGCCCGGCACCAGTGCTTTGGc                                                                                   <  2:612456/150‑1 (MQ=255)
                                                                 cAGTGCCAGCAGCTCAAGGCGCTCGCGGGTGTCGCCACAGGTGCGCCAAGGTTCATTGCTCATGGCCTGCAGCTGTTCAGGGCGTATGCCCTGCCACGGCTGCCCAAGGTCGCAGTCGACGTAGTCCCAGCCCGGCACCAGTGCTTTGGc                                                                                   <  1:62482/150‑1 (MQ=255)
                                                                      ccAGCAGCTCAAGGCGCTCGCGGGTGTCGCCACAGGTGCGCCAAGGTTCATTGCTCATGGCCTGCAGCTGTTCAGGGCGTATGCCCTGCCACGGCTGCCCAAGGTCGCAGTCGAGCGGGTCCCAGCCCGGCACCAGTGCTTTGGCCAAAg                                                                              >  1:1701529/1‑150 (MQ=255)
                                                                       cagcagCTCAAGGCGCTCGCGGGTGTCGCCACAGGTGCGCCAAGGTTCATTGCTCATGGCCTGCAGCTGTTCAGGGCGTATGCCCTGCCACGGCTGCCCAAGGTCGCAGTCGAGCGGGTCCCAGccc                                                                                                    <  2:1818844/127‑1 (MQ=255)
                                                                       cagcagCTCAAGGCGCTCGCGGGTGTCGCCACAGGTGCGCCAAGGTTCATTGCTCATGGCCTGCAGCTGTTCAGGGCGTATGCCCTGCCACGGCTGCCCAAGGTCGCAGTCGAGCGGGTCCCAGCCCGGCACCAGTGCTTTGGCCAAAGc                                                                             <  2:2346880/150‑1 (MQ=255)
                                                                        agcagcTCAAGGCGCTCGCGGGTGTCGCCACAGGTGCGCCAAGGTTCATTGCTCATGGCCTGCAGCTGTTCAGGGCGTATGCCCTGCCACGGCTGCCCAAGGTCGCAGTCGAGCGGGTCCCAGCCCGGCACCAGTGCTTTgg                                                                                    >  2:1333045/1‑142 (MQ=255)
                                                                        agcagcTCAAGGCGCTCGCGGGTGTCGCCACAGGTGCGCCAAGGTTCATTGCTCATGGCCTGCAGCTGTTCAGGGCGTATGCCCTGCCACGGCTGCCCAAGGTCGCAGTCGAGCGGGTCCCAGCCCGGCACCAGTGCTTTGGCCAAAGcc                                                                            >  2:160528/1‑150 (MQ=255)
                                                                        agcagcTCAAGGCGCTCGCGGGTGTCGCCACAGGTGCGCCAAGGTTCATTGCTCATGGCATGCAGCTGTTCAGGGCGTATGCCCTGCCACGGCTGCCCAAGGTCGCAGTCGAGCGGGTCCCAGCCCGGCACCAGTGCTTTgg                                                                                    <  1:1333045/142‑1 (MQ=255)
                                                                          cagcTCAAGGCGCTCGCGGGTGTCGCCACAGGTGCGCCAAGGTTCATTGCTCATGGCCTGCAGCTGTTCAGGGCGTATGCCCTGCCACGGCTGCCCAAGGTCGCCGCCGAGCGGGTCCCAGCCCGGCACCAGTGCTTTGGCCAAAGCCCg                                                                          <  2:1264895/150‑1 (MQ=255)
                                                                          cagcTCAAGGCGCTCGCGGGTGTCGCCACAGGTGCGCCAAGGTTCATTGCTCATGGCCTGCAGCTGTTCAGGGCGTATGCCCTGCCACGGCTGCCCAAGGTCGCAGTCGAGCGGGTCCCAGCCCGGCACCAGTGCTTTGGCCAAAGCCCg                                                                          <  1:1422332/150‑1 (MQ=255)
                                                                                aaGGCGCTCGCGGGTGTCGCCACAGGTGCGCCAAGGTTCATTGCTCATGGCCTGCAGCTGTTCAGGGCGTATGCCCTGCCACGGCTGCCCAAGGTCGCAGTCGAGCGGGTCCCAGCCCGGCACCAGTGCTTTGGCCAAAGCCCGcagcag                                                                   >  1:2261554/1‑150 (MQ=255)
                                                                                aaGGCGCTCGCGGGTGTCGCCACAGGTGCGCCAAGGTTCATTGCTCATGGCCTGCAGCTGTTCAGGGCGTATGCCCTGCCACGGCTGCCCAAGGTCGCAGTCGAGCGGGTCCCAGCCCGGCACCAGTGCTTTGGCCAAAGCCCGcagcag                                                                   >  2:2138120/1‑150 (MQ=255)
                                                                                 aGGCGCTCGCGGGTGTCGCCACAGGTGCGCCAAGGTTCATTGCTCATGGCCTGCAGCTGTTCAGGGCGTATGCCCTGCCACGGCTGCCCAAGGTCGCAGTCGAGc                                                                                                                <  1:812770/105‑1 (MQ=255)
                                                                                 aGGCGCTCGCGGGTGTCGCCACAGGTGCGCCAAGGTTCATTGCTCATGGCCTGCAGCTGTTCAGGGCGTATGCCCTGCCACGGCTGCCCAAGGTCGCAGTCGAGCGGGTCCCAGCCCGGCACCAGTGCTTTGGCCAAAGCCCGCAGCAgg                                                                  >  1:451379/1‑150 (MQ=255)
                                                                                  ggCGCTCGCGGGTGTCGCCACAGGTGCGCCAAGGTTCATTGCTCATGGCCTGCAGCTGTTCAGGGCGTATGCCCTGCCACGGCTGCCCAAGGTCGCAGTCGAGCGGGTCCCAGCCCGGCACCAGTGCTTTGTCCAAAGCCCGCAGCAGGc                                                                 >  1:1944190/1‑150 (MQ=255)
                                                                                   gcgcTCGCGGGTGTCGCCACAGGTGCGCCAAGGTTCATTGCTCATGGCCTGCAGCTGTTCAGGGCGTATGCCCTGCCACGGCTGCCCAAGGTCGCAGTCGAGCGGGTCCCAGCCCGGCACCAGTGCTTTGGCCAAAGCCCGCAGCAGGCt                                                                <  2:154393/150‑1 (MQ=255)
                                                                                    cgcTCGCGGGTGTCGCCACAGGTGCGCCAAGGTTCATTGCTCATGGCCTGCAGCTGTTCAGGGCGTATGCCCTGCCACGGCTGCCCAAGGTCGCAGTCGAGCGGGTCCCAGCCCGGCACCAGTGCTTTGGCCAAAGCCCGCAGCAGGCTg                                                               <  2:1336307/150‑1 (MQ=255)
                                                                                       tCGCGGGTGTCGCCACAGGTGCGCCAAGGTTCATTGCTCATGGCCTGCAGCTGTTCAGGGCGTATGCCCTGCCACGGCTGCCCAAGGTCGCAGTCGAGCGGGTCCCAGCCCGGCACCAGTGCTTTGGCCAAAGCCCGCAGCAGGCTGGCa                                                            >  1:1039297/1‑150 (MQ=255)
                                                                                         gcgGGTGTCGCCACAGGTGCGCCAAGGTTCATTGCTCATGGCCTGCAGCTGTTCAGGGCATACGCCCTgaa                                                                                                                                          >  1:156743/1‑69 (MQ=255)
                                                                                         gcgGGTGTCGCCACAGGTGCGCCAAGGTTCATTGCTCATGGCCTGCAGCTGTTCAGGGCATACGCCCTgaa                                                                                                                                          <  2:156743/71‑3 (MQ=255)
                                                                                                tCGCCACAGGTGCGCCAAGGTTCATTGCTCATGGCCTGCAGCTGTTCAGGGCGTATGCCCTGCCACGGCTGCCCAAGGTCGCAGTCGAGCGGGTCCCAGCCCGGCACCAGTGCTTTGGCCAAAGCCCGCAGCAGGCTGGCATTACCCCCa                                                   >  1:1346676/1‑150 (MQ=255)
                                                                                                     acaGGTGCGCCAAGGTTCATTGCTCATGGCCTGCAGCTGTTCAGGGCGTATGCCCTGCCACGGCTGCCCAAGGTCGCAGTCGAGCGGGTCCCAGCCCGGCACCAGTGCTTTGGCCAAAGCCCGCAGCAGGCTGGCATTACCCCCACGAcc                                              <  2:451379/150‑1 (MQ=255)
                                                                                                      caGGTGCGCCAAGGTTCATTGCTCATGGCCTGCAGCTGTTCAGGGCGTATGCCCTGCCACGGCTGCCCAAGGTCGCAGTCGAGCGGGTCCCAGCCCGGCACCAGTGCTTTGGCCAAAGCCCGCAGCAGGCTGGCATTACCCCCACGACCg                                             >  1:1337555/1‑150 (MQ=255)
                                                                                                       aGGTGCGCCAAGGTTCATTGCTCATGGCCTGCAGCTGTTCAGGGCGTATGCCCTGCCACGGCTGCCCAAGGTCGCAGTCGAGCGGGTCCCAGCCCGGCACCAGTGCTTTg                                                                                     <  1:1932773/110‑1 (MQ=255)
                                                                                                       aGGTGCGCCAAGGTTCATTGCTCATGGCCTGCAGCTGTTCAGGGCGTATGCCCTGCCACGGCTGCCCAAGGTCGCAGTCGAGCGGGTCCCAGCCCGGCACCAGTGCTTTg                                                                                     >  2:1932773/1‑110 (MQ=255)
                                                                                                       aGGTGCGCCAAGGTTCATTGCTCATGGCCTGCAGCGGTTCAGGGCGTATGCCCTGCCACGGCTGCCCAAGGTCGCAGTCGAGCGGGTCCCAGCCCGGCACCAGTGCTTTGGCCAAAGCCCGCAGCAGGCTGGCATTACCCCCACGACCGt                                            >  2:426701/1‑150 (MQ=255)
                                                                                                         gTGCGCCAAGGTTCATTGCTCATGGCCTGCAGCTGTTCAGGGCGTATGCCCTGCCACGGCTGCCCAAGGTCGCAGTCGAGCGGGTCCCAGCCCGGCACCAGTGCTTTGGCCAAAGCCCGCAGCAGGCTGGCATTACCCCCACGACCGTCg                                          >  1:467717/1‑150 (MQ=255)
                                                                                                             gcCAAGGTTCATTGCTCATGGCCTGCAGCTGTTCAGGGCGTATGCCCTGCCACGGCTGCCCAAGGTCTCAGTCGAGCGGGTCCCAGCCCGGCACCAGTGCTTTGGCCAAAGCCCGCAGCAGGCTGGCATTACCCCCACGACCGTCGCCAc                                      >  2:71430/1‑150 (MQ=255)
                                                                                                              ccAAGGTTCATTGCTCATGGCCTGCAGCTGTTCAGGGCGTATGCCCTGCCACGGCTGCCCAAGGTCGCAGTCGAGCGGGTCCCAGCCCGGCACCAGTGCTTTGGCCAAAGCCCGCAGCAGGCTGGCATTACCCCCACGACCGTCGCCACg                                     >  1:2023309/1‑150 (MQ=255)
                                                                                                                aaGGTTCATTGCTCATGGCCTGCAGCTGTTCAGGGCGTATGCCCTGCCACGGCTGCCCAAGGTCg                                                                                                                         <  2:181030/65‑1 (MQ=255)
                                                                                                                aaGGTTCATTGCTCATGGCCTGCAGCTGTTCAGGGCGTATGCCCTGCCACGGCTGCCCAAGGTCg                                                                                                                         >  1:181030/1‑65 (MQ=255)
                                                                                                                   gTTCATTGCTCATGGCCTGCAGCTGTTCAGGGCGTATGCCCTGCCACGGCTGCCCAAGGTCGCAGTCGAGCGGGTCCCAGCCCGGCACCAGTGCTTTGGCCAAAGCCCGCAGCAGGCTGGCATTACCCCCACGACCGTCGCCACGCTCGa                                >  1:973340/1‑150 (MQ=255)
                                                                                                                    ttCATTGCTCATGGCCTGCAGCTGTTCAGGGCGTATGCCCTGCCACGGCTGCCCAAGGTCGCAGTCGAGCGGGTCCCAGCCCGGCACCAGTGCTTTGGCCAAAGCCCGCAGCAGGCTGGCATTACCCCCACGACCGTCGCCACGCTCGAc                               <  2:1944190/150‑1 (MQ=255)
                                                                                                                    gttgttGCTCATGGCCTGCAGCTGTTCAGGGCGTATGCCCTGCCACGGCTGCCCAAGGTCGCAGTCGAGCGGGTCCCAGCCCGGCACCAGTGCTTTGGCCAAAGCCCGCAGCAGGCTGGCATTACCCCCACGACCGTCGCCACGCTCGAc                               >  2:1619157/5‑150 (MQ=255)
                                                                                                                       aTTGCTCATGGCCTGCAGCTGTTCAGGGCGTATGCCCTGCCACGGCTGCCCAAGGTCGCAGTCGAGCGGGTCCCAGCCCGGCACCAGTGCTTTGGCCAAAGCCCGCAGCAGGCTGGCATTACCCCCACGa                                                <  1:1520723/130‑1 (MQ=255)
                                                                                                                       aTTGCTCATGGCCTGCAGCTGTTCAGGGCGTATGCCCTGCCACGGCTGCCCAAGGTCGCAGTCGAGCGGGTCCCAGCCCGGCACCAGTGCTTTGGCCAAAGCCCGCAGCAGGCTGGCATTACCCCCACGa                                                >  2:1520723/1‑130 (MQ=255)
                                                                                                                        ttGCTAATGGCCTGCAGCTGTTCAGGGCGTATGCCCTGCCACGGCTGCCCAAGGTCGCAGTCGAGCGGGTCCCAGCCCGGCACCAGTGCTTTGGCCAAAGCCCGCAGCAGGCTGGCATTACCCCCACGACCGTCGCCACGCTCGACCCGc                           >  2:2441629/1‑150 (MQ=255)
                                                                                                                          gCTCATGGCCTGCAGCTGTTCAGGGCGTATGCCCTGCCACGGCTGCCCAAGGTCGCAGTCGAGCGGGTCCCAGCCCGGCACCAGTGCTTTGGCCAAAGCCCGCAGCAGGCTGGCATTACCCCCACGACCGTCGCCACGCTCGACCCGcaa                         <  2:1130073/150‑1 (MQ=255)
                                                                                                                           cTCATGGCCTGCAGCTGTTCAGGGCGTATGCCCTGCCACGGCTGCCCAAGGTCGCAGTCGAGCGGGTCCCAGCCCGGCACCAGTGCTTTGGCCAAAGCCCGCAGCAAGGCCGCATTACCCCCACGACCGTCGCCACGCTCGACCCGcaac                        >  1:23902/1‑150 (MQ=255)
                                                                                                                            tCATGGCCTGCAGCTGTTCAGGGCGTATGCCCTGCCACGGCTGCCCAAGGTCGCAGTCGAGCGGGTCCCAGCCCGGCACCAGTGCTTTGGCCAAAGCCCGCAGCAGGCTGGCATTACCCCCACGACCGTCGCCACGCTCGACCCGcaaca                       >  1:941013/1‑150 (MQ=255)
                                                                                                                            tCATGGCCTGCAGCTGTTCAGGGCGTATGCCCTGCCACGGCTGCCCAAGGTCGCAGTCGAGCGGGTCCCAGCCCGGCACCAGTGCTTTGGCCAAAGCCCGCAGCAGGCTGGCATTACCCCCACGACCGTCGCCACGCTCGACCCGcaaca                       >  1:494008/1‑150 (MQ=255)
                                                                                                                              aTGGCCTGCAGCTGTTCAGGGCGTATGCCCTGCCACGGCTGCCCAAGGTCGCAGTCGAGCGGGTCCCAGCCCGGCACCAGTGCTTTGGCCAAAGCCCGCAGCAGGCTgg                                                              <  1:2335003/109‑1 (MQ=255)
                                                                                                                              aTGGCCTGCAGCTGTTCAGGGCGTATGCCCTGCCACGGCTGCCCAAGGTCGCAGTCGAGCGGGTCCCAGCCCGGCACCAGTGCTTTGGCCAAAGCCCGCAGCAGGCTGGCATTACCCCCACGACCGTCGCCACGCTCGACCCGCaacaac                     <  1:1619157/150‑1 (MQ=255)
                                                                                                                               tGGCCTGCAGCTGTTCAGGGCGTATGCCCTGCCACGGCTGCCCAAGGTCGCAGTCGAGCGGGTCCCAGCCCGGCACCAGTGCTTTGGCCAAAGCCCGCAGCAGGCTGGCATTACCCCCACGACCGTCGCCACGCTCGACCCGCAACAACg                    >  1:771521/1‑150 (MQ=255)
                                                                                                                                   cTGCAGCTGTTCAGGGCGTATGCCCTGCCACGGCTGCCCAAGGTCGCAGTCGAGCGGGTCCCAGCCCGGCACCAGTGCTTTGGCCAAAGCCCGCAGCAGGCTGGCATTAGCCCCACGACCGTCGCCACGCTCGACCCGCAACAACGCcag                >  2:1484694/1‑150 (MQ=255)
                                                                                                                                    tGCAGCTGTTCAGGGCGTATGCCCTGCCACGGCTGCCCACGGTCGCAGTCGAGCGGGTCCCAGCCCGGCACCAGTGCTTTGGCCAAAGCCCGCAGCAGGCTGGCATTACCCCCACGACCGTCGCCACGCTCGACCCGCAACAACGcgacc               >  2:894679/1‑146 (MQ=255)
                                                                                                                                     gCAGCTGTTCAGGGCGTATGCCCTGCCACGGCTGCCCAAGGTCGCAGTCGAGCGGGTCCCAGCCCGGCACCAGTGCTTTGGCCAAAGCCCGCAGCAGGCTGGCATTACCCCCACGACCGTCGCCACGCTCGACCCGCAACAACGCcagca              <  1:265497/150‑1 (MQ=255)
                                                                                                                                            ttCAGGGCGTATGCCCTGCCACGGCTGCCCAAGGTCGCAGTCGAGCGGGTCCCAGCCCGGCACCAGTGCTTTGGCCAAAGCCCGCAGCAGGCTGGCATTACCCCCACGACCGTCGCCACGCTCGACCCGCAACAACGCCAGCAGCGTATc       >  2:219889/1‑150 (MQ=255)
                                                                                                                                              cAGGGCGTATGCCCTGCCACGGCTGCCCAAGGTCGCAGTCGAGCGGGTCCCAGCCCGGCACCAGTGCTTTGGCCAAAGCCCGCAGCAGGCTGGCATTACCCCCACGACCGTCGCCACGCTCGACCCGCAACAACGCCAGCAGCGTATCCa     <  2:2378387/150‑1 (MQ=255)
                                                                                                                                               aGGGCGTATGCCCTGCCACGGCTGCCCAAGGTCGCAGTCGAGCGGGTCCCAGCCCGGCACCAGTGCTTTGGCCAAAGCCCGCAGCAGGCTGGCATTACCCCCACGACCGTCGCCACGCTCGACCCGCAACAACGCCAGCAGCGTATCCAg    >  1:2011314/1‑150 (MQ=255)
                                                                                                                                                gggCGTATGCCCTGCCACGGCTGCCCAAGGTCGCAGTCGAGCGGGTCCCAGCCCGGCACCAGTGCTTTGGCCAAAGCCCGCAGCAGGCTGGCATTACCCCCACGACCGTCGCCACGCTCGACCCGCAACAACGCCAGCAGCGTATCCAgg   >  2:1620108/1‑150 (MQ=255)
                                                                                                                                                gggCGTATGCCCTGCCACGGCTGCCCAAGGTCGCAGTCGAGCGGGTCCCAGCCCGGCACCAGTGCTTTGGCCAAAGCCGAAGAGAAGCTGGCATTACCCCCACGACCGTCGCCACGCTCGACCCGCAACAACGCCAGCAGCGTATCCAgg   >  1:413075/1‑150 (MQ=255)
                                                                                                                                                 ggCGTATGCCCTGCCACGGCTGCCCAAGGTCGCAGTCGAGCGGGTCCCAGCCCGGCACCAGTGCTTTGGCCAAAGCCCGCAGCAAGCTGGCATTACCCCCACGACCGTCGCCACGCTCGACCCGCAACAACGCCAGCAGCGTATCCAGGc  >  2:2033844/1‑150 (MQ=255)
                                                                                                                                                    gTATGCCCTGCCACGGCTGCCCAAGGTCGCAGTCGAGCGGGTCCCAGCCCGGCACCAGTGCTTTGGCCaaa                                                                               >  1:835246/1‑71 (MQ=255)
                                                                                                                                                    gTATGCCCTGCCACGGCTGCCCAAGGTCGCAGTCGAGCGGGTCCCAGCCCGGCACCAGTGCTTTGGCCaaa                                                                               <  2:835246/71‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
ACCGGTTGCCACTCACTCAGGTCGGGTAACTGTGCCGTGCCCTCCAATGCCTGCCCGATCACCTGCAGTGCCAGCAGCTCAAGGCGCTCGCGGGTGTCGCCACAGGTGCGCCAAGGTTCATTGCTCATGGCCTGCAGCTGTTCAGGGCGTATGCCCTGCCACGGCTGCCCAAGGTCGCAGTCGAGCGGGTCCCAGCCCGGCACCAGTGCTTTGGCCAAAGCCCGCAGCAGGCTGGCATTACCCCCACGACCGTCGCCACGCTCGACCCGCAACAACGCCAGCAGCGTATCCAGGC  >  NC_002947/3978267‑3978561

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