Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 875,393 C→A 7.8% S20* (TCA→TAA)  yliF → putative membrane‑anchored diguanylate cyclase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913875,3930CA 7.8% 203.0 / 5.5 72S20* (TCA→TAA) yliFputative membrane‑anchored diguanylate cyclase
Reads supporting (aligned to +/- strand):  ref base C (40/26);  new base A (2/4);  total (42/30)
Fisher's exact test for biased strand distribution p-value = 2.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.51e-01

TTCAGGCTCAGGGGTATCTGTACCAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAATTGATCTTATCAACACGCTAAATCC  >  NC_000913/875248‑875540
                                                                                                                                                 |                                                                                                                                                   
ttCAGGCTCAGGGGTATCTGTACCAGAAACCCATGCCATTCTCCGACTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTcag                                                                                                                                                   >  2:32286/1‑148 (MQ=255)
     gCTCAGGGGTATCTGTACCAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAGTTTACTTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTg                                                                                                                                             >  2:324120/1‑149 (MQ=255)
         aGGGGTATCTGTACCAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTca                                                                                                                                                    >  3:308189/1‑138 (MQ=255)
         aGGGGTATCTGTACCAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTca                                                                                                                                                    <  4:308189/138‑1 (MQ=255)
            ggTATCTGTACCAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCgg                                                                                                                                      <  1:100340/149‑1 (MQ=255)
                tcggTACCAGAAACCCATGCCATTCTCCGCCTGGTATAAAAGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATAAATAAGTTCGTACTTACAGTAAGTCTGCTGATTTTTATAATGATTTCAGCAGTTGCCTGCGGGATc                                                                                                                                  >  8:9780/4‑149 (MQ=255)
                tCTGTACCAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATc                                                                                                                                  >  8:9782/1‑149 (MQ=255)
                tCTGTACCAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATc                                                                                                                                  >  5:198892/1‑149 (MQ=255)
                  tGTACCAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTa                                                                                                                                >  5:331030/1‑149 (MQ=255)
                      ccAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTTTTTTTAATTATTCTGTTTCCTGCGTGATCTacac                                                                                                                             <  2:133199/148‑1 (MQ=255)
                      ccAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTacac                                                                                                                             <  1:317059/148‑1 (MQ=255)
                      ccAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACt                                                                                                                            >  2:236842/1‑149 (MQ=255)
                       cAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACt                                                                                                                            >  4:78005/1‑148 (MQ=255)
                       cAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACt                                                                                                                            >  4:210772/1‑148 (MQ=255)
                       cAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACt                                                                                                                            >  1:279937/1‑148 (MQ=255)
                          aaaCCCATGCCATTCTCCGCCTGGGGTAAAAGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCaaa                                                                                                                        >  2:95265/1‑149 (MQ=255)
                              ccATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTAAGCAGTTGCCTGCGGGATCTACACTCAAATgg                                                                                                                     >  6:158525/1‑148 (MQ=255)
                               cATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAGTTAACGTATGACCAGAAGCAATAAGATCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTAAGCAGTTCCCTGCGGGATCTACACTCAAATGGt                                                                                                                    >  3:203800/1‑148 (MQ=255)
                                    cATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTAAAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGAGGGATCTACACACAAATGGTAAAgg                                                                                                               >  3:268686/1‑148 (MQ=255)
                                     aTTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGac                                                                                                             >  4:209516/1‑148 (MQ=255)
                                         tcCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACggg                                                                                                         >  2:152769/1‑149 (MQ=255)
                                           cGCCTGGGATAAAAGTGGAAAATTATTAAAAGAGTAGTTTACGTCTGTCAAGCATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTTTTATTATTTATTCTGTTTCCTGCGTGATCTACACTCAAATGGTAAAGGAACGGgt                                                                                                        <  6:55748/148‑1 (MQ=255)
                                           cGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGTTTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGgt                                                                                                        <  5:8439/148‑1 (MQ=255)
                                           cGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGGCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGgt                                                                                                        <  1:236842/148‑1 (MQ=255)
                                            gCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGgtgt                                                                                                      <  7:336426/149‑1 (MQ=255)
                                                 ggATAAAAGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATTATTATTTCATCTGTTGCCTGCGTGATCTACACTCAAATGGTAAAGGAACGGGTGtata                                                                                                   <  4:188119/147‑1 (MQ=255)
                                                    tAAAAGTGGAAACTTAGTAAAAGAGTAGTTTACGTATGTCCATCATCAATAATTTCTTACTTACAGGCATGCTGCTGATTTTTATTATTATTTATTCTGTTTCCTGCGTGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTg                                                                                              <  4:203800/149‑1 (MQ=255)
                                                    tAAAAGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGCATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATCCTTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGCTAAAGGAACGGGTGTATAGCCt                                                                                               >  4:212630/1‑148 (MQ=255)
                                                      aaaGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTa                                                                                                                                >  6:17292/1‑113 (MQ=255)
                                                      aaaGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTa                                                                                                                                <  5:17292/113‑1 (MQ=255)
                                                      aaaGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGa                                                                                             <  1:32286/148‑1 (MQ=255)
                                                       aaGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGaaa                                                                                           >  2:318700/1‑149 (MQ=255)
                                                         gTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACa                                                                                         <  1:152769/149‑1 (MQ=255)
                                                            gAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTc                                                                                      >  8:267726/1‑149 (MQ=255)
                                                              aaaTTAGTAAAAGAGTAGTTTACGTATGTCCATAATCAATAAGTTCGTACTTAAAGTCAGTCTGCTGATTTTTATCATGATTTCAGCATTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCg                                                                                    <  5:208021/149‑1 (MQ=255)
                                                              aaaTTAGTAAAAGAGAAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCg                                                                                    <  7:86420/149‑1 (MQ=255)
                                                                    gTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATTATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTAtt                                                                               <  6:198892/148‑1 (MQ=255)
                                                                       aaaGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGACTGAAACAGTCCGTTATTGATa                                                                           >  1:32302/1‑149 (MQ=255)
                                                                         agagTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACt                                                                         >  7:258707/1‑149 (MQ=255)
                                                                                tttACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTg                                                                  <  3:210772/149‑1 (MQ=255)
                                                                                            cAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGGGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAata                                                      >  6:250618/1‑149 (MQ=255)
                                                                                             aGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTAATTTATACTGCTTTTGCGGTGGCAAatat                                                     >  6:294320/1‑149 (MQ=255)
                                                                                              gAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTAtt                                                                               <  1:333227/122‑1 (MQ=255)
                                                                                              gAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTAtt                                                                               >  2:333227/1‑122 (MQ=255)
                                                                                               aaTCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGCACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTg                                                   <  1:318700/149‑1 (MQ=255)
                                                                                                aTCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGc                                                  >  6:240973/1‑149 (MQ=255)
                                                                                                 tCAATAAGTTCGTACTTACAGTCAGTCTGCTTATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCt                                                 <  5:158525/149‑1 (MQ=255)
                                                                                                 tCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCCGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCt                                                 <  3:272762/149‑1 (MQ=255)
                                                                                                 tCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCt                                                 <  8:258707/149‑1 (MQ=255)
                                                                                                 tCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCt                                                 >  3:65400/1‑149 (MQ=255)
                                                                                                 tCAATAACTTCGTACTTACAGTCAGTCTGCTGATTTTTTTTTTTATTTAATCTGTTTCCTGCGTGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCt                                                 <  4:202703/149‑1 (MQ=255)
                                                                                                     tAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGaa                                              >  2:41563/1‑148 (MQ=255)
                                                                                                       aGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAatat                                           <  4:65400/149‑1 (MQ=255)
                                                                                                            gTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGATTTTGAGGTGGCAAATATTGCTGAATATCGGCg                                      >  6:192248/1‑149 (MQ=255)
                                                                                                                  aCAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCg                                 <  1:95265/148‑1 (MQ=255)
                                                                                                                      tcagtcTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATTGCCTGAAACAGTCCGTTATTGCTACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGc                             >  5:133128/1‑148 (MQ=255)
                                                                                                                        agtcTGCTGATTTTTATCATGATTTCATCTGTTTCCTGCGTGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGACACTGCTTTTGCTGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAAt                          <  6:331030/149‑1 (MQ=255)
                                                                                                                          tcTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAAtt                         >  3:62255/1‑148 (MQ=255)
                                                                                                                            tgctgATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAATTGAt                      >  3:243228/1‑149 (MQ=255)
                                                                                                                              ctgATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAATTGATCt                    >  1:80568/1‑149 (MQ=255)
                                                                                                                              ctgATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAATTGATCt                    >  7:283651/1‑149 (MQ=255)
                                                                                                                               tgATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAATTGATCtt                   >  5:83155/1‑149 (MQ=255)
                                                                                                                               tgATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAATTGATCtt                   >  7:113069/1‑149 (MQ=255)
                                                                                                                                 aTTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAATTGATCTTa                  >  2:79458/1‑148 (MQ=255)
                                                                                                                                 aTTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAATTGATCTTa                  >  8:38326/1‑148 (MQ=255)
                                                                                                                                  tttttATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGTGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAATTGATCTTAt                 <  1:79458/148‑1 (MQ=255)
                                                                                                                                     ttATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAATTGATCTTATCaa              <  7:9782/148‑1 (MQ=255)
                                                                                                                                     ttATCATGATTTCAGCAGTTGACGTCGGGATCTACCCTCAAATTGTATAGGAACTTGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGAGTTTGCAAATATTGCTGAATATCTGCGTAGCGTGGCAATTGATCTTCTCaa              <  7:9780/148‑1 (MQ=255)
                                                                                                                                      tATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAATTGATCTTATCAac             <  7:38326/148‑1 (MQ=255)
                                                                                                                                       aTCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAATTGATCTTATCAaca            <  1:41563/148‑1 (MQ=255)
                                                                                                                                        tCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAATTGATCTTATCAACACg          >  4:142174/1‑149 (MQ=255)
                                                                                                                                         cATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAATTGATCTTATCAACACGc         >  6:85811/1‑149 (MQ=255)
                                                                                                                                         cATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAATTGATCTTATCAACACGc         >  3:52295/1‑149 (MQ=255)
                                                                                                                                            gATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAATTGATCTTATCAACACGCTa       >  4:101726/1‑148 (MQ=255)
                                                                                                                                             aTTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAATTGATCTTATCAACACGCTaaa     >  1:177930/1‑149 (MQ=255)
                                                                                                                                                 cagcagTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAATTGATCTTATCAACACGCTAAATcc  >  7:97841/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
TTCAGGCTCAGGGGTATCTGTACCAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAATTGATCTTATCAACACGCTAAATCC  >  NC_000913/875248‑875540

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: