Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,981,619 G→T 5.3% I326I (ATC→ATA cobN ← aerobic cobaltochelatase subunit CobN

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,981,6190GT 5.3% 718.9 / 19.6 208I326I (ATC→ATAcobNaerobic cobaltochelatase subunit CobN
Reads supporting (aligned to +/- strand):  ref base G (92/105);  new base T (3/8);  total (95/113)
Fisher's exact test for biased strand distribution p-value = 2.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 212 total aligned reads displayed.

TGGGCATTGGGCAGGCGGGCCAGCTCGACCCAGCGGCGGGCCAGCTCGGCGACAAAGTCCATGCGCTCGGGGTGGGCGCGGTAGCAGACAACGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCTGCAGCACCGGGATGTCGCGGCGAAACGGGCGCAGGTT  >  NC_002947/3981470‑3981768
                                                                                                                                                     |                                                                                                                                                     
  ggCATTGGGCAGGCGGGCCAGCTCGACCCAGCGGCGGGCCAGCTCGGCGACAAAGTCCATGCGCTCGGGGTGGGCGCGGTAGCAGACAACGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGAt                                                                                                                                                     >  1:3355386/1‑150 (MQ=255)
   gCATTGGGCAGGCGGGCCAGCTCGACCCAGCGGCGGGCCAGCTCGGCGACAAAGTCCATGCGCTCGGGGTGGGCGCGGTAGCAGACAACGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATa                                                                                                                                                    <  1:792915/150‑1 (MQ=255)
     aTTGGGCAGGCGGGCCAGCTCGACCCAGCGGCGGGCCAGCTCGGCGACAAAGTCCATGCGCTCGGGGTGGGCGCGGTAGCAGACAACGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAAt                                                                                                                                                  <  2:3338021/150‑1 (MQ=255)
         ggCAGGCGGGCCAGCTCGACCCAGCGGCGGGCCAGCTCGGCGACAAAGTCCATGCGCTCGGGGTGGGCGCGGTAGCAGACAACGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATgcgc                                                                                                                                              <  1:1967907/150‑1 (MQ=255)
           cAGGCGGGCCAGCTCGACCCAGCGGCGGGCCAGCTCGGCGACAAAGTCCATGCGCTCGGGGTGGGCGCGGTAGCAGACAACGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGccc                                                                                                                                            <  2:1454835/150‑1 (MQ=255)
            aGGCGGGCCAGCTCGACCCAGCGGCGGGCCAGCTCGGCGACAAAGTCCATGCGCTCGGGGTGGGCGCGGTAGCAGACAACGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCATTATCACGCGCCCg                                                                                                                                           <  2:2104585/150‑1 (MQ=255)
               cGGGCCAGCTCGACCCAGCGGCGGGCCAGCTCGGCGACAAAGTCCATGCGCTCGGGGTGGGCGCGGTAGCAGACAACGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTcc                                                                                                                                        <  1:1875385/150‑1 (MQ=255)
               cGGGCCAGCTCGACCCAGCGGCGGGCCAGCTCGGCGACAAAGTCCATGCGCTCGGGGTGGGCGCGGTAGCAGACAACGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCATTATCACGCGCCCggtc                                                                                                                                        <  2:1074263/150‑4 (MQ=255)
               cGGGCCAGCTCGACCCAGCGGCGGGCCAGCTCGGCGACAAAGTCCATGCGCTCGGGGTGGGCGCGGTAGCAGACAACGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCATTATCACGCGCCCggtc                                                                                                                                        <  2:1604181/150‑4 (MQ=255)
               cGGGCCAGCTCGACCCAGCGGCGGGCCAGCTCGGCGACAAAGTCCATGCGCTCGGGGTGGGCGCGGTAGCAGACAACGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCATTATCACGCGCCCggtc                                                                                                                                        <  2:3879688/150‑4 (MQ=255)
                gggCCAGCTCGACCCAGCGGCGGGCCAGCTCGGCGACAAAGTCCATGCGCTCGGGGTGGGCGCGGTAGCAGACAACGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCa                                                                                                                                       <  1:1319598/150‑1 (MQ=255)
                  gCCAGCACGCCCCAGCGGCGGGCCAGCTCGGCGACACCGTCCATGCGCTCGGGGTGGGCGCGGTAGCAGACACCGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGGACGGGCGCATTATCACGCGCCCGGTCAGc                                                                                                                                     <  2:3783983/150‑1 (MQ=255)
                   ccAGCTCGACCCAGCGGCGGGCCAGCTCGGCGACAAAGTCCATGCGCTCGGGGTGGGCGCGGTAGCAGACAACGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCt                                                                                                                                    >  1:3100348/1‑150 (MQ=255)
                   ccAGCTCGACCCAGCGGCGGGCCAGCTCGGCGACAAAGTCCATGCGCTCGGGGTGGGCGCGGTAGCAGACAACGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGGCCTTGTAGCTGACCGGGCGCGTGTTAATGCGCCCGTCCAGCt                                                                                                                                    >  2:3662473/1‑150 (MQ=255)
                    cAGCTCGACCCAGCGGCGGGCCAGCTCGGCGACAAAGTCCATGCGCTCGGGGTGGGCGCGGTAGCAGACAACGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTTCAGCtt                                                                                                                                   <  1:2870107/150‑2 (MQ=255)
                    cAGCTCGACCCAGCGGCGGGCCAGCTCGGCGACAAAGTCCATGCGCTCGGGGTGGGCGCGGTAGCAGACAACGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTTCAGCtt                                                                                                                                   <  1:2694210/150‑2 (MQ=255)
                       cTCGACCCAGCGGCGGGCCAGCTCGGCGACAAAGTCCATGCGCTCGGGGTGGGCGCGGTAGCAGACAACGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGGACGGGCGCATTATCACGCGCCCGGTCAGCttcaa                                                                                                                                <  2:1535197/150‑5 (MQ=38)
                       cTCGACCCAGCGGCGGGCCAGCTCGGCGACAAAGTCCATGCGCTCGGGGTGGGCGCGGTAGCAGACAACGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCgg                                                                                                                                <  1:3694708/150‑1 (MQ=255)
                        tCGACCCAGCGGCGGGCCAGCTCGGCGACAAAGTCCATGCGCTCGGGGTGGGCGCGGTAGCAGACAACGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGgc                                                                                                                               >  2:2932266/1‑150 (MQ=255)
                         cGACCCAGCGGCGGGCCAGCTCGGCGACAAAGTCCATGCGCTCGGGGTGGGCGCGGTAGCAGACAACGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGgca                                                                                                                              >  1:2902120/1‑150 (MQ=255)
                           aCCCAGCGGCGGGCCAGCTCGGCGACAAAGTCCATGCGCTCGGGGTGGGCGCGGTAGCAGACAACGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGgcaat                                                                                                                            <  2:4095480/150‑1 (MQ=255)
                            cccAGCGGCGGGCCAGCTCGGCGACAAAGTCCATGCGCTCGGGGTGGGCGCGGTAGCAGACAACGTCGGACCGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGgcaatg                                                                                                                           >  1:4320986/1‑150 (MQ=255)
                              cAGCGGCGGGCCAGCTCGGCGACAAAGTCCATGCGCTCGGGGTGGGCGCGGTAGCAGACTACGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGgcaatgca                                                                                                                         <  2:3124631/150‑1 (MQ=255)
                               aGCGGCGGGCCAGCTCGGCGACAAAGTCCATGCGCTCGGGGTGGGCGCGGTAGCAGACAACGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGgcaatgcaa                                                                                                                        <  2:23522/150‑1 (MQ=255)
                                  ggcggGCCAGCTCGGCGACCAAGTCCATGCGCTCGGGGTGGGCGCGGTAGCAGACAACGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAAtgt                                                                                                                     <  1:655586/150‑1 (MQ=255)
                                    cggGCCAGCTCGGCGACAAAGTCCATGCGCTCGGGGTGGGCGCGGTAGCAGACAACGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGc                                                                                                                   <  1:2249032/150‑1 (MQ=255)
                                      ggCCAGATCGGCGACAAAGTCCATGCGCTCGGGGTTGGAGCGGTAGAATACAAAGTCTTAGTTGCTGCGGGCTCTGGGGCAGGAGATTTCCTTGAAGCTGAGAGGTCGCGTGATAATGGGACCGTCGAGCTTCGGGGATGGTATGTGGAt                                                                                                                 >  2:69986/1‑150 (MQ=255)
                                       gCCAGCTCGGCGACAAAGTCCATGCGCTCGGGGTGGGCGCGGTAGCAGACAACGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATg                                                                                                                >  1:1633774/1‑150 (MQ=255)
                                         cAGCTCGGCGACAAAGTCCATGCGCTCGGGGTGGGCGCGGTAGCAGACAACGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGc                                                                                                              >  2:1400526/1‑150 (MQ=255)
                                           gCTCGGCGACAAAGTCCATGCGCTCGGGGTGGGCGCGGTAGCAGACAACGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCCATGCACATTGCATTGccg                                                                                                            <  1:486361/150‑2 (MQ=255)
                                           gCTCGGCGACAAAGTCCATGCGCTCGGGGTGGGCGCGGTAGCAGACAACGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCATGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCCATGCACATTGCATTGccg                                                                                                            <  2:1660425/150‑2 (MQ=255)
                                            cTCGGCGACAAAGTCCATGCGCTCGGGGTGGGCGCGGTAGCAGACAACGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCCATGCACATTGCATTGccgg                                                                                                           <  2:1942891/150‑3 (MQ=255)
                                              cGGCGACAAAGTCCATGCGCTCGGGGTGGGCGCGGTAGCAGACAACGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGt                                                                                                         <  2:1349779/150‑1 (MQ=255)
                                               ggCGACAAAGTCCATGCGCTCGGGGTGGGCGCGGTAGCAGACAACGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGACCAGGTc                                                                                                        >  1:1198230/1‑150 (MQ=255)
                                                 cGACAAAGTCCATGCGCTCGGGGTGGGCGCGGTAGCAGACAACGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTcgc                                                                                                      >  2:2479214/1‑150 (MQ=255)
                                                   aCAAAGTCCATGCGCTCGGGGTGGGCGCGGTAGCAGACAACGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCTGGgcgcgc                                                                                                    <  1:4312081/150‑1 (MQ=255)
                                                   aCAAAGTCCATGCGCTCGGGGTGGGCGCGGTAGCAGACAACGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTcgcgc                                                                                                    <  1:435661/150‑1 (MQ=255)
                                                   aCAAAGTCCATGCGCTCGGGGTGGGCGCGGTAGCAGACAACGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTcgcgc                                                                                                    >  1:255473/1‑150 (MQ=255)
                                                     aaaGTCCATGCGCTCGGGGTGGGCGCGGTAGCAGACAACGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTcgcgcgc                                                                                                  >  1:307360/1‑150 (MQ=255)
                                                          ccATGCGCTCGGGGTGGGCGCGGTAGCAGACAACGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCa                                                                                             <  2:2838250/150‑1 (MQ=255)
                                                            aTGCGCTCGGGGTGGGCGCGGTAGCAGACAACGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAgg                                                                                           <  1:252324/150‑1 (MQ=255)
                                                             tGCGCTCGGGGTGGGCGCGGTAGCAGACAACGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGc                                                                                          >  2:1235452/1‑150 (MQ=255)
                                                              gcgcTCGGGGTGGGCGCGGTAGCAGACAACGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGcc                                                                                         <  2:3022219/150‑1 (MQ=255)
                                                                gcTCGGGGTGGGCGCGGTAGCAGACAACGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCTGGGCGCGCGCGACCTGGccat                                                                                       <  2:2643405/150‑3 (MQ=255)
                                                                 cTCGGGGTGGGCGCGGTAGCAGACAACGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTg                                                                                      >  1:1163551/1‑150 (MQ=255)
                                                                  tCGGGGTGGGCGCGGTAGCAGACAACGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGc                                                                                     >  1:4044018/1‑150 (MQ=255)
                                                                     gggTGGGCGCGGTAGCAGACAACGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCg                                                                                  <  2:845031/150‑1 (MQ=255)
                                                                     gggTGGGCGCGGTAGCAGACAACGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCg                                                                                  <  1:3885840/150‑1 (MQ=255)
                                                                       gTGGGCGCGGTAGCAGACAACGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGc                                                                                     <  1:3747426/145‑1 (MQ=255)
                                                                       gTGGGCGCGGTAGCAGACAACGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGc                                                                                     >  2:3747426/1‑145 (MQ=255)
                                                                       gTGGGCGCGGTAGCAGACAACGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCt                                                                                <  2:3828830/150‑1 (MQ=255)
                                                                        ggggCGCGGTAGCAGACAACGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCACTGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTg                                                                               <  1:1334409/149‑1 (MQ=255)
                                                                         gggCGCGGTAGCAGACAACGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATCCTGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTgg                                                                              <  2:4143216/150‑1 (MQ=255)
                                                                         gggCGCGGTAGCAGACAACGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTgg                                                                              <  2:1716668/150‑1 (MQ=255)
                                                                              cgGTAGCAGACAACGTCGGACTGGCTGTGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGccgaa                                                                         <  2:1162195/150‑5 (MQ=255)
                                                                                 tAGCAGACAACGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCccag                                                                      <  2:1473266/150‑1 (MQ=255)
                                                                                  aGCAGACAACGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCccagc                                                                     <  1:3864836/150‑1 (MQ=255)
                                                                                    cAGACAACGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCccagcca                                                                   <  2:2271128/150‑1 (MQ=255)
                                                                                    cAGACAACGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCccagcca                                                                   <  2:1843639/150‑1 (MQ=255)
                                                                                    cAGACAACGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCccagcca                                                                   <  1:2443158/150‑1 (MQ=255)
                                                                                      gACAACGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTACATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAgg                                                                 >  1:1536834/1‑150 (MQ=255)
                                                                                        cAACGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCt                                                               >  2:3237189/1‑150 (MQ=255)
                                                                                        cAACGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCt                                                               <  2:2975554/150‑1 (MQ=255)
                                                                                          aCGTCGGACTTGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAAGGCAATGTGCATGGCCAGGTCGCGCGCTCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTgg                                                             <  1:445604/150‑1 (MQ=255)
                                                                                            gTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGtt                                                           >  2:589303/1‑150 (MQ=255)
                                                                                             tCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATACTGCGCCCGTCCAGCTCCGGCCATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCGTCCCAGCctggctgggaa                                                          <  2:2447997/150‑4 (MQ=255)
                                                                                               ggACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTc                                                        <  1:3570990/150‑1 (MQ=255)
                                                                                                gACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTcc                                                       <  2:1197009/150‑1 (MQ=255)
                                                                                                gACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTcc                                                       <  1:4494937/150‑1 (MQ=255)
                                                                                                   tGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTgc                                                    <  2:3696414/150‑1 (MQ=255)
                                                                                                     gCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTgcgc                                                  >  1:2733789/1‑150 (MQ=255)
                                                                                                      cTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGca                                                 >  2:23857/1‑150 (MQ=255)
                                                                                                       tGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGcac                                                >  2:1569249/1‑150 (MQ=255)
                                                                                                         cgctcgctGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAg                                              <  1:4417301/150‑1 (MQ=255)
                                                                                                                tGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCt                                       >  2:225647/1‑150 (MQ=255)
                                                                                                                  cgcCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCTgc                                     <  2:2684152/150‑1 (MQ=255)
                                                                                                                   gcCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCTgca                                    >  2:1055846/1‑150 (MQ=255)
                                                                                                                   gcCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCTgca                                    <  1:1789556/150‑1 (MQ=255)
                                                                                                                   gcCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCTgca                                    <  2:3187314/150‑1 (MQ=255)
                                                                                                                     cAGGCCCTGTCCTTGAAGCTCACCGGGCGCGTGAGAATGCGCCCGTCCAGCTCCGGCACTGCAATGTGCATGGCCCGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGGGGTCCTGCGCCCCGATGGCCGgccgc                                  >  2:2008751/1‑150 (MQ=255)
                                                                                                                          cATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCTGCAGCACCgg                             >  1:4220555/1‑150 (MQ=255)
                                                                                                                          cATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCTGCAGCACCgg                             >  2:351197/1‑150 (MQ=255)
                                                                                                                          cATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCTGCAGCACCgg                             <  2:645880/150‑1 (MQ=255)
                                                                                                                          cATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCTGCAGCACCgg                             <  1:3237189/150‑1 (MQ=255)
                                                                                                                            tGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCACTGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCTGCAGCACCGGGa                           >  2:4271435/1‑150 (MQ=255)
                                                                                                                            tGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCTGCAGCACCGGGa                           >  1:1271498/1‑150 (MQ=255)
                                                                                                                            tGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTAGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCTGCAGCACCGGGa                           >  1:3250119/1‑150 (MQ=255)
                                                                                                                             gTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCTGCAGCACCGGGAt                          <  2:2225645/150‑1 (MQ=255)
                                                                                                                              tCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCTGCAGCACCGGGATg                         >  2:4354422/1‑150 (MQ=255)
                                                                                                                               ccTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCTGCAGCACCGGGATGt                        >  1:1121955/1‑150 (MQ=255)
                                                                                                                               ccTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCTGCAGCACCGGGATGt                        >  1:3716059/1‑150 (MQ=255)
                                                                                                                                cTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCTGCAGCACCGGGATGTc                       >  2:3571402/1‑150 (MQ=255)
                                                                                                                                cTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCTGCAGCACCGGGATGTc                       <  1:2115013/150‑1 (MQ=255)
                                                                                                                                cTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCATGCGCACAGATGGCCTGCCGCACCGGGATGTc                       >  1:1396033/1‑150 (MQ=255)
                                                                                                                                 ttGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCTGCAGCACCGGGATGTcg                      >  2:3861487/1‑150 (MQ=255)
                                                                                                                                 ttGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCTGCAGCACCGGGATGTcg                      >  1:2577640/1‑150 (MQ=255)
                                                                                                                                        tggacGGGCGCATTATCACGCGCCCGGTCAGCTTCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCTGCAGCACCGGGATGTCGCGGCGaa               >  2:443168/5‑150 (MQ=255)
                                                                                                                                        tGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCTGCAGCACCGGGATGTCGCGGCGaa               >  1:2999557/1‑150 (MQ=255)
                                                                                                                                          aCCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCTGCAGCACCGGGATGTCGCGGCGAAAc             <  1:590572/150‑1 (MQ=255)
                                                                                                                                                     gATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCTGCAGCACCGGGATGTCGCGGCGAAACGGGCGCAGGtt  >  1:466582/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TGGGCATTGGGCAGGCGGGCCAGCTCGACCCAGCGGCGGGCCAGCTCGGCGACAAAGTCCATGCGCTCGGGGTGGGCGCGGTAGCAGACAACGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCTGCAGCACCGGGATGTCGCGGCGAAACGGGCGCAGGTT  >  NC_002947/3981470‑3981768

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