Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,508,623 C→T 7.9% A68T (GCC→ACC)  cysG ← multifunctional uroporphyrinogen‑III C‑methyltransferase/precorrin‑2 dehydrogenase/sirohydrochlorin ferrochelatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,508,6230CT 7.9% 234.8 / 7.6 84A68T (GCC→ACC) cysGmultifunctional uroporphyrinogen‑III C‑methyltransferase/precorrin‑2 dehydrogenase/sirohydrochlorin ferrochelatase
Reads supporting (aligned to +/- strand):  ref base C (45/32);  new base T (4/3);  total (49/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.71e-01

ACGATCGCCGGGAAGATCACCGTACACAGGGCCGGCGCATCCACCACATTGACCGGCAGGCTAAGGGCCTGCGCATCGGCCGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCTTCCCGGGCCAACGCGGCCAACTGGCCGTCGACGTCCGGCGCCACCACGCGCAGCACGCCGCCGGCATCGGCCAGCAGGCGCGCCTTGCGCAAGGCGATCTCACCGCCA  >  NC_002947/4508475‑4508768
                                                                                                                                                    |                                                                                                                                                 
aCGATCGCCGGGAAGATCACCGTACACAGGGCCGGCGCATCCACCACATTGACCGGCAGGCTAAGGGCCTGCGCATCGGCCGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCt                                                                                                                                                  >  1:1704607/1‑150 (MQ=255)
  gATCGCCGGGAAGATCACCGTACACAGGGCCGGCGCATCCACCACATTGACCGGCAGGCTAAGGGCCTGCGCATCGGCCGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGc                                                                                                                                                <  2:5079/150‑1 (MQ=255)
       ccGGGAAGATCACCGTACACAGGGCCGGCGCATCCACCACATTGACCGGCAGGCTAAGGGCCTGCGCATCGGCCGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGAt                                                                                                                                           >  1:1388708/1‑150 (MQ=255)
         gggAAGATCACCGTACACAGGGCCGGCGCATCCACCACATTGACCGGCAGGCTAAGGGCCTGCGCATCGGCCGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCCGCCTGATAg                                                                                                                                         >  1:1131674/1‑150 (MQ=255)
         gggAAGATCACCGTACACAGGGCCGGCGCATCCACCACATTGACCGGCAGGCTAAGGGCCTGCGCATCGGCCGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCAACCAGGTCGGCTGCCTGATAg                                                                                                                                         <  1:554685/150‑1 (MQ=255)
                 cACCGTACACAGGGCCGGCGCATCCACCACATTGACCGGCAGGCTAAGGGCCTGCGCATCGGCCGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCAc                                                                                                                                 <  1:1029465/150‑1 (MQ=255)
                   ccGTACACAGGGCCGGCGCATCCACCACATTGACCGGCAGGCTAAGGGCCTGCGCATCGGCCGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCa                                                                                                                               >  2:1220737/1‑150 (MQ=255)
                   ccGTACACAGGGCCGGCGCATCCACCACATTGACCGGCAGGCTAAGGGCCTGCGCATCGGCCGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCa                                                                                                                               >  1:219649/1‑150 (MQ=255)
                    cGTACACAGGGCCGGCGCATCCACCACATTGACCGGCAGGCTAAGGGCCTGCGCATCGGCCGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAg                                                                                                                              <  1:1708683/150‑1 (MQ=255)
                       acacAGGGCCGGCGCATCCACCACATTGACCGGCAGGCTAAGGGCCTGCGCATCGGCCGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCAc                                                                                                                           <  1:547993/150‑1 (MQ=255)
                         acaGGGCCGGCGCATCCACCACATTGACCGGCAGGCTAAGGGCCTGCGCATCGGCCGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCt                                                                                                                         >  2:1246525/1‑150 (MQ=255)
                             ggCCGGCGCATCCACCACATTGACCGGCAGGCTAAGGGCCTGCGCATCGGCCGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCAcc                                                                                                                     <  1:868505/150‑1 (MQ=255)
                                  gcgcATCCACCACATTGACCGGCAGGCTAAGGGCCTGCGCATCGACAGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCAc                                                                                                                >  2:1212847/1‑150 (MQ=255)
                                          accacATTGACCGGCAGGCTAAGGGCCTGCGCATCGGCCGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGc                                                                                                                                   >  1:1700259/1‑123 (MQ=255)
                                          accacATTGACCGGCAGGCTAAGGGCCTGCGCATCGGCCGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCTTCCCgg                                                                                                        >  1:708905/1‑150 (MQ=255)
                                          accacATTGACCGGCAGGCTAAGGGCCTGCGCATCGGCCGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCTTCCCgg                                                                                                        >  2:1242892/1‑150 (MQ=255)
                                          accacATTGACCGGCAGGCTAAGGGCCTGCGCATCGGCCGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCTTCCCgg                                                                                                        >  2:1432018/1‑150 (MQ=255)
                                          accacATTGACCGGCAGGCTAAGGGCCTGCGCATCGGCCGACACCTGGGCATTCAGCCCAGGGTCGGCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGc                                                                                                                                   <  2:1700259/123‑1 (MQ=255)
                                           ccacATTGACCGGCAGGCTAAGGGCCTGCGCATCGGCCGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCTTCCCggg                                                                                                       >  1:807809/1‑150 (MQ=255)
                                               aTTGACCGGCAGGCTAAGGGCCTGCGCATCGGCCGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCTTCCCGGGCCaa                                                                                                   >  1:1029220/1‑150 (MQ=255)
                                                    ccGGCAGGCTAAGGGCCTGCGCATCGGCCGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCTTCCCGGGCCAACGCgg                                                                                              <  2:619832/150‑1 (MQ=255)
                                                       gCAGGCTAAGGGCCTGCGCATCGGCCGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCTTCCCGGGCCAACGCGGCCa                                                                                           >  1:840754/1‑150 (MQ=255)
                                                          ggCTAAGGGCCTGCGCATCGGCCGACACCTGGGCATTCCGCCCAGCGCAGTCCGTCGCCCCGATCTCCCGCCGGCGCCTGCCCCCGGTGGTTGCCTGATAGCCTCTTACTAGCACCTCACCGCCACCTTCCCGGGCCAACGCGGCCAACt                                                                                        <  2:1160948/150‑1 (MQ=255)
                                                           gaaaaGGGCCTGCGCATCGGCCGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCTTCCCGGGCCAACGCGGCCAACTg                                                                                       >  1:1316035/4‑150 (MQ=255)
                                                             tAAGGGCCTGCGCATCGGCCGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCTTCCCGGGCCAACGCGGCCAACTGGc                                                                                     >  1:92270/1‑150 (MQ=255)
                                                              aaGGGCCTGCGCATCGGCCGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCTTCCCGGGCCAACGCGGCCAACTGGcc                                                                                    <  2:1704607/150‑1 (MQ=255)
                                                              aaGGGCCTGCGCATCGGCCGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCTTCCCGGGCCAACGCGGCCAACTGGcc                                                                                    <  1:939554/150‑1 (MQ=255)
                                                                        gcATCGGCCGACACCTGCGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCCCACCGCCACCTTCCCGCGCCAACGCGGCCAACTGGCCGTCGACGTcc                                                                          >  2:809096/1‑150 (MQ=255)
                                                                         cATCGGCCGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCTTCCCGGGCCAACGCGGCCAACTGGCCGTCGACGTACg                                                                         >  2:1631818/1‑150 (MQ=255)
                                                                           tCGGCCGACACCTGGGCATGCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCTTCCCGGGCCAACGCGGCCAACTGGCCGTCGACGTCCGgc                                                                       <  2:104770/150‑1 (MQ=255)
                                                                               ccGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCAGCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCAcc                                                                                                                     >  1:1145753/1‑100 (MQ=255)
                                                                               ccGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCAGCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCAcc                                                                                                                     <  2:1145753/100‑1 (MQ=255)
                                                                               ccGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCTTCCCGGGCCAACGCGGCCAACTGGCCGTCGACGTCCGGCGcca                                                                   >  1:671052/1‑150 (MQ=255)
                                                                                cGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCTTCCCGGGCCAACGCGACCAACTGGCCGTCGa                                                                               <  2:1520806/137‑1 (MQ=255)
                                                                                cGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCTTCCCGGGCCAACGCGACCAACTGGCCGTCGa                                                                               >  1:1520806/1‑137 (MQ=255)
                                                                                      cTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCa                                                                                                                                      >  1:1812095/1‑76 (MQ=255)
                                                                                      cTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCa                                                                                                                                      <  2:1812095/76‑1 (MQ=255)
                                                                                      cTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCTTCCCGGGCCAACGCGGCCAACTGGCCGTCGACGTCCGGCGCCACCAcgcg                                                            >  2:875250/1‑150 (MQ=255)
                                                                                       tGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCTTCCCGGGCCAACGCGGCCAACTGGCCGTCGACGTCCGGCGCCACCACgcgc                                                           >  1:874862/1‑150 (MQ=255)
                                                                                       tGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGACTGATAGCCACGCACCAGCACCTCACCGCCACCTTCCCGGGCCAACGCGGCCAACTGGCCGTCGACGTCCGGCGCCACCACgcgc                                                           >  1:1115586/1‑150 (MQ=255)
                                                                                        gggCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCTTCCCGGGCCAACGCGGCCAACTGGCCGTCGACGTCCGGCGCCACCACGCgca                                                          >  2:436367/1‑150 (MQ=255)
                                                                                        gggCATTCAGCCCAGGGCCGTCAGTCGCCGCGATCACCAGCCGGCCCCGGATCAGGTCGGCTGCCTGAAAGCCACGCACCAGCACCTCACCGCCACCTTCCCGGGCCAACGCGGCCAACTGGCCGGCGACGTCCGGCGCCACCACgcgcc                                                          >  2:1708379/1‑149 (MQ=255)
                                                                                              tCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCTTCCCGGGCCAACGCGGCCAACTGGCCGTCGACGTCCGGCGCCACCACGCGCAGCAcgc                                                    >  1:1524116/1‑150 (MQ=255)
                                                                                              tCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCTTCCCGGGCCAACGCGGCCAACTGGCCGTCGACGTCCGGCGCCACCACGCGCAGCAcgc                                                    <  1:436367/150‑1 (MQ=255)
                                                                                                 gCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCTTCCCGGGCCAACGCGGCCAACTGGCCGTCGACGTCCGGCGCCACCACGCgca                                                          <  2:476891/141‑1 (MQ=255)
                                                                                                 gCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCTTCCCGGGCCAACGCGGCCAACTGGCCGTCGACGTCCGGCGCCACCACGCgca                                                          >  1:476891/1‑141 (MQ=255)
                                                                                                   ccAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCTTCCCGGGCCAACGCGGCCAACTGGCCGTCGACGTCCGGCGCCACCACGCGCAGCACGccgccg                                               >  1:1678581/1‑150 (MQ=255)
                                                                                                       ggTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCTTCCCGGGCCAACGCGGCCAACTg                                                                                       <  1:286259/106‑1 (MQ=255)
                                                                                                       ggTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCTTCCCGGGCCAACGCGGCCAACTg                                                                                       >  2:286259/1‑106 (MQ=255)
                                                                                                        gtcgtcAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCTTCCCGGGCCAACGCGGCCAACTGGCCGTCGACGTCCGGCGCCACCACGCGCAGCACGCCGCCGGCATc                                          <  2:1502004/150‑1 (MQ=255)
                                                                                                            tcAGTCGCCGCGATCACAAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCTTCCCGGGCCAACGCGGCCAACTGGCCGTCGACGTCCGGCGCCACCACGCGCAGCACGCCGCCGGCATCGGcc                                      >  2:1551722/1‑150 (MQ=255)
                                                                                                             cAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCTTCCCGGGCCAACGCGGCCAACTGGCCGTCGACGTCCGGCGCcaccac                                                               >  2:1879010/1‑124 (MQ=255)
                                                                                                             cAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCTTCCCGGGCCAACGCGGCCAACTGGCCGTCGACGTCCGGCGCcaccac                                                               <  1:1879010/124‑1 (MQ=255)
                                                                                                                     gcgATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCTTCCCGGGCCAACGCGGCCAACTGGCCGTCGACGTCCGGCGCCACCACGCGCAGCACGCCGCCGGCATCGGCCAGCAGgcgc                             <  2:1029220/150‑1 (MQ=255)
                                                                                                                       gATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCTTCCCGGGCCAACGCGGCCAACTGGCCGTCGACGTCCGGCGCCACCACGCGCAGCACGCCGCCGGCATCGGCCAGCAGgcgcgc                           >  2:1224694/1‑150 (MQ=255)
                                                                                                                         tCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCTTCCCGGGCCAAc                                                                                                  <  2:1739303/77‑1 (MQ=255)
                                                                                                                         tCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCTTCCCGGGCCAAc                                                                                                  >  1:1739303/1‑77 (MQ=255)
                                                                                                                           aCCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCTTCCCGGGCCAACGCGGCCAACTGGCCGTCGACGTCCGGCGCCACCACGCGCAGCACGCCGCCGGCATCGGCCAGCAGGCGCGCCTTg                       <  1:241663/150‑1 (MQ=255)
                                                                                                                           aCCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCTTCCCGGGCCAACGCGGCCAACTGGCCGTCGACGTCCGGCGCCACCACGCGCAGCACGCCGCCGGCATCGGCCAGCAGGCGCGCCTTg                       <  1:1432018/150‑1 (MQ=255)
                                                                                                                             cAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCTTCCCGGGCCAACGCGGCCAACTGGCCGTCGACGTCCGGCGCCACCACGCGCAGCACGCCGCCGGCATCGGCCAGCAGGCGCGCCTTgcg                     >  1:136864/1‑150 (MQ=255)
                                                                                                                              aGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCTTCCCGGGCCAACGCGGCCAACTGGCCGTCGACGTCCGGCGCCACCACGCGCAGCACGCCGCCGGCATCGGCCAGCAGGCGCGCCTTgcgc                    >  1:95164/1‑150 (MQ=255)
                                                                                                                                ccGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCTTCCCGGGCCAACGCGGCCAACTGGCCGTCGACGTCCGGCGCCACCACGCGCAGCACGCCGCCGGCATCGGCCAGCAGGCGCGCCTTGCGCaa                  <  1:77943/150‑1 (MQ=255)
                                                                                                                                 cGGCAACCGACCAGGTCGGCTGTCTGATAGCCACGCACCAGCACCTCACCGCCACCTTCCCGGGCCAACGCGGCCAACTGGCCGTCGACGTCCGGCGCCACCACGCGCAGCACGCCGCCGGCATCGGCCAGCAGGCGCGCCTTGCGCAAg                 <  1:1604346/150‑1 (MQ=255)
                                                                                                                                  ggCAGCCGACCTGGTCGGTTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCTTCCCGGGCCAACGCGGCCAACTGGCCGTCGAc                                                                              >  1:82452/1‑88 (MQ=39)
                                                                                                                                  ggCAGCCGACCTGGTCGGTTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCTTCCCGGGCCAACGCGGCCAACTGGCCGTCGAc                                                                              <  2:82452/88‑1 (MQ=39)
                                                                                                                                  ggCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCTTCCCGGGCCAACGCGGCCAACTGGCCGTCGACGTCCGGCGCCACCACGCGCAGCACGCCGCCGGCATCGGCCAGCAGGCGCGCCTTGCGCAAgg                <  1:690343/150‑1 (MQ=255)
                                                                                                                                  ggCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCTTCCCGGGCCAACGCGGCCAACTGGCCGTCGACGTCCGGCGCCACCACGCGCAGCACGCCGCCGGCATCGGCCAGAAGGCGCGCCTTGCGCAAgg                >  1:932912/1‑150 (MQ=255)
                                                                                                                                   gCAGCCGACCTGGTCGGTTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCTTCCCGGGCCAACGCGGCCAACTGGCCGTCGACGTCCGGCGCCACCACGCGCAGCACGCCGCCGGCATCGGCCAGCAGGCGCGCCTTGCGCAAGGc               >  2:15267/1‑150 (MQ=255)
                                                                                                                                   gCAGCCGACCTGGTCGGTTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCTTCCCGGGCCAACGCGGCCAACTGGCCGTCGACGTCCGGCGCCAAAGCGCGCAGCACGCCGCCGGCATCGGCCAGCAGGCGCGCCTTGCGCAAGGc               >  1:592864/1‑150 (MQ=255)
                                                                                                                                   gCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCTTCCCGGGCCAACGCGGCCAACTGGCCGTCGACGTCCGgc                                                                       <  1:1579589/94‑1 (MQ=255)
                                                                                                                                   gCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCTTCCCGGGCCAACGCGGCCAACTGGCCGTCGACGTCCGgc                                                                       >  2:1579589/1‑94 (MQ=255)
                                                                                                                                     agccGACCTGGTCGGTTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCTTCCCGGGCCAACGCGGCCAACTGGCCGTCGACGTCCGGCGcca                                                                   <  2:637611/94‑1 (MQ=255)
                                                                                                                                     agccGACCTGGTCGGTTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCTTCCCGGGCCAACGCGGCCAACTGGCCGTCGACGTCCGGCGcca                                                                   >  1:637611/3‑96 (MQ=255)
                                                                                                                                     aaCCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCTTCCCGGGCCAACGCGGCCAACTGGCGGGCTGGAAGAGGCTGATTCGCGCTTCACACTCCTCCGGCACCGGCCAGCAGGCGCGCCTTGCGCAAGGCGa             >  1:696940/1‑150 (MQ=255)
                                                                                                                                     aaCCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCTTCCCGGGCCAACGCGGCCAACTGGCCGTCGACGTCCGGCGCCACCACGCGCAGCACGCCGCCGGCATCGGCCAGCAGGCGCGCCTTGCGCAAGGCgc             <  1:1160468/150‑2 (MQ=255)
                                                                                                                                      aCCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCTTCCCGGGCCAACGCGGCCAACTGGCCGTCGACGTCCGGCGCCACCACGCGCAGCACGCCGCCGGCATCGGCCAGCAGGCGCGCCTTGCGCAAGGCGAt            >  2:710435/1‑150 (MQ=255)
                                                                                                                                      aCCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCTTCCCGGGCCAACGCGGCCAACTGGCCGATGTCAACCGTCTCTACCACGCGGAGCCCGCCGAAGGCTTCATGGCTCGGTCGCGCCGTGCGCAAGGCGAt            >  1:975597/1‑150 (MQ=255)
                                                                                                                                       ccGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCTACCGGGGCCAACGCTGCCAACTGGCCGTCGACGTCCGGCGCCACCACGCGCAGCACGCCGCCGGCATCGGCCAGCAGGCGCGCCTTGCGCAAGGCgcgc           <  2:731202/150‑4 (MQ=255)
                                                                                                                                       ccGACCAGGGCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCTTCCCGGGCCAACGCGGCCAACTGGCCGTCGACGTCCGGCGCCACCACGCGCAGCACGCCGCCGGCATCGGCCAGCAGGCGCGCCTTGCGCAAGGCGAtc           >  1:85919/1‑150 (MQ=255)
                                                                                                                                         gACCTGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCTTCCCGGGCCAACGCGGCCAACTGGCCGTCGACGt                                                                            >  2:324170/1‑83 (MQ=255)
                                                                                                                                         gACCTGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCTTCCCGGGCCAACGCGGCCAACTGGCCGTCGACGt                                                                            <  1:324170/83‑1 (MQ=255)
                                                                                                                                         gACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCTTCCCGGGCCAACGCGGCCAACTGGCCGTCGACGTCCGGCGCCACCACGCGCAGCACGCCGCCGGCATCGGCCAGCAGGCGCGCCTTGCGCAAGGCGAtctc         >  2:651563/1‑150 (MQ=255)
                                                                                                                                           ccAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCTTCCCGGGCCAACGCGGCCAACTGGCCGTCGACGTCCGGCGCCACCACGCGCAGCACGCCGCCGGCATCGGCCAGCAGGCGCGCCTTGCGCAAGGCGATCTCAc       <  1:1212847/150‑1 (MQ=255)
                                                                                                                                            cAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCTTCCCGGGCCAACGCGGCCAACTGGCCGTCGACGTCCGGCGCCACCACGCGCAGCACGCCGCCGGCATCGGCCAGCAGGCGCGCCTTGCGCAAGGCGATCTCAcc      <  1:651563/150‑1 (MQ=255)
                                                                                                                                                gcggcggCCGGATAGCCACGCACCAGCACCTCACCGCCCCCTTCCCGGGCCAACGCGGCCAACTGGCCGTCGACGCCCGGCGCCACCACGCGCAGCACGCCGCCGGCATCGGCCAGCAGGCGCGCCTTGCGCAAGGCGATCTCACCGCCa  <  2:696940/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
ACGATCGCCGGGAAGATCACCGTACACAGGGCCGGCGCATCCACCACATTGACCGGCAGGCTAAGGGCCTGCGCATCGGCCGACACCTGGGCATTCAGCCCAGGGTCGTCAGTCGCCGCGATCACCAGCCGGCAACCGACCAGGTCGGCTGCCTGATAGCCACGCACCAGCACCTCACCGCCACCTTCCCGGGCCAACGCGGCCAACTGGCCGTCGACGTCCGGCGCCACCACGCGCAGCACGCCGCCGGCATCGGCCAGCAGGCGCGCCTTGCGCAAGGCGATCTCACCGCCA  >  NC_002947/4508475‑4508768

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: