Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,898,635 T→C 6.4% V70A (GTC→GCC)  pdeD → putative c‑di‑GMP phosphodiesterase PdeD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,898,6350TC 6.4% 281.5 / 6.7 94V70A (GTC→GCC) pdeDputative c‑di‑GMP phosphodiesterase PdeD
Reads supporting (aligned to +/- strand):  ref base T (44/44);  new base C (3/3);  total (47/47)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.65e-01

ATTTGCGCACTCGTTACGCTCGCTTCGACCCTGAGCGTGCGATTTATTTCACAGCGTAACTTAAATCAACAACGGGTAGTACAATTCGCCAATCACGCTGTAGAGGAATTAGATAAAGTACTGCTTCCCCTACAGGCAGGTAGCGAAGTCTTGCTTCCGCTGATTGGTCTGCCCTGCTCTGTCGCCCATTTGCCATTACGTAAACAGGCGGCAAAACTCCAAACTGTGCGATCCATTGGCCTGGTGCAAGACGGCACACTTTATTGCTCCAGCATTTTTGGTTATCGCAATGTGCCCG  >  NC_000913/1898487‑1898784
                                                                                                                                                    |                                                                                                                                                     
aTTTGCGCACTCGTTACGCTCGCTTCGACCCTGAGCGTGCGATTTATTTCACAGCGTAACTTAAATCAACAACGGGTAGTACAATTCGCCAATCACGCTGTAGAGGAATTAGATAAAGTACTGCTTCCCCTACAGGCAGGTAGCGAAGTc                                                                                                                                                      >  2:608529/1‑150 (MQ=255)
 tttGCGCACTCGTTACGCTCGCTTCGACCCTGAGCGTGCGATTTATTTCACAGCGTAACTTAAATCAACAACGGGTAGTACAATTCGCCAATCCCGCTGTAGAGGAATTAGATAAAGTACTGCTTCCCCTACAGGCAGGTAGCGAAGTCt                                                                                                                                                     >  1:927576/1‑150 (MQ=255)
  ttGCGCACTCGTTACGCTCGCTTCGACCCTGAGCGTGCGATTTATTTCACAGCGTAACTTAAATCAACAACGGGTAGTACAATTCGCCAATCACGCTGTAGAGGAATTAGATAAAGTACTGCTTCCCCTACAGGCAGGTAGCGAAGTCtt                                                                                                                                                    <  1:25170/150‑1 (MQ=255)
  ttGCGCACTCGTTACGCTCGCTTCGACCCTGAGCGTGCGATTTATTTCACAGCGTAACTTAAATCAACAACGGGTAGTACAATTCGCCAATCACGCTGTAGAGGAATTAGATAAAGTACTGCTTCCCCTACAGGCAGGTAGCGAAGTCtt                                                                                                                                                    <  2:89370/150‑1 (MQ=255)
   tGCGCACTCGTTACGCTCGCTTCGACCCTGAGCGTGCGATTTATTTCACAGCGTAACTTAAATCAACAACGGGTAGAACAATTCGCCAATCACGCTGTAGAGGAATTAGATAAAGTACTGCTTCCCCTACAGGCAGGTAGCGAAGTCTTg                                                                                                                                                   <  1:1649607/150‑1 (MQ=255)
    gcgcACTCGTTACGCTCGCTTCGACCCTGAGCGTGCGATTTATTTCACAGCGTAACTTAAATCAACAACGGGTAGTACAATTCGCCAATCACGCTGTAGAGGAATTAGATAAAGTACTGCTTCCCCTACAGGCAGGTAGCGAAGTCTTGc                                                                                                                                                  >  2:284080/1‑150 (MQ=255)
     cgcACTCGTTACGCTCGCTTCGACCCTGAGCGTGCGATTTATTTCACAGCGTAACTTAAATCAACAACGGGTAGTACAATTCGCCAATCACGCTGTAGAGGAATTAGATAAAGTACTGCTTCCCCTACAGGCAGGTAGCGAAGTCTTGCt                                                                                                                                                 >  2:1634164/1‑150 (MQ=255)
      gcACTCGTTACGCTCGCTTCGACCCTGAGCGTGCGATTTATTTCACAGCGTAACTTAAATCAACAACGGGTAGTACAATTCGCCAATCACGCTGTAGAGGAATTAGATAAAGTACTGCTTCCCCTACAGGCAGGTAGCGAAGTCTTGCtt                                                                                                                                                >  2:711814/1‑150 (MQ=255)
       cACTCGTTACGCTCGCTTCGACCCTGAGCGTGCGATTTATTTCACAGCGTAACTTAAATCAACAACGGGTAGTACAATTCGCCAATCACGCTGTAGAGGAATTAGATAAAGTACTGCTTCCCCTACAGGCAGGTAGCGAAGTCTTGCTTc                                                                                                                                               >  1:553079/1‑150 (MQ=255)
          tCGTTACGCTCGCTTCGACCCTGAGCGTGCGATTTATTTCACAGCGTAACTTAAATCAACAACGGGTAGTACAATTCGCCAATCACGCTGTAGAGGAATTAGATAAAGTACTGCTTCCCCTACAGGCAGGTAGCGAAGTCTTGCTTCCGc                                                                                                                                            <  2:1180564/150‑1 (MQ=255)
          tCGTTACGCTCGCTTCGACCCTGAGCGTGCGATTTATTTCACAGCGTAACTTAAATCAACAACGGGTAGTACAATTCGCCAATCACGCTGTAGAGGAATTAGATAAAGTACTGCTTCCCCTACAGGCAGGTAGCGAAGTCTTGCTTCCGc                                                                                                                                            <  1:1434544/150‑1 (MQ=255)
           cGTTACGCTCGCTTCGACCCTGAGCGTGCGATTTATTTCACAGCGTAACTTAAATCAACAACGGGTAGTACAATTCGCCAATCACGCTGTAGAGGAATTAGATAAAGTACTGCTTCCCCTACAGGCAGGTAGCGAAGTCTTGCTTCCGCt                                                                                                                                           >  2:1578870/1‑150 (MQ=255)
              tACGCTCGCTTCGACCCTGAGCGTGCGATTTATTTCACAGCGTAACTTAAATCAACAACGGGTAGTACAATTCGCCAATCACGCTGTAGAGGAATTAGATAAAGTACTGCTTCCCCTACAGGCAGGTAGCGAAGTCTTGCTTCCGCTGAt                                                                                                                                        >  1:1590508/1‑150 (MQ=255)
                  ctcgctTCGACCCTGAGCGTGCGATTTATTTCACAGCGTAACTTAAATCAACAACGGGTAGTACAATTCGCCAATCACGCTGTAGAGGAATTAGATAAAGTACTGCTTCCCCTACAGGCAGGTAGCGAAGTCTTGCTTCCGCTGATTGGt                                                                                                                                    >  1:784567/1‑150 (MQ=255)
                     gctTCGACCCTGAGCGTGCGATTTATTTCACAGCGTAACTTAAATCAACAACGGGTAGTACAATTCGCCAATCACGCTGTAGAGGAATTAGATAAAGTACTGCTTCCCCTACAGGCAGGTAGCGAAGTCTTGCTTCCGCTGATTGGTCTg                                                                                                                                 >  1:284332/1‑150 (MQ=255)
                      ctTCGACCCTGAGCGTGCGATTTATTTCACAGCGTAACTTAAATCAACAACGGGTAGTACAATTCGCCAATCACGCTGTAGAGGAATTAGATAAAGTACTGCTTCCCCTACAGGCAGGTAGCGAAGTCTTGCTTCCGCTGATTGGTCTGc                                                                                                                                >  1:1953707/1‑150 (MQ=255)
                      ctTCGACCCTGAGCGTGCGATTTATTTCACAGCGTAACTTAAATCAACAACGGGTAGTACAATTCGCCAATCACGCTGTAGAGGAATTAGATAAAGTACTGCTTCCCCTACAGGCAGGTAGCGAAGTCTTGCTTCCGCTGATTGGTCTGc                                                                                                                                <  1:1634164/150‑1 (MQ=255)
                        tCGACCCTGAGCGTGCGATTTATTTCACAGCGTAACTTAAATCAACAACGGGTAGTACAATTCGCCAATCACGCTGTAGAGGAATTAGATAAAGTACTGCTTCCCCTACAGGCAGGTAGCGAAGTCTTGCTTCCGCTGATTGGTCTGccc                                                                                                                              <  1:352823/150‑1 (MQ=255)
                        tCGACCCTGAGCGTGCGATTTATTTCACAGCGTAACTTAAATCAACAACGGGTAGTACAATTCGCCAATCACGCTGTAGAGGAATTAGATAAAGTACTGCTTCCCCTACAGGCAGGTAGCGAAGTCTTGCTTCCGCTGATTGGTCTGccc                                                                                                                              <  2:284770/150‑1 (MQ=255)
                           aCCCTGAGCGTGCGATTTATTTCACAGCGTAACTTAAATCAACAACGGGTAGTACAATTCGCCAATCACGCTGTAGAGGAATTAGATAAAGTACTGCTTCCCCTACAGGCAGGTAGCGAAGTCTTGCTTCCGCTGATTGGTCTGCCCTGc                                                                                                                           <  1:783516/150‑1 (MQ=255)
                           aCCCTGAGCGTGCGATTTATTTCACAGCGTAACTTAAATCAACAACGGGTAGTACAATTCGCCAATCACGCTGTAGAGGAATTAGATAAAGTACTGCTTCCCCTACAGGCAGGTAGCGAAGTCTTGCTTCCGCTGATTGGTCTGCCCTGc                                                                                                                           >  2:1699971/1‑150 (MQ=255)
                             ccTGAGCGTGCGATTTATTTCACAGCGTAACTTAAATCAACACCGGGTAGTACAATTCGCCACTCACGCTGTAGAGGAATTAGATAAAGTACTGCTTCCCCTACAGGCAGGTAGCGAAGTCTTGCTTCCGCTGATTGGTCTGCCCTGctc                                                                                                                         <  1:1746115/150‑1 (MQ=255)
                               tGAGCGTGCGATTTATTTCACAGCGTAACTTAAATCAACAACGGGTAGTACAATTCGCCAATCACGCTGTAGAGGAATTAGATAAAGTACTGCTTCCCCTACAGGCAGGTAGCGAAGTCTTGCTTCCGCTGATTGGTCTGCCCTGCTCTg                                                                                                                       >  2:205404/1‑150 (MQ=255)
                                      gCGATTTATTTCACAGCGTAACTTAAATCAACAACGGGTAGTACAATTCGCCAATCACGCTGTAGAGGAATTAGATAAAGTACTGCTTCCCCTACAGGCAGGTAGCGAAGTCTTGCTTCCGCTGATTGGTCTGCCCTGCTCTGTCGCCCa                                                                                                                >  1:619211/1‑150 (MQ=255)
                                        gATTTATTTCACAGCGTAACTTAAATCAACAACGGGTAGTACAATTCGCCAATCACGCTGTAGAGGAATTAGATAAAGTACTGCTTCCCCTACAGGCAGGTAGCGAAGTCTTGCTTCCGCTGATTGGTCTGCCCTGCTCTGTCGCCCAtt                                                                                                              >  1:1710357/1‑150 (MQ=255)
                                           ttatttCACAGCGTAACTTAAATCAACAACGGGTAGTACAATTCGCCAATCACGCTGTAGAGGAATTAGATAAAGTACTGCTTCCCCTACAGGCAGGTAGCGAAGTCTTGCTTCCGCTGATTGGTCTGCCCTGCTCTGTCGCCCAtttcc                                                                                                           <  1:1114499/150‑3 (MQ=255)
                                           ttatttCACAGCGTAACTTAAATCAACAACGGGTAGTACAATTCGCCAATCACGCTGTAGAGGAATTAGATAAAGTACTGCTTCCCCTACAGGCAGGTAGCGAAGTCTTGCTTCCGCTGATTGGTCTGCCCTGCTCTGTCGCCCATTTGc                                                                                                           <  2:1922032/150‑1 (MQ=255)
                                             atttCACAGCGTAACTTAAATCAACAACGGGTAGTACAATTCGCCAATCACGCTGTAGAGGAATTAGATAAAGTACTGCTTCCCCTACAGGCAGGTAGCGAAGTCTTGCTTCCGCTGATTGGTCTGCCCTGCTCTGTCGCCCATTTGCCa                                                                                                         <  2:1184520/150‑1 (MQ=255)
                                               ttCACAGCGTAACTTAAATCAACAACGGGTAGTACAATTCGCCAATCACGCTGTAGAGGAATTAGATAAAGTACTGCTTCCCCTACAGGCAGGTAGCGAAGTCTTGCTTCCGCTGATTGGTCTGCCCTGCTCTGTCGCCCATTTGCCAtt                                                                                                       <  2:1476039/150‑1 (MQ=255)
                                               ttCACAGCGTAACTTAAATCAACAACGGGTAGTACAATTCGCCAATCACGCTGTAGAGGAATTAGATAAAGTACTGCTTCCCCTACAGGCAGGTAGCGAAGTCTTGCTTCCGCTGATTGGTCTGCCCTGCTCTGTCGCCCATTTGCCAtt                                                                                                       <  1:1735202/150‑1 (MQ=255)
                                                 cacaGCGTAACTTAAATCAACAACGGGTAGTACAATTCGCCAATCACGCTGTAGAGGAATTAGATAAAGTACTGCTTCCCCTACAGGCAGGTAGCGAAGTCTTGCTTCCGCTGATTGGTCTGCCCTGCTCTGTCGCCCATTTGCCATTAc                                                                                                     <  2:1388412/150‑1 (MQ=255)
                                                   cggCGTAACTTAAAGCAACAACGGGTAGTACAATTCGGCACTCACGCTGTAGAGGAAGTAGATAAAGTACTGCTTCCCCTACAGGCAGGTAGCGCAGTCTTGCTTCAGCTGATTGGTCTGCCCTGCTCTGTCGCCCATTTGCCAATAcct                                                                                                   >  2:140482/3‑148 (MQ=255)
                                                      cGTAACTTAAATCAACAACGGGTAGTACAATTCGCCAATCACGCTGTAGAGGAATTAGATAAAGTACTGCTTCCGCTACAGGCAGGTAGCGAAGTCTTGCTTCCGCTGATTGGTCTGCCCTGCTCTGTCGCCCATTTGCCATTACGTaaa                                                                                                >  2:1385044/1‑150 (MQ=255)
                                                       gTAACTTAAATCAACAACGGGTAGTACAATTCGCCAATCACGCTGTAGAGGAATTAGATAAAGTACTGCTTCCCCTACAGGCAGGTAGCGAAGTCTTGCTTCCGCTGATTGGTCTGCCCTGCTCTGTCGCCCATTTGCCATTACGTaatc                                                                                               <  2:3397/150‑3 (MQ=255)
                                                       gTAACTTAAATCAACAACGGGTAGTACAATTCGCCAATCACGCTGTAGAGGAATTAGATAAAGTACTGCTTCCCCTACAGGCAGGTAGCGAAGTCTTGCTTCCGCTGATTGGTCTGCCCTGCTCTGTCGCCCATTTGCCATTACGTAAAc                                                                                               <  2:1953707/150‑1 (MQ=255)
                                                           cTTAAATCAACAACGGGTAGTACAATTCGCCAATCACGCTGTAGAGGAATTAGATAAAGTACTGCTTCCCCTACAGGCAGGTAGCGAAGTCTTGCTTCCGCTGATTGGTCTGCCCTGCTCTGTCGCCCATTTGCCATTACGTAAACAggc                                                                                           >  2:1252908/1‑150 (MQ=255)
                                                              aaaTCAACAACGGGTAGTACAATTCGCCAATCACGCTGTAGAGGAATTAGATAAAGTACTGCTTCCCCTACAGGCAGGTAGCGAAGTCTTGCTTCCGCTGATTGGTCTGCCCTGCTCTGTCGCCCATTTGCCATTACGTAAACAggcggc                                                                                        >  2:347137/1‑150 (MQ=255)
                                                               aaTCAACAACGGGTAGTACAATTCGCCAATCACGCTGTAGAGGAATTAGATAAAGTACTGCTTCCCCTACAGGCAGGTAGCGAAGTCTTGCTTCCGCTGATTGGTCTGCCCTGCTCTGTCGCCCATTTGCCATTACGTAAACAGGCGGCa                                                                                       >  1:1604291/1‑150 (MQ=255)
                                                               aaTCAACAACGGGTAGTACAATTCGCCAATCACGCTGTAGAGGAATTAGATAAAGTACTGCTTCCCCTACAGGCAGGTAGCGAAGTCTTGCTTCCGCTGATTGGTCTGCCCTGCTCTGTCGCCCATTTGCCATTACGTAAACAGGCGGCa                                                                                       >  1:1597395/1‑150 (MQ=255)
                                                                aTCAACAACGGGTAGTACAATTCGCCAATCACGCTGTAGAGGAATTAGATAAAGTACTGCTTCCCCTACAGGCAGGTAGCGAAGCCTTGCTTCCGCTGATTGGTCTGCCCTGCTCTGTCGCCCATTTGCCATTACGTAAACAGGCGGCaa                                                                                      >  2:121546/1‑150 (MQ=255)
                                                                aTCAACAACGGGTAGTACAATTCGCCAATCACGCTGTAGAGGAATTAGATAAAGTACTGCTTCCCCTACAGGCAGGTAGCGAAGCCTTGCTTCCGCTGATTGGTCTGCCCTGCTCTGTCGCCCATTTGCCATTACGTAAACAGGCGGCaa                                                                                      >  2:114373/1‑150 (MQ=255)
                                                                aTCAACAACGGGTAGTACAATTCGCCAATCACGCTGTAGAGGAATTAGATAAAGTACTGCTTCCCCTACAGGCAGGTAGCGAAGCCTTGCTTCCGCTGATTGGTCTGCCCTGCTCTGTCGCCCATTTGCCATTACGTAAACAGGCGGCaa                                                                                      >  2:1717639/1‑150 (MQ=255)
                                                                       acGGGTAGTACAATTCGCCAATCACGCTGTAGAGGAATTAGATAAAGTACTGCTTCCCCTACAGGCAGGTAGCGAAGTCTTGCTTCCGCTGATTGGTCTGCCCTGCTCTGTCGCCCATTTGCCATTACGTAAACAGGCGGCAAAACTCCa                                                                               >  2:165347/1‑150 (MQ=255)
                                                                       acGGGTAGTACAATTCGCCAATCACGCTGTAGAGGAATTAGATAAAGTACTGCTTCCCCTACAGGCAGGTAGCGAAGTCTTGCTTCCGCTGATTGGTCTGCCCTGCTCTGTCGCCCATTTGCCATTACGTAAACAGGCGGCAAAACTCCa                                                                               >  2:362398/1‑150 (MQ=255)
                                                                       acGGGTAGTACAATTCGCCAATCACGCTGTAGAGGAATTAGATAAAGTACTGCTTCCCCTACAGGCAGGTAGCGAAGTCTTGCTTCCGCTGATTGGTCTGCCCTGCTCTGTCGCCCATTTGCCATTACGTAAACAGGCGGCAAAACTCCa                                                                               <  2:553079/150‑1 (MQ=255)
                                                                         gggTAGTACAATTCGCCAATCACGCTGTAGAGGAATAAGATAAAGTACTGCTTCCCCTACAGGCAGGTAGCGAAGTCTTGCTTCCGCTGATTGGTCTGCCCTGCTCTGTCGCCCATTTGCCATTACGTAAACAGGCGGCAAAACTCCaaa                                                                             <  2:1938327/150‑1 (MQ=255)
                                                                          ggTAGTACAATTCGCCAATCACGCTGTAGAGGAATTAGATAAAGTACTGCTTCCCCTACAGGCAGGTAGCGAAGTCTTGCTTCCGCTGATTGGTCTGCCCTGCTCTGTCGCCCATTTGCCATTACGTAAACAGGCGGCAAAACTCCAAAc                                                                            <  1:211781/150‑1 (MQ=255)
                                                                              gtaCAATTCGCCAATCACGCTGTAGAGGAATTAGATAAAGTACTGCTTCCCCTACAGGCAGGTAGCGAAGTCTTGCTTCCGCTGATTGGTCTGCCCTGCTCTGTCGCCCATTTGCCATTACGTAAACAGGCGGCAAAACTCCAAACtgtg                                                                        <  2:57009/150‑1 (MQ=255)
                                                                                 cAATTCGCCAATCACGCTGTAGAGGAATTAGATAAAGTACTGCTTCCCCTACAGGCAGGTAGCGAAGTCTTGCTTCCGCTGATTGGTCTGCCCTGCTCTGTCGCCCATTTGCCATTACGTAAACAGGCGGCAAAACTCCAAACTGTGCGa                                                                     <  1:1255363/150‑1 (MQ=255)
                                                                                  aaTTCGCCAATCACGCTGTAGAGGAATTAGATAAAGTACTGCTTCCCCTACAGGCAGGTAGCGAAGTCTTGCTTCCGCTGATTGGTCTGCCCTGCTCTGTCGCCCATTTGCCATTACGTAAACAGGCGGCAAAACTCCAAACTGTGCGAt                                                                    <  1:726327/150‑1 (MQ=255)
                                                                                   aTTCGCCAATCACGCTGTAGAGGAATTAGATAAAGTACTGCTTCCCCTACAGGCAGGTAGCGAAGTCTTGCTTCCGCTGATTGGTCTGCCCTGCTCTGTCGCCCATTTGCCATTACGTAAACAGGCGGCAAAACTCCAAACTGTGCGATc                                                                   <  1:1699971/150‑1 (MQ=255)
                                                                                    ttCGCCAATCACGCTGTAGAGGAATTAGATAAAGTACTGCTTCCCCTACAGGCAGGTAGCGAAGTCTTGCTTCCGCTGATTGGTCTGCCCTGCTCTGTCGCCCATTTGCCATTACGTAAACAGGCGGCAAAACTCCAAACTGTGCGATcc                                                                  <  2:598285/150‑1 (MQ=255)
                                                                                      cGCCAATCACGCTGTAGAGGAATTAGATAAAGTACTGCTTCCCCTACAGGCAGGTAGCGAAGTCTTGCTTCCGCTGATTGGTCTGCCCTGCTCTGTCGCCCATTTGCCATTACGTAAACAGGCGGCAAAACTCCAAACTGTGCGATCCAt                                                                >  2:1706195/1‑150 (MQ=255)
                                                                                        ccAATCACGCTGTAGAGGAATTAGATAAAGTACTGCTTCCCCTACAGGCAGGTAGCGAAGTCTTGCTTCCGCTGATTGGTCTGCCCTGCTCTGTCGCCCATTTGCCATTACGTAAACAGGCGGCAAAACTCCAAACTGTGCGATCCATTg                                                              <  1:165347/150‑1 (MQ=255)
                                                                                             cACGCTGTAGAGGAATTAGATAAAGTACTGCTTCCCCTACAGGCAGGTAGCGAAGTCTTGCTTCCGCTGATTGGTCTGCCCTGCTCTGTCGCCCATTTGCCATTACGTAAACAGGCGGCAAAACTCCAAACTGTGCGATCCATTGGCCTg                                                         >  1:1353190/1‑150 (MQ=255)
                                                                                               cGCTGTAGAGGAATTAGATAAAGTACTGCTTCCCCTACAGGCAGGTAGCGAAGTCTTGCTTCCGCTGATTGGTCTGCCCTGCTCTGTCGCCCATTTGCCATTACGTAAACAGGCGGCAAAACTCCAAACTGTGCGATCCATTGGcc                                                           >  2:1903803/1‑146 (MQ=255)
                                                                                               cGCTGTAGAGGAATTAGATAAAGTACTGCTTCCCCTACAGGCAGGTAGCGAAGTCTTGCTTCCGCTGATTGGTCTGCCCTGCTCTGTCGCCCATTTGCCATTACGTAAACAGGCGGCAAAACTCCAAACTGTGCGATCCATTGGcc                                                           <  1:1903803/146‑1 (MQ=255)
                                                                                                gCTGTAGAGGAATTAGATAAAGTACTGCTTCCCCTACAGGCAGGTAGCGAAGTCTTGCTTCCGCTGATTGGTCTGCCCTGCTCTGTCGCCCATTTGCCATTACGTAAACAGGCGGCAAAACTCCAAACTGTGCGATCCATTGGCCTGGTg                                                      <  1:1310452/150‑1 (MQ=255)
                                                                                                 cTGTAGAGGAATTAGATAAAGTACTGCTTCCCCTACAGGCAGGTAGCGAAGTCTTGCTTCCGCTGATTGGTCTGCCCTGCTCTGTCGCCCATTTGCCATTACGTAAACAGGCGGCAAAACTCCAAACTGTGCGATCCATTGGCCTGGTGc                                                     <  2:443169/150‑1 (MQ=255)
                                                                                                  tGTAGAGGAATTAGATAAAGTACTGCTTCCCCTACAGGCAGGTAGCGAAGTCTTGCTTCCGCTGATTGGTCTGCCCTGCTCTGTCGCCCATTTGCCATTACGTAAACAGGCGGCAAAACTCCAAACTGTGCGATCCATTGGCCTGGTGCa                                                    <  1:711814/150‑1 (MQ=255)
                                                                                                        ggAATTAGATAAAGTACTGCTTCCCCTACAGGCAGGTAGCGAAGTCTTGCTTCCGCTGATTGGTCTGCCCTGCTCTGTCGCCCATTTGCCATTACGTAAACAGGCGGCAAAACTCCAAACTGTGCGATCCATTGGCCTGGTGCAAGACgg                                              >  1:133334/1‑150 (MQ=255)
                                                                                                        ggAATTAGATAAAGTACTGCTTCCCCTACAGGCAGGTAGCGAAGTCTTGCTTCCGCTGATTGGTCTGCCCTGCTCTGTCGCCCATTTGCCATTACGTAAACAGGCGGCAAAACTCCAAACTGTGCGATCCATTGGCCTGGTGCAAGACgg                                              >  2:204779/1‑150 (MQ=255)
                                                                                                         gAATTAGATAAAGTACTGCTTCCCCTACAGGCAGGTAGCGAAGTCTTGCTTCCGCTGATTGGTCTGCCCTGCTCTGTCGCCCATTTGCCATTACGTAAACAGGCGGCAAAACTCCAAACTGTGCGATCCATTGGCCTGGTGCAAGACGGc                                             <  1:1256201/150‑1 (MQ=255)
                                                                                                             tAGATAAAGTACTGCTTCCCCTACAGGCAGGTAGCGAAGTCTTGCTTCCGCTGATTGGTCTGCCCTGCTCTGTCGCCCATTTGCCATTACGTAAACAGGCGGCAAAACTCCAAACTGTGCGATCCATTGGCCTGGTGCAAGACGGCacac                                         >  1:1558642/1‑150 (MQ=255)
                                                                                                               gATAAAGTACTGCTTCCCCTACAGGCAGGTAGCGAAGTCTTGCTTCCGCTGATTGGTCTGCCCTGct                                                                                                                          <  2:1496808/67‑1 (MQ=255)
                                                                                                               gATAAAGTACTGCTTCCCCTACAGGCAGGTAGCGAAGTCTTGCTTCCGCTGATTGGTCTGCCCTGct                                                                                                                          >  1:1496808/1‑67 (MQ=255)
                                                                                                                aTAAAGTACTGCTTCCCCTACAGGCAGGTAGCGAAGCCTTGCTTCCGCTGATTGGTCTGCCCTGCTCTGTCGCCCATTTGCCATTACGTAAACAGGCGGCAAAACTCCAAACTGTGCGATCCATTGGCCTGGTGCAAGACGGCACACttt                                      <  1:114373/150‑1 (MQ=255)
                                                                                                                aTAAAGTACTGCTTCCCCTACAGGCAGGTAGCGAAGCCTTGCTTCCGCTGATTGGTCTGCCCTGCTCTGTCGCCCATTTGCCATTACGTAAACAGGCGGCAAAACTCCAAACTGTGCGATCCATTGGCCTGGTGCAAGACGGCACACttt                                      <  1:121546/150‑1 (MQ=255)
                                                                                                                aTAAAGTACTGCTTCCCCTACAGGCAGGTAGCGAAGCCTTGCTTCCGCTGATTGGTCTGCCCTGCTCTGTCGCCCATTTGCCATTACGTAAACAGGCGGCAAAACTCCAAACTGTGCGATCCATTGGCCTGGTGCAAGACGGCACACttt                                      <  1:1717639/150‑1 (MQ=255)
                                                                                                                 tAAAGTACTGCTTCCCCTACAGGCAGGTAGCGAAGTCTTGCTTCCGCTGATTGGTCTGCCCTGCTCTGTCGCCCATTTGCCATTACGTAAACAGGCGGCAAAACTCCAAACTGTGCGATCCATTGGCCTGGTGCAAGACGGCACACTTTa                                     >  2:1656648/1‑150 (MQ=255)
                                                                                                                  aaaGTACTGCTTCCCCTACAGGCACGTAGCGAACTCTAGCTTCCTCTCATTGGTCTGCACTGCTCTGTCGCCCATTTGCCATTAAGTAAACAGGCGGCAAAACTCCAAACAGTGCCATCCATTGGCCTGGTGCAAGACGTCACACTTTAt                                    >  2:874623/1‑150 (MQ=255)
                                                                                                                    aGTACTGCTTCCCCTACAGGCAGGTAGCGAAGTCTTGCTTCCGCTGATTGGTCTGCCCTGCTCTGTCGCCCATTTGCCATTACGTAAACAGGCGGCAAAACTCCAAACTGTGCGATCCATTGGCCTGGTGCAAGACGGCACACTTTATTg                                  <  1:284080/150‑1 (MQ=255)
                                                                                                                    aGTACTGCTTCCCCTACAGGCAGGTAGCGAAGTCTTGCTTCCGCTGATTCGTCTGCCCTGCTCTGTCGCCCATTTGCCATTACGTAAACAGGCGGCAAAACTCCAAACTGTGCGATCCATTGGCCTGGTGCAAGACGGCACACTTTATTg                                  >  2:992166/1‑150 (MQ=255)
                                                                                                                     gTACTGCTTCCCCTACAGGCAGGTAGCGAAGTCTTGCTTCCGCTGATTGGTCTGCCCTGCTCTGTCGCCCATTTGCCATTACGTAAACAGGCGGCAAAACTCCAAACTGTGCGATCCATTGGCCTGGTGCAAGACGGCACACTTTATTGc                                 >  1:627858/1‑150 (MQ=255)
                                                                                                                     gTACTGCTTCCCCTACAGGCAGGTAGCGAAGTCTTGCTTCCGCTGATTGGTCTGCCCTGCTCTGTCGCCCATTTGCCATTACGTAAACAGGCGGCAAAACTCCAAACTGTGCGATCCATTGGCCTGGTGCAAGACGGCACACTTTATTGc                                 >  2:1783999/1‑150 (MQ=255)
                                                                                                                        cTGCTTCCCCTACAGGCAGGTAGCGAAGTCTTGCTTCCGCTGATTGGTCTGCCCTGCTCTGTCGCCCATTTGCCATTACGTAAACAGGCGGCAAAACTCCAAACTGTGCGATCCATTGGCCTGGTg                                                      >  2:1952221/1‑126 (MQ=255)
                                                                                                                        cTGCTTCCCCTACAGGCAGGTAGCGAAGTCTTGCTTCCGCTGATTGGTCTGCCCTGCTCTGTCGCCCATTTGCCATTACGTAAACAGGCGGCAAAACTCCAAACTGTGCGATCCATTGGCCTGGTg                                                      <  1:1952221/126‑1 (MQ=255)
                                                                                                                          gCTTCCCCTACAGGCAGGTAGCGAAGTCTTGCTTCCGCTGATTGGTCTGCCCTGCTCTGTCGCCCATTTGCCATTACGTAAACAGGCGGCAAAACTCCAAACTGTGCGATc                                                                   <  1:630163/111‑1 (MQ=255)
                                                                                                                          gCTTCCCCTACAGGCAGGTAGCGAAGTCTTGCTTCCGCTGATTGGTCTGCCCTGCTCTGTCGCCCATTTGCCATTACGTAAACAGGCGGCAAAACTCCAAACTGTGCGATc                                                                   >  2:630163/1‑111 (MQ=255)
                                                                                                                            ttCCCCTACAGGCAGGTAGCGAAGTCTTGCTTCCGCTGATTGGTCTGCCCTGCTCTGTCGCCCATTTGCCATTACGTAAACAGGCGGCAAAACTCCAAACTGTGCGATCCATTGGCCTGGTGCAAGACGGCACACTTTATTGCTCCAGCa                          <  2:374195/150‑1 (MQ=255)
                                                                                                                             tCCCCTACAGGCAGGTAGCGAAGTCTTGCTTCCGCTGATTGGTCTGCCCTGCTCTGTCGCCCATTTGCCATTACGTAAACAGGCGGCAAAACTCCAAACTGTGCGATCCATTGGCCTGGTGCAAGACTGCACACTTTATTGCTCCAGCAt                         <  1:1755218/150‑1 (MQ=255)
                                                                                                                                  tACAGGCAGGTAGCGAAGTCTTGCTTCCGCTGATTGGTCTGCCCTGCTCTGTCGCCCATTTGCCATTACGTAAACAGGCGGCAAAACTCCAAACTGTGCGATCCATTGGCCTGGTGCAAGACGGCACACTTTATTGCTCCAGCATTTTTg                    >  2:1636075/1‑150 (MQ=255)
                                                                                                                                  tACAGGCAGGTAGCGAAGTCTTGCTTCCGCTGATTGGTCTGCCCTGCTCTGTCGCCCATTTGCCATTACGTAAACAGGCGGCAAAACTCCAAACTGTGCGATCCATTGGCCTGGTGCAAGACGGCACACTTTATTGCTCCAGCATTTTTg                    <  1:1636075/150‑1 (MQ=255)
                                                                                                                                   aCAGGCAGGTAGCGAAGTCTTGCTTCCGCTGATTGGTCTGCCCTGCTCTGTCGCCCATTTGCCATTACGTAAACAGGCGGCAAAACTCCAAACTGTGCGATCCATTGGCCTGGTGCAAGACGGCACACTTTATTGCTCCAGCATTTTTgg                   <  1:347137/150‑1 (MQ=255)
                                                                                                                                    caggcaggTAGCGAAGTCTTGCTTCCGCTGATTGGTCTGCCCTGCTCTGTCGCCCATTTGCCATTACGTAAACAGGCGGCAAAACTCCAAACTGTGCGATCCATTGGCCTGGTGCAAGACGGCACACTTTATTGCTCCAGCATTTTTGGt                  >  2:1662398/1‑150 (MQ=255)
                                                                                                                                    caggcaggTAGCGAAGTCTTGCTTCCGCTGATTGGTCTGCCCTGCTCTGTCGCCCATTTGCCATTACGTAAACAGGCGGCAAAACTCCAAACTGTGCGATCCATTGGCCTGGTGCAAGACGGCACACTTTATTGCTCCAGCATTTTTGGt                  <  2:695653/150‑1 (MQ=255)
                                                                                                                                      ggcaggTAGCGAAGTCTTGCTTCCGCTGATTGGTCTGCCCTGCTCTGTCGCCCATTTGCCATTACGTAAACAGGCGGCAAAACTCCAAACTGTGCTATCCATTGGCCTGGTGCAAGACGGCACACTTTATTGCTCCAGCATTTTTGGTTa                >  1:1986727/1‑150 (MQ=255)
                                                                                                                                          ggTAGCGAAGTCTTGCTTCCGCTGATTGGTCTGCCCTGCTCTGTCGCCCATTTGCCATTACGTAAACAGGCGGCAAAACTCCAAACTGTGCGATCCATTGGCCTGGTGCAAGACGGCACACTTTATTGCTCCAGCATTTTTGGTTATCGc            >  2:1456125/1‑150 (MQ=255)
                                                                                                                                          ggTAGCGAAGTCTTGCTTCCGCTGATTGGTCTGCCCTGCTCTGTCGCCCATTTGCCATTACGTAAACAGGCGGCAAAACTCCAAACTGTGCGATCCATTGGCCTGGTGCAAGACGGCACACTTTATTGCTCCAGCATTTTTGGTTATCGc            <  1:362398/150‑1 (MQ=255)
                                                                                                                                           gTAGCGAAGTCTTGCTTCCGCTGATTGGTCTGCCCTGCTCTGTCGCCCATTTGCCATTACGTAAACAGGCGGCAAAACTCCAAACTGTGCGATCCATTGGCCTGGTGCAAGACGGCACACTTTATTGCTCCAGCATTTTTGGTTATCGCa           <  1:205404/150‑1 (MQ=255)
                                                                                                                                           gTAGCGAAGTCTTGCTTCCGCTGATTGGTCTGCCCTGCTCTGTCGCCCATTTGCCATTACGTAAACAGGCGGCAAAACTCCAAACTGTGCGATCCATTGGCCTGGTGCAAGACGGCACACTTTATTGCTCCAGCATTTTTGGTTATCGCa           >  2:1124222/1‑150 (MQ=255)
                                                                                                                                               cGAAGTCTTGCTTCCGCTGATTGGTCTGCCCTGCTCTGTCGCCCATTTGCCATTACGTAAACAGGCGGCAAAACTCCAAACTGTGCGATCCATTGGCCTGGTGCAAGACGGCACACTTTATTGCCCCAGCATTTTTGGTTATCGCAAtgt       >  1:1217923/1‑150 (MQ=255)
                                                                                                                                                  aGTCTTGCTTCCGCTGATTGGTCTGCCCTGCTCTGTCGCCCATTTGCCATTACGTAAACAGGCGGCAAAACTCCAAACTGTGCGATCCATTGGCCTGGTGCAAGACGGCACACTTTATTGCTCCAGCATTTTTGGTTATCGCAATGTGcc    >  1:1333138/1‑150 (MQ=255)
                                                                                                                                                  aGTCTTGCTTCCGCTGATTGGTCTGCCCTGCTCTGTCGCCCATTTGCCATTACGTAAACAGGCGGCAAAACTCCAAACTGTGCGATCCATTGGCCTGGTGCAAGACGGCACACTTTATTGCTCCAGCATTTTTGGTTATCGCAATGTGcc    <  2:627858/150‑1 (MQ=255)
                                                                                                                                                    tCTTGCTTCCGCTGATTGGTCTGCCCTGCTCTGTCGCCCATTTGCCATTACGTAAACAGGCGGCAAAACTCCAAACTGTGCGATCCATTGGCCTGGTGCAAGACGGCACACTTTATTGCTCCAGCATTTTTGGTTATCGCAATGTGCCcg  >  1:788134/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
ATTTGCGCACTCGTTACGCTCGCTTCGACCCTGAGCGTGCGATTTATTTCACAGCGTAACTTAAATCAACAACGGGTAGTACAATTCGCCAATCACGCTGTAGAGGAATTAGATAAAGTACTGCTTCCCCTACAGGCAGGTAGCGAAGTCTTGCTTCCGCTGATTGGTCTGCCCTGCTCTGTCGCCCATTTGCCATTACGTAAACAGGCGGCAAAACTCCAAACTGTGCGATCCATTGGCCTGGTGCAAGACGGCACACTTTATTGCTCCAGCATTTTTGGTTATCGCAATGTGCCCG  >  NC_000913/1898487‑1898784

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