Predicted mutation
evidence seq id position mutation annotation gene description
RA CP006584 2,265,347 G→A intergenic (‑100/‑74) LY180_10855 ← / → LY180_10860 precorrin‑8W decarboxylase/terminase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0065842,265,3470GA96.3% 40.8 / NA 27intergenic (‑100/‑74)LY180_10855/LY180_10860precorrin‑8W decarboxylase/terminase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (14/12);  minor base C (0/1);  total (14/13)
Fisher's exact test for biased strand distribution p-value = 4.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGGTGTCACCCTCAACGCCGATACGAAAGAATTTTGAGACTTTTTTTGCCATTGTCAGGAGTCCTGAATAGTGATTAGAGGAGTCACATGTCGGCATCAGTTTCCCGACGATGCGCATCCTCCGCCATCAGTCCCGGATGGCTTATCACTGGCACAACAGCACCTTAGCGAATCGCGGGGCGCGACTCAGTAGCCTTGCCATGGATTCATCACGGCGAGGTATTCATGACCATCACCACAGACACCACTCTTTTACACGACCCGCGTCGTCAGGCGGCGCTGCTGTACTGGCAGGGATTTTCCGTGCCGCAGATTGCCGCCATGTTGCAGATGAAACGCCCGACGGTGCAGAGCTGGAAACAGCGCGACGGCTGGGACAGCGTTGCCCCCATCAGCCGTGTCGAAATGAGTCTGGAAGCGCGGCTGACCCAGCTCATCATCA  >  CP006584/2265196‑2265637
                                                                                                                                                       |                                                                                                                                                                                                                                                                                                  
aGGTGTCACCCTCAACGCCGATACGAAAGAATTTTGAGACTTTTTTTGCCATTGTCAGGAGTCCTGAATAGTGATTAGAGGAGTCACATGTCGGCATCAGTTTCCCGACGATGCGCATCCTCCGCCATCAGTCCCGGATGGCTTATCACTGACACAACAGCACCTTAGCGAATCGCGGGGCGCGACTCAGTAGCCTTGCCGTGTATTCATCACGGCGAGGTATTCATGACCATCACCACAGACACCACTCTTTTACACGACCCGCGTCGTCAGGCGGCGCTGCTGTACTGGCAGGGAt                                                                                                                                                  <  1:40795/298‑1 (MQ=11)
        cccTCAACGCCGATACGAAAGAATTTTGAGACTTTTTTTGCCATTGTCAGGAGTCCTGAATAGTGATTAGAGGAGTCACATGTCGGCATCAGTTTCCCGACGATGCGCATCCTCCGCCATCAGTCCCGGATGGCTTATCACTGACACAACAGCACCTTAGCGAATCGCGGGGCGCGACTCAGTAGCCTTGCCGTGTATTCATCACGGCGAGGTATTCATGACCATCACCACAGACACCACTCTTTTACACGACCCGCGTCGTCAGGCGGCGCTGCTGTACTGGCAGGGATTTTCCGTg                                                                                                                                          >  1:480955/1‑298 (MQ=14)
          cTCAACGCCAGTATAAATCAAGTTTGGCACTGTTTTTGCCATTGTCAGGAGTCGTGAATGGGGATTGGAGGAGTCACATGCCGGCAGCAGTTTACCGACGATGCGCATACTCCGCCATCGGTCCCGGAGGGCTTATCACTGACACAACAGCACCTTAGCGAATCGCGGGGCTCGACTCAGTAGCCTTGCCGTGTATTCATCACGGCGAGGTATTCATGACCATCACCACATACACCACTCTTTTACACGACCCGCGTCGTCAGGCGGCGCTGCTGTACTGGCAGGGATTTTCCGTGcc                                                                                                                                        <  2:166949/298‑1 (MQ=11)
           tCAACGCCGATACGAAAGAATTTTGAGACTTTTTTTGCCATTGTCAGGAGTCCTGAATAGTGATTAGAGGAGTCACATGTCGGCATCAGTTTCCCGACGATGCGCATCCTCCGCCATCAGTCCCGGATGGCTTATCACTGACACAACAGCACCTTAGCGAATCGCGGGGCGCGACTCAGTAGCCTTGCCGTGTATTCATCACGGCGAGGTATTCATGACCATCACCACAGACACCACTCTTTTACACGACCCGCGTCGTCAGGCGGCGCTGCTGTACTGGCAGGGATTTTCCGTGCCg                                                                                                                                       >  1:574104/1‑298 (MQ=14)
                  cGATACGAAAGAATTTTGAGACTTTTTTTGCCATTGTCAGGAGTCCTGAATAGTGATTAGAGGAGTCACATGTCGGCATCAGTTTCCCGACGATGCGCATCCTCCGCCATCAGTCCCGGATGGCTTATCACTGACACAACAGCACCTTAGCGAATCGCGGGGCGCGACTCAGTAGCCTTGCCGTGTATTCATCACGGCGAGGTATTCATGACCATCACCACAGACCCCACTCTTTTACACGACCCCCGTCGTCAGGCGGCCCCGCCGTACTGGCCGGGATTTTCCGTGCCGCCGATTg                                                                                                                                >  2:489694/1‑298 (MQ=11)
                   gATACGAAAGAATTTTGAGACTTTTTTTGCCATTGTCAGGAGTCCTGAATAGTGATTAGAGGAGTCACATGTCGGCATCAGTTTCCCGACGATGCGCATCCTCCGCCATCAGTCCCGGATGGCTTATCACTGACACAACAGCACCTTAGCGAATCGCGGGGCGCGACTCAGTAGCCTTGCCGTGTATTCATCACGGCGAGGTATTCATGACCATCACCACAGACACCACTCTTTTACACGACCCGCGTCGTCAGGCGGCGCTGCTGTACTGGCAGGGAGTTTCCGTGCCGCAGATTgc                                                                                                                               >  1:246032/1‑298 (MQ=14)
                                  tGAGACTTTTTTTGCCATTGTCAGGAGTCCTGAATAGTGATTAGAGGAGTCACATGTCGGCATCAGTTTCCCGACGATGCGCATCCTCCGCCATCAGTCCCGGATGGCTTATCACTGACACAACAGCACCTTAGCGAATCGCGGGGCGCGACTCAGTAGCCTTGCCGTGTATTCATCACGGCGAGGTATTCATGACCATCACCACAGACACCACTCTTTTACACGACCCGCGTCGTCAGGCGGCGCTGCTGTACTGGCAGGGATTTTCCGTGCCGCAGATTGCCGCCATGTTGCAGAt                                                                                                                <  2:160484/298‑1 (MQ=14)
                                                                  aaTAGTGATTAGAGGAGTCACATGTCGGCATCAGTTTCCCGACGATGCGCATCCTCCGCCATCAGTCCCGGATGGCTTATCACTGACACAACAGCACCTTAGCGAATCGCGGGGCGCGACTCAGTAGCCTTGCCGTGTATTCATCACGGCGAGGTATTCATGACCATCACCACAGACACCACTCTTTTACACGACCCGCGTCGTCAGGCGGCGCTGCTGTACTGGCAGGGATTTTCCGTGCCGCAGATTGCCGCCATGTTGCAGATGAAACGCCCGACGGTGCAGAGCTGGAAACAgc                                                                                <  1:415674/298‑1 (MQ=14)
                                                                           gaggagagTCACATGTCGGCATCAGTTTCCCGACTGTGCGCATCCCCCGCCTCCCGTCCCGGTTGGTTTTGCACTGAAACAACAAGACCTTAGCGAAGCGAGGGGCGCGACTCAGTAGCCTTGCCGTGTATTCATCACTGCGAGGTATTCATGTCCATCACCACAGTCCCCACTCTTTTGCACGACCCGCGTCGTCAGGCGGCGCTGCTGTACTGGCAGGGATTTTCCGTGCCGCAGATTGCCGCCATGTTGCAGATGAAACGCCCGACGGTGCAGAGCTGGAAACAGCGCGACGGCt                                                                       <  2:347342/293‑1 (MQ=11)
                                                                             gaggagTCACATGTCGGCATCAGTTTCCCGACGATGCGCATCCTCCGCCATCAGTCCCGGATGGCTTATCACTGACACAACAGCACCTTAGCGAATCGCGGGGCGCGACTCAGTAGCCTTGCCGTGTATTCATCACGGCGAGGTATTCATGACCATCACCACAGACACCACTCTTTTACACGACCCGCGTCGTCAGGCGGCGCTGCTGTACTGGCAGGGATTTTCCGTGCCGCAGATTGCCGCCATGTTGCAGATGAAACGCCCGACGGTGCAGAGCTGGAAACAGCGCGACGGCTgg                                                                     <  1:489694/298‑1 (MQ=14)
                                                                                   tCACATGTCGGCATCAGTTTCCCGACGATGCGCATCCTCCGCCATCAGTCCCGGATGGCTTATCACTGACACAACAGCACCTTAGCGAATCGCGGGGCGCGACTCAGTAGCCTTGCCGTGTATTCATCACGGCGAGGTATTCATGACCATCACCACAGACACCACTCTTTTACACGACCCGCGTCGTCAGGCGGCGCTGCTGTACTGGCAGGGATTTTCCGTGCCGCAGATTGCCGCCATGTTGCAGATGAAACGCCCGACGGTGCAGAGCTGGAAACAGCGCGACGGCTGGGACAGc                                                               <  1:577887/298‑1 (MQ=14)
                                                                                     ccagttCGGCATCGTTTTCCCGCCGTGTGGCTTCCACCGCCATCAGGCCCCGGTGCGCTAACACTGCTCCAACTGCACCTTCGCGAATGGGGGGCGGCGCCTCAGTGGACGTGCCGTGTGTTCAGCGCGGCGGGGTATTCATGACCATCAGCACAGACGCCACTCTTTTACGCGACCCGCGTCGTCAGGCGGCGCTGCTGTACTGGCAGGGATTTTCCGTGCCGCCGATTGCCGCCATGTTGCAGATGAAACGCCCGACGGTGCAGAGCTGGAAACAGCGCGACGGCTGGGACAGCGt                                                             <  2:170869/293‑1 (MQ=2)
                                                                                     acaTGTCGGCATCAGTTTCCCGACGATGCGCATCCTCCGCCATCAGTCCCGGATGGCTTATCACTGACACAACAGCACCTTAGCGAATCGCGGGGCGCGACTCAGTAGCCTTGCCGTGTATTCATCACGGCGAGGTATTCATGACCATCACCACAGACACCACTCTTTTACACGACCCGCGTCGTCAGGCGGCGCTGCTGTACTGGCAGGGATTTTCCGTGCCGCAGATTGCCGCCATGTTGCAGATGAAACGCCCGACGGTGCAGAGCTGGAAACAGCGCGACGGCTGGGACAGCGt                                                             >  1:505910/1‑298 (MQ=14)
                                                                                          tCGGCATCAGTTTCCCGACGATGCGCATCCTCCGCCATCAGTCCCGGATGGCTTATCACTGACACAACAGCACCTTAGCGAATCGCGGGGCGCGACTCAGTAGCCTTGCCGTGTATTCATCACGGCGAGGTATTCATGACCATCACCACAGACACCACTCTTTTACACGACCCGCGTCGTCAGGCGGCGCTGCTGTACTGGCAGGGATTTTCCGTGCCGCAGATTGCCGCCCCGTTGCAGATGAAACGCCCGACGGGGCACGGCAGGAAACAGCCCGCCGGCTGGGCCAGCGTCGccc                                                        >  2:356453/1‑298 (MQ=11)
                                                                                           ccccctCAGCTTCCCCCCGCTCCCCATCCCCCGCCCTCCGTCCCGCCTTGCTTATCCCTTACACAACCCCACCTTCGCGAAGCGCTGGGCGCGACTCAGTAGCCTTGCCGTGTATTCATCACGGCGAGGTATTCATGACCATCACCACAGACACCACTCTTTTACACGACCCGCGTCGTCAGGCGGCGCTGCTGTACTGGCAGGGATTTTCCGTGCCGCAGATTGCCGCCATGTTGCAGATGAAACGCCCGACGGTGCAGAGCTGGAAACAGCGCGACGGCTGGGACAGCGTTGcccc                                                       <  2:73976/293‑1 (MQ=11)
                                                                                               aTCAGTTTCCCAACGATGCGCATCCACCGCCATCAGTCCCGGATGGCTTCGCCCTGACACAACAGTACCTTAGCGAATCGCGGGGCGCGCCTCAGTAGCCTTGCCGTGTATTCAGCACGGCGAGGTATTCATGACCATCACCACAGACACCACTCTTTTACACGACCCGCGTCGTCAGGCGGCGCTGCTGTACTGGCAGGGATTTTCCGTGCCGCAGATTGCCGCCATGTTGCAGATGAAACGCCCGACGGTGCAGAGCTGGAAACAGCGCGACGGCTGGGACAGCGTTGCCCCCATc                                                   <  2:352353/298‑1 (MQ=11)
                                                                                                tCAGTTTCCCGACGATGCGCATCCTCCGCCATCAGTCCCGGATGGCTTATCACTGACACAACAGCACCTTAGCGAATCGCGGGGCGCGACTCAGTAGCCTTGCCGTGTATTCATCACGGCGAGGTATTCATGACCATCACCACAGACACCACTCTTTTACACGACCCGCGTCGTCAGGCGGCGCTGCTGTACTGGCAGGGATTTTCCGTGCCGCAGATTGCCGCCATGTTGCAGATGAAACGCCCGACGGTGCAGAGCTGGAAACAGCGCGACGGCTGGGACAGCGTTGCCCCCATCa                                                  >  1:309095/1‑298 (MQ=14)
                                                                                                   gTTTCCCGACGATGCGCATCCTCCGCCATCAGTCCCGGATGGCTTATCACTGACACAACAGCACCTTAGCGAATCGCGGGGCGCGACTCAGTAGCCTTGCCGTGTATTCATCACGGCGAGGTATTCATGACCATCACCACAGACACCACTCTTTTACACGACCCGCGTCGTCAGGCGGCGCTGCTGTACTGGCAGGGATTTTCCGTGCCGCAGATTGCCGCCATGTTGCAGCTGAAACGCCCGACGGCGCAGAGCTGGAAACAGCCCGCCGGGCGGGCCAGCGTCGgcccccgccgcc                                               >  2:43886/1‑291 (MQ=11)
                                                                                                    tttACCGCCGATGCGCATCCTCCGCCATCAGGCCCGGATGGCTTATCACTGACACAACAGCACCTTAGCGAATCGCGGGGCGCGACTCAGTAGCCTTGCCGTGTATTCATCACGGCGAGGTATTCATGACCATCACCACAGACACCACTCTTTTACACGACCCGCGTCGTCAGGCGGCGCTGCTGTACTGGCAGGGATTTTCCGTGCCGCAGATTGCCGCCATGTTGCAGATGAAACGCCCGACGGTGCAGAGCTGGAAACAGCGCGACGGCTGGGACAGCGTTGCCCCCATCAGCCg                                              <  1:73494/298‑1 (MQ=9)
                                                                                                          gacgaTGCGCATCCTCCGCCATCAGTCCCGGATGGCTTATCACTGACACAACAGCACCTTAGCGAATCGCGGGGCGCGACTCAGTAGCCTTGCCGTGTATTCATCACGGCGAGGTATTCATGACCATCACCACAGACACCACTCTTTTACACGACCCGCGTCGTCAGGCGGCGCTGCTGTACTGGCAGGGATTTTCCGTGCCGCAGATTGCCGCCATGTTGCAGATGAAACGCCCGACGGTGCATAGCTGGAAACAGCGCGACGGCTGGGACAGCGTTGCCCCCATCAGCCGTGTCGa                                        >  1:453439/1‑298 (MQ=14)
                                                                                                                   cATCCGCGGCCATCAGTCCCGGATGGCTTATCACTGACCCAACAGCACCTTAGCGAATCGCGGGGCGCGACTCAGTAGCCTTGCCGTGTATTCATCACGGCGAGGTATTCATGACCATCACCACAGACACCACTCTTTTACACGACCCGCGTCGTCAGGCGGCGCTGCTGTACTGGCAGGGATTTTCCGTGCCGCAGATTGCCGCCATGTTGCAGATGAAACGCCCGACGGTGCAGAGCTGGAAACAGCGCGACGGCTGGGACAGCGTTGCCCCCATCAGCCGTGTCGAAATGAGTCt                               <  2:480955/298‑1 (MQ=9)
                                                                                                                              tccAGCCCCGACTGGCTTAGCCCTGACCCAACAGCCCGTTAGGGAAGCGCGGGGCGCGACGCAGTAGCCTTGCCGTGTAGTCATCACGGCGAGGTTTTCATGACCATCACAACAGACCCCACACTTTTACACGACCCGCGTCGTCAGGCGGCGCTGCTGTACTGGCAGGGATTTTCCGTGCCGCAGATTGCCGCCATGTTGCAGATGAAACGCCCGACGGTGCAGAGCTGGAAACAGCGCGACGGCTGGGACAGCGTTGCCCCCATCAGCCGTGTCGAAATGAGTCTGGAAGCGCGGc                    <  2:201307/296‑1 (MQ=11)
                                                                                                                               tCAGTCCCGGATGGCTTATCACTGACACAACAGCACCTTAGCGAATCGCGGGGCGCGACTCAGTAGCCTTGCCGTGTATTCATCACGGCGAGGTATTCATGACCATCACCACAGACACCACTCTTTTACACGACCCGCGTCGTCAGGCGGCGCTGCTGTACTGGCAGGGATTTTCCGTGCCGCAGATTGCCGCCATGTTGCAGATGAAACGCCCGACGGTGCAGAGCTGGAAACAGCGCGACGGCTGGGACAGCGTTGCCCCCATCAGCCGTGTCGAAATGAGTCTGGAAGCGCGGCt                   >  1:492782/1‑298 (MQ=14)
                                                                                                                                 aGTCCCGGATGGCTTATCACTGACACAACAGCACCTTAGCGAATCGCGGGGCGCGACTCAGTAGCCTTGCCGTGTATTCATCACGGCGAGGTATTCATGACCATCACCACAGACACCACTCTTTTACACGACCCGCGTCGTCAGGCGGCGCTGCTGTACTGGCAGGGATTTTCCGTGCCGCAGATTGCCGCCATGTTGCAGATGAAACGCCCGACGGTGCAGAGCTGGAAACAGCGCGACGGCTGGGACAGCGTTGCCCCCATCAGCCGTGTCGAAATGAGTCTGGAAGCGCGGCtag                 >  1:201307/1‑296 (MQ=14)
                                                                                                                                              ttATCACTGACACAACAGCACCTTAGCGAATCGCGGGGCGCGACTCAGTAGCCTTGCCGTGTATTCATCACGGCGAGGTATTCATGACCATCACCACAGACACCACTCTTTTACACGACCCGCGTCGTCAGGCGGCGCTGCTGTACTGGCAGGGATTTTCCGTGCCGCAGATTGCCGCCATGTTGCAGATGAAACGCCCGACGGTGCAGAGCTGGAAACAGCGCGACGGCTGGGACAGCGTTGCCCCCATCAGCCGTGTCGAAATGAGTCTGGAAGCGCGGCTGACCCAGCccatcat    >  1:472627/1‑298 (MQ=14)
                                                                                                                                               tATCACTGACACAACAGCACCTTAGCGAATCGCGGGGCGCGACTCAGTAGCCTTGCCGTGTATTCATCACGGCGAGGTATTCATGACCATCACCACAGACACCACTCTTTTACACGACCCGCGTCGTCAGGCGGCGCTGCTGTACTGGCAGGGATTTTCCGTGCCGCAGATTGCCGCCATGTTGCAGATGAAACGCCCGACGGTGCAGAGCTGGAAACAGCGCGACGGCTGGGACAGCGTTGCCCCCCTCAGCCGTCGCGAAATGAGCTTGGAAAGCCGGCCGACCCAGccgctcctc   >  2:279495/1‑290 (MQ=11)
                                                                                                                                                aTCACTGACACAACAGCACCTTAGCGAATCGCGGGGCGCGACTCAGTAGCCTTGCCGTGTATTCATCACGGCGAGGTATTCATGACCATCACCACAGACACCACTCTTTTACACGACCCGCGTCGTCAGGCGGCGCTGCTGTACTGGCAGGGATTTTCCGTGCCGCAGATTGCCGCCATGTTGCAGATGAAACGCCCGACGGTGCAGAGCTGGAAACAGCGCGACGGCTGGGACAGCGTTGCCCCCATCAGCCGTGTCGAAATGAGTCTGGAAGCGCGGCTGACCCAGCtcatcatca  >  1:395058/1‑298 (MQ=14)
                                                                                                                                                       |                                                                                                                                                                                                                                                                                                  
AGGTGTCACCCTCAACGCCGATACGAAAGAATTTTGAGACTTTTTTTGCCATTGTCAGGAGTCCTGAATAGTGATTAGAGGAGTCACATGTCGGCATCAGTTTCCCGACGATGCGCATCCTCCGCCATCAGTCCCGGATGGCTTATCACTGGCACAACAGCACCTTAGCGAATCGCGGGGCGCGACTCAGTAGCCTTGCCATGGATTCATCACGGCGAGGTATTCATGACCATCACCACAGACACCACTCTTTTACACGACCCGCGTCGTCAGGCGGCGCTGCTGTACTGGCAGGGATTTTCCGTGCCGCAGATTGCCGCCATGTTGCAGATGAAACGCCCGACGGTGCAGAGCTGGAAACAGCGCGACGGCTGGGACAGCGTTGCCCCCATCAGCCGTGTCGAAATGAGTCTGGAAGCGCGGCTGACCCAGCTCATCATCA  >  CP006584/2265196‑2265637

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 17 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: