Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,407,926 T→A 5.2% intergenic (‑131/‑124) ydaM ← / → ydaN diguanylate cyclase, csgD regulator/putative Zn(II) transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,407,9260TA 5.2% 268.7 / 2.4 91intergenic (‑131/‑124)ydaM/ydaNdiguanylate cyclase, csgD regulator/putative Zn(II) transporter
Reads supporting (aligned to +/- strand):  ref base T (50/36);  new base A (3/2);  total (53/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGAAGTTGTGCGTAATCATTGAGATCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGATCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCATTTTGCTTTTTGTATTCAATTTAGACAGAATTTTATTAATCATTTCAGGGTAATGGGGTGATGAGATGTTGCGTAACAGGGCCAGAAGGCTAGACTACAAAATAATGCGTTGATGATGGAGGCACTGTGGAAGCGATTAAGGGATCGGACGT  >  NC_000913/1407777‑1408075
                                                                                                                                                     |                                                                                                                                                     
tGAAGTTGTGCGTAATCATTGAGATCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGATCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCAttt                                                                                                                                                       >  1:631912/1‑150 (MQ=255)
 gAAGTTGTGCGTAATCCTTGAGATCCCGTTATTTATATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGATCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCAtttt                                                                                                                                                      <  2:174160/150‑1 (MQ=255)
   aGTTGTGCGTAATCATTGAGATCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGATCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCATTTTGc                                                                                                                                                    >  1:718802/1‑150 (MQ=255)
    gTTGTGCGTAATCATTGAGATCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGATCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCATTTTGCt                                                                                                                                                   >  1:586231/1‑150 (MQ=255)
      tgtgCGTAATCATTGAGATCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGATCAGATTATTCTGATAATTATAGAGGCTATTTAACAATTCACACGGTGTGTAGACGGAAAAGCATTTTGCttt                                                                                                                                                 >  2:311071/1‑150 (MQ=255)
      tgtgCGTAATCATTGAGATCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGATCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCATTTTGCttt                                                                                                                                                 >  1:1502204/1‑150 (MQ=255)
           gTAATCATTGAGATCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGATCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCAAAATGCttttccgt                                                                                                                                            <  2:391930/150‑5 (MQ=255)
             aaTCATTGAGATCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGATCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCATTTTGCTTTTTGTAtt                                                                                                                                          <  2:613503/150‑1 (MQ=255)
              aTCATTGAGATCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGATCAGATAATTCTGATAATGTTAGACTCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCATTTTGCTTTTTGTATTc                                                                                                                                         <  2:1197203/150‑1 (MQ=255)
              aTCATTGAGATCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGACACGTCATATGATCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCATTTTGCTTTTTGTATTc                                                                                                                                         <  1:931210/150‑1 (MQ=255)
                cATTGAGATCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGATCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCATTTTGCTTTTTGTATTCaa                                                                                                                                       >  2:1328225/1‑150 (MQ=255)
                cATTGAGATCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGATCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCATTTTGCTTTTTGTATTCaa                                                                                                                                       >  2:1566776/1‑150 (MQ=255)
                 aTTGAGATCCCGTTATTTTTATCTATTTTTGTTGCTATCCGATCGCAAAATAGCCACGTCATATGATCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGTAAAGCATTTTGCTTTTTGTATTCAAt                                                                                                                                      <  2:956034/150‑1 (MQ=255)
                    gagaTCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGATCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCATTTTGCTTTTTGTATTCAAtttt                                                                                                                                   >  2:837615/1‑149 (MQ=255)
                           cGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGATCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGTAAAAGCATTTTGCTTTTTGTATTCAATTTAGaatgaa                                                                                                                            >  2:340046/1‑145 (MQ=255)
                           cGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGATCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCATTTTGCTTTTTGTATTCAATTTAGACAGaa                                                                                                                            >  2:409932/1‑150 (MQ=255)
                                    tATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGATCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCATTTTGCTTTTTGTATTCAATTTAGACAGAATTTTATTaa                                                                                                                   >  1:530850/1‑150 (MQ=255)
                                    tATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGATCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCATTTTGCTTTTTGTATTCAATTTAGACAGAATTTTATTaa                                                                                                                   >  1:1524650/1‑150 (MQ=255)
                                     aTCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGATCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCATTTTGCTTTTTGTATTCAATTTAGACAGAATTTTATTAAt                                                                                                                  <  2:1189628/150‑1 (MQ=255)
                                     aTCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGATCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCATTTTGCTTTTTGTATTCAATTTAGACAGAATTTTATTAAt                                                                                                                  >  2:690568/1‑150 (MQ=255)
                                     aTCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGATCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCATTTTGCTTTTTGTATTCAATTTAGACAGAATTTTATTAAt                                                                                                                  >  2:538225/1‑150 (MQ=255)
                                     aTCAATTTTTGTAGCTATCCGATCGCAAAAAAGCCACGTCATATGATCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCATTTTCCTATTTGTATTCAAATTAGACAGAATTTTATTAAt                                                                                                                  >  1:1555652/1‑150 (MQ=255)
                                      tCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGATCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCATTTTGCTTTTTGTATTCAATTTAGACAGAATTTTATTAATc                                                                                                                 >  1:840615/1‑150 (MQ=255)
                                      tCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGATCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCATTTTGCTTTTTGTATTCAATTTAGACAGAATTTTATTAATc                                                                                                                 <  1:1702605/150‑1 (MQ=255)
                                         aTTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGATCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCATTTTGCTTTTTGTATTCAATTTAGACAGAATTTTATTAATCAtt                                                                                                              <  2:1668770/150‑1 (MQ=255)
                                         aTTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGATCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCATTTTGCTTTTTGTATTCAATTTAGACAGAATTTTATTAATCAtt                                                                                                              <  2:1411470/150‑1 (MQ=255)
                                         aTTTTTGTTGCTACCCGATCGCAAAAAAGCCACGTCATATGATCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCATTTTGCTTTTTGTATTCAATTTAGACAGAATTTTATTAATCAtt                                                                                                              <  2:1020015/150‑1 (MQ=255)
                                               gttgCTATCCGATCGCAAAAAAGCCACGTCATATGATCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCATTTTGCTTTTTGTATTCAATTTAGACAGAATTTTATTAATCATTTCAggg                                                                                                        >  2:933426/1‑150 (MQ=255)
                                                ttgCTATCCGATCGCAAAAAAGCCACGTCATATGATCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTAAACGGAAAAGCATTTTGCTTTTTGTATTCAATTTAGACAGAATTTTATTAATCATTTCAGGGt                                                                                                       >  1:1074057/1‑150 (MQ=255)
                                                  gCTATCCGATCGCAAAAAAGCCACGTCATATGATCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCATTTTGCTTTTTGTATTCAATTTAGACAGAATTTTATTAATCATTTCAGGGTaa                                                                                                     >  2:882237/1‑150 (MQ=255)
                                                   cTATCCGATCGCAAAAAAGCCACGTCATATGATCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCATTTTGCTTTTTGTATTCAATTTACTCAGAATTTTATTAATCATTTCAGGGTaaa                                                                                                    >  1:1458013/1‑149 (MQ=255)
                                                    tATCCGATCGCAAAAAAGCCACGTCATATGATCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCATTTTGCTTTTTGTATTCAATTTAGACAGAATTTTATTAATCATTTCAGGGTAATg                                                                                                   >  1:1307933/1‑150 (MQ=255)
                                                    tATCCGATCGCAAAAAAGCCACGTCATATGATCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCATTTTGCTTTTTGTATTCAATTTAGACAGAATTTTATTAATCATTTCAGGGTAATg                                                                                                   >  1:616819/1‑150 (MQ=255)
                                                      tCCGATCGCAAAAAAGCCACGTCATATGATCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCATTTTGCTTTTTGTATTCAATTTAGACAGAATTTTATTAATCATTTCAGGGTAATggg                                                                                                 <  2:800682/150‑1 (MQ=255)
                                                       ccGATCGCAAAAAAGCCACGTCATATGATCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCATTTTGCTTTTTGTATTCAATTTAGACAGAATTTTATTAATCATTTCAGGGTAATgggg                                                                                                <  2:1228345/150‑1 (MQ=255)
                                                         gATCGCAAAAAAGCCACGTCATATGATCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCATTTTGCTTTTTGTATTCAATTTAGACAGAATTTTATTAATCATTTCAGGGTAATGGGGtg                                                                                              >  2:1101560/1‑150 (MQ=255)
                                                            cGCAAAAAAGCCACGTCATAAGATCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCATTTTGCTTTTTGTATTCAATTTAGACAGAATTTTATTAATCATTTCAGGGTAATGGGGtgatg                                                                                           <  2:19937/150‑1 (MQ=255)
                                                                aaaaaGCCACGTCATATGATCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCATTTTGCTTTTTGTATTCAATTTAGACAGAATTTTATTAAACATTTCAGGGTAATGGGGTGATGAGAt                                                                                       >  2:137866/1‑150 (MQ=255)
                                                                  aaaGCCACGTCATATGATCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCATTTTGCTTTTTGTATTCAATTTAGACAGAATTTTATTAATCATTTCAGTGTAATGGGGTGATGAGATGt                                                                                     >  2:677675/1‑150 (MQ=255)
                                                                     gCCACGTCATATGATCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCATTTTGCTTTTTGTATTCAATTTAGACAGAATTTTATTAATCATTTCAGGGTAATGGGGTGATGAGATGTTGc                                                                                  >  1:1434900/1‑150 (MQ=255)
                                                                      ccACGTCATATGATCAGATAATTCTGATAATGATAGACGCTATTTAACACTACACACGGTTTGTATACGGAAAAGCATTTTGCTTTTTGTATTCAATTTAGACAGAATTTTATTAATCATTTCAGGGTAATGGGGTGATGAGATGTTGCg                                                                                 <  2:485451/150‑1 (MQ=255)
                                                                       cACGTCATATGATCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCATTTTGCTTTTTGTATTCAATTTAGACAGAATTTTATTAATCATTTCAGGGTAATGGGGTGATGAGATGTTGCGt                                                                                <  1:1181186/150‑1 (MQ=255)
                                                                       cACGTCATATGATCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCATTTTGCTTTTTGTATTCAATTTAGACAGAATTTTATTAATCATTTCAGGGTAATGGGGTGATGAGATGTTGCGt                                                                                <  2:636987/150‑1 (MQ=255)
                                                                       cACGTCATATGATCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCATATTGATTATTGTAATCAATTTAGACAGAATTTTATTAATCATTTCAGGGTAATGGGTTGATGATATGTTTCGt                                                                                >  1:281752/1‑150 (MQ=255)
                                                                          gTCATATGATCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCATTTTGCTTTTTGTATTCAATTTAGACAGAATTTTATTAATCATTTCAGGGTAATGGGGTGATGAGATGTTGCGTAAc                                                                             >  2:1723395/1‑150 (MQ=255)
                                                                           tCATATGATCAGATAATTCTGATAATGATAGACTCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCATTTTGCTTTTTGTATTCAATTTAGACAGAATTTTATTAATCATTTCAGGGTAATGGGGTGATGAGATGTTGCGTAACa                                                                            <  2:170319/150‑1 (MQ=255)
                                                                              tatGATCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCATTTTGCTTTTTGTATTCAATTTAGACAGAATTTTATTAATCATTTCAGGGTAATGGGGTGATGAGATGTTGCGTAACAggg                                                                         >  2:1557604/1‑150 (MQ=255)
                                                                               atGATCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCATTTTGCTTTTTGTATTCAATTTAGACAGAATTTTATTAATCATTTCAGGGTAATGGGGTGATGAGATGTTGCGTAACAGGGc                                                                        <  1:716391/150‑1 (MQ=255)
                                                                                 gATCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCATTTTGCTTTTTGTATTCAATTTAGACAGAATTTTATTAATCATTTCAGGGTAATGGGGTGATGAGATGTTGCGTAACAGGGCCa                                                                      >  2:1350577/1‑150 (MQ=255)
                                                                                  aTCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCATTTTGCTTTTTGTATTCAATTTAGACAGAATTTTATTAATCATTTCAGGGTAATGGGGTGATGAGATGTTGCGTAACAGGGCCAg                                                                     <  1:133225/150‑1 (MQ=255)
                                                                                   tCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCATTTTGCTTTTTGTATTCAATTTAGACAGAATTTTATTAATCATTTCAGGGTAATGGGGTGATGAGATGTTGCGTAACAGGGCCAGa                                                                    >  2:1530646/1‑150 (MQ=255)
                                                                                     aGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCATTTTGCTTTTTGTATTCAATTTAGACAGATTTTTATTAATCATTTCAGGGTAATGGGGTGATGAGATGTTGCGTAACAGGGCCAGAAg                                                                  >  2:277175/1‑150 (MQ=255)
                                                                                      gATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCATTTTGCTTTTTGTATTCAATTTAGACAGAATTTTATTAATCATTTCAGGGTAATGGGGTGATGAGATGTTGCGTAACAGGGCCAGAAgg                                                                 >  2:967176/1‑150 (MQ=255)
                                                                                      gATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCATTTTGCTTTTTGTATTCAATTTAGACAGAATTTTATTAATCATTTCAGGGTAATGGGGTGATGAGATGTTGCGTAACAGGGCCAGAAgg                                                                 >  1:795145/1‑150 (MQ=255)
                                                                                       aTAATTCTGATAATGATAGAAGCTATTTAACACTTCACACGGTTTGTATAAGGAAAAGCATTTTGCTTTTTGTATTCAATTTAGCCAGAATTTTATTAATCATTTCAGGGTAATGGGGTGATGAGATGTTGCGTAACAGGGCCAGAAGGc                                                                <  1:284116/150‑1 (MQ=255)
                                                                                        tAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCATTTTGCTTTTTGTATTCAATTTAGACAGAATTTTATTAATCATTTCAGTGTAATGGGGTGATGAGATGTTGCGTAACAGGGCCAGAATGCt                                                               <  1:837615/150‑1 (MQ=255)
                                                                                             cTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCATTTTGCTTTTTGTATTCAATTTAGACAGAATTTTATTAATCATTTCAGGGTAATGGTGTGATGAGATGTTGCGTAACAGGGCCAGAAGGCTAGACt                                                          >  1:1109088/1‑150 (MQ=255)
                                                                                             cTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCATTTTGCTTTTTGTATTCAATTTAGACAGAATTTTATTAATCATTTCAGGGTAATGGGGTGATGAGATGTTGCGTAACAGGGCCAGAAGGCTAGACt                                                          >  1:343091/1‑150 (MQ=255)
                                                                                                 tAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCATTTTGCTTTTTGTATTCAATTTAGACAGAATTTTATTAATCATTTCAGGGTAATGGGGTGATGAGATGTTGCGTAACAGGGCCAGAAGGCTAGACTACaa                                                      <  1:1256582/150‑1 (MQ=255)
                                                                                                   aTGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCATTTTGCTTTTTGTATTCAATTTAGACAGAATTTTATTAATCATTTCAGGGTAATGGGGTGATGAGATGTTGCGTAACAGGGCCAGAAGGCTAGACTACaaaa                                                    >  1:1382611/1‑150 (MQ=255)
                                                                                                    tGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCATTTTGCTTTTTGTATTCAATTTAGACAGAATTTTATTAATCATTTCAGGGTAATGGGGTGATGAGATGTTGCGTAACAGGGCCAGAAGGCTAGACTACAAaat                                                   <  1:1328225/150‑1 (MQ=255)
                                                                                                      aTAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCATTTTGCTTTTTGTATTCAATTTAGACAGAATTTTATTAATCATTTCAGGGTAATGGGGTGATGAGATGTTGCGTAACAGGGCCAGAAGGCTAGACTACAAaataa                                                 <  2:1205529/150‑1 (MQ=255)
                                                                                                            gCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCATTTTGCTTTTTGTATTCAATTTAGACAGAATTTTATTAATCATTTCAGGGTAATGGGGTGATGAGATGTTGCGTAACAGGGCCAGAAGGCTAGACTACAAAATAATGCGtt                                           >  1:1256920/1‑150 (MQ=255)
                                                                                                             cTATTTAACACTTCACACGGTTTGTATACGGAAAAGCATTTTGCTTTATGTATTCAATTTAGACAGAATTTTATTAATCATTTCAGGGTAATGGGGTGATGAGATGTTGCGTAACAGGGCCAGAAGGCTAGACTACAAAATAATGCGttt                                          >  1:383698/1‑149 (MQ=255)
                                                                                                                tttATCACTTCACACGGTTTGTATACGGAAAAGCATTTTGCTTTTTGTATTCAATTTAGACAGAATTTTATTAATCATTTCAGGGTAATGGGGTGATGAGATGTTGCGTAACAGGGCCAGAAGGCTAGACTACAAAATAATGCGTtgatg                                       <  1:135274/150‑1 (MQ=255)
                                                                                                                tttAACACTTCACACGGTTTGTATATGGAAAAGGATTATGCTTTTTGTAGTCAATTTAGACAGAATTTTATTAATCATTTCAGGGTAATGGGGTGATGAGATGTTGCGTAACAGGGCCAGAAGGCTAGACTACAAAATAATGCGTtgatg                                       <  1:84937/150‑1 (MQ=255)
                                                                                                                tttAACACTTCACACGGTTTGTATACGGAAAAGCATTTTGCTTTTTGTATTCAATTTAGACAGAATTTTATTAATCATTTCAGGGTAATGGGGTGATGAGATGTTGCGTAACAGGGCCAGAAGGCTAGACTACAAAATAATGCGTtgatg                                       >  2:1727470/1‑150 (MQ=255)
                                                                                                                 ttAACACTTCACACGGTTTGTATACGGAAAAGCATTTTGCTTTTTGTATTCAATTTAGACAGAATTTTATTAATCATTTCAGGGTAATGGGGTGATGAGATGTTGCGTAACAGGGCCAGAAGGCTAGACTACAAAATAATGCGTtgatga                                      <  1:203951/150‑1 (MQ=255)
                                                                                                                          cacaCGGTTTGTATACGGAAAAGCATTTTGCTTTTTGTATTCAATTTAGACAGAATTTTATTAATCATTTCAGGGTAATTGGGTGATGAGATGTTGCGTAACAGTGCCAGAAGGCTAGACTACAAAATAATGCTTTGATGATGGAGGCAc                             >  1:303958/1‑150 (MQ=255)
                                                                                                                          cacaCGGTTTGTATACGGAAAAGCATTTTGCTTTTTGTATTCAATTTAGACAGAATTTTATTAATCATTTCAGGGTAATGGGGTGATGAGATGTTGCGTAACAGGGCCAGAAGGCTAGACTACAAAATAATGCTTTGATGATGGAGGCaa                             >  1:1293039/1‑149 (MQ=255)
                                                                                                                          cacaCGGTTTGTATACGGAAAAGCATTTTGCTTTTTGTATTCAATTTAGACAGAATTTTATTAATCATTTCAGGGTAATGGGGTGATGAGATGTTGCGTAACAGGGCCAGAAGGCTAGACTACAAAATAATGCGTTGATGATGGAGGCAc                             <  2:572864/150‑1 (MQ=255)
                                                                                                                           acacGGTTTGTATACGGAAAAGCATTTTGCTTTTTGTATTCAATTTAGACAGAATTTTATTAATCATTTCAGGGTAATGGGGTGATGAGATGTTGCGTAACAGGGCCAGAAGGCTAGACTACAAAATAATGCGTTGATGATGGAGGCACt                            <  2:586231/150‑1 (MQ=255)
                                                                                                                            cacGGTTTGTATACGGAAAAGCATTTTGCTTTTTGTATTCAATTTAGACAGAATTTTATTAATCATTTCAGGGTAATGGGGTGATGAGATGTTGCGTAACAGGGCCATAAGGCTAGACTACAAAATAATGCGTTGATGATGGAGGCACtg                           >  2:954677/1‑150 (MQ=255)
                                                                                                                             acGGTTTGTATACGGAAAAGCATTTTGCTTTTTGTATTCAATTTAGACAGAATTTTATTAATCATTTCAGGGTAATGGGGTGATGAGATGTTGCGTAACAGGGCCAGAAGGCTAGACTACAAAATAATGCGTTGATGATGGAGGCACtgt                          <  2:1488007/150‑1 (MQ=255)
                                                                                                                              cGGTTTGTATACGGAAAAGCATTTTGCTTTTTGTATTCAATTTAGACAGAATTTTATTAATCATTTCAGGGTAATGGGGTGATGAGATGTTGCGTAACAGGGCCAGAAGGCTAGACTACAAAATAATGCGTTGATGATGGAGGCACtgtg                         >  1:1559319/1‑150 (MQ=255)
                                                                                                                              cGGTTTGTATACGGAAAAGCATTTTGCTTTTTGTATTCAATTTAGACAGAATTTTATTAATCATTTCAGGGTAATGGGGTGATGAGATGTTGCGTAACAGGGCCAGAAGGCTAGACTACAAAATAATGCGTTGATGATGGAGGCACtgtg                         >  2:1662496/1‑150 (MQ=255)
                                                                                                                                gTTTGTATACGGTAAAGCATTTTTCTTTTTGTATTCAATTTAGAAAGAATTTTATTAATAATTTCAGGGTAATTGGGTGATGAGATGTTGCGTAACAGGGCCAGAAGGCTAGACTACAAAATAATGCGTTGATGATGGAGGCACTGTGGa                       <  1:107423/150‑1 (MQ=255)
                                                                                                                                 tttGTATACGGAAAAGCATTTTGCTTTTTGTATTCAATTTAGACAGAATTTTATTAATCATTTCAGGGTAATGGGGTGATGAGATGTTGCGTAACAGGGCCTGAAGGCTAGACTACAAAATAATGCGTTGATGATGGAGGCACTGTGGaa                      >  2:1266957/1‑150 (MQ=255)
                                                                                                                                 tttGTATACGGAAAAGCATTTTGCTTTTTGTATTCAATTTAGACAGAATTTTATTAATCATTTCAGGGTAATGGGGTGATGAGATGTTGCGTAACAGGGCCAGAAGGCTAGACTACAAAATAATGCGTTGATGATGGAGGCACTGTGGaa                      <  1:882237/150‑1 (MQ=255)
                                                                                                                                  ttGTATACGGAAAAGCATTTTGCTTTTTGTATTCAATTTAGACAGAATTTTATTAATCATTTCAGGGTAATGGGGTGATGAGATGTTGCGTAACAGGGCCAGAAGGCTAGACTACAAAATAATGCGTTGATGATGGAGGCACTGTGGAAg                     <  1:1336876/150‑1 (MQ=255)
                                                                                                                                   tGTATACGGAAAAGCATTTTGCTTTTTGTATTCAATTTAGACAGAATTTTATTAATCATTTCAGGGTAATGGGGTGATGAGATGTTGCGTAACAGGGCCAGAAGGCTAGACTACAAAATAATGCGTTGATGATGGAGGCACTGTGGAAGc                    <  2:840615/150‑1 (MQ=255)
                                                                                                                                     tataCGGAAAAGCATTTTGCTTTTTGTATTCAATTTAGACAGAATTTTATTAATCATTTCAGGGTAATGGGGTGATGAGATGTTGCGTAACAGGGCCAGAAGGCTAGACTACAAAATAATGCGTTGATGATGGAGGCACTGTGGAAGCGa                  >  2:324224/1‑150 (MQ=255)
                                                                                                                                      ataCGGAAAAGCATTTTGCTTTTTGTATTCAATTTAGACAGAATTTTATTAATCATTTCAGGGTAATGGGGTGATGAGATGTTGCGTAACAGGGCCAGAAGGCTAGACTACAAAATAATGCGTTGATGATGGAGGCACTGTGGAAGCGAt                 >  1:567005/1‑150 (MQ=255)
                                                                                                                                      ataCGGAAAAGCATTTTGCTTTTTGTATTCAATTTAGACAGAATTTTATTAATCATTTCAGGGTAATGGGGTGATGAGATGTTGCGTAACAGGGCCAGAAGGCTAGACTACAAAATAATGCGTTGATGATGGAGGCACTGTGGAAGCGAt                 >  1:733342/1‑150 (MQ=255)
                                                                                                                                       taCGGAAAAGCATTTTGCTTTTTGTATTCAATTTAGACAGAATTTTATTAATCATTTCAGGGTAATGGGGTGATGAGATGTTGCGTAACAGGGCCAGAAGGCTAGACTACAAAATAATGCGTTGATGATGGAGGCACTGTGGAAGCGAtt                <  1:1595432/150‑1 (MQ=255)
                                                                                                                                       atacaaaaaGCAAAATGCTTTTTGTTTTCAATTTAGACAGAATTTTATTAATCATTTCAGGGTAATGGGGTGATGAGATGTTGCGTAACAGGGCAAGAAGGCTAGACTACAAAATAATGCGTTGATGATGGAGGCACTGTGGAAGCGAtt                >  1:888057/6‑150 (MQ=255)
                                                                                                                                        tacaaaaaGCAAAATGCTTTTTGTATTCAATTTAGACAGAATTTTATTAATCATTTCAGGGTAATGGGGTGATGAGATGTTGCGTAACAGGGCCAGAAGGCTAGACTACAAAATAATGCGTTGATGATGGAGGCACTGTGGAAGCGATTa               >  1:1172485/5‑150 (MQ=255)
                                                                                                                                         cGGAAAAGCATTTTGCTTTTTGTATTCAATTTAGACAGAATTTTATTAATCATTTCAGGGTAATGGGGTGATGAGATGTTGCGTAACAGGGCCAGAAGGCTAGACTACAAAATAATGCGTTGATGATGGAGGCACTGTGGAAGCGATTaa              <  2:1502204/150‑1 (MQ=255)
                                                                                                                                            aaaaGCATTTTGCTTTTTGTATTCAATTTAGACAGAATTTTATTAATCATTTCAGGGTAATGGGGTGATGAGATGTTGCGTAACAGGGCCAGAAGGCTAGACTACAAAATAATGCGTTGATGATGGAGGCACTGTGGAAGCGATTAAggg           <  2:1307933/150‑1 (MQ=255)
                                                                                                                                            aaaaGCAAAATGCTTTTTGTATTCAATTTAGACAGAATTTTATTAATCATTTCAGGGTAATGGGGTGATGAGATGTTGCGTAACAGGGCCAGAAGGCTAGACTACAAAATAATGCGTTGATGATGGAGGCACTGTGGAAGCGATTAAggg           >  2:532990/1‑150 (MQ=255)
                                                                                                                                              aaGCATTTTGCTTTTTGTATTCAATTTAGACAGAATTTTATTAATCATTTCAGGGTAATGGGGTGATGAGATGTTGCGTAACAGGGCCAGAACGCTAGACTACAAAATAATGCGTTGATGATGGAGGCACTGTGGAAGCGATTAAGGGAt         >  2:1706772/1‑150 (MQ=255)
                                                                                                                                                  aTTTTGCTTTTTGTATTCAATTTAGACAGAATTTTATTAATCATTTCAGGGTAATGGGGTGATGAGATGTTGCGTAACAGGGCCAGAAGGCTAGACTACAAAATAATGCGTTGATGATGGAGGCACTGTGGAAGCGATTAAGGGATCGGa     <  1:652871/150‑1 (MQ=255)
                                                                                                                                                  aTTTTGCTTTTTGTATTCAATTTAGACAGAAATTTATTAATCATTTCAGGGTAATGGGGTGATGAGATGTTGCGTAACAGGGCAAGAAGGCAAGACTACAAAATAATCAGATGATGATAGAGGCACTGAGGAAGCTAATAAGGGATCTGa     >  1:123103/1‑150 (MQ=255)
                                                                                                                                                     ttGCTTTTTGTATTCAATTTAGACAGAATTTTATTAATCATTTCAGGGTAATGGGGTGATGAGATGTTGCGTAACAGGGCCAGAAGGCTAGACTACAAAATAATGCGTTGATGATGGAGGCACTGTGGAAGCGATTAAGGGATCGGACGt  <  2:1074057/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TGAAGTTGTGCGTAATCATTGAGATCCCGTTATTTTTATCAATTTTTGTTGCTATCCGATCGCAAAAAAGCCACGTCATATGATCAGATAATTCTGATAATGATAGACGCTATTTAACACTTCACACGGTTTGTATACGGAAAAGCATTTTGCTTTTTGTATTCAATTTAGACAGAATTTTATTAATCATTTCAGGGTAATGGGGTGATGAGATGTTGCGTAACAGGGCCAGAAGGCTAGACTACAAAATAATGCGTTGATGATGGAGGCACTGTGGAAGCGATTAAGGGATCGGACGT  >  NC_000913/1407777‑1408075

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