Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 915,468 C→T A194T (GCT→ACT)  aqpZ ← water channel AqpZ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913915,4680CT98.3% 195.9 / ‑4.5 58A194T (GCT→ACT) aqpZwater channel AqpZ
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (31/26);  total (31/27)
Fisher's exact test for biased strand distribution p-value = 4.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

ACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGCACCACCCAGAAGAACCACAGTTGTTCTAATGCCCAGCCGCCCTGGAAGATAGCAACCGCGGTGCTGCGCGCCGGGTTAACAGAAGTGTTAGTCACCGGAATACTAATTAAGTGAATCAGGGTTAAGGCCAGACCAATAGCGATCGGCGCAAAACCTGCCGGCGCGAATTTGTCGGTTGCGCCGTGGAT  >  NC_000913/915331‑915603
                                                                                                                                         |                                                                                                                                       
aCGCCAGGCCTTCTGCGTTAATTAATAACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGCACCACCCAGAAGAACCACAGTTGTTCTAATGCCCAGCCGCCCTGGAAGATAGTaa                                                                                                                                       >  1:685533/1‑140 (MQ=255)
 cGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGCACCACCCAGAAGAACCACAGTTGTTCTAATGCCCAGCCGCCCTGGAAGATAGTAAc                                                                                                                                      >  2:898401/1‑140 (MQ=255)
       gCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGCACCACCCAGAAGAACCACAGTTGTTCTAATGCCCAGCCGCCCTGGAAGATAGTAACCGCGGt                                                                                                                                >  1:639330/1‑140 (MQ=255)
       gCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGCACCACCCAGAAGAACCACAGTTGTTCTAATGCCCAGCCGCCCTGGAAGATAGTAACCGCGGt                                                                                                                                >  1:602155/1‑140 (MQ=255)
         ttgcTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGCACCACCCAGAAGAACCACAGTTGTTCTAATGCCCAGCCGCCCTGGAAGATAGTAACCGCGGtgc                                                                                                                              <  2:411881/137‑1 (MQ=255)
            cTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGCACCACCCAGAAGAACCACAGTTGTTCTAATGCCCAGCCGCCCTGGAAGATAGTAACCGCGGtgctgc                                                                                                                           >  2:62227/1‑140 (MQ=255)
             tGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGCACCACCCAGAAGAACCACAGTTGTTCTAATGCCCAGCCGCCCTGGAAGATAGTAACCGCGGTGCTgcg                                                                                                                          >  2:880864/1‑140 (MQ=255)
              gCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGCACCACCCAGAAGAACCACAGTTGTTCTAATGCCCAGCCGCCCTGGAAGATAGTAACCGCGGTGCTgcgc                                                                                                                         <  1:137351/140‑1 (MQ=255)
                gTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGCACCACCCAGAAGAACCACAGTTGTTCTAATGCCCAGCCGCCCTGGAAGATAGTAACCGCGGTGCTgcgcgc                                                                                                                       <  1:546000/140‑1 (MQ=255)
                         tCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGCACCACCCAGAAGAACCACAGTTGTTCTAATGCCCAGCCGCCCTGGAAGATAGTAACCGCGGTGCTGCGCGCCGGGTTAAc                                                                                                              >  2:837216/1‑140 (MQ=255)
                          cACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGCACCACCCAGAAGAACCACAGTTGTTCTAATGCCCAGCCGCCCTGGAAGATAGTAACCGCGGTGCTGCGCGCCGGGTTAACa                                                                                                             <  2:463806/140‑1 (MQ=255)
                           aCGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGCACCACCCAGAAGAACCACAGTTGTTCTAATGCCCAGCCGCCCTGGAAGATAGTAACCGCGGTGCTGCGCGCCGGGTTAACAg                                                                                                            <  1:898401/140‑1 (MQ=255)
                              cTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGCACCACCCAGAAGAACCACAGTTGTTCTAATGCCCAGCCGCCCTGGAAGATAGTAACCGCGGTGCTGCGCGCCGGGTTAACAGAAg                                                                                                         >  2:322767/1‑140 (MQ=255)
                               ttttCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGCACCACCCAGAAGAACCACAGTTGTTCTAATGCCCAGCCGCCCTGGAAGATAGTAACCGCGGTGCTGCGCGCCGGGTTAACAGAAgt                                                                                                        >  1:851452/1‑140 (MQ=255)
                                 ttCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGCACCACCCAGAAGAACCACAGTTGTTCTAATGCCCAGCCGCCCTGGAAGATAGTAACCGCGGTGCTGCGCGCCGGGTTAACAGAAgtgt                                                                                                      <  2:434877/140‑1 (MQ=255)
                                     agcagGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGCACCACCCAGAAGAACCACAGTTGTTCTAATGCCCAGCCGCCCTGGAAGATAGTAACCGCGGTGCTGCGCGCCGGGTTAACAGAAGTGTTAGt                                                                                                  <  1:880864/140‑1 (MQ=255)
                                               ggTAAATCAGACCACCGATAATGCCGCCGACAATTGGCACCACCCAGAAGAACCACAGTTGTTCTAATGCCCAGCCGCCCTGGAAGATAGTAACCGCGGTGCTGCGCGCCGGGTTAACAGAAGTGTTAGTCACCGGAATa                                                                                        <  1:62227/140‑1 (MQ=255)
                                                    aTCAGACCACCGATAATGCCGCCGACAATTGGCACCACCCAGAAGAACCACAGTTGTTCTAATGCCCAGCCGCCCTGGAAGATAGTAACCGCGGTGCTGCGCGCCGGGTTAACAGAAGTGTTAGTCACCg                                                                                             <  1:184157/130‑1 (MQ=255)
                                                    aTCAGACCACCGATAATGCCGCCGACAATTGGCACCACCCAGAAGAACCACAGTTGTTCTAATGCCCAGCCGCCCTGGAAGATAGTAACCGCGGTGCTGCGCGCCGGGTTAACAGAAGTGTTAGTCACCg                                                                                             >  2:184157/1‑130 (MQ=255)
                                                       aGACCACCGATAATGCCGCCGACAATTGGCACCACCCAGAAGAACCACAGTTGTTCTAATGCCCAGCCGCCCTGGAAGATAGTAACCGCGGTGCTGCGCGCCGGGTTAACAGAAGTGTTAGTCACCGGAATACTAATTaa                                                                                <  2:433455/140‑1 (MQ=255)
                                                       aGACCACCGATAATGCCGCCGACAATTGGCACCACCCAGAAGAACCACAGTTGTTCTAATGCCCAGCCGCCCTGGAAGATAGTAACCGCGGTGCTGCGCGCCGGGTTAACAGAAGTGTTAGTCACCGGAATACTAATTaa                                                                                <  2:685533/140‑1 (MQ=255)
                                                       aGACCACCGATAATGCCGCCGACAATTGGCACCACCCAGAAGAACCACAGTTGTTCTAATGCCCAGCCGCCCTGGAAGATAGTAACCGCGGTGCTGCGCGCCGGGTTAACAGAAGTGTTAGTCACCGGAATACTAATTaa                                                                                <  2:522023/140‑1 (MQ=255)
                                                       aGACCACCGATAATGCCGCCGACAATTGGCACCACCCAGAAGAACCACAGTTGTTCTAATGCCCAGCCGCCCTGGAAGATAGTAACCGCGGTGCTGCGCGCCGGGTTAACAGAAGTGTTAGTCACCGGAATACTAATTaa                                                                                <  2:16116/140‑1 (MQ=255)
                                                        gACCACCGATAATGCCGCCGACAATTGGCACCACCCAGAAGAACCACAGTTGTTCTAATGCCCAGCCGCCCTGGAAGATAGTAACCGCGGTGCTGCGCGCCGGGTTAACAGAAGTGTTAGTCACCGGAATACTAATTAAg                                                                               >  1:409422/1‑140 (MQ=255)
                                                             ccGATAATGCCGCCGACAATTGGCACCACCCAGAAGAACCACAGTTGTTCTAATGCCCAGCCGCCCTGGAAGATAGTAACCGCGGTGCTGCGCGCCGGGTTAACAGAAGTGTTAGTCACCGGAATACTAATTAAGTGAAt                                                                          <  2:284755/140‑1 (MQ=255)
                                                              cGATAATGCCGCCGACAATTGGCACCACCCAGAAGAACCACAGTTGTTCTAATGCCCAGCCGCCCTGGAAGATAGTAACCGCGGTGCTGCGCGCCGGGTTAACAGAAGTGTTAGTCACCGGAATACTAATTAAGTGAATc                                                                         <  2:163619/140‑1 (MQ=255)
                                                               gATAATGCCGCCGACAATTGGCACCACCCAGAAGAACCACAGTTGTTCTAATGCCCAGCCGCCCTGGAAGATAGTAACCGCGGTGCTGCGCGCCGGGTTAACAGAAGTGTTAGTCACCGGAATACTAATTAAGTGAATCa                                                                        <  2:426735/140‑1 (MQ=255)
                                                                aTAATGCCGCCGACAATTGGCACCACCCAGAAGAACCACAGTTGTTCTAATGCCCAGCCGCCCTGGAAGATAGTAACCGCGGTGCTGCGCGCCGGGTTAACAGAAGTGTTAGTCACCGGAATACTAATTAAGTGAATCAg                                                                       >  2:368555/1‑140 (MQ=255)
                                                                 tAATGCCGCCGACAATTGGCACCACCCAGAAGAACCACAGTTGTTCTAATGCCCAGCCGCCCTGGAAGATAGTAACCGCGGTGCTGCGCGCCGGGTTAACAGAAGTGTTAGTCACCGGAATACTAATTAAGTGAATCa                                                                        >  1:843324/1‑138 (MQ=255)
                                                                 tAATGCCGCCGACAATTGGCACCACCCAGAAGAACCACAGTTGTTCTAATGCCCAGCCGCCCTGGAAGATAGTAACCGCGGTGCTGCGCGCCGGGTTAACAGAAGTGTTAGTCACCGGAATACTAATTAAGTGAATCa                                                                        <  2:843324/138‑1 (MQ=255)
                                                                    tGCCGCCGACAATTGGCACCACCCAGAAGAACCACAGTTGTTCTAATGCCCAGCCGCCCTGGAAGATAGTAACCGCGGTGCTGCGCGCCGGGTTAACAGAAGTGTTAGTCACCGGAATACTAATTAAGTGAATCAGGGtt                                                                   <  2:639330/140‑1 (MQ=255)
                                                                    tGCCGCCGACAATTGGCACCACCCAGAAGAACCACAGTTGTTCTAATGCCCAGCCGCCCTGGAAGATAGTAACCGCGGTGCTGCGCGCCGGGTTAACAGAAGTGTTAGTCACCGGAATACTAATTAAGTGAATCAGGGtt                                                                   <  2:602155/140‑1 (MQ=255)
                                                                        gccgACAATTGGCACCACCCAGAAGAACCACAGTTGTTCTAATGCCCAGCCGCCCTGGAAGATAGTAACCGCGGTGCTGCGCGCCGGGTTAACAGAAGTGTTAGTCACCGGAATACTAATTAAGTGAATCAGGGTTAAgg                                                               >  2:63996/1‑140 (MQ=255)
                                                                        gccgACAATTGGCACCACCCAGAAGAACCACAGTTGTTCTAATGCCCAGCCGCCCTGGAAGATAGTAACCGCGGTGCTGCGCGCCGGGCTAACAGAAGTGTTAGTCACCGGAATACTAATTAAGTGAATCAGGGTTAAgg                                                               >  2:292324/1‑140 (MQ=255)
                                                                         ccgACAATTGGCACCACCCAGAAGAACCACAGTTGTTCTAATGCCCAGCCGCCCTGGAAGATAGCAACCGCGGTGCTGCGCGCCGGGTTAACAGAAGTGTTAGTCACCGGAATACTAATTAAGTGAATCAGGGTTAAGGc                                                              <  1:830789/140‑1 (MQ=255)
                                                                                        accCAGAAGAACCACAGTTGTTCTAATGCCCAGCCGCCCTGGAAGATAGTAACCGCGGTGCTGCGCGCCGGGTTAACAGAAGTGTTAGTCACCGGAATACTAATTAAGTGAATCAGGGTTAAGGCCAGACCAATAGCGAt                                               >  2:51544/1‑140 (MQ=255)
                                                                                         cccAGAAGAACCACAGTTGTTCTAATGCCCAGCCGCCCTGGAAGATAGTAACCGCGGTGCTGCGCGCCGGGTTAACAGAAGTGTTAGTCACCGGAATACTAATTAAGTGAATCAGGGTTAAGGCCAGACCAATAGCGATc                                              <  1:322767/140‑1 (MQ=255)
                                                                                         cccAGAAGAACCACAGTTGTTCTAATGCCCAGCCGCCCTGGAAGATAGTAACCGCGGTGCTGCGCGCCGGGTTAACAGAAGTGTTAGTCACCGGAATACTAATTAAGTGAATCAGGGTTAAGGCCAGACCAATAGCGATc                                              <  2:209898/140‑1 (MQ=255)
                                                                                          ccAGAAGAACCACAGTTGTTCTAATGCCCAGCCGCCCTGGAAGATAGTAACCGCGGTGCTGCGCGCAGGGTTAACAGAAGTGTTAGTCACCGGAATACTAATTAAGTGAATCAGGGTTAAGGCCAGACCAATAGCGATCg                                             >  1:178001/1‑140 (MQ=255)
                                                                                              aagaaCCACAGTTGTTCTAATGCCCAGCCGCCCTGGAAGATAGTAACCGCGGTGCTGCGCGCCGGGTTAACAGAAGTGTTAGTCACCGGAATACTAATTAAGTGAATCAGGGTTAAGGCCAGACCAATAGCGATCGgcgc                                         <  2:752541/140‑1 (MQ=255)
                                                                                               agaaCCACAGTTGTTCTAATGCCCAGCCGCCCTGGAAGATAGTAACCGCGGTGCTGCGCGCCGGGTTAACAGAAGTGTTAGTCACCGGAATACTAATTAAGTGAATCAGGGTTAAGGCCAGACCAATAGCGATCGGCGCa                                        <  1:63996/140‑1 (MQ=255)
                                                                                                   ccACAGTTGTTCTAATGCCCAGCCGCCCTGGAAGATAGTAACCGCGGTGCTGCGCGCCGGGTTAACAGAAGTGTTAGTCACCGGAATACTAATTAAGTGAATCAGGGTTAAGGCCAGACCAATAGCGATCGGCGCAAAAc                                    >  2:243097/1‑140 (MQ=255)
                                                                                                   ccACAGTTGTTCTAATGCCCAGCCGCCCTGGAAGATAGTAACCGCGATGCTGCGCGCCGGGTTAACAGAAGTGTTAGTCACCGGAATACTAATTAAGTGAATCAGGGTTAAGGCCAGACCAATAGCGATCGGCGCAAAAc                                    >  1:413881/1‑140 (MQ=255)
                                                                                                    cacaGTTGTTCTAATGCCCAGCCGCCCTGGAAGATAGTAACCGCGGTGCTGCGCGCCGGGTTAACAGAAGTGTTAGTCACCGGAATACTAATTAAGTGAATCAGGGTTAAGGCCAGACCAATAGCGATCGGCGCAAAAcc                                   >  1:8634/1‑140 (MQ=255)
                                                                                                    cacaGTTGTTCTAATGCCCAGCCGCCCTGGAAGATAGTAACCGCGGTGCTGCGCGCCGGGTTAACAGAAGTGTTAGTCACCGGAATACTAATTAAGTGAATCAGGGTTAAGGCCAGACCAATAGCGATCGGCGCAAAAcc                                   >  1:756244/1‑140 (MQ=255)
                                                                                                           gttCTAATGCCCAGCCGCCCTGGAAGATAGTAACCGCGGTGCTGCGCGCCGGGTTAACAGAAGTGTTAGTCACCGGAATACTAATTAAGTGAATCAGGGTTAAGGCCAGACCAATAGCGATCGGCGCAAAACCTGCCGgc                            >  1:744302/1‑140 (MQ=255)
                                                                                                                 aTGCCCAGCCGCCCTGGAAGATAGTAACCGCGGTGCTGCGCGCCGGGTTAACAGAAGTGTTAGTCACCGGAATACTAATTAAGTGAATCAGGGTTAAGGCCAGACCAATAGCGATCGGCGCAAAACCTGCCGGCGCGAAt                      >  2:471688/1‑140 (MQ=255)
                                                                                                                 aTGCCCAGCCGCCCTGGAAGATAGTAACCGCGGTGCTGCGCGCCGGGTTAACAGAAGTGTTAGTCACCGGAATACTAATTAAGTGAATCAGGGTTAAGGCCAGACCAATAGCGATCGGCGCAAAACCTGCCGGCGCGAAt                      >  2:633835/1‑140 (MQ=255)
                                                                                                                   gCCCAGCCGCCCTGGAAGATAGTAACCGCGGTGCTGCGCGCCGGGTTAACAGAAGTGTTAGTCACCGGAATACTAATTAAGTGAATCAGGGTTAAGGCCAGACCAATAGCGATCGGCGCAAAACCTGCCGGCGCGAAttt                    >  2:417736/1‑140 (MQ=255)
                                                                                                                      cAGCCGCCCTGGAAGATAGTAACCGCGGTGCTGCGCGCCGGGTTAACAGAAGTGTTAGTCACCGGAATACTAATTAAGTGAATCAGGGTTAAGGCCAGACCAATAGCGATCGGCGCAAAACCTGCCGGCGCGAATTTGTc                 >  2:721913/1‑140 (MQ=255)
                                                                                                                      cAGCCGCCCTGGAAGATAGTAACCGCGGTGCTGCGCGCCGGGTTAACAGAAGTGTTAGTCACCGGAATACTAATTAAGTGAATCAGGGTTAAGGCCAGACCAATAGCGATCGGCGCAAAACCTGCCGGCGCGAATTTGTc                 >  1:534374/1‑140 (MQ=255)
                                                                                                                      cAGCCGCCCTGGAAGATAGTAACCGCGGTGCTGCGCGCCGGGTTAACAGAAGTGTTAGTCACCGGAATACTAATTAAGTGAATCAGGGTTAAGGCCAGACCAATAGCGATCGGCGCAAAACCTGCCGGCGCGAATTTGTc                 >  1:145761/1‑140 (MQ=255)
                                                                                                                      cAGCCGCCCTGGAAGATAGTAACCGCGGTACTGCGCGCCGGGTTAACAGAAGTGTTAGTCACCGGAATACTAATTAAGTGAATCAGGGTTAAGGCCAGACCAATAGCGATCGGCGCAAAACCTGCCGGCGCGAATTTGTc                 >  1:507326/1‑140 (MQ=255)
                                                                                                                       aGCCGCCCTGGAAGATAGTAACCGCGGTGCTGCGCGCCGGGTTAACAGAAGTGTTAGTCACCGGAATACTAATTAAGTGAATCAGGGTTAAGGCCAGACCAATAGCGATCGGCGCAAAACCTGCCGGCGCGAATTTGTCg                >  1:402213/1‑140 (MQ=255)
                                                                                                                         ccgccCTGGAAGATAGTAACCGCGGTGCTGCGCGCCGGGTTAACAGAAGTGTTAGTCACCGGAATACTAATTAAGTGAATCAGGGTTAAGGCCAGACCAATAGCGATCGGCGCAAAACCTGCCGGCGCGAATTTGTCGGt              <  2:409422/140‑1 (MQ=255)
                                                                                                                                 gAAGATAGTAACCGCGGTGCTGCGCGCCGGGTTAACAGAAGTGTTAGTCACCGGAATACTAATTAAGTGAATCAGGGTTAAGGCCAGACCAATAGCGATCGGCGCAAAACCTGCCGGCGCGAATTTGTCGGTTGCGCCGt      >  2:872732/1‑140 (MQ=255)
                                                                                                                                     aTAGTAACCGCGGTGCTGCGCGCCGGGTTAACAGAAGTGTTAGTCACCGGAATACTAATTAAGTGAATCAGGGTTAAGGCCAGACCAATAGCGATCGGCGCAAAACCTGCCGGCGCGAATTTGTCGGTTGCGCCGTGGAt  <  2:413881/140‑1 (MQ=255)
                                                                                                                                      tAGTAACCGCGGTGCTGCGCGCCGGGTTAACAGAAGTGTTAGTCACCGGAATACTAATTAAGTGAATCAGGGTTAAGGCCAGACCAATAGCGATCGGCGCAAAACCTGCCGGCGCGAATTTGTCGGTTGCGCCGt      <  1:799307/135‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                       
ACGCCAGGCCTTCTGCGTTAATTAATCACGCTTTTCCAGCAGGGTCCGGTAAATCAGACCACCGATAATGCCGCCGACAATTGGCACCACCCAGAAGAACCACAGTTGTTCTAATGCCCAGCCGCCCTGGAAGATAGCAACCGCGGTGCTGCGCGCCGGGTTAACAGAAGTGTTAGTCACCGGAATACTAATTAAGTGAATCAGGGTTAAGGCCAGACCAATAGCGATCGGCGCAAAACCTGCCGGCGCGAATTTGTCGGTTGCGCCGTGGAT  >  NC_000913/915331‑915603

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑