Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,294,001 T→C 5.3% Q16R (CAG→CGG)  PP_3763 ← precorrin‑6A synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,294,0010TC 5.3% 275.0 / 3.9 94Q16R (CAG→CGG) PP_3763precorrin‑6A synthase
Reads supporting (aligned to +/- strand):  ref base T (33/56);  new base C (2/3);  total (35/59)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.08e-01

CCACCGGCACCAGGCGATAACCACCCTCGGGGCGGTAGCGCTGCAGAATGGCCTTGCGCAGGCGCACCAGTTCATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAACCCGCGACCATTGTCGGGCGGTTACCGGTATTTGGCAAAGCAGCTGCTTAGAAGATCGAGAAGTTGTAGCTGACAATGACCC  >  NC_002947/4293852‑4294149
                                                                                                                                                     |                                                                                                                                                    
ccaccGGCACCAGGCGATAACCACCCTCGGGGCGGTAGCGCTGCAGAATGGCCTTGCGCAGGCGCACCAGTTCATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCt                                                                                                                                                      >  1:437065/1‑150 (MQ=255)
  accGGCACCAGGCGATAACCACCCTCGGGGCGGTAGCGCTGCAGAATGGCCTTGCGCAGGCGCACCAGTTCATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTgg                                                                                                                                                    >  1:1618611/1‑150 (MQ=255)
   ccGGCACCAGGCGATAACCACCCTCGGGGCGGTAGCGCTGCAGAATGGCCTTGCGCAGGCGCACCAGTTCATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGgc                                                                                                                                                   >  2:1751875/1‑150 (MQ=255)
      gCACCAGGCGATAACCACCCTCGGGGCGGTAGCGCTGCAGAATGGCCTTGCGCAGGCGCACCAGTTCATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGcgcg                                                                                                                                                <  1:1469327/150‑1 (MQ=255)
       cACCAGGCGATAACCACCCTCGGGGCGGTAGCGCTGCAGAATGGCCTTGCGCAGGCGCACCAGTTCATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCgg                                                                                                                                               <  2:1618611/150‑1 (MQ=255)
       cACCAGGCGATAACCACCCTCGGGGCGGTAGCGCTGCAGAATGGCCTTGCGCAGGCGCACCAGTTCATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCgg                                                                                                                                               <  1:1390443/150‑1 (MQ=255)
         ccAGGCGATAACCACCCTCGGGGCGGTAGCGCTGCAGAATGGCCTTGCGCAGGCGCACCAGTTCATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGGt                                                                                                                                             >  2:994323/1‑150 (MQ=255)
          cAGGCGATAACCACCCTCGGGGCGGTAGCGCTGCAGAATGGCCTTGCGCAGGCGCACCAGTTCATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCgg                                                                                                                                               <  2:904437/147‑1 (MQ=255)
          cAGGCGATAACCACCCTCGGGGCGGTAGCGCTGCAGAATGGCCTTGCGCAGGCGCACCAGTTCATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCgg                                                                                                                                               >  1:904437/1‑147 (MQ=255)
          cAGGCGATAACCACCCTCGGGGCGGTAGCGCTGCAGAATGGCCTTGCGCAGGCGCACCAGTTCATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGGTc                                                                                                                                            <  1:1002780/150‑1 (MQ=255)
            ggCGATAACCACCCTCGGGGCGGTAGCGCTGCAGAATGGCCTTGCGCAGGCGCACCAGTTCATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGGTCAc                                                                                                                                          <  2:130915/150‑1 (MQ=255)
             gCGATAACCACCCTCGGGGCGGTAGCGCTGCAGAATGGCCTTGCGCAGGCGCACCAGTTCATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGGTCAcc                                                                                                                                         <  2:1653970/150‑1 (MQ=255)
             gCGATAACCACCCTCGGGGCGGTAGCGCTGCAGAATGGCCTTGCGCAGGCGCACCAGTTCATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGGTCAcc                                                                                                                                         <  2:1457035/150‑1 (MQ=255)
               gATAACCACCCCCGGGGCGGTAGCGCTGCAGAATGGCCTTGCGCAGGCGCACCAGTTCATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGGTCACCcg                                                                                                                                       >  2:1209236/1‑150 (MQ=255)
                  aaCCACCCTCGGGGCGGTAGCGCTGCAGAATGGCCTTGCGCAGGCGCACCAGTTCATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGGTCACCCgcgc                                                                                                                                    <  1:269972/150‑1 (MQ=255)
                  aaCCACCCTCGGGGCGGTAGCGCTGCAGAATGGCCTTGCGCAGGCGCACCAGTTCATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGGTCACCCgcgc                                                                                                                                    <  1:757356/150‑1 (MQ=255)
                  aaCCACCCTCGGGGCGGTAGCGCTGAAGTAGTGCCTTGCGCAGGCGCACCAGTTCATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGGTCCCCCgcgc                                                                                                                                    <  1:972588/150‑1 (MQ=255)
                   accaccCTCGGGGCGGTAGCGCTGCAGAATGGCCTTGCGCAGGCGCACCAGTTCATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGGTCACCCGCGcc                                                                                                                                   <  1:253413/150‑1 (MQ=255)
                    ccaccCTCGGGGCGGTAGCGCTGCAGAATGGCCTTGCGCAGGCGCACCAGTTCATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGGTGACCCGCGcca                                                                                                                                  <  2:982692/150‑2 (MQ=255)
                     caccCTCGGGGCGGTAGCGCTGCAGAATGGCCTTGCGCAGGCGCACCAGTTCATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGGTCACCCGCGCCGa                                                                                                                                 >  2:1104188/1‑150 (MQ=255)
                     caccCTCGGGGCGGTAGCGCTGCAGAATGGCCTTGCGCAGGCGCACCAGTTCATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGGTCACCCGCGCCGa                                                                                                                                 >  2:1185264/1‑150 (MQ=255)
                     caccCTCGGGGCGGTAGCGCTGCAGAATGGCCTTGCGCAGGCGCACCAGTTCATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGGTCACCCGCGCCGa                                                                                                                                 >  1:389140/1‑150 (MQ=255)
                     caccCTCGGGGCGGTAGCGCTGCAGAATGGCCTTGCGCAGGCGCACCAGTTCATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCCGGCGCGGGTGACCCGCGccag                                                                                                                                 <  1:815963/150‑3 (MQ=255)
                       cccTCGGGGCGGTAGCGCTGCAGAATGGCCTTGCGCAGGCGCACCAGTTCATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGGTCACCCGCGCCGATa                                                                                                                               <  1:1174179/150‑1 (MQ=255)
                        ccTCGGGGCGGTAGCGCTGCAGAATGGCCTTGCGCAGGCGCACCAGTTCATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGGCCGCTTCATAGGTGACCTGGCGCGGGTCACCCGCGCCGATAc                                                                                                                              <  1:1927490/150‑1 (MQ=255)
                        ccTCGGGGCGGTAGCGCTGCAGAATGGCCTTGCGCAGGCGCACCAGTTCATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATAGACCGCTTCATAGGTGACCTGGCGCGGGTCACCCGCGCCGATAc                                                                                                                              <  1:345549/150‑1 (MQ=255)
                         cTCGGGGCGGTAGCGCTGCAGAATGGCCTTGCGCAGGCGCACCAGTTCATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGGTCACCCGCGCCGATAcc                                                                                                                             <  1:1860453/150‑1 (MQ=255)
                           cGGGGCGGTAGCGCTGCAGAATGGCCTTGCGCAGGCGCACCAGTTCATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGTTCACCCGCGCCGATACCGa                                                                                                                           <  2:437065/150‑1 (MQ=255)
                            ggggCGGTAGCGCTGCAGAATGGCCTTGCGCAGGCGCACCAGTTCATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGAt                                                                                                                          <  1:729705/150‑1 (MQ=255)
                            ggggCGGTAGCGCTGCAGAATGGCCTTGCGCAGGCGCACCAGTTCATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGAATCGACCGCTTCATAGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGAt                                                                                                                          <  1:542106/150‑1 (MQ=255)
                               gCGGTAGCGCTGCAGAATGGCCTTGCGCAGGCGCACCAGTTCATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATcaa                                                                                                                       <  2:389140/150‑1 (MQ=255)
                                 ggTAGCGCTGCAGAATGGCCTTGCGCAGGCGCACCAGTTCATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCACGGGGCACCCGCGCGGATACCGATcaaca                                                                                                                     >  2:507325/1‑150 (MQ=255)
                                 ggTAGCGCTGCAGAATGGCCTTGCGCAGGCGCACCAGTTCATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACCCGTGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGGTCACCCGCGCCGATCCCGATcaaca                                                                                                                     >  1:127739/1‑150 (MQ=255)
                                  gTAGCGCTGCAGAATGGCCTTGCGCAGGCGCACCAGTTCATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATcaa                                                                                                                       <  1:1784831/147‑1 (MQ=255)
                                  gTAGCGCTGCAGAATGGCCTTGCGCAGGCGCACCAGTTCATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATcaa                                                                                                                       >  2:1784831/1‑147 (MQ=255)
                                   tAGCGCTGCAGAATGGCCTTGCGCAGGCGCACCAGTTCATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCaacaac                                                                                                                   <  1:1780715/150‑1 (MQ=255)
                                   tAGCGCTGCAGAATGGCCTTGCGCAGGCGCACCAGTTCATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCaacaac                                                                                                                   <  2:1545049/150‑1 (MQ=255)
                                    aGCGCTGCAGAATGGCATTGCAGAGGCGCACCAGTTCATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAacaaca                                                                                                                  <  1:592179/150‑1 (MQ=255)
                                     gcgcTGCAGAATGGCCTTGCGCAGGCGCACCAGTTCATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAg                                                                                                                 <  1:651425/150‑1 (MQ=255)
                                       gcTGCAGAATGGCCTTGCGCAGGCGCACCAGTTCATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGt                                                                                                               <  2:1432727/150‑1 (MQ=255)
                                       gcTGCAGAATGGCCTTGCGCAGGCGCACCAGTTCATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGt                                                                                                               <  2:980641/150‑1 (MQ=255)
                                        cTGCAGAATGGCCTTGCGCAGGCGCACCAGTTCATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTc                                                                                                              >  2:1240835/1‑150 (MQ=255)
                                           cAGAATGGCCTTGCGCAGGCGCACCAGTTCATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCttt                                                                                                           <  1:1185264/150‑1 (MQ=255)
                                           cAGAATGGCCTTGCGCAGGCGCACCAGTTCATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCttt                                                                                                           <  1:309403/150‑1 (MQ=255)
                                           cAGAATGGCCTTGCGCAGGCGCACCAGTTCATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCttt                                                                                                           <  1:1751875/150‑1 (MQ=255)
                                             gAATGGCCTTGCGCAGGCGCACCAGTTCATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGAt                                                                                                                          <  2:666142/133‑1 (MQ=255)
                                             gAATGGCCTTGCGCAGGCGCACCAGTTCATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGAt                                                                                                                          >  1:666142/1‑133 (MQ=255)
                                                   ccTTGCGCAGGCGCACCAGTTCATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATtgccc                                                                                                   <  1:1280636/150‑1 (MQ=255)
                                                    cTTGCGCAGGCGCACCAGTTCATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCCGGCGCGGGtga                                                                                                                                           <  1:1632995/109‑3 (MQ=255)
                                                    cTTGCGCAGGCGCACCAGTTCATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCCGGCGCGGGtga                                                                                                                                           >  2:1632995/1‑107 (MQ=255)
                                                          cAGGCGCACCAGTTCATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCtt                                                                                            <  1:1700303/150‑1 (MQ=255)
                                                             gcgcACCAGTTCATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGc                                                                                         <  2:292695/150‑1 (MQ=255)
                                                               gcACCAGTTCATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTg                                                                                                                                                     >  1:41620/1‑88 (MQ=255)
                                                               gcACCAGTTCATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTg                                                                                                                                                     <  2:41620/88‑1 (MQ=255)
                                                               gcACCAGTTCATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCgg                                                                                       >  1:1127600/1‑150 (MQ=255)
                                                               gcACCAGTTCATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCgg                                                                                       >  1:1513132/1‑150 (MQ=255)
                                                                cACCAGTTCATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTc                                                                                                          <  2:311967/130‑1 (MQ=255)
                                                                cACCAGTTCATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTc                                                                                                          >  1:311967/1‑130 (MQ=255)
                                                                   cAGTTCATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAAc                                                                                   >  1:1840453/1‑150 (MQ=255)
                                                                   cAGTTCATCTTTGTCGTCGCCCTTGGCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAAc                                                                                   <  2:722301/150‑1 (MQ=255)
                                                                      ttCATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTAGCGGGAACCcg                                                                                >  1:397947/1‑150 (MQ=255)
                                                                         aTCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAACCCGCGa                                                                             >  2:1590851/1‑150 (MQ=255)
                                                                         aTCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAACCCGCGa                                                                             <  1:1240835/150‑1 (MQ=255)
                                                                             ttGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAACCCGCGACCAt                                                                         >  1:1128812/1‑150 (MQ=255)
                                                                               gtcgtcGCCCTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAACCCGCGACCATTg                                                                       <  2:1128812/150‑1 (MQ=255)
                                                                                tcgtcgCCCTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAACCCGCGACCATTGt                                                                      >  2:115517/1‑150 (MQ=255)
                                                                                      cccTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAACCCGCGACCATTGTCGGGCg                                                                >  2:946730/1‑150 (MQ=255)
                                                                                       ccTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAACCCGCGACAATGGTCgcgggt                                                               <  1:32625/150‑6 (MQ=255)
                                                                                         ttGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCCTTGCGGGAACCCGCGACCATTGTCGGGCGGtt                                                             >  2:246556/1‑150 (MQ=255)
                                                                                              gAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGAt                                                                                                                          >  1:1483979/1‑84 (MQ=255)
                                                                                              gAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGAt                                                                                                                          <  2:1483979/84‑1 (MQ=255)
                                                                                                   cGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAACCCGCGACCATTGTCGGGCGGTTACCGGTAttt                                                   <  1:111919/150‑1 (MQ=255)
                                                                                                    gaagaaCACGCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCCGGCGCGGGTGACCCGCGccag                                                                                                                                 <  2:1519313/71‑3 (MQ=38)
                                                                                                    gaagaaCACGCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCCGGCGCGGGTGACCCGCGccag                                                                                                                                 >  1:1519313/1‑69 (MQ=39)
                                                                                                         acacGCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAACCCGCGACCATTGTCGGGCGGTTACCGGTATTTGGCaaa                                             <  1:1806730/150‑1 (MQ=255)
                                                                                                                 gCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAACCCGCGACCATTGTCGGGCGGTTACCGGTATTTGGCAAAGCAgctgc                                     >  1:1613679/1‑150 (MQ=255)
                                                                                                                   aCTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAACCCGCGACCATTGTCGGGCGGTTACCGGTATTTGGCAAAGCAGCTGCtt                                   >  2:1092584/1‑150 (MQ=255)
                                                                                                                   aCTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAACCCGCGACCATTGTCGGGCGGTTACCGGTATTTGGCAAAGCAGCTGCtt                                   <  2:453325/150‑1 (MQ=255)
                                                                                                                    cTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAACCCGCGACCATTGTCGGGCGGTTACCGGTATTTGGCAAAGCAGCTGCTTa                                  <  2:1513132/150‑1 (MQ=255)
                                                                                                                         cagcgcagcgACCGCTTCATATGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACCGGTCTTTCATTGCGCTGCCCTTTGCGGGAACCCGCGACCATTGTCGGGCGGTTACCGGTATTTGGCAAAGCAGCTGCTTagaaga                             <  2:896959/150‑1 (MQ=255)
                                                                                                                            cgcaTCGACCGCTTCATAGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAACCCGCGACCATTGTCGGGCGGTTACCGGTATTTGGCAAAGCAGCTGCTTAGAAGATCg                          <  2:666910/150‑1 (MQ=255)
                                                                                                                             gcaTCGACCGCTTCATAGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAACCCGCGACCATTGTCGGGCGGTTACCg                                                         >  1:1209280/1‑118 (MQ=255)
                                                                                                                             gcaTCGACCGCTTCATAGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAACCCGCGACCATTGTCGGGCGGTTACCg                                                         <  2:1209280/118‑1 (MQ=255)
                                                                                                                              caTCGACCGCTTCATAGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAACCCGCGACCATTGTCGGGCGGTTACCGGTATTTGGCAAAGCAGCTGCTTAGAAGATCgag                        <  2:1840453/150‑1 (MQ=255)
                                                                                                                               aTCGACCGCTTAATAGGTGACATGGCGCGGGTCACCCGCGACGATACCGATCAACAACATGTCTTTCATTGCCCTGCCCTTTGCGGGAACCCGCGACCATTGTCGGGCGGTTACCGGTATTTGGCAAAGCAGCTGCTTAGAAGATCgaga                       <  2:803368/150‑1 (MQ=255)
                                                                                                                                tCGACCGCTTCATAGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAACCCGCGACCATTGTCGGGCGGTTACCGGTATTTGGCAAAGCAGCTGCTTAGAAGATCGAGaa                      <  1:994323/150‑1 (MQ=255)
                                                                                                                                 cGACCGCTTCATAGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAACCCGCGACCATTGTCGGGCGGTTACCGGTATTTGGCAAAGCAGCTGCTTAGAAGATCGAGAAg                     <  1:1837252/150‑1 (MQ=255)
                                                                                                                                    ccGCTTCATAGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAACCCGCGACCATTGTCGGGCGGTTACCGGTATTTGGCAAAGCAGCTGCTTAGAAGATCGAGAAGTTg                  <  1:115517/150‑1 (MQ=255)
                                                                                                                                           aTAGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAACCCGCGACCATTGTCGGGCGGTTACCGGTATTTGGCAAAGCAGCTGCTTAGAAGATCGAGAAGTTGTAGCTGa           >  2:1004503/1‑150 (MQ=255)
                                                                                                                                              ggTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAACCCGCGACCATTGTCGGGCGGTTACCGGTATTTGGCAAAGCAGCTGCtt                                   <  1:1347076/123‑1 (MQ=255)
                                                                                                                                              ggTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAACCCGCGACCATTGTCGGGCGGTTACCGGTATTTGGCAAAGCAGCTGCtt                                   >  2:1347076/1‑123 (MQ=255)
                                                                                                                                               gTGACCTGGCGCGGGTCACCCGCGCCGATACCGCTCAACAGCAGGTCATTCATGGCCCTGCCCTTTGCGGGAAACTGCGTTCATACGCGGGCGGTTACCGGTATTTGGCAAAGCAGCTGCTTAGAAGATCGAGAAGTAGGTGCTGACAAt       >  1:925258/1‑150 (MQ=255)
                                                                                                                                               gTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAACCCGCGACCATTGTCGGGCGGTTACCGGTATTTGGCAAAGCAGCTGCTTAGAAGATCGAGAAGTTGTAGCTGACAAt       >  1:1080779/1‑150 (MQ=255)
                                                                                                                                                  aCCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCCGCCCTTTGCGGGAACCCGCGACCATTCTCGGGCGGTTACCGGTATTTGGCAAAGCATCTGCTTAGAAGATCGAGAAGTTGTAGCTGACTATGAc    <  2:113944/150‑1 (MQ=255)
                                                                                                                                                    cTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAACCCGCGACCATTGTCGGGCGGTTACCGGTATTTGGCAAAGCAGCTGCTTAGAAGATCGAGAAGTTGTAGCTGACAATGAccc  >  1:1745222/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CCACCGGCACCAGGCGATAACCACCCTCGGGGCGGTAGCGCTGCAGAATGGCCTTGCGCAGGCGCACCAGTTCATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGGTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAACCCGCGACCATTGTCGGGCGGTTACCGGTATTTGGCAAAGCAGCTGCTTAGAAGATCGAGAAGTTGTAGCTGACAATGACCC  >  NC_002947/4293852‑4294149

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: ‑