Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,492,598 A→G 6.0% V489A (GTG→GCG)  cobJ ← precorrin‑3B C17‑methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,492,5980AG 6.0% 240.1 / 3.7 83V489A (GTG→GCG) cobJprecorrin‑3B C17‑methyltransferase
Reads supporting (aligned to +/- strand):  ref base A (38/40);  new base G (3/2);  total (41/42)
Fisher's exact test for biased strand distribution p-value = 6.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.06e-01

ACCATGGTGCGCATGTCGACCATTTCTGGCAGCAGTTCGGCCAAGGTGACTACCTTCAGCGTCTGCCCCGGCCTACCGATGTCACGGCCAAGGACTACAGGCGTGTTCCCGTCACGATGTCGGCGCACAACCTCCAGGGCCACACCGAGCTGGTGCGGCCTGGCCTTGGAGATCGGGTTGTAGAAGGCCAGCACCAGGTCGGCCTGGGCAGCCAGGTCCAGGCGCTTTTCAATGATCGACCAGGGCTTGAGGTTGTCCGACAGCGACATCACGCAGAAGTCATGCCCCAG  >  NC_002947/5492457‑5492746
                                                                                                                                             |                                                                                                                                                     
tgcaccgctggggcatgacttctgcgtgatgtcgctgtcggacaacctcaagccctggtcgatcattgaaaagcgcctggacctggctgcccaggccgacctggtgctggccttctacaacccgatctccaaggccaGGCCGCACCGAGc                                                                                                                                               >  2:481462‑M1/138‑150 (MQ=255)
 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                                                                                                                   <  1:1392302/121‑1 (MQ=255)
                                                        cAGCGTCTGCCCCGGCCTACCGATGTCACGGCCAAGGACTACAGGCGTGTTCCCGTCACGATGTCGGCGCACAACCTCCAGGGCCACACCGAGCTGGTGCGGCCTGGCCTTGGAGATCggg                                                                                                                   >  2:1392302/1‑121 (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                                                                           >  2:812642/1‑150 (MQ=255)
                                                                   ccGGCCTACCGATGTCACGGCCAAGGACTACAGGCGTGTTCCCGTCACGATGTCGGCGCACAACCTCCAGGGCCACACCGAGCTGGTGCGGCCTGGCCTTGGAGATCGGGTTGTAGAAGGCCAGCACCAGGTCGGCCTGGGCAGCCAGGt                                                                           <  2:664382/150‑1 (MQ=255)
                                                                   ccGGCCTACCGATGTCACGGCCAAGGACTACAGGCGTGTTCCCGTCACGATGTCGGCGCACAACCTCCAGGGCCACACCGAGCTGGTGCGGCCTGGCCTTGGAGATCGGGTTGTAGAAGGCCAGCACCAGGTCGGCCTGGGCAGCCAGGt                                                                           <  1:1402649/150‑1 (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:837427/1‑150 (MQ=255)
                                                                            cGATGTCACGGCCAAGGACTACAGGCGTGTTCCCGTCACGATGTCGGCGCACAACCTCCAGGGCCACACCGAGCTGGTGCGGCCTGGCCTTGGAGATCGGGTTGTAGAAGGCCAGCACCAGGTCGGCCTGGGCAGCCAGGTCCAGGCGCt                                                                  >  1:111115/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                                                            >  2:1700541/1‑132 (MQ=255)
                                                                                                    gCGTGTTCCCGTCACGATGTCGGCGCACAACCTCCAGGGCCACACCGAGCTGGTGCGGCCTGGCCTTGGAGATCGGGTTGTAGAAGGCCAGCACCAGGTCGGCCTGGGCAGCCAGGTCCAGGCGCTTTTCaa                                                            <  1:1700541/132‑1 (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)
                                                                                                                                             |                                                                                                                                                     
ACCATGGTGCGCATGTCGACCATTTCTGGCAGCAGTTCGGCCAAGGTGACTACCTTCAGCGTCTGCCCCGGCCTACCGATGTCACGGCCAAGGACTACAGGCGTGTTCCCGTCACGATGTCGGCGCACAACCTCCAGGGCCACACCGAGCTGGTGCGGCCTGGCCTTGGAGATCGGGTTGTAGAAGGCCAGCACCAGGTCGGCCTGGGCAGCCAGGTCCAGGCGCTTTTCAATGATCGACCAGGGCTTGAGGTTGTCCGACAGCGACATCACGCAGAAGTCATGCCCCAG  >  NC_002947/5492457‑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: