Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 875,401 G→T 13.8% A23S (GCC→TCC)  yliF → putative membrane‑anchored diguanylate cyclase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913875,4010GT13.8% 156.5 / 4.4 67A23S (GCC→TCC) yliFputative membrane‑anchored diguanylate cyclase
Reads supporting (aligned to +/- strand):  ref base G (39/17);  new base T (2/9);  total (41/26)
Fisher's exact test for biased strand distribution p-value = 2.17e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCTCAGGGGTATCTGTACCAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAATTGATCTTATCAACACGCTAAATCCCACGGAGGAA  >  NC_000913/875252‑875550
                                                                                                                                                     |                                                                                                                                                     
ggCTCAGGGGTATCTGTACCAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTg                                                                                                                                                       <  1:1228871/150‑1 (MQ=255)
ggCTCAGGGGTATCTGTACCAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTg                                                                                                                                                       <  1:1277616/150‑1 (MQ=255)
  cTCAGGGGTAGCTGTTCCATATACCCATGCCATTCTCCGCCTGGGATTAAATTTGAAAATTAGTAAAAGAGTAGTTTTCGTATGTCCATAATCAATAAGTTCGTACTTACAGTCAGTGTGCTGATTTTTATCATGTTTTCAGCAGTTGcc                                                                                                                                                     <  1:1103068/150‑1 (MQ=255)
       ggttATCTGTACCAGAATGCCATGTCTTTATCCGCATGTGATAAAAGTGGAAAATTAGTATAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTTCTTAGATTGAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCgg                                                                                                                                                <  1:826523/147‑1 (MQ=255)
        ggTATCTGTACCAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCggg                                                                                                                                               <  1:1264741/150‑1 (MQ=255)
         gTATCTGTACCAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGa                                                                                                                                              >  2:1242848/1‑150 (MQ=255)
          tATCTGTACCAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAGATTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCGGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGAt                                                                                                                                             <  2:819159/150‑1 (MQ=255)
           aTCTGTACCAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATc                                                                                                                                            <  1:1274777/150‑1 (MQ=255)
             cTGTACCAGAAACCCATGCCATTCTCCGACTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTAGTACTTACAGTCAGTCTTCAGATTTTTATAATGATTTCAGCAGTTGCCTGAGGGAtata                                                                                                                                          >  2:962765/1‑150 (MQ=255)
                tACCAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATTATTATTTATTCTGTTTCCTGCGTGATCTacac                                                                                                                                       <  2:317557/150‑1 (MQ=255)
                    aGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGAAGTTTCCTGCGGGATCTACACTCaa                                                                                                                                   >  1:630699/1‑150 (MQ=255)
                           cATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTaa                                                                                                                            >  2:786630/1‑150 (MQ=255)
                                cATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGaa                                                                                                                       >  1:1330756/1‑150 (MQ=255)
                                        gCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGta                                                                                                               <  1:113015/150‑1 (MQ=255)
                                        gCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTGCGTACTCACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGta                                                                                                               <  1:834427/150‑1 (MQ=255)
                                          cTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGtata                                                                                                             >  1:790416/1‑150 (MQ=255)
                                          cTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGtata                                                                                                             >  2:1060888/1‑150 (MQ=255)
                                              taTAAAAGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGACCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTTTTTTTATTTATTCTGTTTCCTGCGTGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCt                                                                                                         <  2:294703/149‑1 (MQ=255)
                                              gATAAAAGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCt                                                                                                         >  2:1277691/1‑150 (MQ=255)
                                              gATAAAAGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCt                                                                                                         >  2:706320/1‑150 (MQ=255)
                                               aTAAAAGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCtt                                                                                                        >  2:494211/1‑149 (MQ=255)
                                                aaaaaGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAAGCAATAAGTTCGTACTTACATTCAGTCTGCTGATTTTTTTTTTTATTGATTCTGTTTCCTGCGTGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGa                                                                                                       <  2:388166/149‑1 (MQ=255)
                                                  aaaGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTAAAATTGTAAAGGAACGGGTGTATAGCCTGaaa                                                                                                     >  1:795641/1‑150 (MQ=255)
                                                    aGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGTTCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACa                                                                                                   >  1:701942/1‑150 (MQ=255)
                                                    aGTGGAAAATTAGTAAAAGAGTAGTTTAAGTATGTCCAGAATCAATAAGTTAGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTAAGAAGTTGCGTGCGGGATCTACACTCAATTGGTAAAGGAACGGGTGTATTGCCTGAAACa                                                                                                   >  1:485946/1‑150 (MQ=255)
                                                       ggAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTc                                                                                                >  1:424800/1‑150 (MQ=255)
                                                       ggAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTc                                                                                                <  1:588190/150‑1 (MQ=255)
                                                       ggAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAAAGGGTGTATAGCCTGAAACAGTc                                                                                                >  1:634764/1‑150 (MQ=255)
                                                       ggAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTCCTGATTTTTATCATCATTTCATCAGTTTCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGCTGAATATCCAGAAACAGTc                                                                                                >  2:536204/1‑150 (MQ=255)
                                                           aaTTAGTAAAAGAGAAGTTTAAGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGATGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATTCTAAAGTAACGGCTATATAGCATGAAACAGTCCGta                                                                                            >  2:998653/1‑149 (MQ=255)
                                                             ttAGTTAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCTTGATTTAAGCTGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTAt                                                                                          <  2:1323255/150‑1 (MQ=255)
                                                                 tAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGAt                                                                                      <  2:1330756/150‑1 (MQ=255)
                                                                    aaGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCTTACTTACAGTCATTCTGCTTATTTTTTTTTTTATTTATTCTGTGTCCTGCGTGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACt                                                                                   <  2:568946/150‑1 (MQ=255)
                                                                      gagTAGTTTACGTATTTCCAGAATCTATATGTTAGTACTTTGAGTCAGTCTGCTGATTTTTTTTTTTATTTATTCTGTTTCCTGCGTGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGc                                                                                 <  2:152441/150‑1 (MQ=255)
                                                                             ttACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACCGGTGTATAGCCTGAAACAGTCCGTAATTGATACTGCTTTTGCg                                                                          >  2:970741/1‑150 (MQ=255)
                                                                               aCGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATTGTAAAGTAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGt                                                                        <  1:970741/150‑1 (MQ=255)
                                                                                cGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTg                                                                       >  2:990146/1‑150 (MQ=255)
                                                                                cGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTAGCGGTg                                                                       >  2:1188307/1‑150 (MQ=255)
                                                                                  tATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGc                                                                     >  1:1257216/1‑150 (MQ=255)
                                                                                     gTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCtaa                                                                  >  1:787477/1‑150 (MQ=255)
                                                                                      tCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCCGGATCTACACTCAAATGGTTAAGGAACGGGTGTATAGCCTGAAAAAGACCGTTATTGATACTGATTTACCGGTGGCAAat                                                                 >  1:1009696/1‑150 (MQ=255)
                                                                                      gacAGAATCAATCAGTTCGTACTTACAGTCATTCTGCTGATTTTTATCATTATTTCACCAGTTTCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAat                                                                 <  2:245158/148‑1 (MQ=255)
                                                                                       ccATAATCAATAAGTTCGTACTTACAGTCATTCTTATTTTTTTTTTTTTTATTGATTCTGTGTCCTTAGTGATGTACACTCAAATGGTAAAGGATCGGGTGTATAGCCTGAAACAGTCCGTTATTGATACGGCTTTTGCGGTGGCaaaga                                                                <  2:1199680/150‑3 (MQ=255)
                                                                                                   aGTTCGTACTTACAGTCAGTCTGCTGATTTTTTTTATTATTTATTCAGTTTCCTGCGTGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATc                                                    <  2:782855/150‑1 (MQ=255)
                                                                                                    gTTCGTACTTTCAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAACGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCg                                                   <  1:962765/150‑1 (MQ=255)
                                                                                                       cGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCAGAATATCGGCg                                                >  1:552033/1‑150 (MQ=255)
                                                                                                           cTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGATGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGAGGTGGCAAATAATCCTCAATATCGGCGTAga                                            >  1:1246759/1‑149 (MQ=255)
                                                                                                             tACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGt                                          >  2:79913/1‑150 (MQ=255)
                                                                                                               cagtcagtCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAACGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGTTGGCAAATATTGCTGAATATCCGCTTAGCGTgg                                        >  2:166616/1‑150 (MQ=255)
                                                                                                                      tcTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAATTGa                                 <  1:706320/150‑1 (MQ=255)
                                                                                                                        tgatgATTTTTTTTTTGATTTATTCTGTTTCCTGCGTGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAAGATGGCTGAAGATCGGCGTAGCGTGGCAATTGATc                               <  2:1210490/147‑1 (MQ=255)
                                                                                                                            gATTTTTATCATGATTTCCGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCTAATATTGCTGAATATCGGCGTAGCGTGGCAATTGATCTTAt                           >  1:1305196/1‑150 (MQ=255)
                                                                                                                             aTTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAATTGATCTTATc                          <  1:1242848/150‑1 (MQ=255)
                                                                                                                              tttttATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAATTGATCTTATCa                         <  1:511216/150‑1 (MQ=255)
                                                                                                                                tttATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAATTGATCTTATCAac                       <  1:865637/150‑1 (MQ=255)
                                                                                                                                 ttATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAATTGATCTTATCAaca                      >  1:640681/1‑150 (MQ=255)
                                                                                                                                 ttATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAATTGATCTTATCAaca                      >  1:142345/1‑150 (MQ=255)
                                                                                                                                 ttATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGGAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTAGTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAATTGATCTTATCAaca                      <  1:559395/150‑1 (MQ=255)
                                                                                                                                  tATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAATTGATCTTATCAacac                     >  1:784941/1‑150 (MQ=255)
                                                                                                                                   aTCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAACTGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATAGTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCATTTGATCTTATCAACACg                    >  1:1109433/1‑150 (MQ=255)
                                                                                                                                     cATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAATTGATCTTATCAACACGCt                  >  1:149072/1‑150 (MQ=255)
                                                                                                                                        gATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGACGGGGTGTTTAGCCTGAAACAGTCCGTTATTGTTACTGCTTTTGCGGTGGCAATTATTGCTGAATATCGGCGAAGCGTGGCAATTGATCTTATCAACCCGCTaaa               >  1:318571/1‑150 (MQ=255)
                                                                                                                                         aTTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATAATGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGAAATTGATCTTATCAACACGCTAAAt              >  1:967287/1‑150 (MQ=255)
                                                                                                                                             cagcagTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAATTGATCTTATCAACACGCTAAATCCCa          >  1:667044/1‑150 (MQ=255)
                                                                                                                                             cagcagTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAATTGATCTTATCAACACGCTAAATCCCa          >  1:841825/1‑150 (MQ=255)
                                                                                                                                             cagcagTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTAGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAATTGATCTTATCAACACGCTAAATCCCa          >  2:751029/1‑150 (MQ=255)
                                                                                                                                             cagcagTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGACTGAAACAGTCCGATATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGTCAATTGATCTTATCAACACGCTAAATccaa          >  1:200423/1‑148 (MQ=255)
                                                                                                                                               gcagTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAATTGATCTTATCAACACGCTAAATCCCatg        >  2:706873/1‑148 (MQ=255)
                                                                                                                                                    tGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCCTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAATTGATCTTATCAACACCCTAAATCCCACgcagca   >  2:559947/1‑148 (MQ=255)
                                                                                                                                                     gCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAATTGATCTTATCAACACGCTAAATCCCACGGAGGaa  >  2:406494/1‑150 (MQ=255)
                                                                                                                                                     gCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAATTGATCTTATCAACACGCTAAATCCCACGGAGGaa  >  1:142990/1‑150 (MQ=255)
                                                                                                                                                     gCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGCCAAATATTTCTGAATATCGGCGTAGCGTGGCAATTAATCTTATCAACACGCTAAATCCAACGCAGGaa  >  2:139179/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GGCTCAGGGGTATCTGTACCAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAATTGATCTTATCAACACGCTAAATCCCACGGAGGAA  >  NC_000913/875252‑875550

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