Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,028,159 T→A 16.9% intergenic (‑142/+29) yedQ ← / ← yodC putative membrane‑anchored diguanylate cyclase/uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,028,1590TA16.9% 352.8 / 57.5 141intergenic (‑142/+29)yedQ/yodCputative membrane‑anchored diguanylate cyclase/uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base T (47/70);  new base A (12/12);  total (59/82)
Fisher's exact test for biased strand distribution p-value = 4.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.76e-01

Only 100 of 151 total aligned reads displayed.

CTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTG  >  NC_000913/2028060‑2028257
                                                                                                   |                                                                                                  
cTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGAtgt                                                                                                    <  2:1402893/100‑1 (MQ=255)
cTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGAtgt                                                                                                    <  1:4062508/100‑1 (MQ=255)
 tGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTAATTGCAGGCATTAGCACACCAGGCCGGAtgtg                                                                                                   >  1:3193565/1‑100 (MQ=255)
   gCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGAt                                                                                                 >  2:454210/1‑100 (MQ=255)
   gCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGAt                                                                                                 <  1:2249429/100‑1 (MQ=255)
        aGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTA‑‑ACATCAGGCCGGATCACATCCGGCCt                                                                                          <  1:2919938/100‑1 (MQ=37)
         gCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGcc                                                                                           <  2:1103235/100‑1 (MQ=255)
         gCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGcc                                                                                           >  1:850305/1‑100 (MQ=255)
         gCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGcc                                                                                           <  2:469494/100‑1 (MQ=255)
         gCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGcc                                                                                           >  2:3400597/1‑100 (MQ=255)
          cGTACCGCAAGACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACATCAGGCCGGATCACATCCGGCCt                                                                                          <  1:2888271/100‑1 (MQ=38)
          cGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACATCAGGCCGGATCACATCCGGCCt                                                                                          <  2:3589106/100‑1 (MQ=39)
          cGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACATCAGGCCGGATCACATCCGGCCt                                                                                          <  2:923847/100‑1 (MQ=38)
          cGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGcc                                                                                           <  2:2056270/99‑1 (MQ=255)
          cGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATCACATCCGGCCt                                                                                          <  1:32525/100‑1 (MQ=255)
           gTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTg                                                                                         <  2:2745086/100‑1 (MQ=255)
           gTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATCACATCCGGCCTg                                                                                         <  2:2907171/100‑1 (MQ=255)
            tACCGCAAAACAGCGATAATTTTCTAATGTTGGCCTGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGcc                                                                                           <  1:3969567/97‑1 (MQ=255)
            tACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGa                                                                                        >  2:4294212/1‑100 (MQ=255)
            tACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATCACATCCGGCCTgg                                                                                        <  1:3745537/100‑2 (MQ=255)
            tACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATCACATCCGGCCTgg                                                                                        <  1:1815667/100‑2 (MQ=255)
             aCCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGAt                                                                                       <  2:2069878/100‑1 (MQ=255)
             aCCGCAAAACAGCGATAATTTTAGAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACATCAGGCCGGATCACATCCGGCCt                                                                                          <  1:632582/97‑1 (MQ=38)
              ccGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATg                                                                                      >  1:2919886/1‑100 (MQ=255)
              ccGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATg                                                                                      <  2:3171410/100‑1 (MQ=255)
              ccGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATg                                                                                      <  2:70306/100‑1 (MQ=255)
               cGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGt                                                                                     <  2:224409/100‑1 (MQ=255)
               cGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGt                                                                                     >  2:2140215/1‑100 (MQ=255)
               cGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATCACATCCGGCCTGgtgt                                                                                     <  2:1429841/100‑1 (MQ=38)
                gCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACATCAGGCCGGATCACATCCGGCCTGgtgtg                                                                                    <  2:1982297/100‑2 (MQ=37)
                gCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGtt                                                                                    >  2:340426/1‑100 (MQ=255)
                 cAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGACCCGGCCTGATGTTa                                                                                   <  2:1791397/100‑1 (MQ=255)
                  aaaaCAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTaa                                                                                  <  1:1805997/100‑1 (MQ=255)
                  aaaaCAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTaa                                                                                  >  1:569435/1‑100 (MQ=255)
                   aaaCAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAAt                                                                                 <  2:1117267/100‑1 (MQ=255)
                    aaCAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATg                                                                                <  2:4298545/100‑1 (MQ=255)
                    aaCAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATg                                                                                <  2:798670/100‑1 (MQ=255)
                    aaCAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATg                                                                                <  1:253995/100‑1 (MQ=255)
                    aaCAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATg                                                                                <  2:3018558/100‑1 (MQ=255)
                    aaCAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATg                                                                                <  1:2026787/100‑1 (MQ=255)
                        gCGATAATTTTATAATGGTTGCATGATAAGACGGGACTATTATTCATTGCAGGCATTAGCACACCAGGCCTGATGTGATCCGGCCTGATGTTAATGAt                                                                              <  1:3695044/98‑1 (MQ=255)
                          gATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATAccc                                                                          >  1:3826794/1‑100 (MQ=255)
                          gATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATAccc                                                                          <  1:2479988/100‑1 (MQ=255)
                            tAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCtt                                                                        <  2:1959528/100‑1 (MQ=255)
                               ttttATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCTGGCATTAGGACACCAGGCAGGATGAGATCCGGCATGATGTTAATGATACCCTTTca                                                                     >  1:2972817/1‑100 (MQ=255)
                               ttttATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTca                                                                     >  2:2828282/1‑100 (MQ=255)
                                tttATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTcac                                                                    <  2:237171/100‑1 (MQ=255)
                                  tataATGGTGGCATGATAAGACGCGACTATCATTCACTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCacac                                                                  <  2:3097810/100‑1 (MQ=255)
                                   ataATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCGTTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACAct                                                                 <  1:2828265/100‑1 (MQ=255)
                                          tggCGTGATAAGACCCGACTATCATGCATTGCAGGCATTAGCACACCAGGCCGTATGTGATCCAGCCTGATGTTAATGATACCCTTTCACACTTCTTCCg                                                          <  1:2124193/100‑1 (MQ=255)
                                          tggCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTAAATGATACCCTTTCACACTTCTTCCg                                                          <  1:572111/100‑1 (MQ=255)
                                           ggCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGc                                                         <  2:3252674/100‑1 (MQ=255)
                                            gCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCa                                                        <  1:4241634/100‑1 (MQ=255)
                                               tGATAAGACGCGACTATCATTCATGGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGaa                                                     <  1:4115649/100‑1 (MQ=255)
                                                gATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAAc                                                    <  2:4403499/100‑1 (MQ=255)
                                                gATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAAc                                                    <  2:2499764/100‑1 (MQ=255)
                                                 aTAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACg                                                   >  1:2552417/1‑100 (MQ=255)
                                                     gACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGActc                                               >  2:1877365/1‑100 (MQ=255)
                                                      aCGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGActcc                                              <  1:3788056/100‑1 (MQ=255)
                                                      aCGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGActcc                                              <  1:2533993/100‑1 (MQ=255)
                                                      aCGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGActcc                                              <  1:1877358/100‑1 (MQ=255)
                                                           aCTAGCATTCATTGCAGGCATTAGCACGCCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGActcccctcc                                         <  1:895007/100‑1 (MQ=255)
                                                                cattcattGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTccccc                                      >  1:1927463/1‑98 (MQ=255)
                                                                 attcattGCAGGCATTAGCACACCAGGACGGATGTGATGCGGCCTGATGTTAATGATACCATTTCACAATTCTTCCGCACCACGACTCCCCTCCCCCg                                     >  1:4065372/1‑98 (MQ=255)
                                                                  ttcattGCAGGCATTAGCACACCAGGCAGGATGTGATCCGGCCTCATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGaa                                  <  1:646080/100‑1 (MQ=255)
                                                                   tcattGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAAc                                 <  1:3569225/100‑1 (MQ=255)
                                                                    cattGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACa                                >  1:2328208/1‑100 (MQ=255)
                                                                     attGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACaa                               >  2:582795/1‑100 (MQ=255)
                                                                     attGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACaa                               <  2:4408384/100‑1 (MQ=255)
                                                                       tGCGGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGc                             >  2:1147997/1‑100 (MQ=255)
                                                                        gCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCt                            <  1:4317083/100‑1 (MQ=255)
                                                                         cAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGcac                           >  2:3603374/1‑98 (MQ=255)
                                                                          aGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCg                          <  1:1729995/100‑1 (MQ=255)
                                                                             cATTAGCACACCATGCCGGATGTGATCCGGCCTGATGTAAATGATACCCGTCCACACTTCTTCCGCAGAACGACTCCCATCCCCCGGAACAAGCTCgg                         <  1:1180236/98‑1 (MQ=255)
                                                                              tcattaacatcAGGCCGGATCACATCCGGCCTGGTGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGttt                      >  2:1609825/11‑100 (MQ=25)
                                                                                tAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACTCTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCa                    <  2:2226848/100‑1 (MQ=255)
                                                                                tAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCa                    >  2:576908/1‑100 (MQ=255)
                                                                                tAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCa                    <  2:4043564/100‑1 (MQ=255)
                                                                                attaacatcAGGCCGGATCACATCCGGCCTGGTGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCa                    >  1:201597/9‑100 (MQ=25)
                                                                                 ttaacatcAGGCCGGATCACATCCGGCCTGGTGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCAt                   >  2:657241/8‑100 (MQ=25)
                                                                                 aGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCAt                   >  2:254639/1‑100 (MQ=255)
                                                                                  gCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTc                     >  1:2940273/1‑97 (MQ=255)
                                                                                  gCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATg                  >  2:3775018/1‑100 (MQ=255)
                                                                                  gCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATg                  >  1:2510603/1‑100 (MQ=255)
                                                                                  gCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATg                  <  1:784684/100‑1 (MQ=255)
                                                                                    acacCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGaa                <  1:1570837/100‑1 (MQ=255)
                                                                                     catcAGGCCGGATCACATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGaa                >  2:1191285/4‑99 (MQ=38)
                                                                                      acCAGGCCGTATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGaaaa              <  1:4297539/100‑1 (MQ=255)
                                                                                      acCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGaaaa              <  1:3430833/100‑1 (MQ=255)
                                                                                        cacgCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGc            <  1:2849247/97‑1 (MQ=255)
                                                                                        cAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGc            >  2:2850132/1‑100 (MQ=255)
                                                                                         aGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCt           >  2:2914889/1‑100 (MQ=255)
                                                                                          ggCCGGATCACATCCGGCCTGGTGTTAATGATACCCTTTCACCCTTCTTCCGCAGAACGACTCCCCTCACCCGGAACAAGCTCGGTTTCATGAAAAGCtt          >  1:501782/1‑100 (MQ=255)
                                                                                           gCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTc         <  1:216906/100‑1 (MQ=255)
                                                                                            ccGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTcc        >  1:1196221/1‑100 (MQ=255)
                                                                                             cGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTccc       <  1:3813859/100‑1 (MQ=255)
                                                                                               gATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGc     >  2:2500947/1‑100 (MQ=255)
                                                                                                  gtgATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTg  <  1:523355/100‑1 (MQ=255)
                                                                                                  gtgATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTg  <  2:792040/100‑1 (MQ=255)
                                                                                                  gtgATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTg  <  2:1608098/100‑1 (MQ=255)
                                                                                                   |                                                                                                  
CTGGCGATAGCGTACCGCAAAACAGCGATAATTTTATAATGGTGGCATGATAAGACGCGACTATCATTCATTGCAGGCATTAGCACACCAGGCCGGATGTGATCCGGCCTGATGTTAATGATACCCTTTCACACTTCTTCCGCAGAACGACTCCCCTCCCCCGGAACAAGCTCGGTTTCATGAAAAGCTTCCCGCTTG  >  NC_000913/2028060‑2028257

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: