Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,145,086 T→C 5.7% T607T (ACT→ACC dgcE → putative diguanylate cyclase DgcE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,145,0860TC 5.7% 217.0 / 2.7 70T607T (ACT→ACCdgcEputative diguanylate cyclase DgcE
Reads supporting (aligned to +/- strand):  ref base T (30/36);  new base C (2/2);  total (32/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 = 8.69e-01

ATTACGCTTGATTCCATCGGTGAAGCCGTGGTCTGTATTGATATGGCGATGAAAATTACCTTTATGAATCCAGTGGCGGAGAAGATGAGCGGCTGGACGCAGGAAGAAGCGTTAGGTGTTCCGCTCCTGACGGTGTTGCATATTACTTTTGGCGACAACGGACCATTAATGGAGAACATTTACAGTGCCGACACCTCACGTTCCGCGATCGAGCAAGATGTGGTGTTGCACTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGC  >  NC_000913/2144940‑2145233
                                                                                                                                                  |                                                                                                                                                   
aTTACGCTTGATTCCATCGGTGAAGCCGTGGTCTGTATTGATATGGCGATGAAAATTACCTTTATGAATCCAGTGGCGGAGAAGATGAGCGGCTGGACGCAGGAAGAAGCGTTAGGTGTTCCGCTCCTGACGGTGTTGCATATTACtttt                                                                                                                                                  >  2:1719491/1‑150 (MQ=255)
 ttACGCTTGATTCCATCGGTGAAGCCGTGGTCTGTATTGATATGGCGATGAAAATTACCTTTATGAATCCAGTGGCGGAGAAGATGAGCGGCTGGACGCAGGAAGAAGCGTTAGGTGTTCCGCTCCTGACGGTGTTGCATATTACTTTTg                                                                                                                                                 <  1:1234/150‑1 (MQ=255)
   aCGCTTGATTCCATCGGTGAAGCCGTGGTCTGTATTGATATGGCGATGAAAATTACCTTTATGAATCCAGTGGCGGAGAAGATGAGCGGCTGGACGCAGGAAGAAGCGTTAGGTGTTCCGCTCCTGACGGTGTTGCATATTACTTTTGGc                                                                                                                                               <  1:1150219/150‑1 (MQ=255)
        tGATTCCATCGGTGAAGCCGTGGTCTGTATTGATATGGCGATGAAAATTACCTTTATGAATCCAGTGGCGGAGAAGATGAGCGGCTGGACGCAGGAAGAAGCGTTAGGTGTTCCGCTCCTGACGGTGTTGCATATTACTTTTGGCGACaa                                                                                                                                          >  1:603468/1‑150 (MQ=255)
         caTTCCATCGGTGAAGCCTTGGTCTTTATTGATATGGCGATGAAAATTACCTTTATGAATCCAGTGGCGGAGAAGATGAGCGGCTGGACGCAGGAAGAAGCGTTAGGTGTTCCGCTCCTGACTGTGTTGCATATTACTTTTGGCGACAAc                                                                                                                                         <  2:223822/149‑1 (MQ=255)
          aTTCCATCGGTGAAGCCGTGGTCTGTATTGATATGGCGATGAAAATTACCTTTATGAATCCAGTGGCGGAGAAGATGAGCGGCTGGACGCAGGAAGAAGCGTTAGGTGTTCCGCTCCTGACGGTGTTGCATATTACTTTTGGCGACAACg                                                                                                                                        >  1:1618184/1‑150 (MQ=255)
           ttCCATCGGTGAAGCCGTGGTCTGTATTGATATGGCGATGAAAATTACCTTTATGAATCCAGTGGCGGAGAAGATGAGCGGCTGGACGCAGGAAGAAGCGTTAGGTGTTCCGCTCCTGACGGTGTTGCATATTACTTTTGGCGACAACgg                                                                                                                                       <  2:603468/150‑1 (MQ=255)
             ccATCGGTGAAGCCGTGGTCTGTATTGATATGGCGATGAAAATTACCTTTATGAATCCAGTGGCGGAGAAGATGAGCGGCTGGACGCAGGAAGAAGCGTTAGGTGTTCCGCTCCTGACGGTGTTGCATATTACTTTTGGCGACAACGGAc                                                                                                                                     >  2:313235/1‑150 (MQ=255)
              cATCGGTGAAGCCGTGGTCTGTATTGATATGGCGATGAAAATTACCTTTATGAATCCAGTGGCGGAGAAGATGAGCGGCTGGACGCAGGAAGAAGCGTTAGGTGTTCCGCTCCTGACGGTGTTGCATATTACTTTTGGCGACAACGGAcc                                                                                                                                    <  2:1751871/150‑1 (MQ=255)
              cATCGGTGAAGCCGTGGTCTGTATTGATATGGCGATGAAAATTACCTTTATGAATCCAGTGGCGGAGAAGATGAGCGGCTGGACGCAGGAAGAAGCGTTAGGTGTTCCGCTCCTGACGGTGTTGCATATTACTTTTGGCGACAACGGAcc                                                                                                                                    <  1:1721734/150‑1 (MQ=255)
                     gAAGCCGTGGTCTGTATTGATATGGCGATGAAAATTACCTTTATGAATCCAGTGGCGGAGAAGATGAGCGGCTGGACGCAGGAAGAAGCGTTATGTGTTCCGCTCCTGACGGTGTTGCATATTACTTTTGGCGACAACGGACCATTAATg                                                                                                                             <  2:112232/150‑1 (MQ=255)
                      aaGCCGTGGTCTGTATTGATATGGCGATGAAAATTACCTTTATGAATCCAGTGGCGGAGAAGATGAGCGGCTGGACGCAGGAAGAAGCGTTAGGTGTTCCGCTCCTGACGGTGTTGCATATTACTTTTGGCGACAACGGACCATTAATgg                                                                                                                            <  2:753261/150‑1 (MQ=255)
                                cTGTATTGATATGGCGATGAAAATTACCTTTATGAATCCAGTGGCGGAGAAGATGAGCGGCTGGACGCAGGAAGAAGCGTTAGGTGTTCCGCTCCTGACGGTGTTGCATATTACTTTTGGCGACAACGGACCATTAATGGAGAACATTTa                                                                                                                  <  2:1580656/150‑1 (MQ=255)
                                 tGTATTGATATGGCGATGAAAATTACCTTTATGAATCCAGTGGCGGAGAAGATGAGCGGCTGGACGCAGGAAGAAGCGTTAGGTGTTCAGCTCCTGACGGTGTTGCATATTACTTTTGGCGACAACGGACCATTAATGGAGAACATTTAc                                                                                                                 >  1:658271/1‑150 (MQ=255)
                                  gTATTGATATGGCGATGAAAATTACCTTTATGAATCCAGTGGCGGAGAAGATGAGCGGCTGGACGCAGGAAGAAGCGTTAGGTGTTCCGCTCCTGACGGTGTTGCATATTACTTTTGGCGACAACGGACCATTAATGGAGAACATTTACa                                                                                                                <  1:461410/150‑1 (MQ=255)
                                    atagATATGGCGATGAAAATTACCTTTATGAATCCAGTGGCGGAGAAGATGAGCGGCTGGACGCAGGAAGAAGCGTTAGGTGTTCCGCTCCTGACGGTGTTGCATATTACTTTTGGCGACAACGGACCATTAATGGAGAACATTTAc                                                                                                                 >  1:1274321/4‑147 (MQ=255)
                                    atagATATGGCGATGAAAATTACCTTTATGAATCCAGTGGCGGAGAAGATGAGCGGCTGGACGCAGGAAGAAGCGTTAGGTGTTCCGCTCCTGACGGTGTTGCATATTACTTTTGGCGACAACGGACCATTAATGGAGAACATTTAc                                                                                                                 <  2:1274321/144‑1 (MQ=255)
                                        atatGGCGATGAAAATTACCTTTATGAATCCAGTGGCGGAGAAGATGAGCGGCTGGACGCAGGAAGAAGCGTTAGGTGTTCCGCTCCTGACGGTGTTGCATATTACTTTTGGCGACAACGGACCATTAATGGAGAACATTTACAGTGCCg                                                                                                          <  2:77992/150‑1 (MQ=255)
                                         tatGGCGATGAAAATTACCTTTATGAATCCAGTGGCGGAGAAGATGAGCGGCTGGACGCAGGAAGAAGCGTTAGGTGTTCCGCTCCTGACGGTGTTGCATATTACTTTTGGCGACAACGGACCATTAATGGAGAACATTTACAGTGCCGa                                                                                                         >  1:1353709/1‑150 (MQ=255)
                                          atGGCGATGAAAATTACCTTTATGAATCCAGTGGCGGAGAAGATGAGCGGCTGGACGCAGGAAGAAGCGTTAGGTGTTCCGCTCCTGACGGTGTTGCATATTACTTTTGGCGACAACGGACCATTAATGGAGAACATTTACAGTGCCGac                                                                                                        <  1:1200907/150‑1 (MQ=255)
                                          atGGCGATGAAAATTACCTTTATGAATCCAGTGGCGGAGAAGATGAGCGGCTGGACGCAGGAAGAAGCGTTAGGTGTTCCGCTCCTGACGGTGTTGCATATTACTTTTGGCGACAACGGACCATTAATGGAGAACATTTACAGTGCCGac                                                                                                        <  2:658271/150‑1 (MQ=255)
                                          atGGCGATGAAAATTACCTTTATGAATCCAGTGGCGGAGAAGATGAGCGGCTGGACGCAGGAAGAAGCGTTAGGTGTTCCGCTCCTGACGGTGTTGCATATTACTTTTGGCGACAACGGACCATTAATGGAGAACATTTACAGTGCCGac                                                                                                        <  2:1190574/150‑1 (MQ=255)
                                          atGGCGATGAAAATTACCTTTATGAATCCAGTGGCGGAGAAGATGAGCGGCTGGACGCAGGAAGAAGCGTTAGGTGTTCCGCTCCTGACGGTGTTGCATATTACTTTTGGCGACAACGGACCATTAATGGAGAACATTTACAGTGCCGac                                                                                                        <  1:1708032/150‑1 (MQ=255)
                                          atGGCGATGAAAATTACCTTTATGAATCCAGTGGCGGAGAAGATGAGCGGCTGGACGCAGGAAGAAGCGTTAGGTGTTCCGCTCCTGACGGTGTTGCATATTACTTTTGGCGACAACGGACCATTAATGGAGAACATTTACAGTGCCGac                                                                                                        <  2:1353709/150‑1 (MQ=255)
                                             gCGATGAAAATTACCTTTATGAATCCAGTGGCGGAGAAGATGAGCGGCTGGACGCAGGAAGAAGCGTTAGGTGTTCCGCTCCTGACGGTGTTGCATATTACTTTTGGCGACAACGGACCATTAATGGAGAACATTTACAGTGCCGACAcc                                                                                                     <  2:1169503/150‑1 (MQ=255)
                                               gATGAAAATTACCTTTATGAATCCAGTGGCGGAGAAGATGAGCGGCTGGACGCAGGAAGAAGCGTTAGGTGTTCCGCTCCTGACGGTGTTGCATATTACTTTTGGCGACAACGGACCATTAATGGAGAACATTTACAGTGCCGACACCTc                                                                                                   >  2:123808/1‑150 (MQ=255)
                                                    aaaTTACCTTTATGAATCCAGTGGCGGAGAAGATGAGCGGCTGGACGCAGGAAGAAGCGTTAGGTGTTCCGCTCCTGACGGTGTTGCATATTACTTTTGGCGACAACGGACCATTAATGGAGAACATTTACAGTGCCGACACCTCACGtt                                                                                              >  1:939259/1‑150 (MQ=255)
                                                         aCCTTTATGAATCCAGTGGCGGAGAAGATGAGCGGCTGGACGCAGGAAGAAGCGTTAGGTGTTCCGCTCCTGACGGTGTTGCATATTACTTTTGGCGACAACGGACCATTAATGGAGAACATTTACAGTGCCGACACCTCACGTTCcgcg                                                                                         <  2:1618184/150‑1 (MQ=255)
                                                          ccTTTATGAATCCAGTGGCGGAGAAGATGAGCGGCTGGACGCAGGAAGAAGCGTTAGGTGTTCCGCTCCTGACGGTGTTGCATATTACTTTTGGCGACAACGGACCATTAATGGAGAACATTTACAGTGCCGACACCTCACGTTCCGCGa                                                                                        >  1:434285/1‑150 (MQ=255)
                                                            tttATGAATCCAGTGGCGGAGAAGATGAGCGGCTGGACGCAGGAAGAAGCGTTAGGTGTTCCGCTCCTGACGGTGTTGCATATTACTTTTGGCGACAACGGACCATTAATGGAGAACATTTACAGTGCCGACACCTCACGTTCCGCGATc                                                                                      >  1:980649/1‑150 (MQ=255)
                                                             ttATGAATCCAGTGGCGGAGAAGATGAGCGGCTGGACGCAGGAAGAAGCGTTAGGTGTTCCGCTCCTGACGGTGTTGCATATTACTTTTGGCGACAACGGACCATTAATGGAGAACATTTACAGTGCCGACACCTCACGTTCCGCGATCg                                                                                     <  2:1217726/150‑1 (MQ=255)
                                                              tATGAATCCAGTGGCGGAGAAGATGAGCGGCTGGACGCAGGAAGAAGCGTTAGGTGTTCCGCTCCTGACGGTGTTGCATATTACTTTTGGCGACAACGGACCATTAATGGAGAACATTTACAGTGCCGACACCTCACGTTCCGCGATCGa                                                                                    >  1:1824748/1‑150 (MQ=255)
                                                                tGAATCCAGTGGCGGAGAAGATGAGCGGCTGGACGCAGGAAGAAGCCTTAGGTGTTCCGCTCCTGACGCTGTTGCATATTACTATTGGCGACAACGGACCATTAATGGAGAACATTTACAGTGACGACACCTAACGTTCCGCGATCGAga                                                                                  >  2:652527/1‑149 (MQ=255)
                                                                 gAATCCAGTGGCGGAGAAGATGAGCGGCTGGACGCAGGAAGAAGCGTTAGGTGTTCCGCTCCTGACGGTGTTGCATATTACTTTTGGCGAGAACGGACCATTAATGGAGAACATTTACAGTGCAGACACCTCACGTTCAGGGATCGAGCa                                                                                 <  2:374251/150‑1 (MQ=255)
                                                                 gAATCCAGTGGCGGAGAAGATGAGCGGCTGGACGCAGGAAGAAGCGTTAGGTGTTCCGCTCCTGACGGTGTTGCATATTACTTTTGGCGACAACGGACCATTAATGGAGAACATTTACAGTGCCGACACCTCACGTTCCGCGATCGAGCa                                                                                 >  1:1940145/1‑150 (MQ=255)
                                                                   aTCCAGTGGCGGAGAAGATGAGCGGCTGGACGCAGGAAGAAGCGTTAGGTGTTCCGCTCCTGACGGTGTTGCATATTACTTTTGGCGACAACGGACCATTAATGGAGAACATTTACAGTGCCGACACCTCACGTTCCGCGATCGAGCAAg                                                                               <  1:123808/150‑1 (MQ=255)
                                                                      cAGTGGCGGAGAAGATGAGCGGCTGGACGCAGGAAGAAGCGTTAGGTGTTCCGCTCCTGACGGTGTTGCATATTACTTTTGGCGACAACGGACCATTAATGGAGAACATTTACAGTGCCGACACCTCACGTTCCGCGATCGAGCAAGAtg                                                                            >  1:642234/1‑150 (MQ=255)
                                                                      cAGTGGCGGAGAAGAGGAGCGGCTGGACGCAGGAAGAAGCGTTAGGTGTTCCGCTCCTGACGGTGTTGCATATTACTTTTGGCGACAACGGACCATTAATGGAGAACATTTACAGTGCCGACACCTCACGTTCCGCGATCGAGCAAGAtg                                                                            <  2:1596339/150‑1 (MQ=255)
                                                                       aGTGGCGGAGAAGATGAGCGGCTGGACGCAGGAAGAAGCGTTAGGTGTTCCGCTCCTGACGGTGTTGCATATTACTTTTGGCGACAACGGACCATTAATGGAGAACATTTACAGTGCCGACACCTCACGTTCCGCGATCGAGCAAGAtgt                                                                           <  2:1543737/150‑1 (MQ=255)
                                                                                 aagaTGAGCGGCTGGACGCAGGAAGAAGCGTTAGGTGTTCCGCTCCTGACGGTGTTGCATATTACTTTTGTCGACAACGGACCATTAATGGAGAACATTTACAGGGCCGACACCTCACGTTCCGCGATCGAGCAATATGTGGTGTTGCAc                                                                 <  1:825858/150‑1 (MQ=255)
                                                                                 aagaTGAGCGGCTGGACGCAGGAAGAAGCGTTAGGTGTTCCGCTCCTGACGGTGTTGCATATTACTTTTGGCGACAACGGACCATTAATGGAGAACATTTACAGTGCCGACACCTCACGTTCCGCGATCGAGCAAGATGTGGTGTTGCAc                                                                 <  2:1789467/150‑1 (MQ=255)
                                                                                   taTGAGCGGCTGGACGCAGGAAGAAGCGTTAGGTGTTCCGCTCCTGACGGTGTTGCATATTACTTTTGGCGACAACGGACCATTAATGGAGAACATTTACAGTGCCGACACCTCACGTTCCGCGATCGAGCAAGATGTGGTGTTGCACTg                                                               <  1:1719491/149‑1 (MQ=255)
                                                                                   gaTGAGCGGCTGGACGCAGGAAGAAGCGTTAGGTGTTCCGCTCCTGACGGTGTTGCATATTACTTTTGGCGACAACGGACCATTAATGGAGAACATTTACAGTGCCGACACCTCACGTTCCGCGATCGAGCAAGATGTGGTGTTGCACTg                                                               <  2:1940145/150‑1 (MQ=255)
                                                                                             tGGACGCAGGAAGAAGCGTTAGGTGTTCCGCTCCTGACGGTGTTGCATATTACTTTTGGCGACAACGGACCATTAATGGAGAACATTTACAGTGCCGACACCTCACGTTCCGCGATCGAGCAAGATGTGGTGTTGCACTGTCGGAgcggc                                                     >  1:1713678/1‑150 (MQ=255)
                                                                                                  gCAGGAAGAAGCGTTAGGTGTTCCGCTCCTGACGGTGTTGCATATTACTTTTGGCGACAACGGACCATTAATGGAGAACATTTACAGTGCCGACACCTCACGTTCCGCGATCGAGCAAGATGTGGTGTTGCACTGTCGGAGCGGCGGCAg                                                <  1:1103469/150‑1 (MQ=255)
                                                                                                   cAGGAAGAAGCGTTAGGTGTTCCGCTCCTGACGGTGTTGCATATTACTTTTGGCGACAACGGTCCATTAATGGAGAACATTTACAGTGCCGACACCTCACGTTCCGCGATCGAGCAAGATGTGGTGTTGCACTGTCGGAGCGGCGGCAGt                                               >  2:125134/1‑150 (MQ=255)
                                                                                                   cAGGAAGAAGCGTTAGGTGTTCCGCTCCTGACGGTGTTGCATATTACTTTTGGCGACAACGGACCATTAATGGAGAACATTTACAGTGCCGACACCTCACGTTCCGCGATCGAGCAAGATGTGGTGTTGCACTGTCGGAGCGGCGGCAGt                                               >  1:1871740/1‑150 (MQ=255)
                                                                                                           agCGTTAGGTGTTCCGCTCCTGACGGTGTTGCATATTACTTTTGGCGACAACGGACCATTAATGGAGAACATTTACAGTGCCGACACCTCACGTTCCGCGATCGAGCAAGATGTGGTGTTGCACTGTCGGAGCGGCGGCAGTTACGACGt                                       >  2:1458544/1‑150 (MQ=255)
                                                                                                                    tgtTCCGCTCCTGACGGTGTTGCATATTACTTTTGGCGACAACGGACCATTAATGGAGAACATTTACAGTGCCGACACCTCACGTTCCGCGATCGAGCAAGATGTGGTGTTGCACTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAg                              >  2:1494177/1‑150 (MQ=255)
                                                                                                                    tgtTCCGCTCCTGACGGTGTTGCATATTACTTTTGGCGACAACGGACCATTAATGGAGAACATTTACAGTGCCGACACCTCACGTTCCGCGATCGAGCAAGATGTGGTGTTGCACTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAg                              <  1:1073785/150‑1 (MQ=255)
                                                                                                                          gCTCCTGACGGTGTTGCATATTACTTTTGGCGACAACGGACCATTAATGGAGAACATTTACAGTGCCGACACCTCACGTTCCGCGATCGAGCAAGATGTGGTGTTGCACTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTAc                        >  2:1019567/1‑150 (MQ=255)
                                                                                                                            tCCTGACGGTGTTGCATATTACTTTTGGCGACAACGGACCATTAATGGAGAACATTTACAGTGCCGACACCTCACGTTCCGCGATCGAGCAAGATGTGGTGTTGCACTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGc                      <  2:939259/150‑1 (MQ=255)
                                                                                                                              cTGACGGTGTTGCATATTACTTTTGGCGACAACGGACCATTAATGGAGAACATTTACAGTGCCGACACCTCACGTTCCGCGATCGAGCAAGATGTGGTGTTGCACTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAg                              <  1:776927/140‑1 (MQ=255)
                                                                                                                              cTGACGGTGTTGCATATTACTTTTGGCGACAACGGACCATTAATGGAGAACATTTACAGTGCCGACACCTCACGTTCCGCGATCGAGCAAGATGTGGTGTTGCACTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAg                              >  2:776927/1‑140 (MQ=255)
                                                                                                                               tGACGGTGTTGCATATTACCTTTGGCGACAACGGACCATTAATGGAGAACATTTACAGTGCCGACACCTCACGTTCCGCGATCGAGCAAGATGTGGTGTTGCACTGTCGGAGCGGCGGCAGTTACGACGTGc                                     <  1:383761/132‑1 (MQ=255)
                                                                                                                               tGACGGTGTTGCATATTACCTTTGGCGACAACGGACCATTAATGGAGAACATTTACAGTGCCGACACCTCACGTTCCGCGATCGAGCAAGATGTGGTGTTGCACTGTCGGAGCGGCGGCAGTTACGACGTGc                                     <  1:376677/132‑1 (MQ=255)
                                                                                                                               tGACGGTGTTGCATATTACCTTTGGCGACAACGGACCATTAATGGAGAACATTTACAGTGCCGACACCTCACGTTCCGCGATCGAGCAAGATGTGGTGTTGCACTGTCGGAGCGGCGGCAGTTACGACGTGc                                     >  2:376677/1‑132 (MQ=255)
                                                                                                                               tGACGGTGTTGCATATTACCTTTGGCGACAACGGACCATTAATGGAGAACATTTACAGTGCCGACACCTCACGTTCCGCGATCGAGCAAGATGTGGTGTTGCACTGTCGGAGCGGCGGCAGTTACGACGTGc                                     >  2:383761/1‑132 (MQ=255)
                                                                                                                                  cGGTGTTGCATATTACTTTTGGCGACAACGGACCATTAATGGAGAACATTTACAGTGCCGACACCTCACGTTCCGCGATCGAGCAAGATGTGGTGTTGCACTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTaa                >  2:1797462/1‑150 (MQ=255)
                                                                                                                                     tgtTGCATATTACTTTTGGCGACAACGGACCATTAATGGAGAACATTTACAGTGCCGACACCTCACGTTCCGCGATCGAGCAAGATGTGGTGTTGCACTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGcc                     >  2:1318788/1‑142 (MQ=255)
                                                                                                                                     tgtTGCATATTACTTTTGGCGACAACGGACCATTAATGGAGAACATTTACAGTGCCGACACCTCACGTTCCGCGATCGAGCAAGATGTGGTGTTGCACTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGcc                     <  1:1318788/142‑1 (MQ=255)
                                                                                                                                     tgtTGCATATTACTTTTGGCGACAACGGACCATTAATGGAGAACATTTACAGTGCCGACACCTCACGTTCCGCGATCGAGCAAGATGTGGTGTTGCACTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTa             <  1:1068006/150‑1 (MQ=255)
                                                                                                                                        tGCATATTACTTTTGGCGACAACGGACCATTAATGGAGAACATTTACAGTGCCGACACCTCACGTTCCGCGATCGAGCAAGATGTGGTGTTGCACTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTActc          >  1:1840194/1‑150 (MQ=255)
                                                                                                                                        tGCATATTACTTTTGGCGACAACGGACCATTAATGGAGAACATTTACAGTGCCGACACCTCACGTTCCGCGATCGAGCAAGATGTGGTGTTGCACTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTActc          >  1:1786176/1‑150 (MQ=255)
                                                                                                                                         gCATATTACTTTTGGCGACAACGGACCATTAATGGAGAACATTTACAGTGCCGACACCTCACGTTCCGCGATCGAGCAAGATGTGGTGTTGCACTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTActct         <  1:980121/150‑1 (MQ=255)
                                                                                                                                         gCATATTACTTTTGGCGACAACGGACCATTAATGGAGAACATTTACAGTGCCGACACCTCACGTTCCGCGATCGAGCAAGATGTGGTGTTGCACTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTActct         <  1:313235/150‑1 (MQ=255)
                                                                                                                                         gCATATTACTTTTGGCGACAACGGACCATTAATGGAGAACATTTACAGTGCCGACACCTCACGTTCCGCGATCGAGCAAGATGTGGTGTTGCACTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTActct         >  1:1865125/1‑150 (MQ=255)
                                                                                                                                         gCATATTACTTTTGGCGACAACGGACCATTAATGGAGAACATTTACAGTGCCGACACCTCACGTTCCGCGATCGAGCAAGATGTGGTGTTGCACTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTActct         >  2:980121/1‑150 (MQ=255)
                                                                                                                                           atatTACTTTTGGCGACAACGGACCATTAATGGAGAACATTTACAGTGCCGACACCTCACGTTCCGCGATCGAGCAAGATGTGGTGTTGCACTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTgg       >  1:1022409/1‑150 (MQ=255)
                                                                                                                                               tACTTTTGGCGACAACGGACCATTAATGGAGAACATTTACAGTGCCGACACCTCACGTTCCGCGATCGAGCAAGATGTGGTGTTGCACTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACgg   >  2:570329/1‑150 (MQ=255)
                                                                                                                                                aCTTTTGGCGACAACGGACCATTAATGGAGAACATTTACAGTGCCGACACCTCACGTTCCGCGATCGAGCAAGATGTGGTGTTGCACTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGgc  >  1:1715202/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
ATTACGCTTGATTCCATCGGTGAAGCCGTGGTCTGTATTGATATGGCGATGAAAATTACCTTTATGAATCCAGTGGCGGAGAAGATGAGCGGCTGGACGCAGGAAGAAGCGTTAGGTGTTCCGCTCCTGACGGTGTTGCATATTACTTTTGGCGACAACGGACCATTAATGGAGAACATTTACAGTGCCGACACCTCACGTTCCGCGATCGAGCAAGATGTGGTGTTGCACTGTCGGAGCGGCGGCAGTTACGACGTGCATTACAGTATTACGCCGTTAAGTACTCTGGACGGC  >  NC_000913/2144940‑2145233

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: