Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,764,908 G→C 5.4% A46G (GCC→GGC)  PP_2418 ← cobalamin ABC transporter substrate‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,764,9080GC 5.4% 289.8 / 7.0 92A46G (GCC→GGC) PP_2418cobalamin ABC transporter substrate‑binding protein
Reads supporting (aligned to +/- strand):  ref base G (42/45);  new base C (3/2);  total (45/47)
Fisher's exact test for biased strand distribution p-value = 6.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.97e-01

CGCCAGCTTGTCACCCAGCCCCAGGGCATAGAGCAGCTCGGTACCGGCCTGGCCGATGGTGACGGCGCGCTGCGGCGCCTGGGCAAAAACCTGCGGCTTGCCGCAGTTGTCGACCGTCAGCGGGTAGTGGGTGGCAGCGGCCTGGGCCAGGCCAGTGAGGCTTAGGCCCGCGAGCAGGGTGGCGAAACGGAGCAGCATGCTGGGCGATCTCCTATCGGACAGTTCTGGTGAGGAACGCACGGACAGGCATCGCCGAGCCCCAGCGGGTGGCTGGGGCAACACGGATGGCCAATCC  >  NC_002947/2764763‑2765057
                                                                                                                                                 |                                                                                                                                                      
cGCCAGCTTGTCACCCAGCCCCAGGGCATAGAGCAGCTCGGTACCGGCCTGGCCGATGGTGACGGCGCGCTGCGGCGCCTGGGCAAAAACCTGCGGCTTGCCGCAGTTGTCGACCGTCAGCGGGTAGTGGGTGGCAGCGGCCTGggccag                                                                                                                                                    >  2:359220/1‑150 (MQ=255)
cGCCAGCTTGTCACCCAGCCCCAGGGCATAGAGCAGCTCGGTACCGGCCTGGCCGATGGTGACGGCGCGCTGCGGCGCCTGGGCAAAAACCTGCGGCTTGCCGCAGTTGTCGACCGTCAGCGGGTAGTGGGTGGCAGCGGCCTGggccag                                                                                                                                                    >  1:673066/1‑150 (MQ=255)
cGCCAGCTTGTCACCCAGCCCCAGGGCATAGAGCAGCTCGGTACCGGCCTGGCCGATGGCGACGGCGCGCTGCGGCGCCTGGGCAAAAACCTGCGGCTTGCCGCAGTTGTCGACCGTCAGCGGGTAGTGGGTGGCAGCGGCCTGggccag                                                                                                                                                    >  1:1407137/1‑150 (MQ=255)
 gCCAGCTTGTCACCCAGCCCCAGGGCATAGAGCAGCTCGGTACCGGCCTGGCCGATGGTGACGGCGCGCTGCGGCGCCTGGGCAAAAACCTGCGGCTTGCCGCAGTTGTCGACCGTCAGCGGGTAGTGGGTGGCAGCGGCCTGggccagg                                                                                                                                                   <  2:1865233/150‑1 (MQ=255)
      cTTGTCACCCAGCCCCAGGGCATAGAGCAGCTCGGTACCGGCCTGGCCGATGGTGACGGCGCGCTGCGGCGCCTGGGCAAAAACCTGCGGCTTGCCGCAGTTGTCGACCGTCAGCGGGTAGTGGGTGGCAGCGGCCTGGGCCAGGCCAGt                                                                                                                                              <  2:1693662/150‑1 (MQ=255)
         gTCACCCAGCCCCAGGGCATAGAGCAGCTCGGTACCGGCCTGGCCGATGGTGACGGCGCGCTGCGGCGCCTGGGCAAAAACCTGCGGCTTGCCGCAGTTGTCGACCGTCAGCGGGTAGTGGGTGGCAGCGGCCTGGGCCAGGCCAGTGAg                                                                                                                                           >  1:1128403/1‑150 (MQ=255)
         gTCACCCAGCCCCAGGGCATAGAGCAGCTCGGTACCGGCCTGGCCGATGGTGACGGCGCGCTGCGGCGCCTGGGCAAAAACCTGCGGCTTGCCGCAGTTGTCGACCGTCAGCGGGTAGTGGGTGGCAGCGGCCTGGGCCAGGCCAGTGAg                                                                                                                                           <  1:652756/150‑1 (MQ=255)
              ccAGCCCCAGGGCATAGAGCAGCTCGGTACCGGCCTGGCCGATGGTGACGGCGCGCTGCGGCGCCTGGGCAAAAACCTGCGGCTTGCCGCAGTTGTCGACCGTCAGCGGGTAGTGGGTGGCAGCGGCCTGgcccaggccgct                                                                                                                                              >  2:88061/1‑139 (MQ=255)
              ccAGCCCCAGGGCATAGAGCAGCTCGGTACCGGCCTGGCCGATGGTGACGGCGCGCTGCGGCGCCTGGGCAAAAACCTGCGGCTTGCCGCAGTTGTCGACCGTCAGCGGGTAGTGGGTGGCAGCGGCCTGgcccaggccgct                                                                                                                                              <  1:88061/142‑4 (MQ=255)
              ccAGCCCCAGGGCATAGAGCAGCTCGGTACCGGCCTGGCCGATGGTGACGGCGCGCTGCGGCGCCTGGGCAAAAACCTGCGGCTTGCCGCAGTTGTCGACCGTCAGCGGGTAGTGGGTGGCAGCGGCCTGGGCCAGGCCAGTGAGGCTTa                                                                                                                                      >  1:826814/1‑150 (MQ=255)
              ccAGCCCCAGGGCATAGAGCAGCTCGGTACCGGCCTGGCCGATGGTGACGGCGCGCTGCGGCGCCTGGGCAAAAACCTGCGGCTTGCCGCAGTTGTCGACCGTCAGCGGGTAGTGGGTGGCAGCGGCCTGGGCCAGGCCAGTGAGGCTTa                                                                                                                                      >  2:1623196/1‑150 (MQ=255)
              ccAGCCCCAGGGCATAGAGCAGCTCGGTACCGGCCTGGCCGATGGTGACGGCGCGCTGCGGCGCCTGGGCAAAAACCTGCGGCTTGCCGCAGTTGTCGACCGTCAGCGGGTAGTGGGTGGCAGCGGCCTGGGCCAGGCCAGTGAGGCTTa                                                                                                                                      >  1:1923462/1‑150 (MQ=255)
              ccAGCCCCAGGGCATAGAGCAGCTCGGTACCGGCCTGGCCGATGGTGACGGCGCGCTGCGGCGCCTGGGCAAAAACCTGCGGCTTGCCGCAGTTGTCGACCGTCAGCGGGTAGTGGGTGGCAGCGGCCTGGGCCAGGCCAGTGAGGCTTa                                                                                                                                      >  2:350486/1‑150 (MQ=255)
              ccAGCCCCAGGGCATAGAGCAGCTCGGTACCGGCCTGGCCGATGGTGACGCCGCGCTGCGGCGCCTGGGCAAAAACCTGCGGCTTGCCGCAGTTGTCGACCGTCAGCGGGTAGTGGGTGGCAGCGGCCTGGGCCAGGCCAGTGAGGCTTa                                                                                                                                      >  2:1054628/1‑150 (MQ=255)
               cAGCCCCAGGGCATAGAGCAGCTCGGTACCGGCCTGGCCGATGGTGACGGCGCGCTGCGGCGCCTGGGCAAAAACCTGCGGCTTGCCGCAGTTGTCGACCGTCAGCGGGTAGTGGGTGGCAGCGGCCTGGGCCAGGCCAGTGAGGCTTAg                                                                                                                                     <  1:1512390/150‑1 (MQ=255)
               cAGCCCCAGGGCATAGAGCAGCTCGGTACCGGCCTGGCCGATGGTGACGGCGCGCTGCGGCGCCTGGGCAAAAACCTGCGGCTTGCCGCAGTTGTCGACCGTCAGCGGGTAGTGGGTGGCAGCGGCCTGGGCCAGGCCAGTGAGGCTTAg                                                                                                                                     <  2:113074/150‑1 (MQ=255)
               cAGCCCCAGGGCATAGAGCAGCTCGGTACCGGCCTGGCAGATGGTGACGGCGCGCTGCGGCGCCTGGGCAAAAACCTGCGGCTTGCCGCAGTTGTCGACCGTCAGCGGGTAGTGGGTGGCAGCGGCCTGGGCCAGGCCAGTGAGGCTTAg                                                                                                                                     <  1:313316/150‑1 (MQ=255)
                   cccAGGGCATAGAGCAGCTCGGTACCGGCCTGGCCGATGGTGACGGCGCGCTGCGGCGCCTGGGCAAAAACCTGCGGCTTGCCGCAGTTGTCGACCGTCAGCGGGTAGTGGGTGGCAGCGGCCTGGGCCAGGCCAGTGAGGCTTAGGccc                                                                                                                                 >  1:1642325/1‑150 (MQ=255)
                     cAGGGCATAGAGCAGCTCGGTACCGGCCTGGCCGATGGTGACGGCGCGCTGCGGCGCCTGGGCAAAAACCTGCGGCTTGCCGCAGTTGTCGACCGTCAGCGGGTAGTGGGTGGCAGCGGCCTGGGCCAGGCCAGTGAGGCTTAg                                                                                                                                     >  1:1188753/1‑144 (MQ=255)
                     cAGGGCATAGAGCAGCTCGGTACCGGCCTGGCCGATGGTGACGGCGCGCTGCGGCGCCTGGGCAAAAACCTGCGGCTTGCCGCAGTTGTCGACCGTCAGCGGGTAGTGGGTGGCAGCGGCCTGGGCCAGGCCAGTGAGGCNTAg                                                                                                                                     <  2:1188753/144‑1 (MQ=255)
                          cATAGAGCAGCTCGGTACCGGCCTGGCCGATGGTGACGGCGCGCTGCGGCGCCTGGGCAAAAACCTGCGGCTTGCCGCAGTTGTCGACCGTCAGCGGGTAGTGGGTGGCAGCGGCCTGGGCCAGGCCAGTGAGGCTTAGGCCCGCGAGCa                                                                                                                         <  2:1726993/150‑1 (MQ=255)
                            tAGAGCAGCTCGGTACCTGCCTGGCCGATGGTGACGGCGCGCTGCGGCGCCTGGGCAAAAACCTGCGGCTTGCCGCAGTTGTCGACCGTCAGCGGGTAGTGGGTGGCAGCGGCCTGGGCCAGGCCAGTGAGGCTTAGGCCCGCGAGCAgg                                                                                                                       <  1:673625/150‑1 (MQ=255)
                               agcagcTCGGTACCGGCCTGGCCGATGGTGACGGCGCGCTGCGGCGCCTGGGCAAAAACCTGCGGCTTGCCGCAGTTGTCGACCGTCAGCGGGTAGTGGGTGGCAGCGGCCTGGGCCAGGCCAGTGAGGCTTAGGCCCGCGAGCAGGGTg                                                                                                                    <  1:1424393/150‑1 (MQ=255)
                               agcagcTCGGTACCGGCCTGGCCGATGGTGACGGCGCGCTGCGGCGCCTGGGCAAAAACCTGCGGCTTGCCGCAGTTGTCGACCGTCAGCGGGTAGTGGGTGGCAGCGGCCTGGGCCAGGCCAGTGAGGCTTAGGCCCGCGAGCAGGGTg                                                                                                                    <  1:732459/150‑1 (MQ=255)
                               agcagcTCGGTACCGGCCTGGCCGATGGTGACGGCGCGCTGCGGCGCCTGGGCAAAAACCTGCGGCTTGCCGCAGTTGTCGACCGTCAGCGGGTAGTGGGTGGCAGCGGCCTGGGCCAGGCCAGTGAGGCTTAGGCCCGCGAGCAGGGTg                                                                                                                    <  1:592781/150‑1 (MQ=255)
                                        gTACCGGCCTGGCCGATGGTGACGGCGCGCTGCGGCGCCTGGGCAAAAACCTGCGGCTTGCCGCAGTTGTCGACCGTCAGCGGGTAGTGGGTGGCAGCGGCCTGGGCCAGGCCAGTGAGGCTTAGGCCCGCGAGCAGGGTGGCGAAACgg                                                                                                           <  1:175222/150‑1 (MQ=255)
                                         tACCGGCCTGGCCGATGGTGACGGCGCGCTGCGGCGCCTGGGCAAAAACCTGCGGCTTGCCGCAGTTGTCGCCCGGCAGCGGGTAGTGGGTGGCAGCGGCCTGGGCCAGGCCAGTGAGGCTTAGGCCCGCGAGCAGGGTGGCGAAACGGa                                                                                                          <  1:124838/150‑1 (MQ=255)
                                         tACCGGCCTGGCCGATGGTGACGGCGCGCTGCGGCGCCTGGGCAAAAACCTGCGGCTTGCCGCAGTTGTCGACCGTCAGCGGGTAGTGGGTGCCAGCGGCCTGGGCCGCGCCAGTGCGGCGCAGGCCCGCGAGCAGTGTGGCGAAACGGa                                                                                                          >  2:414441/1‑150 (MQ=255)
                                          aCCGGCCTGGCCGATGGTGACGGCGCGCTGCGGCGCCTGGGCAAAAACCTGCGGCTTGCCGCAGTTGTCGACCGTCAGCGGGTAGTGGGTGGCAGCGGCCTGGGCCAGGCCAGTGAGGCTTAGGCCCGCGAGCAGGGTGGCGAAACGGag                                                                                                         <  2:1007886/150‑1 (MQ=255)
                                           ccGGCCTGGCCGATGGTGACGGCGCGCTGCGGCGCCTGGGCAAAAACCTGCGGCTTGCCGCAGTTGTCGACCGTCAGCGGGTAGTGGGTGGCAGCGGCCTGGGCCAGGCCAGTGAGGCTTAGGCCCGCGAGCATGGTGGCGAAACGGagc                                                                                                        <  2:11044/150‑1 (MQ=255)
                                            cGGCCTGGCCGATGGTGACGGCGCGCTGCGGCGCCTTGGCAAAAACCTGCGGCTTGCCGCAGTTGTCGACCGTCAGCGGGTAGTGGGTGGCAGCGGCCTGGGCCAGGCCAGTGAGGCTTATGCCCGCGATCAGGGTGGCGAAACTGagca                                                                                                       <  2:1341926/150‑1 (MQ=255)
                                            cGGCCTGGCCGATGGTGACGGCGCGCTGCGGCGCCTGGGCAAAAACCTGCGGCTTGCCGCAGTTGTCGACCGTCAGCGGGAAGTGGATGGCAGCGGCCTGGGCCAGGCCAATGAGGCTTAGGCCCGCGACCAGGGTGGCTAAACAGagcc                                                                                                       >  2:356868/1‑149 (MQ=255)
                                             ggCCTGGCCGATGGTGACGGCGCGCTGCGGCGCCTGGGCAAAAACCTGCGGCTTGCCGCAGTTGTCGACCGTCAGCGGGTAGTGGGTGGCAGCGGCCTGGGCCAGGCCAGTGAGGCTTAGGCCCGCGAGCAGGGTGGCGAAACGGagcag                                                                                                      <  2:1744498/150‑1 (MQ=255)
                                             ggCCTGGCCGATGGTGACGGCGCGCTGCGGCGCCTGGGCAAAAACCTGCGGCTTGCCGCAGTTGTCGACCGTCAGCGGGTAGTGGGTGGCAGCGGCCTGGGCCAGGCCAGTGAGGCTTAGGCCCGCGAGCAGGGTGGCGAAACGGagcag                                                                                                      >  2:1627459/1‑150 (MQ=255)
                                                 tGGCCGATGGCGACGGCGCGCTGCGGCGCCTGGGCAAAAACCTGCGGCTTGCCGCAGTTGTCGACCGTCAGCGGGTAGTGGGTGGCAGCGGCCTGGGCCAGGCCAGTGAGGCTTAGGCCCGCGAGCAGGGTGGCGAAACGGAGCAGCATg                                                                                                  <  2:1407137/150‑1 (MQ=255)
                                                    ccGATGGTGACGGCGCGCTGCGGCGCCTGGGCAAAAACCTGCGGCTTGCCGCAGTTGTCGACCGTCAGCGGGTAGTGGGTGGCAGCGGCCTGGGCCAGGCCAGTGAGGCTTAGGCCCGCGAGCAGGGTGGCGAAACGGAGCAGCATGCTg                                                                                               <  1:190355/150‑1 (MQ=255)
                                                    ccGATGGTGACGGCGCGCTGCGGCGCCTGGGCAAAAACCTGCGGCTTGCCGCAGTTGTCGACCGTCAGCGGGTAGTGGGTGGCAGCGGCCTGGGCCAGGCCAGTGAGGCTTAGGCCCGCGAGCAGGGTGGCGAAACGGAGCAGCATGCTg                                                                                               <  2:1081228/150‑1 (MQ=255)
                                                     cGATGGTGACGGCGCGCTGCGGCGCCTGGGCAAAAACCTGCGGCTTGCCGCAGTTGTCGACCGTCAGCGGGTAGTGGGTGGCAGCGGCCTGGGCCAGGCCAGTGAGGCTTAGGCCCGCGAGCAGGGTGGCGAAACGGAGCAGCATGCTgg                                                                                              <  2:1923462/150‑1 (MQ=255)
                                                     cGATGGTGACGGCGCGCTGCGGCGCCTGGGCAAAAACCTGCGGCTTGCCGCAGTTGTCGACCGTCAGCGGGTAGTGGGTGGCAGCGGCCTGGGCCAGGCCAGTGAGGCTTAGGCCCGCGAGCAGGGTGGCGAAACGGAGCAGCATGCTgg                                                                                              >  2:1246468/1‑150 (MQ=255)
                                                     cGATGGTGACGGCGCGCTGCGGCGCCTGGGCAAAAACCTGCGGCTTGCCGCAGTTGTCGACCGTCAGCGGGTAGTGGGTGGCAGCGGCCTGGGCCAGGCCAGTGAGGCTTAGGCCCGCGAGCAGGGTGGCGAAACGGAGCAGCATGCTgg                                                                                              >  1:1145180/1‑150 (MQ=255)
                                                      gATGGTGACGGCGCGCTGCGGCGCCTGGGCAAAAACCTGCGGCTTGCCGCAGTTGTCGACCGTCAGCGGGTAGTGGGTGGCAGCGGCCTGGGCCAGGCCAGTGAGGCTTAGGCCCGCGAGCAGGGTGGCGAAACGGAGCAGCATGCTggg                                                                                             >  1:340303/1‑150 (MQ=255)
                                                       aTGGTGACGGCGCGCTGCGGCGCCTGGGCAAAAACCTGCGGCTTGCCGCAGTTGTCGACCGTCAGCGGGTAGTGGGTGGCAGCGGCCTGGGCCAGGCCAGTGAGGCTTAGGCCCGCGAGCAGGGTGGCGAAACGGAGCAGCATGCTgtgg                                                                                            <  2:34771/150‑4 (MQ=255)
                                                             aCGGCGCGCTGCGGCGCCGGGGCAAAAACCTGCGGCTTGCCGCAGTTGTCGACCGTCAGCGGGTAGTGGGTGGCAGCGGCCTGGGCCAGGCCAGTGAGGCTTAGGCCCGCGAGCAGGGTGGCGAAACGGAGCAGCATGCTGGGCGAtctc                                                                                      <  2:46156/150‑1 (MQ=255)
                                                               ggCGCGCTGCGGCGCCTGGGCAAAAACCTGCGGCTTGCCGCAGTTGTCGACCGTCAGCGGGAAGTGGGTGGCAGCGGCCTGGGCCACGCACGCCGGCGTCTGCTTGAAAAACAAAGCGGCGAAACGGAGCAGCATGCTGGGCGATCTCCt                                                                                    >  1:1375863/1‑150 (MQ=255)
                                                                 cgcgcTGCGGCGCCTGGGCAAAAACCTGCGGCTTGCCGCAGTTGTCGACCGTCAGCGGGTAGTGGGTGGCAGCGGCCTGGGCCAGGCCAGTGAGGCTTAGGCCCGCGAGCAGGGTGGCGAAACGGAGCAGCACGCTGGGCGATCTCCTAt                                                                                  <  1:1591776/150‑1 (MQ=255)
                                                                    gcTGCGGCGCCTGGGCAAAAACCTGCGGCTTGCCGCAGTTGTCGACCGTCAGCGGGTAGTGGGTGGCAGCGGCCTGGGCCAGGCCAGTGAGGCTTAGGCCCGCGAGCAGGGTGGCGAAACGGAGCAGCATGCTGGGCGATCTCCTATCgg                                                                               >  2:1237749/1‑150 (MQ=255)
                                                                       gcggcgCCTGGGCAAAAACCTGCGGCTTGCCGCAGTTGTCGACCGTCAGCGGGTAGTGGGTGGCAGCGGCCTGGGCCAGGCCAGTGAGGCTTAGGCCCGCGAGCAGGGTGGCGAAACGGAGCAGCATGCTGGGCGATCTCCTATCGGACa                                                                            <  2:1394253/150‑1 (MQ=255)
                                                                       gcggcgCCTGGGCAAAAACCTGCGGCTTGCCGCAGTTGTCGACCGTCAGCGGGTAGTGGGTGGCAGCGGCCTGGGCCAGGCCAGTGAGGCTTAGGCCCGCGAGCAGGGTGGCGAAACGGAGCAGCATGCTGGGCGATCTCCTATCGGACa                                                                            <  2:340303/150‑1 (MQ=255)
                                                                          gcgcCTGGGCAAAAACCTGCGGCTTGCCGCAGTTGTCGACCGTCAGCGGGTAGTGGGTGGCAGCGGCCTGGGCCAGGCCAGt                                                                                                                                              <  1:1508862/82‑1 (MQ=255)
                                                                          gcgcCGGGGCAAAAACCTGCGGCTTGCCGCAGTTGTCGACCGTCAGCGGGTAGTGGGTGGCAGCGGCCTGGGCCAGGCCAGt                                                                                                                                              >  2:1508862/1‑82 (MQ=255)
                                                                             ccTGGGCAAAAACCTGCGGCTTGCCGCAGTTGTCGACCGTCAGCGGGTAGTGGGTGGCAGCGGCCTGGGCCAGGCCAGTGAGGCTTAGGCCCGCGAGCAGGGTGGCGAAACGGAGCAGCATGCTGGGCGATCTCCTATCGGACAGTTCTg                                                                      <  1:1237749/150‑1 (MQ=255)
                                                                                 ggCAAAAACCTGCGGCTTGCCGCAGTTGTCGACCGTCAGCGGGTAGTGGGTGGCAGCGGCCTGGGCCAGGCCAGTGAGGCTTAGGCCCGCGAGCAGGGTGGCGAAACGGAGCAGCATGCTGGGCGATCTCCTATCGGACAGTTCTGGTGa                                                                  >  1:1442764/1‑150 (MQ=255)
                                                                                  gCAAAAACCTGCGGCTTGCCGCAGTTGTCGACCGTCAGCGGGTAGTGGGTGGCAGCGGCCTGGGCCAGGCCAGTGAGGCTTAGGCCCGCGAGCAGGGTGGCGAAACGGAGCAGCATGCTGGGCGATCTCCTATCGGACAGTTCTGGTGAg                                                                 >  2:1527766/1‑150 (MQ=255)
                                                                                  gCAAAAACCTGCGGCTTGCCGCAGTTGTCGACCGTCAGCGGGTAGTGGGTGGCAGCGGCCTGGGCCAGGCCAGTGAGGCTTAGGCCCGCGAGCAGGGTGGCGAAACGGAGCAGCATGCTGGGCGATCTCCTATCGGACAGTTCTGGTGAg                                                                 <  1:1307817/150‑1 (MQ=255)
                                                                                   cAAAAACCTGCGGCTTGCCGCAGTTGTCGACCGTCAGCGGGTAGTGGGTGGCAGCGGCCTGGGCCAGGCCAGTGAGGCTTAGGCCCGCGAGCAGGGTGGCGAAACGGAGCAGCATGCTGGGCGATCTCCTATCGGACAGTTCTGGTGAgg                                                                <  2:571326/150‑1 (MQ=255)
                                                                                     aaaaCCTGCGGCTTGCCGCAGTTGTCGACCGTCAGCGGGTAGTGGGTGGCAGCGGCCTGGGCCAGGCCAGTGAGGCTTAGGCCCGCGAGCAGGGTGGCGAAACGGAGCAGCATGCTgg                                                                                              >  1:1666813/1‑118 (MQ=255)
                                                                                     aaaaCCTGCGGCTTGCCGCAGTTGTCGACCGTCAGCGGGTAGTGGGTGGCAGCGGCCTGGGCCAGGCCAGTGAGGCTTAGGCCCGCGAGCAGGGTGGCGAAACGGAGCAGCATGCTgg                                                                                              <  2:1666813/118‑1 (MQ=255)
                                                                                     aaaaCCTGCGGCTTGCCGCAGTTGTCGACCGTCAGCGGGTAGTGGGTGGCAGCGGCCTGGGCCAGGCCAGTGAGGCTTAGGCCCGCGAGCAGGGTGGCGAAACGGAGCAGCATGCTGGGCGATCTCCTATCGGACAGTTCTGGTGAGGaa                                                              <  2:1854083/150‑1 (MQ=255)
                                                                                      aaaCCTGCGGCTTGCCGCAGTTGTCGACCGTCAGCGGGTAGTGGGTGGCAGCGGCCTGGGCCAGGCCAGTGAGGCTTAGGCCCGCGAGCAGGGTGGCGAAACGGAGCAGCATGCTGGGCGATCTCCTATCGGACAGTTCTGGTGAGGAAc                                                             >  2:605709/1‑150 (MQ=255)
                                                                                      aaaCCTGCGGCTTGCCGCAGTTGTCGACCGTCAGCGGGTAGTGGGTGGCAGCGGCCTGGGCCAGGCCAGTGAGGCTTAGGCCCGCGAGCAGGGTGGCGAAACGGAGCAGCATGCTGGGCGATCTCCTATCGGACAGTTCTGGTGAGGAAc                                                             <  2:1642325/150‑1 (MQ=255)
                                                                                        aacTGCGGCAAGCCGCAGTTGTCGACCGTCAGCGGGTAGTGGGTGGCAGCGGCCTGGGCCAGGCCAGTGAGGCTTAGGCCCGCGAGCAGGGTGGCGAAACGGAGCAGCATGCTGGGCGATCTCCTATCGGACAGTTCTGGTGAGGAACGc                                                           >  2:1741875/3‑150 (MQ=255)
                                                                                         ccTGCGGCTTGCCGCAGTTGTCGACCGTCAGCGGGTAGTGGGTGGCAGCGGCCTGGGCCAGGCCAGTGAGGCTTAGGCCCGCGAGCAGGGTGGCGAAACGGAGCAGCATGCTGGGCGATCTCCTATCGGACAGTTCTGGTGAGGAACGCa                                                          <  2:1145180/150‑1 (MQ=255)
                                                                                         ccTGCGGCTTGCCGCAGTTGTCGACCGTCAGCGGGTAGTGGGTGGCAGCGGCCTGGGCCAGGCCAGTGAGGCTTAGGCCCGCGAGCAGGGTGGCGAAACGGAGCAGCATGCTGGGCGATCTCCTATCGGACAGTTCTGGTGAGGAACGCa                                                          >  2:442563/1‑150 (MQ=255)
                                                                                          cTGCGGCTTGCCGCAGTTGTCGACCGTCAGCGGGTAGTGGGTGGCAGCGGCCTGGGCCAGGCCAGTGAGGCTTAGGCCCGCGAGCAGGGTGGCGAAACGGAGCAGCATGCTGGGCGATCTCCTATCGGACAGTTCTGGTGAGGAACGCAc                                                         >  2:1144710/1‑150 (MQ=255)
                                                                                          cTGCGGCTTGCCGCAGTTGTCGACCGTCAGCGGGTAGTGGGTGGCAGCGGCCTGGGCCAGGCCAGTGAGGCTTAGGCCCGCGAGCAGGGTGGCGAAACGGAGCAGCATGCTGGGCGATCTCCTATCGGACAGTTCTGGTGAGGAACGCAc                                                         >  2:285258/1‑150 (MQ=255)
                                                                                             cGGCTTGCCGCAGTTGTCGACCGTCAGCGGGTAGTGGGTGGCAGCGGCCTGGGCCAGGCCAGTGAGGCTTAGGCCCGCGAGCAGGGTGGCGAAACGGAGCAGCATGCTGGGCGATCTCCTATCGGACAGTTCTGGTGAGGAACGCACGGa                                                      <  1:878832/150‑1 (MQ=255)
                                                                                              ggCTTGCCGCAGTTGTCGACCGTCAGCGGGTAGTGGGTGGCAGCGGCCTGGGCCAGGCCAGTGAGGCTTAGGCCCGCGAGCAGGGTGGCGAAACGGAGCAGCATGCTGGGCGATCTCCTATCGGACAGTTCTGGTGAGGAACGCACGGAc                                                     <  1:1044617/150‑1 (MQ=255)
                                                                                              ggCTTGCCGCAGTTGTCGACCGTCAGCGGGTAGTGGGTGGCAGCGGCCTGGGCCAGGCCAGTGAGGCTTAGGCCCGCGAGCAGGGTGGCGAAACGGAGCAGCATGCTGGGCGATCTCCTATCGGACAGTTCTGGTGAGGAACGCACGGAc                                                     >  1:1141509/1‑150 (MQ=255)
                                                                                                     cGCAGTTGTCGACCGTCAGCGGGTAGTGGGTGGCAGCGGCCTGGGCCAGGCCAGTGAGGCTTAg                                                                                                                                     <  2:1214063/64‑1 (MQ=255)
                                                                                                     cGCAGTTGTCGACCGTCAGCGGGTAGTGGGTGGCAGCGGCCTGGGCCAGGCCAGTGAGGCTTAg                                                                                                                                     >  1:1214063/1‑64 (MQ=255)
                                                                                                          ttGTCGACCGTCAGCGGGTAGTGGGTGGCAGCGGCCTGGGCCAGGCCAGTGAGGCTTAGGCCCGCGAGCAGGGTGGCGAAACGGAGCAGCATGCTGGGCGATCTCCTATCGGACAGTTCTGGTGAGGAACGCACGGACAGGCATCGCCGa                                         <  1:605709/150‑1 (MQ=255)
                                                                                                               gACCGTCAGCGGGTAGTGGGTGGCAGCGGCCTGGGCCAGGCCAGTGAGGCTTAGGCCCGCGAGCAGGGTGGCGAAACGGAGCAGCATGCTGGGCGATCTCCTATCGGACAGTTCTGGTGAGGAACGCACGGACAGGCATCGCCGAGcccc                                    >  1:1838460/1‑150 (MQ=255)
                                                                                                               gACCGTCAGCGGGTAGTGGGTGGCAGCGGCCTGGGCCAGGCCAGTGAGGCTTAGGCCCGCGAGCAGGGTGGCGAAACGGAGCAGCATGCTGGGCGATCTCCTATCGGACAGTTCTGGGGAGGAACGCACGGACAGGCATCGCCGAGcaaa                                    >  2:1350185/1‑147 (MQ=255)
                                                                                                                aCCGTCAGCGGGTAGTGGGTGGCAGCGGCCTGGGCCAGGCCAGTGAGGCTTAGGCCCGCGAGCAGGGTGGCGAAACGGAGCAGCATGCTGGGCGATCTCCTATCGGACAGTTCTGGTGAGGAACGCACGGACAGGCAGCGCCGAGccccc                                   >  2:141327/1‑149 (MQ=255)
                                                                                                                       gCGGGTAGTGGGTGGCAGCGGCCTGGGCCAGGCCAGTGAGGCTTAGGCCCGCGAGCAGGGTGGCGAAACGGAGCAGCATGCTGGGCGATCTCCTATCGGACAGTTCTGGTGAGGAACGCACGGACAGGCATCGCCGAGCCCCAGCGGGTg                            >  1:90435/1‑150 (MQ=255)
                                                                                                                       gCGGGTAGTGGGTGGCAGCGGCCTGGGCCAGGCCAGTGAGGCTTAGGCCCGCGAGCAGGGTGGCGAAACGGAGCAGCATGCTGGGCGATCTCCTATCGGACAGTTCTGGTGAGGAACGCACGGACAGGCATCGCCGAGCCCCAGCGGGTg                            >  2:1293069/1‑150 (MQ=255)
                                                                                                                           gTAGTGGGTGGCAGCGGCCTGGGCCAGGCCAGTGAGGCTTAGGCCCGCGAGCAGGGTGGCGAAACGGAGCAGCATGCTGGGCGATCTCCTATCGGACAGTTCTGGTGAGGAACGCACGGACAGGCATCGCCGAGCCCCAGCGGGTGGCTg                        >  2:586886/1‑150 (MQ=255)
                                                                                                                            tAGTGGGTGGCAGCGGCCTGGGCCAGGCCAGTGAGGCTTAGGCCCGCGAGCAGGGTGGCGAAACGGAGCAGCATGCTGGGCGATCTCCTATCGGACAGTTCTGGTGAGGAACGCACGGACAGGCATCGCCGAGCCCCAGCGGGTGGCTgg                       >  1:1527138/1‑150 (MQ=255)
                                                                                                                            tAGTGGGTGGCAGCGGCCTGGGCCAGGCCAGTGAGGCTTAGGCCCGCGAGCAGGGTGGCGAAACGGAGCAGCATGCTGGGCGATCTCCTATCGGACAGTTCTGGTGAGGAACGCACGGACAGGCATCGCCGAGCCCCAGCGGGTGGCTgg                       <  1:608844/150‑1 (MQ=255)
                                                                                                                            tAGTGGGTGGCAGCGGCCTGGGCCAGGCCAGTGAGGCTTAGGCCCGCGAGCAGGGTGGCGAAACGGAGCAGCATGCTGGGCGATCTCCTATCGGACAGTTCTGGTGAGGAACGCACGGACAGGCATCGCCGAGCCCCAGCGGGTGGCTgg                       <  1:359220/150‑1 (MQ=255)
                                                                                                                             aGTGGGTGGCAGCGGCCTGGGCCAGGCCAGTGAGGCTTAGGCCCGCGAGCAGGGTGGCGAAACGGAGCAGCATGCTGGGCGATCTCCTATCGGACAGTTCTGGTGAGGAACGCACGGACAGGCATCGCCGAGCCCCAGCGGGTGGCTggg                      <  2:1845309/150‑1 (MQ=255)
                                                                                                                             aGTGGGTGGCAGCGGCCTGGGCCAGGCCAGTGAGGCTTAGGCCCGCGAGCAGGGTGGCGAAACGGAACAGCCAGCAGTTTGTTCTCCAAACGGTGTTTTCTGGTGAGGAACGCACGGACAGGCATCGCCGAGCCCCAGCGGGTGGCTggg                      <  1:1671711/150‑1 (MQ=255)
                                                                                                                              gtgggtggCAGCGGCCTGGGCCAGGCCAGTGAGGCTTAGGCCCGCGAGCAGGGTGGCGAAACGGAGCAGCATGCTGGGCGATCTCCTATCGGACAGTTCTGGTGAGGAACGCACGGACAGGCATCGCCGAGCCCCAGCGGGTGGCTgggg                     <  1:1627459/150‑1 (MQ=255)
                                                                                                                                 ggtggCAGCGGCCTGGGCCAGGCCAGTGAGGCTTAGGCCCGCGAGCAGGGTGGCGAAACGGAGCAGCATGCTGGGCGATCTCCTATCGGACAGTTCTGGTGAGGAACGCACGGACAGGCATCGCCGAGCCCCAGCGGGTGGCTGGGGCaa                  >  1:137675/1‑150 (MQ=255)
                                                                                                                                    ctcactggCCTGGCCCAGGCCGGTGAGGCTTAGGCCCGCGAGCAGGGTGGCGAAACGGAGCAGCATGCTGGGCGATCTCCTATCGGACAGTTCTGGTGAGGAACGCACGGACAGGCATCGCCGAGCCCCAGCGGGTGGCTGGGGCAacac               >  1:1649251/7‑150 (MQ=255)
                                                                                                                                     tcactggCCTGGCCCAGGCCGGTGAGGCTTAGGCCCGCGAGCAGGGTGGCGAAACGGAGCAGCATGCTGGGCGATCTCCTATCGGACAGTTCTGGTGAGGAACGCACGGACAGGCATCGCCGAGCCCCAGCGGGTGGCTGGGGCAACACg              <  2:1649251/145‑1 (MQ=255)
                                                                                                                                     gCAGCGGCCTGGGCCAGGCCAGTGAGGCTTAGGCCCGCGAGCAGGGTGGCGAAACGGAGCAGCATGCTGGGCGATCTCCTATCGGACAGTTCTGGTGAGGAACGCACGGACAGGCATCGCCGAGCCCCAGCGGGTGGCTGGGGCAACACg              >  2:1481388/1‑150 (MQ=255)
                                                                                                                                      cAGCGGCCTGGGCCAGGCCAGTGAGGCTTAGGCCCGCGAGCAGGGTGGCGAAACGGAGCAGCATGCTGGGCGATCTCCTATCGGACAGTTCTGGTGAGGAACGCACGGACAGGCATCGCCGAGCCCCAGCGGGTGGCTGGGGCAACACgg             >  1:1744761/1‑150 (MQ=255)
                                                                                                                                      cAGCGGCCTGGGCCAGGCCAGTGAGGCTTAGGCCCGCGAGCAGGGTGGCGAAACGGAGCAGCATGCTGGGCGATCTCCTATCGGACAGTTCTGGTGAGGAACGCACGGACAGGCATCGCCGAGCCCCAGCGGGTGGCTGGGGCAACACgg             >  1:1415798/1‑150 (MQ=255)
                                                                                                                                       actggCCTGGCCCAGGCCAGTGAGGCTTAGGCCCGCGAGCAGGGTGGCGAAACGGAGCAGCATGCTGGGCGATCTCCTATCGGACAGTTCTGGTGAGGAACGCACGGACAGGCATCGCCGAGCCCCAGCGGGt                             >  2:729262/4‑133 (MQ=255)
                                                                                                                                       actggCCTGGCCCAGGCCAGTGAGGCTTAGGCCCGCGAGCAGGGTGGCGAAACGGAGCAGCATGCTGGGCGATCTCCTATCGGACAGTTCTGGTGAGGAACGCACGGACAGGCATCGCCGAGCCCCAGCGGGt                             <  1:729262/130‑1 (MQ=255)
                                                                                                                                        gCGGCCTGGGCCAGGCCAGTGAGGCTTAGGCCCGCGAGCAGGGTGGCGAAACGGAGCAGCATGCTGGGCGATCTCCTATCGGACAGTTCTGGTGAGGAACGCACGGACAGGCATCGCCGAGCCCCAGCGGGTGGCTGGGGCAACACGGAt           <  1:586886/150‑1 (MQ=255)
                                                                                                                                        gCGGCCTGGGCCAGGCCAGTGAGGCTTAGGCCCGCGAGCAGGGTGGCGAAACGGAGCAGCATGCTGGGCGATCTCCTATCGGACAGTTCTGGTGAGGAACGCACGGACAGGCATCGCCGAGCCCCAGCGGGTGGCTGGGGCAACACGGAt           >  1:1306888/1‑150 (MQ=255)
                                                                                                                                         tggCCTGGCCCAGGCCAGTGAGGCTTAGGCCAGACCGCATGGTGGCGAAACGGAGCAGCATGCTGGGCGATCTCCTATCGGACAGTTCTGGTGAGGAACGCACGGACAGGCATCGCCGAGCCCCAGCGGGTGGCTGGGGCAACACGGATg          >  1:912618/2‑150 (MQ=255)
                                                                                                                                           gCCTGGGCCAGGCCAGTGAGGCTTAGGCCCGCGAGCAGGGTGGCGAAACGGAGCAGCATGCTGGGCGATCTCCTATCGGACAGTTCTGGTGAGGAACGCACGGACAGGCATCGCCGAGCCCCAGCGGGTGGCTGGGGCAACACGGATGGc        <  1:1144710/150‑1 (MQ=255)
                                                                                                                                             cTGGGCCAGGCCAGTGAGGCTTAGGCCCGCGAGCAGGGTGGCGAAACGGAGCAGCATGCTGGGCGATCTCCTATCGGACAGTTCTGGTGAGGAACGCACGGACAGGCATCGCCGAGCCCCAGCGGGTGGCTGGGGCAACACGGATGGCCa      >  1:784314/1‑150 (MQ=255)
                                                                                                                                             cTGGGCCAGGCCAGTGAGGCTTAGGCCCGCGAGCAGGGTGGCGAAACGGAGCAGCATGCTGGGCGATCTCCTATCGGACAGTTCTGGTGAGGAACGCACGGACAGGCATCGCCGAGCCCCAGCGGGTGGCTGGGGCAACACGGATGGCCa      >  1:1618865/1‑150 (MQ=255)
                                                                                                                                             cTGGGCCAGGCCAGTGAGGCTTAGGCCCGCGAGCAGGGTGGCGAAACGGAGCAGCATGCTGGGCGATCTCCTATCGGACAGTTCTGGTGAGGAACGCACGGACAGGCATCGCCGAGCCCCAGCGGGTGGCTGGGGCAACACGGATGGCCa      >  1:662750/1‑150 (MQ=255)
                                                                                                                                                 gccaggccagTGAGGCTTAGGCCCGCGAGCAGGGTGGCGAAACGGAGCAGCATGCTGGGCGATCTCCTATCGGACAGTTCTGGTGAGGAACGCACGGACAGGCATCGCCGAGCCCCAGCGGGTGGCTGGGGCAACACGGATGGCCAATcc  <  1:350486/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                      
CGCCAGCTTGTCACCCAGCCCCAGGGCATAGAGCAGCTCGGTACCGGCCTGGCCGATGGTGACGGCGCGCTGCGGCGCCTGGGCAAAAACCTGCGGCTTGCCGCAGTTGTCGACCGTCAGCGGGTAGTGGGTGGCAGCGGCCTGGGCCAGGCCAGTGAGGCTTAGGCCCGCGAGCAGGGTGGCGAAACGGAGCAGCATGCTGGGCGATCTCCTATCGGACAGTTCTGGTGAGGAACGCACGGACAGGCATCGCCGAGCCCCAGCGGGTGGCTGGGGCAACACGGATGGCCAATCC  >  NC_002947/2764763‑2765057

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