Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,492,607:1 +C 5.1% coding (1457/1698 nt) cobJ ← precorrin‑3B C17‑methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,492,6071.C 5.1% 318.8 / 10.4 78coding (1457/1698 nt)cobJprecorrin‑3B C17‑methyltransferase
Reads supporting (aligned to +/- strand):  ref base . (35/39);  new base C (2/2);  total (37/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.06e-01

CCATGGTGCGCATGTCGACCATTTCTGGCAGCAGTTCGGCCAAGGTGACTACCTTCAGCGTCTGCCCCGGCCTACCGATGTCACGGCCAAGGACTACAGGCGTGTTCCCGTCACGATGTCGGCGCACAACCTCCAGGGCCACACCGAGCTGGTGCGGCCTGGCCTTGGAGATCGGGTTGTAGAAGGCCAGCACCAGGTCGGCCTGGGCAGCCAGGTCCAGGCGCTTTTCAATGATCGACCAGGGCTTGAGGTTGTCCGACAGCGACATCACGCAGAAGTCATGCCCCAG  >  NC_002947/5492458‑5492746
                                                                                                                                                      |                                                                                                                                           
ccATGGTGCGCATGTCGACCATTTCTGGCAGCAGTTCGGCCAAGGTGACTACCTTCAGCGTCTGCCCCGGCCTACCGATGTCACGGCCAAGGACTACAGGCGTGTTCCCGTCACGATGTCGGCGCACAACCTCCAGGGCCACACCGAGCt                                                                                                                                              >  2:1875662/1‑150 (MQ=255)
ccATGGTGCGCATGTCGACCATTTCTGGCAGCAGTTCGGCCAAGGTGACTACCTTCAGCGTCTGCCCCGGCCTACCGATGTCACGGCCAAGGACTACAGGCGTGTTCCCGTCACGATGTCGGCGCACAACCTCCAGGGCCACACCGAGCt                                                                                                                                              >  1:1558435/1‑150 (MQ=255)
     gTGCGCATGTCGACCATTTCTGGCAGCAGTTCGGCCAAGGTGACTACCTTCAGCGTCTGCCCCGGCCTACCGATGTCACGGCCAAGGACTACAGGCGTGTTCCCGTCACGATGTCGGCGCACAACCTCCAGGGCCACACCGAGCTGGTGc                                                                                                                                        >  1:1334374/1‑150 (MQ=255)
        cgcATGTCGACCATTTCTGGCAGCAGTTCGGCCAAGGTGACTACCTTCAGCGTCTGCCCCGGCCTACCGATGTCACGGCCAAGGACTACAGGCGTGTTCCCGTCACGATGTCGGCGCACAACCTCCAGGGCCACACCGAGCTGGTGCggc                                                                                                                                     <  2:1684729/150‑1 (MQ=255)
            tGTCGACCATTTCTGGCAGCAGTTCGGCCAAGGTGACTACCTTCAGCGTCTGCCCCGGCCTACCGATGTCACGGCCAAGGACTACAGGCGTGTTCCCGTCACGATGTCGGCGCACAACCTCCAGGGCCACACCGAGCTGGTGCggcctgg                                                                                                                                 >  2:432954/1‑150 (MQ=255)
               cGACCATTTCTGGCAGCAGTTCGGCCAAGGTGACTACCTTCAGCGTCTGCCCCGGCCTACCGATGTCACGGCCAAGGACTACAGGCGTGTTCCCGTCACGATGTCGGCGCACAACCTCCAGGGCCACACCGAGCTGGTGCggcctggcct                                                                                                                              <  2:194281/150‑1 (MQ=255)
                gACCATTTCTGGCAGCAGTTCGGCCAAGGTGACTACCTTCAGCGTCTGCCCCGGCCTACCGATGTCACGGCCAAGGACTGCTCGCGTGTTCCCGTCACGATGTCGGCGCACAACCTCCAGGGCCACACCGAGCTGGTGCGGCCTGGCCtt                                                                                                                             >  1:1861385/1‑150 (MQ=255)
                gACCATTTCTGGCAGCAGTTCGGCCAAGGTGACTACCTTCAGCGTCTGCCCCGGCCTACCGATGTCACGGCCAAGGACTACAGGCGTGTTCCCGTCACGATGTCGGCGCACAACCTCCAGGGCCACACCGAGCTg                                                                                                                                            >  2:146916/1‑135 (MQ=255)
                gACCATTTCTGGCAGCAGTTCGGCCAAGGTGACTACCTTCAGCGTCTGCCCCGGCCTACCGATGTCACGGCCAAGGACTACAGGCGTGTTCCCGTCACGATGTCGGCGCACAACCTCCAGGGCCACACCGAGCTg                                                                                                                                            <  1:146916/135‑1 (MQ=255)
                 aCCATTTCTGGCAGCAGTTCGGCCAAGGTGACTACCTTCAGCGTCTGCCCCGGCCTACCGATGTCACGGCCAAGGACTACAGGCGTGTTCCCGTCACGGTGTCGGCGCACAACCTCCAGGGGCACACCGAGCTGGTGCggcctggcactg                                                                                                                            >  1:676378/1‑146 (MQ=255)
                   cATTTCTGGCAGCAGTTCGGCCAAGGTGACTACCTTCAGCGTCTGCCCCGGCCTACCGATGTCACGGCCAAGGACTACAGTCGTGTTCCCGTCACGATGTCGGCGCACAACCTCCAGGGCCACACCGAGCTGGTGCGGCCTGGCCTTGga                                                                                                                          >  2:1402649/1‑150 (MQ=255)
                         tGGCAGCAGTTCGGCCAAGGTGACTACCTTCAGCGTCTGCCCCGGCCTACCGATGTCACGGCCAAGGACTACAGGCGTGTTCCCGTCACGATGTCGGCGCACAACCTCCAGGGCCACACCGAGCTGGTGCGGCCTGGCCTTGGAGATCgg                                                                                                                    >  2:30984/1‑150 (MQ=255)
                          ggCAGCAGTTCGGCCAAGGTGACTACCTTCAGCGTCTGCCCCGGCCTACCGATGTCACGGCCAAGGACTACAGGCGTGTTCCCGTCACGATGTCGGCGCACAACCTCCAGGGCCACACCGAGCTGGTGCGGCCTGGCCTTGGAGATCggg                                                                                                                   <  1:30984/150‑1 (MQ=255)
                           gcagcaGTTCGGCCAAGGTGACTACCTTCAGCGTCTGCCCCGGCCTACCGATGTCACGGCCAAGGACTACAGGCGTGTTCCCGTCACGATGTCGGCGCACAACCTCCAGGGCCACACCGAGCTGGTGCGGCCTGGCCTTGGAGATCGGGt                                                                                                                  <  2:801911/150‑1 (MQ=255)
                               cagTTCGGCCAAGGTGACTACCTTCAGCGTCTGCCCCGGCCTACCGATGTCACGGCCAAGGACTACAGGCGTGTTCCCGTCACGATGTCGGCGCACAACCTCCAGGGCCACACCGAGCTGGTGCGGCCTGGCCTTGGAGATCGGGTTGTa                                                                                                              <  2:1558435/150‑1 (MQ=255)
                                  ttCGGCCAAGGTGACTACCTTCAGCGTCTGCCCCGGCCTACCGATGTCACGGCCAAGGACTACAGGCGTGTTCCCGTCACGATGTCGGCGCACAACCTCCAGGGCCACACCGAGCTGGTGCGGCCTGGCCTTGGAGATCGGGTTGTAGaa                                                                                                           <  1:502953/150‑1 (MQ=255)
                                             tGACTACCCTCAGCGCTGGCCCCGGCCTACCGATGTCACGGCCAAGGACTACAGGAGTATCGGCACCCCGATGTCGGCCAACAACCTCCCGGGCCAAAAAGCGCTGGTTTGGCCTGGCCTTGGAGATCGGGTTGTAGAAGGCCAGCACCa                                                                                                <  1:1875662/150‑1 (MQ=255)
                                              gACTACCTTCAGCGTCTGCCCCGGCCTACCGATGTCACGGCCAAGGACTACAGGCGTGTTCCCGTCACGATGTCGGCGCACAACCTCCAGGGCCACACCGAGCTGGTGCGGCCTGGCCTTGGAGATCGGGTTGTAGAAGGCCAGCACCAg                                                                                               >  1:373283/1‑150 (MQ=255)
                                                 tACCTTCAGCGTCTGCCCCGGCCTACCGATGTCACGGCCAAGGACTACAGGCGTGTTCCCGTCACGATGTCGGCGCACAACCTCCAGGGCCACACCGAGCTGGTGCGGCCTGGCCTTGGAGATCGGGTTGTAGAAGGCCAGCACCAGGTc                                                                                            >  2:217130/1‑150 (MQ=255)
                                                     ttCAGCGTCTGCCCCGGCCTACCGATGTCACGGCCAAGGACTACAGGCGTGTTCCCGTCACGATGTCGGCGCACAACCTCCAGGGCCACACCGAGCTGGTGCGGCCTGGCCTTGGAGATCGGGTTGTAGAAGGCCAGCACCAGGTCGGcc                                                                                        >  2:336854/1‑150 (MQ=255)
                                                     ttCAGCGTCTGCCCCGGCCTACCGATGTCACGGCCAAGGACTACAGGCGTGTTCCCCTCACGATGTCGGCGCACAACCTCCAGGGCCACACCGAGCTGGTGCGGCCTGGCCTTGGAGATCGGGTTGTAGAAGGCCAGCACCAGGTCGGcc                                                                                        >  1:1921671/1‑150 (MQ=255)
                                                       cAGCGTCTGCCCCGGCCTACCGATGTCACGGCCAAGGACTACAGGCGTGTTCCCGTCACGATGTCGGCGCACAACCTCCAGGGCCACACCGAGCTGGTGCGGCCTGGCCTTGGAGATCggg                                                                                                                   >  2:1392302/1‑121 (MQ=255)
                                                       cAGCGTCTGCCCCGGCCTACCGATGTCACGGCCAAGGACTACAGGCGTGTTCCCGTCACGATGTCGGCGCACAACCTCCAGGGCCACACCGAGCTGGTGCGGCCTGGCCTTGGAGATCggg                                                                                                                   <  1:1392302/121‑1 (MQ=255)
                                                        aGCGTCTGCCCCGGCCTACCGATGTCACGGCCAAGGACTACAGGCGTGTTCCCGTCACGATGTCGGCGCACAACCTCCAGGGCCACACCGAGCTGGTGCGGCCTGGCCTTGGAGATCGGGTTGTAGAAGGCCAGCACCAGGTCGGCCTgg                                                                                     <  1:357002/150‑1 (MQ=255)
                                                         gCGTCTGCCCCGGCCTACCGATGTCACGGCCAAGGACTACAGGCGTGTTCCCGTCACGATGTCGGCGCACAACCTCCAGGGCCACACCGAGGTGGTGCGGCCTGGCCTTGGAGATCGGGTTGTAGAAGGCCAGCACCAGGTCGGCCTggg                                                                                    >  1:313134/1‑150 (MQ=255)
                                                           gTCTGCCCCGGCCTACCGCTGTCACGGCCAAGGACTACAGGCGTGTTCCCGTCACGATGTCGGCGCACAACCTCCAGGGCCACACCGAGCTGGTGCGGCCTGGCCTTGGAGATCGGGTTGTAGAAGGCCAGCACCAGGTCGGCCTGGGCa                                                                                  <  1:1133253/150‑1 (MQ=255)
                                                           gTCTGCCCCGGCCTACCGATGTCACGGCCAAGGACTACAGGCGTGTTCCCGTCACGATGTCGGCGCACAACCTCCAGGGCCACACCGAGCTGGTGCGGCCTGGCCTTGGAGATCGGGTTGTAGAAGGCCAGCACCAGGTCGGCCTGGGCa                                                                                  <  2:1334374/150‑1 (MQ=255)
                                                           gTCTGCCCCGGCCTACCGATGTCACGGCCAAGGACTACAGGCGTGTTCCCCTCACGATGTCGGCGCACAACCTCCAGGGCCACACCGAGCTGGTGCGGCCTGGCCTTGGAGATCGGGTTGTAGAAGGCCAGCACCAGGTCGGCCTGGGCa                                                                                  <  2:1921671/150‑1 (MQ=255)
                                                            tCTGCCCCGGCCTACCGATGTCACGGCCAAGGACTACAGGCGTTTTCCCGTCACGATGTCGGCGCACAACCTCCAGGGCCACACCGAGCTGGTGCGGCCTGGCCTTGGAGATCGGGTTGTAGAAGGCCAGCACCAGGTCGGCCTGGGCAg                                                                                 >  1:101997/1‑150 (MQ=255)
                                                                ccccGGCCTACCGATGTCACGGCCAAGGACTACAGGCGTGTTCACGTCACGATGTCGGCGCACAACCTCCAGGGCCACACCGAGCTGGTGCGGCCTGGCCTTGGAGATCGGGTTGTAGAAGGCCAGCACCAGGTCGGCCTGGGCAGCCAg                                                                             <  2:397193/150‑1 (MQ=255)
                                                                 cccGGCCTGCCGATGTCACGGCCAAGGACTACAGGCGTGTTCCCGTCACGATGTCGGCGCACAACCTCCAGGGCCACACCGAGCTGGTGCGGCCTGGCCTTGGAGATCGGGTTGTAGAAGGCCAGCACCAGGTCGGCCt                                                                                       >  1:1580244/1‑139 (MQ=255)
                                                                 cccGGCCTGCCGATGTCACGGCCAAGGACTACAGGCGTGTTCCCGTCACGATGTCGGCGCACAACCTCCAGGGCCACACCGAGCTGGTGCGGCCTGGCCTTGGAGATCGGGTTGTAGAAGGCCAGCACCAGGTCGGCCt                                                                                       <  2:1580244/139‑1 (MQ=255)
                                                                 cccGGCCTACCGATGTCACGGCCAAGGACTACAGGCGTGTTCCCGTCACGATGTCGGCGCACAACCTCCAGGGCCACACCGAGCTGGTGCGGCCTGGCCTTGGAGATCGGGTTGTAGAAg                                                                                                          >  2:1022138/1‑120 (MQ=255)
                                                                 cccGGCCTACCGATGTCACGGCCAAGGACTACAGGCGTGTTCCCGTCACGATGTCGGCGCACAACCTCCAGGGCCACACCGAGCTGGTGCGGCCTGGCCTTGGAGATCGGGTTGTAGAAg                                                                                                          <  1:1022138/120‑1 (MQ=255)
                                                                 cccGGCCTACCAACGTCACGGCCAAGGACTACAGGCGTGTTCCCGTCACGATGTCGGCGCACAACCTCCAGGGCCACACCGAGCTGGTGCGGCCTGGCCTTGGAGATCGGGTTGTAGAAGGCCAGCACCAGGTCGGCCTGGGCAGCCAgg                                                                            <  1:390302/150‑1 (MQ=255)
                                                                  ccGGCCTACCGATGTCACGGCCAAGGACTACAGGCGTGTTCCCGTCACGATGTCGGCGCACAACCTCCAGGGCCACACCGAGCTGGTGCGGCCTGGCCTTGGAGATCGGGTTGTAGAAGGCCAGCACCAGGTCGGCCTGGGCAGCCAGGt                                                                           >  1:811974/1‑150 (MQ=255)
                                                                  ccGGCCTACCGATGTCACGGCCAAGGACTACAGGCGTGTTCCCGTCACGATGTCGGCGCACAACCTCCAGGGCCACACCGAGCTGGTGCGGCCTGGCCTTGGAGATCGGGTTGTAGAAGGCCAGCACCAGGTCGGCCTGGGCAGCCAGGt                                                                           <  1:1402649/150‑1 (MQ=255)
                                                                  ccGGCCTACCGATGTCACGGCCAAGGACTACAGGCGTGTTCCCGTCACGATGTCGGCGCACAACCTCCAGGGCCACACCGAGCTGGTGCGGCCTGGCCTTGGAGATCGGGTTGTAGAAGGCCAGCACCAGGTCGGCCTGGGCAGCCAGGt                                                                           <  2:664382/150‑1 (MQ=255)
                                                                  ccGGCCTACCGATGTCACGGCCAAGGACTACAGGCGTGTTCCCGTCACGATGTCGGCGCACAACCTCCAGGGCCACACCGAGCTGGTGCGGCCTGGCCTTGGAGATCGGGTTGTAGAAGGCCAGCACCAGGTCGGCCTGGGCAGCCAGGt                                                                           >  2:812642/1‑150 (MQ=255)
                                                                  ccGGCCTACCGATGTCACGGCCAAGGACTACAGGCATGTTCCCGTCACGATGTCGGCGCACAACCTCCAGGGCCACACCGAGCTGGTGCGGCCTGGCCTTGGAGATCGGGTTGTAGAAGGCCAGCACCAGGTCGGCCTGGGCAGCCAGGt                                                                           >  1:151341/1‑150 (MQ=255)
                                                                    ggCCTACCGAGGTCACGGCCAAGGACTACAGGCGTGTCCCCGTCACGATGTCGGCGCACAACCTCCAGCGCCACACCGAGCTGGTGCGGCCTGGCCTTGGAGATCGGGTTGTAGAAGGCCAGCACCAGCTCGGCCTAGACAGCCAGGTcc                                                                         <  2:1861385/150‑1 (MQ=255)
                                                                      ccTACCGATGTCACGGCCAAGGACTACAGGCGTGTTCCCGTCACGATGTCGGCGCACAACCTCCAGGGCCACACCGAGCTGGTGCGGCCTGGCCTTGGAGATCGGGTTGTAGAAGGCCAGCACCAGGTCGGCCTGGGCAGCCAGGTCCAg                                                                       <  1:851015/150‑1 (MQ=255)
                                                                           cGATGTCACGGCCAAGGACTACAGGCGTGTTCCCGTCACGATGTCGGCGCACAACCTCCAGGGCCACACCGAGCTGGTGCGGCCTGGCCTTGGAGATCGGGTTGTAGAAGGCCAGCACCAGGTCGGCCTGGGCAGCCAGGTCCAGGCGCt                                                                  >  1:111115/1‑150 (MQ=255)
                                                                           cGATGTCACGGCCAAGGACTACAGGCGTGTTCCCGTCACGATGTCGGCGCACAACCTCCAGGGCCACACCGAGCTGGTGCGGCCTGGCCTTGGAGATCGGGTTGTAGAAGGCCAGCACCAGGTCGGCCTGGGCAGCCAGGTCCAGGCGCt                                                                  >  1:837427/1‑150 (MQ=255)
                                                                            gATGTCACGGCCAAGGACTACAGGCGTGTTCCCGTCACGATGTCGGCGCACAACCTCCAGGGCCACACCGAGCTGGTGCGGCCTGGCCTTGGAGATCg                                                                                                                     <  1:121600/98‑1 (MQ=255)
                                                                            gATGTCACGGCCAAGGACTACAGGCGTGTTCCCGTCACGATGTCGGCGCACAACCTCCAGGGCCACACCGAGCTGGTGCGGCCTGGCCTTGGAGATCg                                                                                                                     >  2:121600/1‑98 (MQ=255)
                                                                            gATGTCACGGCCAAGGACTACAGGCGTGTTCCCGTCACGATGTCGGCGCACAACCTCCAGGGCCACACCGAGCTGGTGCGGCCTGGCCTTGGAGATCGGGTTGTAGAAGGCCAGCACCAGGTCGGCCTGGGCAGCCAGGTCCAGGCGCtt                                                                 >  2:353945/1‑150 (MQ=255)
                                                                             aTGTCACGGCCAAGGACTACAGGCGTGTTCCCGTCACGATGTCGGCGCACAACCTCCAGGGCCACACCGAGCTGGTGCGGCCTGGCCTTGGAGATCGGGTTGTAGAAGGCCAGCACCAGGTCGGCCTGGGCAGCCAGGTCCAGGCGCttt                                                                <  1:336854/150‑1 (MQ=255)
                                                                             aTGTCACGGCCAAGGACTACAGGCGTGTTCCCGTCACGATGTCGGCGCACAACCTCCAGGGCCACACCGAGCTGGTGCGGCCTGGCCTTGGAGATCGGGTTGTAGAAGGCCAGCACCAGGTCGGCCTGGGCAGCCAGGTCCAGGCGCttt                                                                <  1:217130/150‑1 (MQ=255)
                                                                                      ccAAGGACTACAGGCGTGTTCCCGTCACGATGTCGGCGCACAACCTCCAGGGCCACACCGAGCTGGTGCGGCCTGGCCTTGGAGATCGGGTTGTAGAAGGCCAGCACCAGGTCGGCCTGGGCAGCCAGGTCCAGGCGCTTTTCAATGATc                                                       <  2:837427/150‑1 (MQ=255)
                                                                                       cAAGGACTACAGGCGTGTTCCCGTCACGATGTCGGCGCACAACCTCCAGGGCCACACCGAGCTGGTGCGGCCTGGCCTTGGAGATCGGGTTGTAGAAGGCCAGCACCAGGTCGGCCTGGGCAGCCAGGTCCAGGCGCTTTTCAATGATCg                                                      >  2:1433197/1‑150 (MQ=255)
                                                                                        gatctcTACAGGCGTGTCCCCGTCACGATGTCGGCGCACAACCTCCAGCACCACACCGAGCTGGTGCGGCCTGGCCTTGGTGATCGGGTTGTAGAAGGCCAGCACCAGGTCGGCCTGGGCAGCCAGGTCCAGGCGCTTTTCAATGATCGa                                                     <  2:980702/145‑1 (MQ=255)
                                                                                        aaGGACTACAGGCGTGTTCCCGTCACGATGTCGGCGCACAACCTCCAGGGCCACACCGAGCTGGTGCGGCCTGGCCTTGGAGATCGGGTTGTAGAAGGCCAGCACCAGGTCGGCCTGGGCAGCCAGGTCCAGGCGCTTTTCAATGATCGa                                                     >  1:130313/1‑150 (MQ=255)
                                                                                         aGGACTACAGGCGTTTTCCCGTCACGATGTCGGCGCACAACCTCCAGGGCCACACCGAGCTGGTGCGGCCTGGCCTTGGAGATCGGGTTGTAGAAGGCCAGCACCAGGTCGGCCTGGGCAGCCAGGTCCAGGCGCTTTTCAATGATCGAc                                                    <  2:101997/150‑1 (MQ=255)
                                                                                         aGGACTACAGGCGTGTTCCCGTCACGATGTCGGCGCACAACCTCCAGGGCCACACCGAGCTGGTGCGGCCTGGCCTTGGAGATCGGGTTGTAGAAGGCCAGCACCAGGTCGGCCTGGGCAGCCAGGTCCAGGCGCTTTTCAATGATCGAc                                                    <  2:130313/150‑1 (MQ=255)
                                                                                         aGGACTACAGGCGTGTTCCCGTCACGATGTCGGCGCACAACCTCCAGGGCCACACCGAGCTGGTGCGGCCTGGCCTTGGAGATCGGGTTGTAGAAGGCCAGCACCAGGTCGGCCTGGGCAGCCAGGTCCAGGCGCTTTTCAATGATCGAc                                                    >  1:1192648/1‑150 (MQ=255)
                                                                                          ggACTACAGGCGTGTTCCCGTCACGATGTCGGCGCACAACCTCCAGGGCCACACCGAGCTGGTGCGGCCTGGCCTTGGAGATCGGGTTGTAGAAGGCCAGCACCAGGTCGGCCTGGGCAGCCAGGTCCAGGCGCTTTTCAATGCTc                                                       >  1:1183500/1‑146 (MQ=255)
                                                                                          ggACTACAGGCGTGTTCCCGTCACGATGTCGGCGCACAACCTCCAGGGCCACACCGAGCTGGTGCGGCCTGGCCTTGGAGATCGGGTTGTAGAAGGCCAGCACCAGGTCGGCCTGGGCAGCCAGGTCCAGGCGCTTTTCAATGATc                                                       <  2:1183500/146‑1 (MQ=255)
                                                                                              tACAGGCGTGTTCCCGTCACGATGTCGGCGCCCAACCTCCAGGGCCACACCGAGCTGGTGCGGCCTGGCCTTGGAGATCGGGTTGTAGAAGGCCAGCACCAGGTCGGCCTGGGCAGCCAGGTCCAGGCGCTTTTCAATGATCGAcagatc                                               >  1:980702/1‑145 (MQ=255)
                                                                                              tACAGGCGTGTTCCCGTCACGATGTCGGCGCACAACCTCCAGGGCCACACCGAGCTGGTGCGGCCTGGCCTTGGAGATCGGGTTGTAGAAGGCCAGCACCAGGTCGGCCTGGGCAGCCAGGTCCAGGCGCTTTTCAATGATCGACCAggg                                               >  1:813440/1‑150 (MQ=255)
                                                                                                 aGGCGTGTTCCCGTCACGATGTCGGCGCACAACCTCCAGGGCCACACCGAGCTGGTGCGGCCTGGCCTTGGAGATCGGGTTGTAGAAGGCCAGCACCAGGTCGGCCTGGGCAGCCAGGTCCAGGCGCTTTTCAATGATCGACCAGGGCtt                                            >  2:957483/1‑150 (MQ=255)
                                                                                                   gCGTGTTCCCGTCACGATGTCGGCGCACAACCTCCAGGGCCACACCGAGCTGGTGCGGCCTGGCCTTGGAGATCGGGTTGTAGAAGGCCAGCACCAGGTCGGCCTGGGCAGCCAGGTCCAGGCGCTTTTCaa                                                            <  1:1700541/132‑1 (MQ=255)
                                                                                                   gCGTGTTCCCGTCACGATGTCGGCGCACAACCTCCAGGGCCACACCGAGCTGGTGCGGCCTGGCCTTGGAGATCGGGTTGTAGAAGGCCAGCACCAGGTCGGCCTGGGCAGCCAGGTCCAGGCGCTTTTCaa                                                            >  2:1700541/1‑132 (MQ=255)
                                                                                                    cGTGTCCCCGTCACGATGTCGGCGCCCAACCTCCAGGGCCACACCCAGCTGGTGCGGCCTGGCCTTGGAGATCGGGTTGGAGAAGGCCAGCACCAGGTCGGCCTGGGCAGCAAGGTCCAGGCGCTTTTCAATGATCGACCAGGGCTTAAg                                         <  1:1715979/150‑1 (MQ=255)
                                                                                                      tgtTCCCGTCACGATGTCGGCGCACAACCTCCAGGGCCACACCGAGCTGGTGCGGCCTGGCCTTGGAGATCGGGTTGTAGAAGGCCAGCACCAGGTCGGCCTGGGCAGCCAGGTCCAGGCGCTTTTCAATGATCGACCAGGGCTTGAGGt                                       <  2:811974/150‑1 (MQ=255)
                                                                                                       gtTCCCGTCACGATGTCGGCGCACAACCTCAAGGAACACACCGAGCTGGTGCGGCCTGGCCTTGGAGATCGGGTTGTAGAAGGCCAGCACCAGGTCGGCCTGGGCAGCCAGGTCCAGGCGCTTTTCAATGATCGACCAGGGCTTGAGGtt                                      <  2:151341/150‑1 (MQ=255)
                                                                                                        ttCCCGTCACGATGTCGGCGCACAACCTCCAGGGCCACACCGAGCTGGTGCGGCCTGGCCTTGGAGATCGGGTTGTAGAAGGCCAGCACCAGGTCGGCCTGGGCAGCCAGGTCCAGGCGCTTTTCAATGATGGACCAGGGCTTGAGGTTg                                     >  2:1786803/1‑150 (MQ=255)
                                                                                                        ttCCCGTCACGATGTCGGCGCACAACCTCCAGGGCCACACCGAGCTGGTGCGGCCTGGCCTTGGAGATCGGGTTGTAGAAGGCCAGCACCAGGTCGGCCTGGGCAGCCAGGTCCAGGCGCTTTTCAATGATCGACCAGGGCTTGAGGTTg                                     >  2:1515123/1‑150 (MQ=255)
                                                                                                        ttCCCGTCACGATGTCGGCGCACAACCTCCAGGGCCACACCGAGCTGGTGCGGCCTGGCCTTGGAGATCGGGTTGTAGAAGGCCAGCACCAGATCGGCCTGGGCAGCCAGGTCCAGGCGCTTTTCAATGATCGACCAGGGCTTGAGGTTg                                     >  2:1308063/1‑150 (MQ=255)
                                                                                                           ccGTCACGATGTCGGCGCACAACCTCCAGGGCCACACCGAGCTGGTGCGGCCTGGCCTTGGAGATCGGGTTGTAGAAGGCCAGCACCAGGTCGGCCTGGGCAGCCAGGTCCAGGCGCTTTTCAATGATCGACCAGGGCTTGAGGTTGTcc                                  <  1:1433197/150‑1 (MQ=255)
                                                                                                                        ggCGCACAACCTCCAGGGCCACACCGAGCTGGTGCGGCCTGGCCTTGGAGATCGGGTTGt                                                                                                               <  1:1061396/60‑1 (MQ=255)
                                                                                                                        ggCGCACAACCTCCAGGGCCACACCGAGCTGGTGCGGCCTGGCCTTGGAGATCGGGTTGt                                                                                                               >  2:1061396/1‑60 (MQ=255)
                                                                                                                             acaACCTCCAGGGCCACACCGAGCTGGTGCGGCCTGGCCTTGGAGATCGGGTTGTAGAAGGCCAGCACCAGGTCGGCCTGGGCAGCCAGGTCCAGGCGCTTTTCAATGATCGACCAGGGCTTGAGGTTGTCCGACAGCGACATCACGCAg                <  1:957483/150‑1 (MQ=255)
                                                                                                                                caaggccaGGCCGCACCAGCTCGGTGTGGCCTGGCCTTGGAGATCGGGTTGTAGAAGGCCAGCACCAGGTCGGCCt                                                                                       <  2:548795‑M1/68‑1 (MQ=255)
                                                                                                                                caaggccaGGCCGCACCAGCTCGGTGTGGCCTGGCCTTGGAGATCGGGTTGTAGAAGGCCAGCACCAGGTCGGCCt                                                                                       >  1:548795‑M1/9‑76 (MQ=255)
                                                                                                                                  aggccaGGCCGCACCAGCTCGGTGTGGCCTGGCCTTGGAGATCGGGTTGTAGAAGGCCAGCACCAg                                                                                               >  2:553947‑M1/7‑66 (MQ=255)
                                                                                                                                  aggccaGGCCGCACCAGCTCGGTGTGGCCTGGCCTTGGAGATCGGGTTGTAGAAGGCCAGCACCAg                                                                                               <  1:553947‑M1/60‑1 (MQ=255)
                                                                                                                                   tCCAGGGCCACACCGAGCTGGTGCGGCCTGGCCTTGGAGATCGGGTTGCAGAAGGCCAGCACCAGGTCGGCCTGGGCAGCCAGGTCCAGGCGCTTTTCAATGATCGACCAGGGCTTGAGGTTGTCCGACAGCGACATCACGCAGAAGTCa          >  1:96901/1‑150 (MQ=255)
                                                                                                                                      aGGGCCACACCGAGCTGGTGCGGCCTGGCCTTGGAGATCGGGTTGTAGAAGGCCAGCACCAGGTCGGCCTGGGCAGCCAGGTCCAGGCGCTTTTCAATGATCGACCAGGGCTTGAGGTTGTCCGACAGCGACATCACGCAGAAGTCATGc       <  1:353945/150‑1 (MQ=255)
                                                                                                                                      aGGGCCACACCGAGCTGGTGCGGCCTGGCCTTGGAGATCGGGTTGTAGAAGGCCAGCACCAGGTCGGCCTGGGCAGCCAGGTCCAGGCGCTTTTCAATGATCGACCAGGGCTTGAGGTTGTCCGACAGCGACATCACGCAGAAGTCATGc       <  2:813440/150‑1 (MQ=255)
                                                                                                                                         GCCACACCGAGCTGGTGCGGCCTGGCCTTGGAGATCGGGTTGTAGAAGGCCAGCACCAGGTCGGCCTGGGCAGCCAGGTCCAGGCGCTTTTCAATGATCGACCAGGGCTTGAGGTTGTCCGACAGCGACATCACGACAAAGTCATGcccc    <  1:1239917‑M1/150‑1 (MQ=255)
                                                                                                                                           CACACCGAGCTGGTGCGGCCTGGCCTTGGAGATCGGGTTGTAGAAGGCCAGCACCAGGTCGGCCTGGGCAGCCAGGTCGGCGCGCTTTTCAATGATCGACCAGGGCTTGAGGTTGTCCGACAGCGACATCACGCAGAAGTCATGCCCCAg  <  1:797668‑M1/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                           
CCATGGTGCGCATGTCGACCATTTCTGGCAGCAGTTCGGCCAAGGTGACTACCTTCAGCGTCTGCCCCGGCCTACCGATGTCACGGCCAAGGACTACAGGCGTGTTCCCGTCACGATGTCGGCGCACAACCTCCAGGGCCACACCGAGCTGGTGCGGCCTGGCCTTGGAGATCGGGTTGTAGAAGGCCAGCACCAGGTCGGCCTGGGCAGCCAGGTCCAGGCGCTTTTCAATGATCGACCAGGGCTTGAGGTTGTCCGACAGCGACATCACGCAGAAGTCATGCCCCAG  >  NC_002947/5492458‑5492746

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: