Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 640,333 A→C W43G (TGG→GGG)  rna ← ribonuclease I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273640,3330AC95.0% 207.8 / 2.7 60W43G (TGG→GGG) rnaribonuclease I
Reads supporting (aligned to +/- strand):  ref base A (1/2);  new base C (27/30);  total (28/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.95e-01

CCGATTTAGGCAATCCTGGCCACAGACCATGTACGGTCAGAAAATCAGCTTTGTTGGTCGTTTCGGTTTGCAGGCGACATTCATCTCGTTCGTTACGATTTCGATCGTGTTGACTCTGGCAAAATCCGGTTTGCCAGGAGAGGGCCAGGACATAGCGATCAAAATCGCCATACTGTTTTGCCTGCAACGCTAAGGCGTTGGCAGAAGAGAAGGGAAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGG  >  NZ_CP009273/640198‑640464
                                                                                                                                       |                                                                                                                                   
ccGATTTAGGCAATCCTGGCCACAGACCATGTACGGTCAGAAAATCAGCTTTGTTGGTCGTTTCGGTTTGCAGGCGACATTCATCTCGTTCGTTACGATTTCGATCGTGTTGACTCTGGCAAAATCCGGTTTGCCCGga                                                                                                                                  <  2:832311/139‑1 (MQ=255)
 cGATTTAGGCAATCCTGGCCACAGACCATGTACGGTCAGAAAATCAGCTTTGTTGGTCGTTTCGGTTTGCAGGCGACATTCATCTCGTTCGTTACGATTTCGATCGTGTTGACTCTGGCAAAATCCGGTTTGCCCGgag                                                                                                                                 <  1:822185/139‑1 (MQ=255)
 cGATTTAGGCAATCCTGGCCACAGACCATGTACGGTCAGAAAATCAGCTTTGTTGGTCGTTTCGGTTTGCAGGCGACATTCATCTCGTTCGTTACGATTTCGATCGTGTTGACTCTGGCAAAATCCGGTTTGCCCGgag                                                                                                                                 <  1:462191/139‑1 (MQ=255)
 cGATTTAGGCAATCCTGGCCACAGACCATGTACGGTCAGAAAATCAGCTTTGTTGGTCGTTTCGGTTTGCAGGCGACATTCATCTCGTTCGTTACGATTTCGATCGTGTTGACTCTGGCAAAATCCGGTTTGCCCGgag                                                                                                                                 <  2:60489/139‑1 (MQ=255)
 cGATTTAGGCAATCCTGGCCACAGACCATGTACGGTCAGAAAATCAGCTTTGTTGGTCGTTTCGGTTTGCAGGCGACATTCATCTCGTTCGTTACGATTTCGATCGTGTTGACTCTGGCAAAATCCGGTTTGCCCGgag                                                                                                                                 <  2:1105476/139‑1 (MQ=255)
 cGATTTAGGCAATCCTGGCCACAGACCATGTACGGTCAGAAAATCAGCTTTGTTGGTCGTTTCGGTTTGCAGGCGACATTCATCTCGTTCGTTACGATTTCGATCGTGTTGACTCTGGCAAAATCCGGTTTGCCAGgag                                                                                                                                 <  1:86611/139‑1 (MQ=255)
       aGGCAATCCTGGCCACAGACCATGTACGGTCAGAAAATCAGCTTTGTTGGTCGTTTCGGTTTGCAGGCGACATTCATCTCGTTCGTTACGATTTCGATCGTGTTGACTCTGGCAAAATCCGGTTTGCCCGGAGAGGGc                                                                                                                            <  2:86065/138‑1 (MQ=255)
             tCCTGGCCACAGACCATGTACGGTCAGAAAATCAGCTTTGTTGGTCGTTTCGGTTTGCAGGCGACATTCATCTCGTTCGTTACGATTTCGATCGTGTTGACTCTGGCAAAATCCGGTTTGCCCGGAGAGGGCCAGGACa                                                                                                                     >  1:1192022/1‑139 (MQ=255)
                  gCCACAGACCATGTACGGTCAGAAAATCAGCTTTGTTGGTCGTTTCGGTTTGCAGGCGACATTCATCTCGTTCGTTACGATTTCGATCGTGTTGACTCTGGCAAAATCCGGTTTGCCCGGAGAGGGCCAGGACagcgat                                                                                                                >  1:1257083/1‑134 (MQ=255)
                  gCCACAGACCATGTACGGTCAGAAAATCAGCTTTGTTGGTCGTTTCGGTTTGCAGGCGACATTCATCTCGTTCGTTACGATTTCGATCGTGTTGACTCTGGCAAAATCCGGTTTGCCCGGAGAGGGCCAGGACATAGc                                                                                                                 >  1:866494/1‑138 (MQ=255)
                  gCCACAGACCATGTACGGTCAGAAAATCAGCTTTGTTGGTCGTTTCGGTTTGCAGGCGACATTCATCTCGTTCGTTACGATTTAGATCGTGTTGACTCTGGCAAAATCCGGTTTGCCCGGAGAGGGCCAGGACATAGc                                                                                                                 >  2:1214714/1‑138 (MQ=255)
                       aGACCATGTACGGTCAGAAAATCAGCTTTGTTGGTCGTTTCGGTTTGCAGGCGACATTCATCTCGTTCGTTACGATTTCGATCGTGTTGACTCTGGCAAAATCCGGTTTGCCCGgag                                                                                                                                 >  2:1145710/1‑117 (MQ=255)
                       aGACCATGTACGGTCAGAAAATCAGCTTTGTTGGTCGTTTCGGTTTGCAGGCGACATTCATCTCGTTCGTTACGATTTCGATCGTGTTGACTCTGGCAAAATCCGGTTTGCCCGgag                                                                                                                                 <  1:1145710/117‑1 (MQ=255)
                       aGACCATGTACGGTCAGAAAATCAGCTTTGTTGGTCGTTTCGGTTTGCAGGCGACATTCATCTCGTTCGTTACGATTTCGATCGTGTTGACTCTGGCAAAATCCGGTTTGCCCGGAGAGGGCCAGGACATAGCg                                                                                                                <  1:98778/134‑1 (MQ=255)
                       aGACCATGTACGGTCAGAAAATCAGCTTTGTTGGTCGTTTCGGTTTGCAGGCGACATTCATCTCGTTCGTTACGATTTCGATCGTGTTGACTCTGGCAAAATCCGGTTTGCCCGGAGAGGGCCAGGACATAGCg                                                                                                                >  2:98778/1‑134 (MQ=255)
                       aGACCATGTACGGTCAGAAAATCAGCTTTGTTGGTCGTTTCGGTTTGCAGGCGACATTCATCTCGTTCGTTACGATTTCGATCGTGTTGACTCTGGCAAAATCCGGTTTGCCCGGAGAGGGCCAGGACATAGCGATCaa                                                                                                           <  2:804368/139‑1 (MQ=255)
                         aCCATGTACGGTCAGAAAATCAGCTTTGTTGGTCGTTTCGGTTTGCAGGCGACATTCATCTCGTTCGTTACGATTTCGATCGTGTTGACTCTGGCAAAATCCGGTTTGCCCGGAGAGGGCCAGGACATAGCGATCaaaa                                                                                                         <  1:715491/139‑1 (MQ=255)
                            aTGTACGGTCAGAAAATCAGCTTTGTTGGTCGTTTCGGTTTGCAGGCGACATTCATCTCGTTCGTTACGATTTCGATCGTGTTGACTCTGGCAAAATCCGGTTTGCCCGGAGAGGGCCAGGACATAGCGATCAAAATCg                                                                                                      >  1:373855/1‑139 (MQ=255)
                            aTGTACGGTCAGAAAATCAGCTTTGTTGGTCGTTTCGGTTTGCAGGCGACATTCATCTCGTTCGTTACGATTTCGATCGTGTTGACTCTGGCAAAATCCGGTTTGCCCGGAGAGGGCCAGGACATAGCGATCAAAATCg                                                                                                      >  1:605915/1‑139 (MQ=255)
                             tGTACGGTCAGAAAATCAGCTTTGTTGGTCGTTTCGGTTTGCAGGCGACATTCATCTCGTTCGTTACGATTTCGATCGTGTTGACTCTGGCAAAATCCGGTTTGCCCGGAGAGGGCCAGGACATAGCGATCAAAATCGc                                                                                                     >  1:908989/1‑139 (MQ=255)
                                    tCAGAAAATCAGCTTTGTTGGTCGTTTCGGTTTGCAGGCGACATTCATCTCGTTCGTTACGATTTCGATCGTGTTGACTCTGGCAAAATCCGGTTTGCCCGGAGAGGGCCAGGACATAGCGATCAAAATCGCCATACTg                                                                                              >  2:179025/1‑139 (MQ=255)
                                    tCAGAAAATCAGCTTTGTTGGTCGTTTCGGTTTGCAGGCGACATTCATCTCGTTCGTTACGATTTCGATCGTGTTGACTCTGGCAAAATCCGGTTTGCCCGGAGAGGGCCAGGACATAGCGATCAAAATCGCCATACTg                                                                                              >  1:501027/1‑139 (MQ=255)
                                    tCAGAAAATCAGCTTTGTTGGTCGTTTCGGTTTGCAGGCGACATTCATCTCGTTCGTTACGATTTCGATCGTGTTGACTCTGGCAAAATCCGGTTTGCCCGGAGAGGGCCAGGACATAGCGATCAAAATCGCCATACTg                                                                                              >  2:820111/1‑139 (MQ=255)
                                    tCAGAAAATCAGCTTTGTTGGTCGTTTCGGTTTGCAGGCGACATTCATCTCGTTCGTTACGATTTCGATCGTGTTGACTCTGGCAAAATCCGGTTTGCCAGGAGAGGGCCAGGACATAGCGATCAAAATCGCCATACTg                                                                                              >  1:1073591/1‑139 (MQ=255)
                                       gAAAATCAGCTTTGTTGGTCGTTTCGGTTTGCAGGCGACATTCATCTCGTTCGTTACGATTTCGATCGTGTTGACTCTGGCAAAATCCGGTTTGCCCGGAGAGGGCCAGGACATAGCg                                                                                                                <  1:118587/118‑1 (MQ=255)
                                       gAAAATCAGCTTTGTTGGTCGTTTCGGTTTGCAGGCGACATTCATCTCGTTCGTTACGATTTCGATCGTGTTGACTCTGGCAAAATCCGGTTTGCCCGGAGAGGGCCAGGACATAGCg                                                                                                                >  2:118587/1‑118 (MQ=255)
                                       gAAAATCAGCTTTGTTGGTCGTTTCGGTTTGCAGGCGACATTCATCTCGTTCGTTACGATTTCGATCGTGTTGACTCTGGAAAAATCCGGTTTGCCCGGAGAGGGCCAGGACATAGCGATCAAAATCGCCATACTGttt                                                                                           >  1:1102851/1‑139 (MQ=255)
                                        aaaaTCAGCTTTGTTGGTCGTTTCGGTTTGCAGGCGACATTCATCTCGTTCGTTACGATTTCGATCGTGTTGACTCTGGCAAAATCCGGTTTGCCCGGAGAGGGCCAGGACATAGCGATCAAAATCGCCATACTGtttt                                                                                          >  2:961367/1‑139 (MQ=255)
                                        aaaaTCAGCTTTGTTGGTCGTTTCGGTTTGCAGGCGACATTCATCTCGTTCGTTACGATTTCGATCGTGTTGACTCTGGCAAAATCCGGTTTGCCCGGAGAGGGCCAGGACATAGCGATCAAAATCGCCATACTGtttt                                                                                          >  1:975808/1‑139 (MQ=255)
                                        aaaaTCAGCTTTGTTGGTCGTGTCGGTTTGCAGGCGACATTCATCTCTTTCGTTACGATTTCGATCGTGTTGACTCTGGCAAAATCCGGTTTGCCCGGAGAGGGCCAGGACATAGCGATCAAAATCGCCATACTGtttt                                                                                          <  2:222433/139‑1 (MQ=255)
                                              aGCTTTGTTGGTCGTTTCGGTTTGCAGGCGACATTCATCTCGTTCGTTACGATTTCGATCGTGTTGACTCTGGCAAAATCCGGTTTGCCAGGAGAGGGCCAGGACATAGCGATCAAAATCGCCATACTGTTTTGCCTg                                                                                     <  1:654993/138‑1 (MQ=255)
                                                     ttgGTCGTTTCGGTTTGCAGGCGACATTCATCTCGTTCGTTACGATTTCGATCGTGTTGACTCTGGCAAAATCCGGTTTGCCCGgag                                                                                                                                 <  2:528621/87‑1 (MQ=255)
                                                     ttgGTCGTTTCGGTTTGCAGGCGACATTCATCTCGTTCGTTACGATTTCGATCGTGTTGACTCTGGCAAAATCCGGTTTGCCCGgag                                                                                                                                 >  1:528621/1‑87 (MQ=255)
                                                                    tGCAGGCGACATTCATCTCGTTCGTTACGATTTCGATCGTGTTGACTCTGGCAAAATCCGGTTTGCCCGGAGAGGGCCAGGACATAGCGATCAAAATCGCCATACTGTTTTGCCTGCAACGCTAAGGCGTTGGCag                                                                 >  2:859171/1‑136 (MQ=255)
                                                                     gCAGGCGACATTCATCTCGTTCGTTACGATTTCGATCGTGTTGACTCTGGCAAAATCCGGTTTGCCCGGAGAGGGCCAGGACATAGCGATCAAAATCGCCATACTGTTTTGCCTGCAACGCTAAGGCGTTGGCagaag                                                              <  1:941768/138‑1 (MQ=255)
                                                                                       gttcgttACGATTTCGATCGTGTTGACTCTGGCAAAATCCGGTTTGCCCGGAGAGGGCCAGGACATAGCGATCAAAATCGCCATACTGTTTTGCCTGCAACGCTAAGGCGt                                                                       >  1:881235/1‑111 (MQ=255)
                                                                                       gttcgttACGATTTCGATCGTGTTGACTCTGGCAAAATCCGGTTTGCCCGGAGAGGGCCAGGACATAGCGATCAAAATCGCCATACTGTTTTGCCTGCAACGCTAAGGCGt                                                                       <  2:881235/111‑1 (MQ=255)
                                                                                        ttcgttACGATTTCGATCGTGTTGACTCTGGCAAAATCCGGTTTGCCCGgag                                                                                                                                 >  1:721469/1‑52 (MQ=255)
                                                                                        ttcgttACGATTTCGATCGTGTTGACTCTGGCAAAATCCGGTTTGCCCGgag                                                                                                                                 <  2:721469/52‑1 (MQ=255)
                                                                                        ttcgttACGATTTCGATCGTGTTGACTCTGGCAAAATCCGGTTTGCCCGGAGAGGGCCAGGACATAGCGATCAAAATCGCCATACTGTTTTGCCTGCAACGCTAAGGCGTTGGCAGAAGAGAAGGGAAGCAGAGAAAcc                                          >  2:452827/1‑139 (MQ=255)
                                                                                          cgttACGATTTCGATCGTGTTGACTCTGGCAAAATCCGGTTTGCCCGGAGAGGGCCAGGACATAGCGATCAAAATCGCCATACTGTTTTGCCTGCAACGCTAAGGCGTTGGCAGAAGAGAAGGGAAGCAGAGAAACcgc                                        >  2:380928/1‑139 (MQ=255)
                                                                                            ttACGATTTCGATCGTGTTGACTCTGGCAAAATCCGGTTTGCCCGGAGAGGGc                                                                                                                            <  2:289421/53‑1 (MQ=255)
                                                                                            ttACGATTTCGATCGTGTTGACTCTGGCAAAATCCGGTTTGCCCGGAGAGGGc                                                                                                                            >  1:289421/1‑53 (MQ=255)
                                                                                              aCGATTTCGATCGTGTTGACTCTGGCAAAATCCGGTTTGCCCGGAGAGGGCCAGGACATAGCGATCAAAATCGCCATACTGTTTTGCCTGCAACGCTAAGGCGTTGGCAGAAGAGAAGGGAAGCAGAGAAACCGCGAg                                     <  1:452827/138‑1 (MQ=255)
                                                                                                       aTCGTGTTGACTCTGGCAAAATCCGGTTTGCCCGGAGAGGGCCAGGACATAGCGATCAAAATCGCCATACTGTTTTGCCTGCAACGCTAAGGCGTTGGCAGAAGAGAAGGGAAGCagag                                               >  2:684290/1‑119 (MQ=255)
                                                                                                       aTCGTGTTGACTCTGGCAAAATCCGGTTTGCCCGGAGAGGGCCAGGACATAGCGATCAAAATCGCCATACTGTTTTGCCTGCAACGCTAAGGCGTTGGCAGAAGAGAAGGGAAGCagag                                               <  1:684290/119‑1 (MQ=255)
                                                                                                       aTCGTGTTGACTCTGGCAAAATCCGGTTTGCCCGGAGAGGGCCAGGACATAGCGATCAAAATCGCCATACTGTTTTGCCTGCAACGCTAAGGCGTTGGCAGAAGAGAAGGGAAGCAGAGAAACCGCGAGCAACgcggcg                           <  2:605915/139‑1 (MQ=255)
                                                                                                        tCGTGTTGACTCTGGCAAAATCCGGTTTGCCCGGAGAGGGCCAGGACATAGCGATCAAAATCGCCATACTGTTTTGCCTGCAACGCTAAGGCGTTGGCAGAAGAGAAGGGAAGCAGAGAAACCGCGAGCAACGCGGCGt                          <  2:724038/139‑1 (MQ=255)
                                                                                                                 ctctGGCAAAATCCGGTTTGCCCGGAGAGGGCCAGGACATAGCGATCAAAATCGCCATACTGTTTTGCCTGCAACGCTAAGGCGTTGGCAGAAGAGAAGGGAAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGa                 <  1:247308/139‑1 (MQ=255)
                                                                                                                 ctctGGCAAAATCCGGTTTGCCCGGAGAGGGCCAGGACATAGCGATCAAAATCGCCATACTGTTTTGCCTGCAACGCTAAGGCGTTGGCAGAAGAGAAGGGAAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGa                 <  1:394829/139‑1 (MQ=255)
                                                                                                                    tGGCAAAATCCGGTTTGCCCGGAGAGGGCCAGGACATAGCGATCAAAATCGCCATACTGTTTTGCCTGCAACGCTAAGGCGTTGGCAGAAGAGAAGGGAAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGaa                <  2:1141281/137‑1 (MQ=255)
                                                                                                                     ggCAAAATCCGGTTTGCCCGGAGAGGGCCAGGACATAGCGATCAAAATCGCCATACTGTTTTGCCTGCAACGCTAAGGCGTTGGCAGAAGAGAAGGGAAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGc             >  1:106481/1‑139 (MQ=255)
                                                                                                                      gCAAAATCCGGTTTGCCCGGAGAGGGCCAGGACATAGCGATCAAAATCGCCATACTGTTTTGCCTGCAACGCt                                                                              >  2:653181/1‑73 (MQ=255)
                                                                                                                      gCAAAATCCGGTTTGCCCGGAGAGGGCCAGGACATAGCGATCAAAATCGCCATACTGTTTTGCCTGCAACGCt                                                                              <  1:653181/73‑1 (MQ=255)
                                                                                                                      gCAAAATCCGGTTTGCCCGGAGAGGGCCAGGACATAGCGATCAAAATCGCCATACTGTTTTGCCTGCAACGCTAAGGCGTTGGCAGAAGAGAAGGGAAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCt            <  1:421298/139‑1 (MQ=255)
                                                                                                                              cGGTTTGCCCGGAGAGGGCCAGGACATAGCGATCAAAATCGCCATACTGTTTTGCCTGCAACGCTAAGGCGTTGGCAGAAGAGAAGGGAAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAAt    <  1:294997/139‑1 (MQ=255)
                                                                                                                              cGGTTTGCCCGGAGAGGGCCAGGACATAGCGATCAAAATCGCCATACTGTTTTGCCTGCAACGCTAAGGCGTTGGCAGAAGAGAAGGGAAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAAt    <  1:951030/139‑1 (MQ=255)
                                                                                                                              cGGTTTGCCCGGAGAGGGCCAGGACATAGCGATCAAAATCGCCATACTGTTTTGCCTGCAACGCTAAGGCGTTGGCAGAAGAGAAGGGAAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAAt    <  1:859171/139‑1 (MQ=255)
                                                                                                                               ggTTTGCCCGGAGAGGGCCAGGACATAGCGATCAAAATCGCCATACTGTTTTGCCTGCAACGCTAAGGCGTTGGCAGAAGAGAAGGGAAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATg   <  1:179025/139‑1 (MQ=255)
                                                                                                                                gTTTGCCCGGAGAGGGCCAGGACATAGCGATCAAAATCGCCATACTGTTTTGCCTGCAACGCTAAGGCGTTGGCAGAAGAGAAGGGAAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATgg  >  2:274159/1‑139 (MQ=255)
                                                                                                                                       |                                                                                                                                   
CCGATTTAGGCAATCCTGGCCACAGACCATGTACGGTCAGAAAATCAGCTTTGTTGGTCGTTTCGGTTTGCAGGCGACATTCATCTCGTTCGTTACGATTTCGATCGTGTTGACTCTGGCAAAATCCGGTTTGCCAGGAGAGGGCCAGGACATAGCGATCAAAATCGCCATACTGTTTTGCCTGCAACGCTAAGGCGTTGGCAGAAGAGAAGGGAAGCAGAGAAACCGCGAGCAACGCGGCGTTACGCCAGAATGCTTTCATAATGG  >  NZ_CP009273/640198‑640464

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑