Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,980,198 T→A 6.6% Q159H (CAA→CAT cobN ← aerobic cobaltochelatase subunit CobN

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,980,1980TA 6.6% 242.8 / 6.7 77Q159H (CAA→CATcobNaerobic cobaltochelatase subunit CobN
Reads supporting (aligned to +/- strand):  ref base T (34/37);  new base A (2/3);  total (37/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.42e-01

CGGTAGAACAGCAGCGGTGCCACTGGCTGCTCGGGGTGCCAGTGTGGGTACCAGTCCTCGAGTGTCGCGCTGCCCTTGGCGGGGTGATAGACCGCCGTGCGCGGCAAAGGCTGCGGCTCGTCCCAGGCATAGTCGCGGCCCAGCCATTGGCTGGCCAGGCAGTTGAACAGGTTGATGGCGTTGGCCTTGCCGCCCTGGCGCAGGTAATGCCAAAGCCGCTCGGCCTGCTCACCGCGGACGTTGCCCAGGCTGGTCAGCTCCGGGTCCGGGCGGTCATCGCCCGGCACCAGAATCAG  >  NC_002947/3980052‑3980347
                                                                                                                                                  |                                                                                                                                                      
cGGTAGAACAGCAGCGGTGCCACTGGCTGCTCGGGGTGCCAGTGTGGGTACCAGTCCTCGAGTGTCGCGCTGCCCTTGGCGGGGTGATAGACCGCCGTGCGCGGCAAAGGCTGCGGCTCGTCCCAGGCATAGTCGCGGCCCAGCCATtgg                                                                                                                                                     <  1:686729/150‑1 (MQ=255)
      aaCAGCAGCGGTGCCACTGGCTGCTCGGGGTGCCAGTGTGGGTACCAGTCCTCGAGTGTCGCGCTGCCCTTGGCGGGGTGATAGACCGCCGTGCGCGGCAAAGGCTGCGGCTCGTCCCAGGCATAGTCGCGGCCCAGCCATTGGCTGGcc                                                                                                                                              >  1:1424937/1‑150 (MQ=255)
      aaCAGCAGCGGTGCCACTGGCTGCTCGGGGTGCCAGTGTGGGTACCAGTCCTCGAGTGTCGCGCTGCCCTTGGCGGGGTGATAGACCGCCGTGCGCGGCAAAGGCTGCGGCTCGTCCCAGGCATAGTCGCGGCCCAGCCAAtggctgggc                                                                                                                                              <  1:1257028/150‑3 (MQ=255)
      aaCAGCAGCGGTGCCACTGGCTGCTCGGGGTGCCAGTGTGGGTACCAGTCCTCGAGTGTCGCGCTGCCATTGGCGGAGTGATAGACCGCCTGGCGCGGCAAAGGCTGCGGCTCGTCCCAGGCATAGTCGCGGCCCAGCCAAtggctgggc                                                                                                                                              <  1:884225/150‑3 (MQ=255)
         agcagcGGTGCCACTGGCTGCTCGGGGTGCCAGTGTGGGTACCAGTCCTCGAGTGTCGCGCTGCCCTTGGCGGGGTGATAGACCGCCGTGCGCGGCAAAGGCTGCGGCTCGTCCCAGGCATAGTCGCGGCCCAGCCATTGGCTGGCCAgg                                                                                                                                           >  2:440061/1‑150 (MQ=255)
           cagcGGTGCCACTGGCTGCTCGGGGTGCCAGTGTGGGTACCAGTCCTCGAGTGTCGCGCTGCCCTTGGCGGGGTGATAGACCGCCGTGCGCGGCAAAGGCTGCGGCTCGTCCCAGGCATAGTCGCGGCCCAGCCATTGGCTGGCCAGGCa                                                                                                                                         <  2:927929/150‑1 (MQ=255)
                gTGCCACTGGCTGCTCGGGGTGCCAGTGTGGGTACCAGTCCTCGAGTGTCGCGCTGCCCTTGGCGGGGTGATAGACCGCCGTGCGCGGCAAAGGCTGCGGCTCGTCCCAGGCATAGTCGCGGCCCAGCCATTGGCTGGCCAGGCAGTTGa                                                                                                                                    >  2:827072/1‑150 (MQ=255)
                   ccACTGGCTGCTCGGGGTGCCAGTGTGGGTACCAGTCCTCGAGTGTCGCGCTGCCCTTGGCGGGGTGATAGACCGCCGTGCGCGGCAAAGGCTGCGGCTCGTCCCAGGCATAGTCGCGGCCCAGCCATTGGCTGGCCAGGCAGTTGAACa                                                                                                                                 >  2:635562/1‑150 (MQ=255)
                       tGGCTGCTCGGGGTGCCAGTGTGGGTACCAGTCCTCGAGTGTCGCGCTGCCCTTGGCGGGGTGATAGACCGCCGTGCGCGGCAAAGGCTGCGGCTCGTCCCAGGCATAGTCGCGGCCCAGCCATTGGCTGGCCAGGCAGTTGAACAg                                                                                                                                >  2:855752/1‑147 (MQ=255)
                       tGGCTGCTCGGGGTGCCAGTGTGGGTACCAGTCCTCGAGTGTCGCGCTGCCCTTGGCGGGGTGATAGACCGCCGTGCGCGGCAAAGGCTGCGGCTCGTCCCAGGCATAGTCGCGGCCCAGCCATTGGCTGGCCAGGCAGTTGAACAg                                                                                                                                <  1:855752/147‑1 (MQ=255)
                           tgctCGGGGTGCCAGTGTGGGTACCAGTCCTCGAGTGTCGCGCTGCCCTTGGCGGGGTGATAGACCGCCGTGCGCGGCAAAGGCTGCGGCTCGTCCCAGGCATAGTCGCGGCCCAGCCATTGGCTGGCCAGGCAGTTGAACAGGTTGATg                                                                                                                         >  1:1294234/1‑150 (MQ=255)
                                ggggTGCCAGTGTGGGTACCAGTCCTCGAGTGTCGCGCTGCCCTTGGCGGGGTGATAGACCGCCGTGCGCGGCAAAGGCTGCGGCTCGTCCCAGGCATAGTCGCGGCCCAGCCATTGGCTGGCCAGGCAGTTGAACAGGTTGATGGCGtt                                                                                                                    <  1:1025585/150‑1 (MQ=255)
                                   gTGCCAGTGTGGGTACCAGTCCTCGAGTGTCGCGCTGCCCTTGGCGGGGTGATAGACCGCCGTGCGCGGCAAAGGCTGCGGCTCGTCCCAGGCATAGTCGCGGCCCAGCCATTGGCTGGCCAGGCCGTTGAACAGGTTGATGGCGTTGGc                                                                                                                 >  1:778083/1‑150 (MQ=255)
                                   gTGCCAGTGTGGGTACCAGTCCTCGAGTGTCGCGCTGCCCTTGGCGGGGTGATAGACCGCCGTGCGCGGCAAAGGCTGCGGCTCGTCCCAGGCATAGTCGCGGCCCAGCCATTGGCTGGCCAGGCAGTTGAACAGGTTGATGGCGTTGGc                                                                                                                 >  2:64509/1‑150 (MQ=255)
                                   gTGCCAGTGTGGGTACCAGTCCTCGAGTGTCGCGCTGCCCTTGGCGGGGTGATAGACCGCCGTGCGCGGCAAAGGCTGCGGCTCGTCCCAGGCATAGTCGCGGCCCAGCCATTGGCTGGCCAGGCAGTTGAACAGGTTGATGGCGTTGGc                                                                                                                 <  1:827072/150‑1 (MQ=255)
                                       cattgtgGGTACCAGTCCTCGAGTGTCGCGCTGCCCTTGGCGGGGTGATAGACCGTCGTGCGCGGCAAAGGCTGCGGCTCGTCCCAGGCATAGTCGCGGCCCAGCCATTGGCTGGCCAGGCAGTTGAACAGGTTGATGGCGTTGGCCTTg                                                                                                             <  1:1427577/147‑1 (MQ=255)
                                        aGTGTGGGTACCAGTCCTCGAGTGTCGCGCTGCCCTTGGCGGGGTGATAGACCGCCGTGCGCGGCAAAGGCTGCGGCTCGTCCCAGGCATAGTCGCGGCCCAGCCATTGGCTGGCCAGGCAGTTGAACAGGTTGATGGCGTTGGCCTTgc                                                                                                            <  1:1318002/150‑1 (MQ=255)
                                        aGTGTGGGTACCAGTCCTCGAGTGTCGCGCTGCCCTTGGCGGGGTGATAGACCGCCGTGCGCGGCAAAGGCTGCGGCTCGTCCCAGGCATAGTCGCGGCCACGCCTCAGGGTGAGCGGTTGCGTTAGCAGGTTGATGGCGTTGGCCTTgc                                                                                                            >  1:1384672/1‑150 (MQ=255)
                                               gTACCAGTCCTCGAGTGTCGCGCTGCCCTTGGCGGGGTGATAGACCGCCGTGCGCGGCAAAGGCTGCGGCTCGTCCCAGGCATAGTCGCGGCCCAGCCATTGGCTGGCCAGGCAGTTGAACAGGTTGATGGCGTTGGCCTTGCCGCCCTg                                                                                                     >  1:1067823/1‑150 (MQ=255)
                                                 aCCAGTCCTCGAGTGTCGCGCTGCCCTTGGCGGGGTGATAGACCGCCGTGCGCGGCAAAGGCTGCGGCTCGTCCCAGGCATAGTCGCGGCCCAGCCATTGGCTGGCCAGGCAGTTGAACAGGTTGATGGCGTTGGCCTTGCCGCCCTGgc                                                                                                   >  2:274220/1‑150 (MQ=255)
                                                    aGTCCTCGAGTGTCGCGCTGCCCTTGGCGGGGTGATAGACCGCCGTGCGCGGCAAAGGCTGCGGCTCGTCCCAGGCATAGTCGCGGCCCAGCCATTGGCTGGCCAGGCAGTTGAACAGGTTGATGGCGTTGGCCTTGCCGCCCTGGCGCa                                                                                                >  2:453702/1‑150 (MQ=255)
                                                    aGTCCTCGAGTGTCGCGCTGCCCTTGGCGGGGTGATAGACCGCCGTGCGCGGCAAAGGCTGCGGCTCGTCCCAGGCATAGTCGCGGCCCAGCCATTGGCTGGCCAGGCAGTTGAACAGGTTGATGGCGTTGGCCTTGCCGCCCTGGCGCa                                                                                                >  1:667692/1‑150 (MQ=255)
                                                     gTCCTCGAGTGTCGCGCTGCCCTTGGCGGGGTGATAGACCGCCGTGCGCGGCAAAGGCTGCGGCTCGTCCCAGGCATAGTCGCGGCCCAGCCATTGGCTGGCCAGGCAGTTGAACAGGTTGATGGCGTTGGCCTTGCCGCCCTGGCGCAg                                                                                               <  2:778083/150‑1 (MQ=255)
                                                      tCCTCGAGTGTCGCGCTGCCCTTGGCGGGGTGATAGACCGCCGTGCGCGGCAAAGGCTGCGGCTCGTCCCAGGCATAGTCGCGGCCCAGCCATTGGCTGGCCAGGCAGTTGAACAGGTTGATGGCGTTGGCCTTGCCGCCCTGGCGCAgg                                                                                              <  1:321628/150‑1 (MQ=255)
                                                        cTCGAGTGTCGCGCTGCCCTTGGCGGGGTGATAGACCGCCGTGCGCGGCAAAGGCTGCGGCTCGTCCCAGGCATAGTCGCGGCCCAGTCATTGGCTGGCCAGGCAGTTGAACAGGTTGATGGCGTTGGCCTTGCCGCCCTGGCGCAGGTa                                                                                            <  1:860940/150‑1 (MQ=255)
                                                        cTCGAGTGTCGCGCTGCCCTTGGCGGGGTGATAGACCGCCGTGCGCGGCAAAGGCTGCGGCTCGTCCCAGGCATAGTCGCGGCCCAGCCATTGGCTGGCCAGGCAGTTGAACAGGTTGATGGCGTTGGCCTTGCCGCCCTGGCGCGGGTa                                                                                            >  1:769071/1‑150 (MQ=255)
                                                         tCGAGTGTCGCGCTGCCCTTGGCGGGGTGATAGACCGCCGTGCGCGGCAAAGGCTGCGGCTCGTCCCAGGCATAGTCGCGGCCCAGCCATTGGCTGGCCAGGCAGTTGaa                                                                                                                                   >  2:297858/1‑110 (MQ=255)
                                                         tCGAGTGTCGCGCTGCCCTTGGCGGGGTGATAGACCGCCGTGCGCGGCAAAGGCTGCGGCTCGTCCCAGGCATAGTCGCGGCCCAGCCATTGGCTGGCCAGGCAGTTGaa                                                                                                                                   <  1:297858/110‑1 (MQ=255)
                                                         tCGAGTGTCGCGCTGCCCTTGGCGGGGTGATAGACCGCCGTGCGCGGCAAAGGCTGCGGCTCGTCCCAGGCATAGTCGCGGCCCAGCCATTGGCTGGCCAGGCAGTTGAACAGGTTGATGGCGTTgg                                                                                                                  <  2:920714/127‑1 (MQ=255)
                                                         tCGAGTGTCGCGCTGCCCTTGGCGGGGTGATAGACCGCCGTGCGCGGCAAAGGCTGCGGCTCGTCCCAGGCATAGTCGCGGCCCAGCCATTGGCTGGCCAGGCAGTTGAACAGGTTGATGGCGTTgg                                                                                                                  >  1:920714/1‑127 (MQ=255)
                                                          cGAGTGTCGCGCTGCCCTTGGCGGGGTGATAGACCGCCGTGCGCGGCAAAGGCTGCGGCTCGTCCCAGGCATAGTCGCGGCCCAGCCATTGGCTGGCCAGGCAGTTGAACAGGTTGATGGCGTTGGCCTTGCCGCCCTGGCGCAGGTAAt                                                                                          <  1:858267/150‑1 (MQ=255)
                                                          cGAGTGTCGCGCTGCCCTTGGCGGGGTGATAGACCGCCGTGCGCGGCAAAGGCTGCGGCTCGTCCCAGGCATAGTCGCGGCCCAGCCATTGGCTGGCCAGGCAGTTGAACAGGTTGATGGCGTTGGCCTTGCCGCCCTGGCGCAGGTAAt                                                                                          <  2:1067823/150‑1 (MQ=255)
                                                          cGAGTGTCGCGCTGCCCTTGGCGGGGTGATAGACCGCCGTGCGCGGCAAAGGCTGCGGCTCGTCCCAGGCATAGTCGCGGCCCAGCCATTGGCTGGCCAGGCAGTTGAACAGGTTGATGGCGTTGGCCTTGCCGCCCTGGCCCAGGTAAt                                                                                          <  1:904957/150‑1 (MQ=255)
                                                           gAGTGTCGCGCTGCCCTTGGCGGGGTGATAGACCGCCGTGCGCGGCAAAGGCTGCGGCTCGTCCCAGGCATAGTCGCGGCCCAGCCATTGGCTGGCCAGGCAGTTGAACAGGTTGATGGCGTTGGCCTTGCCGCCCTGGCGCAGGTAATg                                                                                         <  2:591317/150‑1 (MQ=255)
                                                            aGTGTCGCGCTGCCCTTGGCGGGGTGATAGACCGCCGTGCGCGGCAAAGGCTGCGGCTCGTCCCAGGCATAGTCGCGGCCCAGCCATTGGCTGGCCAGGCAGTTGAACAGGTTGATGGCGTTGGCCTTGCCGCCCTGGCGCAGGTAATGc                                                                                        <  1:453702/150‑1 (MQ=255)
                                                               gtCGCGCTGCCCTTGGCGGGGTGATAGACCGCCGTGCGCGGCAAAGGCTGCGGCTCGTCCCAGGCATAGTCGCGGCCCAGCCATTGGCTGGCCAGGCAGTTGAACAGGTTGATGGCGTTGGCCTTGCCGCCCTGGCGCAGGTAATGCCaa                                                                                     >  1:46862/1‑150 (MQ=255)
                                                                tCGCGCTGCCCTTGGCGGGGTGATAGACCGCCGTGCGCGGCAAAGGCTGCGGCTCGTCCCAGGCATAGTCGCGGCCCAGCCATTGGCTGGCCAGGc                                                                                                                                          <  1:310983/96‑1 (MQ=255)
                                                                tCGCGCTGCCCTTGGCGGGGTGATAGACCGCCGTGCGCGGCAAAGGCTGCGGCTCGTCCCAGGCATAGTCGCGGCCCAGCCATTGGCTGGCCAGGc                                                                                                                                          >  2:310983/1‑96 (MQ=255)
                                                                tCGCGCTGCCCTTGGCGGGGTGATAGACCGCCGTGCGCGGCAAAGGCTGCGGCTCGTCCCAGGCATAGTCGCGGCCCAGCCATTGGCTGGCCAGGCAGTTGAACAGGTTGATGGCGTTGGCCTTGCCGCCCTGGCGCAGGTAATGCCaaa                                                                                    <  2:622195/150‑1 (MQ=255)
                                                                  gcgcTGCCCTTGGCGGGGTGATAGACCGCCGTGCGCGGCAAAGGCTGCGGCTCGTCCCAGGCATAGTCGCGGCCCAGCCATTGGCTGGCCAGGCAGTTGAACAGGTTGATGGCGTTGGCCTTGCCGCCCTGGCGCGGGTAATGCCAAAGc                                                                                  <  2:769071/150‑1 (MQ=255)
                                                                       gCCCTTGGCGGGGTGATAGACCGCCGTGCGCGGCAAAGGCTGCGGCTCGTCCCAGGCATAGTCGCGGCCCAGCCATTGGCTGGCCAGGCAGTTGAACAGGTTGATGGCGTTGGCCTTGCCGCCCTGGCGCAGGTAATGCCAAAGCCGCTc                                                                             <  1:82106/150‑1 (MQ=255)
                                                                        cccTTGGCGGGGTGATAGACCGCCGTGCGCGGCAAAGGCTGCGGCTCGTCCCAGGCATAGTCGCGGCCCAGCCATTGGCTGGCCAGGCAGTTGaa                                                                                                                                   >  1:633192/1‑95 (MQ=255)
                                                                        cccTTGGCGGGGTGATAGACCGCCGTGCGCGGCAAAGGCTGCGGCTCGTCCCAGGCATAGTCGCGGCCCAGCCATTGGCTGGCCAGGCAGTTGaa                                                                                                                                   <  2:633192/95‑1 (MQ=255)
                                                                         ccTTGGCGGGGTGATAGACCGCCGTGCGCGGCAAAGGCTGCGGCTCGTCCCAGGCATAGTCGCGGCCCAGCCATTGGCTGGCCAGGCAGTTGAACAGGTTGATGGCGTTGGCCTTGCCGCCCTGGCGCAgg                                                                                              >  1:1197764/1‑131 (MQ=255)
                                                                         ccTTGGCGGGGTGATAGACCGCCGTGCGCGGCAAAGGCTGCGGCTCGTCCCAGGCATAGTCGCGGCCCAGCCATTGGCTGGCCAGGCAGTTGAACAGGTTGATGGCGTTGGCCTTGCCGCCCTGGCGCAgg                                                                                              <  2:1197764/131‑1 (MQ=255)
                                                                              gCGGGGTGATAGACCGCCGTGCGCGGCAAAGGCTGCGGCTCGTCCCAGGCATAGTCGCGGCCCAGCCATTGGCTGGCCAGGCAGTTGAACAGGTTGATGGCGTTGGCCTTGCCGCCCTGGCGCa                                                                                                <  1:1282636/124‑1 (MQ=255)
                                                                              gCGGGGTGATAGACCGCCGTGCGCGGCAAAGGCTGCGGCTCGTCCCAGGCATAGTCGCGGCCCAGCCATTGGCTGGCCAGGCAGTTGAACAGGTTGATGGCGTTGGCCTTGCCGCCCTGGCGCa                                                                                                >  2:1282636/1‑124 (MQ=255)
                                                                                 gggTGATAGACCGCCGTGCGCGGCAAAGGCTGCGGCTCGTCCCAGGCATAGTCGCGGCCCAGCCATTGGCTGGCCAGGCAGTTGAACAGGTTGATGGCGTTGGCCTTGCCGCCCTGGCGCAGGTAATGCCAAAGCCGCTCGGCCTGCTCa                                                                   >  1:619285/1‑150 (MQ=255)
                                                                                   gTGATAGACCGCCGTGCGCGGCAAAGTCTGCGGCTCGTCCTAGGCATAGTCGCGGCCCAGCCATTGGCTGGCCAGGCAGTTGAACAGGTTGATGGCGTTGGCCTTGCCGCCCTGGCGCAGGTAATGCCAAAGCCGCTCGTCCTGCTCAcc                                                                 <  2:1361631/150‑1 (MQ=255)
                                                                                   gTGATAGACCGCCGTGCGCGGCAAAGGCTGCGGCTCGTCCCAGGCATAGTCGCGGCCCAGCCATTGGCTGGCCAGGCAGTTGAACAGGTTGATGGCGTTGGCCTTGCCGCCCTGGCGCAGGTAATGCCAAAGCCGCTCGGCCTGCTCAcc                                                                 >  2:1045530/1‑150 (MQ=255)
                                                                                   gTGATAGACCGCCGTGCGCGGCAAAGGCTGCGGCTCGTCCCAGGCATAGTCGCGGCCCAGACATTGGCTGGCCAGGCAGTTGAACAGGTTGATGGCGTTGGCCTTGCCGCCCTGGCGCAGGTAATGCCAAAGCCGCTCGGCCTGCTCAcc                                                                 >  2:1411832/1‑150 (MQ=255)
                                                                                     gATAGACCGCCGTGCGCGGCAAAGGCTGCGGCTCGTCCCAGGCATAGTCGCGGCCCAGCCATTGGCTGGCCAGGCAGTTGAACAGGTTGATGGCGTTGGCCTTGCCGCCCTGGCGCAGGTAATGCCAAAGCCGCTCGGCCTGCTCCCcgc                                                               >  1:240087/1‑150 (MQ=255)
                                                                                      aTAGACCGCCGTGCGCGGCAAAGGCTGCGGCTCGTCCCAGGCATAGTCGCGGCCCAGCCATTGGCTGGCCAGGCAGTTGAACAgg                                                                                                                               <  1:327390/85‑1 (MQ=255)
                                                                                      aTAGACCGCCGTGCGCGGCAAAGGCTGCGGCTCGTCCCAGGCATAGTCGCGGCCCAGCCATTGGCTGGCCAGGCAGTTGAACAgg                                                                                                                               >  2:327390/1‑85 (MQ=255)
                                                                                             gccgTGCGCGGCAAAGGCTGCGGCTCGTCCCAGGCATAGTCGCGGCCCAGCCATTGGCTGGCCAGGCAGTTGAACAGGTTGATGGCGTTGGCCTTGCCGCCCTGGCGCAGGTAATGCCAAAGCCGCTCGGCCTGCTCACCGCGGACGTTg                                                       <  1:274220/150‑1 (MQ=255)
                                                                                                 tGCGCGGCAAAGGCTGCGGCTCGTCCCAGGCATAGTCGCGGCCCAGCCATTGGCTGGCCAGGCAGTTGAACAGGTTGATGGCGTTGGCCTTGCCGCCCTGGCGCAGGTAATGCCAAAGCCGCTCGGCCTGCTCACCGCGGACGTTGCCCa                                                   <  2:46862/150‑1 (MQ=255)
                                                                                                    gcgGCAAAGGCTGCGGCTCGTCCCAGGCATAGTCGCGGCCCAGCCATTGGCTGGCCAGGCAGTTGAACAGGTTGATGGCGTTGGCCTTGCCGCCCTGGCGCAGGTAATGCCAAAGCCGCTCGGCCTGCTCACCGCGGACGTTGCCCAGGc                                                <  1:64509/150‑1 (MQ=255)
                                                                                                                     tCGTCCCAGGCATAGTCGCGGCCCAGCCATTGGCTGGCCAGGCAGTTGAACAGGTTGATGGCGTTGGCCTTGCCGCCCTGGCGCAGGTAATGCCaaa                                                                                    >  2:876065/1‑97 (MQ=255)
                                                                                                                     tCGTCCCAGGCATAGTCGCGGCCCAGCCATTGGCTGGCCAGGCAGTTGAACAGGTTGATGGCGTTGGCCTTGCCGCCCTGGCGCAGGTAATGCCaaa                                                                                    <  1:876065/97‑1 (MQ=255)
                                                                                                                       gTCCCAGGCATAGTCGCGGCCCAGCCATTGGCTGGCCAGGCAGTTGAACAGGTTGATGGCGTTGGCCTTGCCGCCCTGGCGCAGGTAATGCCAAAGCCGCTCGGCCTGCTCACCGCGGACGTTGCCCAGGCTGGTCAGCtccgggtccgg                             <  1:1045530/150‑1 (MQ=255)
                                                                                                                                aTAGTCGCGGCCCAGCCATTGGCTGGCCAGGCAGTTGAACAGGTTGATGGCGTTGGCCTTGCCGCCCTGGCGCAGGTAATGCCAAAGCCGCTCGGCCTGCTCACCGCGGACGTTGCCCAg                                                  >  2:850278/1‑120 (MQ=255)
                                                                                                                                aTAGTCGCGGCCCAGCCATTGGCTGGCCAGGCAGTTGAACAGGTTGATGGCGTTGGCCTTGCCGCCCTGGCGCAGGTAATGCCAAAGCCGCTCGGCCTGCTCACCGCGGACGTTGCCCAg                                                  <  1:850278/120‑1 (MQ=255)
                                                                                                                                    tCGCGGCCCAGCCATTGGCTGGCCAGGCAGTTGAACAGGTTGATGGCGTTGGCCTTGCCGCCCTGGCGCAGGTAATGCCAAAGCCGCTCGGCCTGCTCACCGCGGACGTTGCCCAGGCTGGTCAGCtccgg                                   <  1:568040/131‑1 (MQ=255)
                                                                                                                                    tCGCGGCCCAGCCATTGGCTGGCCAGGCAGTTGAACAGGTTGATGGCGTTGGCCTTGCCGCCCTGGCGCAGGTAATGCCAAAGCCGCTCGGCCTGCTCACCGCGGACGTTGCCCAGGCTGGTCAGCtccgg                                   >  2:568040/1‑131 (MQ=255)
                                                                                                                                    tCGCGGCCCAGCCATTGGCTGGCCAGGCAGTTGAACAGGTTGATGGCGTTGGCCTTGCCGCCCTGGCGCAGGTAATGCCAAAGCCGCTCGGCCTGCTCACCGCGGACGTTGCCCAGGCTGGTCAGCTCCGGGTCCGGGCGGTCATCGccc                <  1:1422859/150‑1 (MQ=255)
                                                                                                                                    tCGCGGCCCAGCCATTGGCTGGCCAGGCAGTTGAACAGGTTGATGGCGTTGGCCTTGCCGCCCTGGCGCAGGTAATGCCAAAGCCGCTCGGCCTGCTCACCGCGGACGTTGCCCAGGCTGGTCAGCTCCGGGTCCGGGCGGTCATCGccc                <  1:635562/150‑1 (MQ=255)
                                                                                                                                    tCGCGGCCCAGCCATTGGCTGGCCAGGCAGTTGAACAGGTTGATGGCGTTGGCCTTGCCGCCCTGGCGCAGGTAATGCCAAAGCCGCTCGGCCTGCTCACCGCGGACGTTGCCCAGGCTGGTCAGCTCCGGGTCCGGGCGGTCATCGccc                <  1:832545/150‑1 (MQ=255)
                                                                                                                                      gcgGCCCAGCCATTGGCTGGCCAGGCAGTTGAACAGGTTGATGGCGTTGGCCTTGCCACCCTGGCGCAGGTAATGCCAAAGCCGCTCGGCCTGCTCACCGCGGACGTTGCCCAGGCTGGTCAGCTCCGGGTCCGGGCGGTCATCGCCCgg              >  2:1009387/1‑150 (MQ=255)
                                                                                                                                        tggccAGCCAATGGCTGGGCCAGGCAGTTGAACAGGTTGATGGCGTTGGCCTTGCCGCCCTGGCGCAGGTAATGCCAAAGCCGCTCGGCCTGCTCACCGCGGACGtt                                                        >  1:1426603/4‑107 (MQ=255)
                                                                                                                                        tggccAGCCAATGGCTGGGCCAGGCAGTTGAACAGGTTGATGGCGTTGGCCTTGCCGCCCTGGCGCAGGTAATGCCAAAGCCGCTCGGCCTGCTCACCGCGGACGtt                                                        <  2:1426603/104‑1 (MQ=255)
                                                                                                                                        tggccAGCCAATGGCTGGCCAGGCAGTTGAACAGGTTGATGGCGTTGGCCTTGCCGCCCTGGCGCAGGTAATGCCAAAGCCGCTCGGCCTGCTCACCGCGGACGTTGCCCAGGCTGGTCAGCTCCGGGTCCGGGCGGTCATCGCCCGGCa            >  1:1333549/4‑150 (MQ=255)
                                                                                                                                         gCCCAGCCATTGGCTGGCCAGGCAGTTGAACAGGTTGATGGCGTTGGCCTTGCCGCCCTGGCGCAGGTAATGCCAAAGCCGCTCGGCCTGCTCACCGCGGACGTTCCCCAGGCt                                               >  2:579797/1‑114 (MQ=255)
                                                                                                                                         gCCCAGCCATTGGCTGGCCAGGCAGTTGAACAGGTTGATGGCGTTGGCCTTGCCGCCCTGGCGCAGGTAATGCCAAAGCCGCTCGGCCTGCTCACCGCGGACGTTCCCCAGGCt                                               <  1:579797/114‑1 (MQ=255)
                                                                                                                                                cATTGGCTGGCCAGGCAGTTGAACAGGTTGATGGCGTTGGCCTTGCCGCCCTGGCGCAGGTAATGCCAAAGCCGCTCGGCCTGCTCACCGCGGACGTTGCCCAGGCTGGTCAGCTCCGGGTCCGGGCGGTCATCGCCCGGCACCAGAATc    >  2:1154216/1‑150 (MQ=255)
                                                                                                                                                 aTTGGCTGGCCAGGCAGTTGAACAGGTTGATGGCGTTGGCCTTGCCGCCCTGGCGCAGGTAATGCCAAAGCCGCTCGGCCTGCTCACCGCGGACGTTGCCCAGGCTGGTCAGCTCCGGGTCCGGGCGGTCATCGCCCGGCACCAGAATCa   >  2:1147191/1‑150 (MQ=255)
                                                                                                                                                 aTTGGCTGGCCAGGCAGTTGAACAGGTTGATGGCGTTGGCCTTGCCGCCCTGGCGCAGGTAATGCCAAAGCCGCTCGGCCTGCTCACCGCGGACGTTGCCCAGGCTGGTCAGCTCCGGGTCCGGGCGGTCATCGCCCGGCACCAGAATCa   >  1:549138/1‑150 (MQ=255)
                                                                                                                                                 aTTGGCTGGCCAGGCAGTTGAACAGGTTGATGGCGTTGGCCTTGCCGCCCTGGCGCAGGTAATGCCAAAGCCGCTCGGCCTGCTCACCGCGGACGTTGCCCAGGCTGGTCAGCTCCGGGTCCGGGCGGTCATCGCCCGGCACCAGAATCa   >  2:1361232/1‑150 (MQ=255)
                                                                                                                                                  ttGGCTGGCCAGGCAGTTGAACAGGTTGATGGCGTTGGCCTTGCCGCCCTGGCGCAGGTAATGCCAAAGCCGCTCGGCCTGCTCACCGCGGACGTTGCCCAGGCTGGTCAGCTCCGGGTCCGGGCGGTCATCGCCCGGCACCAGAATCAg  <  2:1424937/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                      
CGGTAGAACAGCAGCGGTGCCACTGGCTGCTCGGGGTGCCAGTGTGGGTACCAGTCCTCGAGTGTCGCGCTGCCCTTGGCGGGGTGATAGACCGCCGTGCGCGGCAAAGGCTGCGGCTCGTCCCAGGCATAGTCGCGGCCCAGCCATTGGCTGGCCAGGCAGTTGAACAGGTTGATGGCGTTGGCCTTGCCGCCCTGGCGCAGGTAATGCCAAAGCCGCTCGGCCTGCTCACCGCGGACGTTGCCCAGGCTGGTCAGCTCCGGGTCCGGGCGGTCATCGCCCGGCACCAGAATCAG  >  NC_002947/3980052‑3980347

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: