Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 875,393 C→A 15.0% 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,3930CA15.0% 208.9 / 20.1 104S20* (TCA→TAA) yliFputative membrane‑anchored diguanylate cyclase
Reads supporting (aligned to +/- strand):  ref base C (50/37);  new base A (2/15);  total (52/52)
Fisher's exact test for biased strand distribution p-value = 9.47e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 107 total aligned reads displayed.

GTTCAGGCTCAGGGGTATCTGTACCAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAATTGATCTTATCAACACGCTAAATCCCA  >  NC_000913/875247‑875542
                                                                                                                                                  |                                                                                                                                                     
gtTCAGGCTCAGGGGTATCTGTACCAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAATTTCGTACTTACATTAAGTCTGCTGATTTTTATCATGATTTcagc                                                                                                                                                    >  1:1641713/1‑150 (MQ=255)
gtTCAGGCTCAGGGGTATCTGTACCAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTcagc                                                                                                                                                    <  1:1390822/150‑1 (MQ=255)
       cTCAGGGGTATCTGTACCAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGcc                                                                                                                                             >  2:509075/1‑150 (MQ=255)
         cAGGGGGATCTGTACCAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTg                                                                                                                                           <  2:796909/150‑1 (MQ=255)
          aGGGGTATCTGTACCAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTTCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGc                                                                                                                                          <  1:597116/150‑1 (MQ=255)
             ggTATCTGTACCAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCggg                                                                                                                                       <  2:1299925/150‑1 (MQ=255)
             ggTATCTGTACCAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCggg                                                                                                                                       <  1:412545/150‑1 (MQ=255)
                aTCAGTACCAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATc                                                                                                                                    <  1:1537091/150‑1 (MQ=255)
                 tCTGTACCAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCt                                                                                                                                   <  2:529108/150‑1 (MQ=255)
                   tGTACCAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTAAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGAGGGATCTaa                                                                                                                                 >  1:284999/1‑149 (MQ=255)
                    gTACCAGAAACCCATGCCATTCTCCGCCAGGGATAAAAGTCCAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTAAAGTCAGTCTGCTGATTTTAATCATGATTTCAGCAGTTGCCTGCGGGATCTaca                                                                                                                                >  2:1191117/1‑150 (MQ=255)
                     tACCAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCACAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTacac                                                                                                                               >  2:1126412/1‑150 (MQ=255)
                      aCCAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGAGCTACACt                                                                                                                              <  1:1612197/150‑1 (MQ=255)
                      aCCAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAAAAAGTTCGTACTTACATTCAGTCTGCTGATTTTTATCATTATTTCAGCAGTTGCCTGCGGGATCTAAACt                                                                                                                              >  1:757996/1‑150 (MQ=255)
                       ccAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTc                                                                                                                             >  1:1131979/1‑150 (MQ=255)
                         aGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCaa                                                                                                                           >  1:1656957/1‑150 (MQ=255)
                             acgcATGCCTTTCTCAGCCTGGGATAAAAGTGGAAAAATAGAAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGATCGTACTTACAGTCTTTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATgg                                                                                                                       <  1:659486/147‑1 (MQ=255)
                             aCCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATgg                                                                                                                       <  1:741199/150‑1 (MQ=255)
                             aCCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGATGATTTTTATCATGATTTCAGCAGTTGCCTGCTGGATCTACACTCAAATgg                                                                                                                       >  1:431675/1‑150 (MQ=255)
                                   gCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAgg                                                                                                                 >  2:443170/1‑150 (MQ=255)
                                   gCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGAGCGTACTTACAGTCATTCTGCTTATTTTTTTTTTTATTTATTCTGTTTCCTGCGTGATCTACACTCAAATGGTAAAgg                                                                                                                 <  2:914374/150‑1 (MQ=255)
                                    ccATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCATTCTGCTGATTTTTATCATGATTTAAGCATTTGCCTGCGGGATCTACACTCAAATGGTAAAGGa                                                                                                                >  1:1167935/1‑150 (MQ=255)
                                    ccATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATTGTAAAGGa                                                                                                                >  2:1212139/1‑150 (MQ=255)
                                       ttCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTACAGAATCAATAAGTGCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACg                                                                                                             <  1:66654/150‑1 (MQ=255)
                                         ctcCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACggg                                                                                                           >  2:1276624/1‑150 (MQ=255)
                                         ctcCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACggg                                                                                                           >  1:70175/1‑150 (MQ=255)
                                           ccGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGgtg                                                                                                         >  1:788009/1‑150 (MQ=255)
                                              ccTGGGATAAAAGTGGAAAATTAGTATAAGAGTAGTTTACGTATGTCCAGAATCAATATGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCGTGCGGGTTCTACACTCTAATGGTAATGGTACGGGTGtat                                                                                                      >  2:1498376/1‑150 (MQ=255)
                                                tGGGATAAAAGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAg                                                                                                    >  1:1447813/1‑150 (MQ=255)
                                                  ggATAAAAGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTCCTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGcc                                                                                                  <  1:992449/150‑1 (MQ=255)
                                                  ggATAAAAGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTTTTATTATTGATTCTGGTTCCTGCGTGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGcc                                                                                                  <  2:1027447/150‑1 (MQ=255)
                                                  ggATAAAAGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGcc                                                                                                  <  1:443170/150‑1 (MQ=255)
                                                   gATAAAAGTGGAAAATTAGTAAAAGAGTAGTTTACTTATGTCCAGAATCAATAAGTTCGTACTAACAGTCAGTCTGCTGATTTTTATCATCATTTCAGCACTTGCCTGCTGGATCTACACTCAAATGGTAAAGCAACGGGTGtataccat                                                                                                 >  2:1583621/1‑146 (MQ=255)
                                                   gATAAAAGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCTTACTTACAGTCAGTCTGCTGATTTTTTTTATTATTTAATCAGTTTCCTGCGTGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCt                                                                                                 <  2:1623381/150‑1 (MQ=255)
                                                   gATAAAAGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCt                                                                                                 >  2:496005/1‑150 (MQ=255)
                                                    aTAAAAGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTTTTTTTATTTATTCTGTTTCCTGCGTGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTg                                                                                                <  2:1375961/150‑1 (MQ=255)
                                                    aTAAAAGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTACAGAATCAATATGTTCGTACTTACAGTCTGTCTGCTGATTTTTTTTTTTATTGATTCTGTTTCCTGCGTGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTg                                                                                                <  2:62608/150‑1 (MQ=255)
                                                     tAAAAGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAATTTCGTACTTACAGTCAGTCTGCTGATTTTTTTATTTATTTAATCTGTTTCCTGCGTGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGa                                                                                               <  2:435308/150‑1 (MQ=255)
                                                        aaGTGGAAAATTACTAAAACAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTTCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAAc                                                                                            <  1:1058408/150‑1 (MQ=255)
                                                          gTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAg                                                                                          >  1:394686/1‑150 (MQ=255)
                                                           tGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGt                                                                                         >  2:1468189/1‑150 (MQ=255)
                                                            ggAAAATTAGTAAAAGAGTAGATTACGTATGCCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTTTTTTTATTGATTCTGTTTCCTGCGTGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTc                                                                                        <  2:1233952/150‑1 (MQ=255)
                                                             taaatTTAGTAAAAGAGTAGTTTACGTATGTCCATAATAAAAAAGTTCGTACTTACAGTCATTATGCTTATTTTTTTTTTTATTTATTCGGTGTCCTGCGTGATCTACACTAAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTcc                                                                                       <  2:623477/149‑1 (MQ=255)
                                                             gAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTcc                                                                                       >  2:1349540/1‑150 (MQ=255)
                                                             gAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTcc                                                                                       >  1:1354167/1‑150 (MQ=255)
                                                              aaaaTTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACTGTCCg                                                                                      >  1:835077/1‑150 (MQ=255)
                                                                  ttAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTAt                                                                                  <  2:1656957/150‑1 (MQ=255)
                                                                  ttAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGATGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTAt                                                                                  >  2:615013/1‑150 (MQ=255)
                                                                   tAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTAtt                                                                                 <  1:1316383/150‑1 (MQ=255)
                                                                    aGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATTATTATTTCTTCTGTTTCCTGCGTGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTg                                                                                <  2:438328/150‑1 (MQ=255)
                                                                      tAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGAt                                                                              <  1:808815/150‑1 (MQ=255)
                                                                       aaaaGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATa                                                                             >  1:721796/1‑150 (MQ=255)
                                                                          cgtgtagtTTACGTATGTCCAGACTCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGAAGTTGCCTGAGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTg                                                                          <  1:1438272/147‑1 (MQ=255)
                                                                          agagTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTg                                                                          <  1:1444932/150‑1 (MQ=255)
                                                                           gagTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAg                                                                                                    <  2:1449733/123‑1 (MQ=255)
                                                                           gagTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAg                                                                                                    >  1:1449733/1‑123 (MQ=255)
                                                                           gagTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGc                                                                         <  1:771775/150‑1 (MQ=255)
                                                                           gagTAGTTTACGTATGTCAAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGc                                                                         >  2:639844/1‑150 (MQ=255)
                                                                                gtTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTg                                                                    >  1:641819/1‑150 (MQ=255)
                                                                                   tACGTATGTCCAGAATCAATAAGTTCGTACTTTCAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCgg                                                                 <  1:1126412/150‑1 (MQ=255)
                                                                                    aCGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGt                                                                >  2:855767/1‑150 (MQ=255)
                                                                                    aCGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGt                                                                <  1:878291/150‑1 (MQ=255)
                                                                                     cGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGgtg                                                                                                         >  2:364960/1‑108 (MQ=255)
                                                                                       tATGTCCTGTATCAATTAGTTCGTACTTAAAGTAAGTCTGCTGTTTTTTTTTTTTATTGATTCTGTGTCCTGCGTGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGc                                                             <  2:1256138/150‑1 (MQ=255)
                                                                                       tATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGc                                                             <  1:417638/150‑1 (MQ=255)
                                                                                       tATGTCAAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTTTTTTTATTGATTCTGTGTCCTGCGTGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGc                                                             <  2:652629/150‑1 (MQ=255)
                                                                                        aTGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTGCGGTGGCaa                                                           >  2:1374122/1‑150 (MQ=255)
                                                                                          gTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCaaa                                                          >  1:1161100/1‑150 (MQ=255)
                                                                                            ccAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAata                                                        >  1:1108497/1‑150 (MQ=255)
                                                                                            ccAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAata                                                        >  1:1697177/1‑150 (MQ=255)
                                                                                             cAGCATCAATCAGTTCGTACTTACATTCAGTCTGCTGATTTTTATCATGATTTCAGCCGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAatat                                                       <  1:1510633/150‑1 (MQ=255)
                                                                                               gAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTg                                                     >  2:1604084/1‑150 (MQ=255)
                                                                                                 aTCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATTATTATTTATGCTGTTTCCTGCGTGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCt                                                   <  2:1652439/150‑1 (MQ=255)
                                                                                                    aaTAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCATGAAACAGTCCGTTATTGAAACTGCTTTTGCGGTGGCAAATATTGCTGaa                                                >  1:93146/1‑150 (MQ=255)
                                                                                                     aTAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAat                                               >  2:413705/1‑150 (MQ=255)
                                                                                                     aTAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAat                                               >  2:1679656/1‑150 (MQ=255)
                                                                                                           tCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGc                                         <  1:509075/150‑1 (MQ=255)
                                                                                                            cGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCg                                        <  1:421295/150‑1 (MQ=255)
                                                                                                                cTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGGATATCGGCGTatc                                    >  1:1086314/1‑148 (MQ=255)
                                                                                                                   aCAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTg                                 <  2:1447813/150‑1 (MQ=255)
                                                                                                                    cagtcagtCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTgg                                >  1:454357/1‑150 (MQ=255)
                                                                                                                       tcagtcTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCaa                             <  1:1588027/150‑1 (MQ=255)
                                                                                                                          gtcTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAATTg                          >  1:1158139/1‑150 (MQ=255)
                                                                                                                             tgttgATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAATTGATc                       <  2:1707395/147‑1 (MQ=255)
                                                                                                                             tgctgATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAATTGATc                       <  1:1468189/150‑1 (MQ=255)
                                                                                                                              gctgATTTTTTTTATGATTGATTCTGTTGCCTGCGTGATCTACACTGAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAATTGATCt                      <  2:1641713/150‑1 (MQ=255)
                                                                                                                              gctgATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAATTGATCt                      <  1:413705/150‑1 (MQ=255)
                                                                                                                               ctgATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAATTGATCtt                     >  2:503172/1‑150 (MQ=255)
                                                                                                                               ctgATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAATTGATCtt                     >  1:1454025/1‑150 (MQ=255)
                                                                                                                                tgATTTATATCAAGATTTAAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGTAACGGGTGTATAGCCTGAAACAGTAAGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCAGAATATCGGCGTAGCGTGGCAATTAATATTa                    >  1:592883/1‑150 (MQ=255)
                                                                                                                                 gATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAATTGATCTTAt                   >  2:627044/1‑150 (MQ=255)
                                                                                                                                   tttttATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAATTGATCTTATCa                 <  2:641819/150‑1 (MQ=255)
                                                                                                                                   tttttATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGAGGTGGCAAATATTGCTCAATATCGGCGTAGCGTGGCAATTGATCTTATCa                 >  2:1442169/1‑150 (MQ=255)
                                                                                                                                        aTCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAATTGATCTTATCAACACg            >  1:1185321/1‑150 (MQ=255)
                                                                                                                                         tCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAATTGATCTTATCAACACGc           >  1:1343208/1‑150 (MQ=255)
                                                                                                                                          cATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAATTGATCTTATCAACACGCt          <  1:754954/150‑1 (MQ=255)
                                                                                                                                          cATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAATTGATCTTATCAACACGCt          >  2:754954/1‑150 (MQ=255)
                                                                                                                                              aTTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATTTTGATGAATATCTGAGTAGCTTGGCAATTTATCTTATCAACACGCTaaaa      >  1:93008/1‑149 (MQ=255)
                                                                                                                                                ttCATCAGTTGCCTGCGGGATATACACTCAACTGGGAAAGGAACGGGTGTATAGCCTGAACCAGTCCGTTAATGATACTGCCTTTGCGGTGGCAAATATTGCTGAATATCGGCGGAGCGTGGCAATTGATCTTATCAACACGCTAAATcc    <  1:1042731/150‑1 (MQ=255)
                                                                                                                                                  cagcagTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTAATGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGAGTGGCAATTGATCTTATCAACACGCTAAATCCCa  >  1:1222665/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
GTTCAGGCTCAGGGGTATCTGTACCAGAAACCCATGCCATTCTCCGCCTGGGATAAAAGTGGAAAATTAGTAAAAGAGTAGTTTACGTATGTCCAGAATCAATAAGTTCGTACTTACAGTCAGTCTGCTGATTTTTATCATGATTTCAGCAGTTGCCTGCGGGATCTACACTCAAATGGTAAAGGAACGGGTGTATAGCCTGAAACAGTCCGTTATTGATACTGCTTTTGCGGTGGCAAATATTGCTGAATATCGGCGTAGCGTGGCAATTGATCTTATCAACACGCTAAATCCCA  >  NC_000913/875247‑875542

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: