Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,518,687 A→C E74A (GAA→GCA)  yfeC → DUF1323 family putative DNA‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,518,6870AC98.1% 737.8 / ‑5.3 214E74A (GAA→GCA) yfeCDUF1323 family putative DNA‑binding protein
Reads supporting (aligned to +/- strand):  ref base A (4/0);  new base C (101/109);  total (105/109)
Fisher's exact test for biased strand distribution p-value = 5.63e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 214 total aligned reads displayed.

GTTCACGTCAATGAACAGGTTCGTGAATATATTCGCAATGCTGAACGTCCAGAAGGTCAGGGAGAAGCGCCTGCCCTTTCTGGTGATGCGCCGCTTGAAGTATTGCTGGTGACACTGGCGAAAGAGATGACGCCAGTTGAACAAAAACAGTTTACATCCTTGCTTCTGCGGGAAGGGATTATCGGATTGTTA  >  NC_000913/2518590‑2518781
                                                                                                 |                                                                                              
gTTCACGTCAATGAACAGGTTCGTGAATATATTCGCAATGCTGAACGTCCAGAAGGTCAGGGAGAAGCGCCTGCCCTTTCTGGTGATGCGCCGCTTGCAg                                                                                              <  1:5004377/100‑1 (MQ=255)
gTTCACGTCAATGAACAGGTTCGTGAATATATTCGCAATGCTGAACGTCCAGAAGGTCAGGGAGAAGCGCCTGCCCTTTCTGGTGATGCGCCGCTTGCAg                                                                                              <  1:273155/100‑1 (MQ=255)
gTTCACGACAATGAACAGGTTCGTGAATATATTCGCAATGCTGAACGTCCAGAAGGTCAGGGAGAAGCGCCTGCCCTTTCTGGTGATGCGCCGCTTGCAg                                                                                              <  2:2043620/100‑1 (MQ=255)
 ttCACGTCAATGAACAGGTTCGTGAATATATTCGCAATGCTGAACGTCCAGAAGGTCAGGGAGAAGCGCCTGCCCTTTCTGGTGATGCGCCGCTTGCAGt                                                                                             <  1:3591366/100‑1 (MQ=255)
 ttCACGTCAATGAACAGGTTCGTGAATATATTCGCAATGCTGAACGTCCAGAAGGTCAGGGAGAAGCGCCTGCCCTTTCTGGTGATGCGCCGCTTGCAGt                                                                                             >  2:4107097/1‑100 (MQ=255)
    aCGTCAATGAACAGGTTCGTGAATATATTCGCAATGCTGAACGTCCAGAAGGTCAGGGAGAAGCGCCTGCCCTTTCTGGTGATGCGCCGCTTGCAGTAtt                                                                                          >  1:5141885/1‑100 (MQ=255)
     cGTCAATGAACAGGTTCGTGAATATATTCGCAATGCTGAACGTCCAGAAGGTCAGGGAGAAGCGCCTGCCCTTTCTGGTGATGCGCCGCTTGCAGTATTg                                                                                         >  1:3023885/1‑100 (MQ=255)
     cGTCAATGAACAGGTTCGTGAATATATTCGCAATGCTGAACGTCCAGAAGGTCAGGGAGAAGCGCCTGCCCTTTCTGGTGATGCGCCGCTTGCAGTATTg                                                                                         <  2:3102966/100‑1 (MQ=255)
     cGTCAATGAACAGGTTCGTGAATATATTCGCAATGCTGAACGTCCAGAAGGTCAGGGAGAAGCGCCTGCCCTTTCTGGTGATGCGCCGCTTGCAGTATTg                                                                                         >  1:4921923/1‑100 (MQ=255)
      gTCAATGAACAGGTTCGTGAATATATTCGCAATGCTGAACGTCCAGGAGGTCAGGGAGAAGCGCCTGCCCTTTCTGGTGATGCGCAGCTTGAAGTGTTGc                                                                                        >  2:2753620/1‑100 (MQ=255)
      gTCAATGAACAGGTTCGTGAATATATTCGCAATGCTGAACGTCCAGAAGGTCAGGGAGAAGCGCCTGCCCTTTCTGGTGATGCGCCGCTTGCAGTATTGc                                                                                        <  1:2880785/100‑1 (MQ=255)
       tCAATGAACAGGTTCGTGAATATATTCGCAATGCTGAACGTCCAGAAGGTCAGGGAGAAGCGCCTGCCCTTTCTGGTGATGCGCCGCTTGCAGTATTGCt                                                                                       >  1:3202220/1‑100 (MQ=255)
       tCAATGAACAGGTTCGTGAATATATTCGCAATGCTGAACGTCCAGAAGGTCAGGGAGAAGCGCCTGCCCTTTCTGGTGATGCGCCGCTTGCAGTATTGCt                                                                                       >  2:3705250/1‑100 (MQ=255)
       tCAATGAACAGGTTCGTGAATATATTCGCAATGCTGAACGTCCAGAAGGTCAGGGAGAAGCGCCTGCCCTTTCTGGTGATGCGCCGCTTGCAGTATTGCt                                                                                       >  2:1376163/1‑100 (MQ=255)
        cAATGAACAGGTTCGTGAATATATTCGCAATGCTGAACGTCCAGAAGGTCAGGGAGAAGCGCCTGCCCTTTCTGGTGATGCGCCGCTTGCAGTATTGCTg                                                                                      <  2:2942942/100‑1 (MQ=255)
           tGAACAGGTTCGTGAATATATTCGCAATGCTGAACGTCCAGAAGGTCAGGGAGAAGCGCCTGCCCTTTCTGGTGATGCGCCGCTTGCAGTATTGCTGGTg                                                                                   <  2:4886266/100‑1 (MQ=255)
           tGAACAGGTTCGTGAATATATTCGCAATGCTGAACGTCCAGAAGGTCAGGGAGAAGCGCCTGCCCTTTCTGGTGATGCGCCGCTTGCAGTATTGCTGGTg                                                                                   <  1:3747418/100‑1 (MQ=255)
             aaCAGGTTCGTGAATATATTCGCAATGCTGAACGTCCAGAAGGTCAGGGAGAAGCGCCTGCCCTTTCTGGTGATGCGCCGCTTGCAGTATTGCTGGTGac                                                                                 >  1:3676178/1‑100 (MQ=255)
             aaCAGGTTCGTGAATATATTCGCAATGCTGAACGTCCAGAAGGTCAGGGAGAAGCGCCTGCCCTTTCTGGTGATGCGCCGCTTGCAGTATTGCTGGTGac                                                                                 <  1:2860984/100‑1 (MQ=255)
             aaCAGGTTCGTGAATATATTCGCAATGCTGAACGTCCAGAAGGTCAGGGAGAAGCGCCTGCCCTTTCTGGTGATGCGCCGCTTGCAGTATTGCTGGTGac                                                                                 <  1:1671703/100‑1 (MQ=255)
               cAGGTTCGTGAATATATTCGCAATGCTGAACGTCCAGAAGGTCAGGGAGAAGCGCCTGCCCTTTCTGGTGATGCGCCGCTTGCAGTATTGCTGGTGacac                                                                               <  2:5197844/100‑1 (MQ=255)
                  gTTCGTGAATATATTCGCAATGCTGAACGTCCAGAAGGTCAGGGAGAAGCGCCTGCCCTTTCTGGTGATGCGCCGCTTGCAGTATTGCTGGTGACACTgg                                                                            <  2:1886287/100‑1 (MQ=255)
                  gTTCGTGAATATATTCGCAATGCTGAACGTCCAGAAGGTCAGGGAGAAGCGCCTGCCCTTTCTGGTGATGCGCCGCTTGCAGTATTGCTGGTGACACTgg                                                                            <  2:629816/100‑1 (MQ=255)
                   ttCGTGAATATATTCGCAATGCTGAACGTCCAGAAGGTCAGGGAGAAGCGCCTGCCCTTTCTGGTGATGCGCCGCTTGCAGTATTGCTGGTGACACTGGc                                                                           <  2:1063659/100‑1 (MQ=255)
                   ttCGTGAATATATTCGCAATGCTGAACGTCCAGAAGGTCAGGGAGAAGCGCCTGCCCTTTCTGGTGATGCGCCGCTTGCAGTATTGCTGGTGACACTGGc                                                                           <  1:5088966/100‑1 (MQ=255)
                     cGTGAATATATTCGCAATGCTGAACGTCCAGAAGGTCAGGGAGAAGCGCCTGCCCTTTCTGGTGATGCGCCGCTTGCAGTATTGCTGGTGACACTGGCGa                                                                         >  2:2614313/1‑100 (MQ=255)
                     cGTGAATATATTCGCAATGCTGAACGTCCAGAAGGTCAGGGAGAAGCGCCTGCCCTTTCTGGTGATGCGCCGCTTGCAGTATTGCTGGTGACACTGGCGa                                                                         >  2:2603123/1‑100 (MQ=255)
                     cGTGAATATATTCGCAATGCTGAACGTACAGAAGGTCAGGGAGAAGCGACTGCCCTTTATGGTGATGAGCCGCTTGCAGTATTGCTGGTGACACTGGCGa                                                                         >  2:576615/1‑100 (MQ=255)
                       tGAATATATTCGCAATGCTGAACGTCCAGAAGGTCAGGGAGAAGCGCCTGCCCTTTCTGGTGATGCGCCGCTTGCAGTATTGCTGGTGACACTGGCGaaa                                                                       <  1:783613/100‑1 (MQ=255)
                       tGAATATATTCGCAATGCTGAACGTCCAGAAGGTCAGGGAGAAGCGCCTGCCCTTTCTGGTGATGCGCCGCTTGCAGTATTGCTGGTGACACTGGCGaaa                                                                       >  2:4863503/1‑100 (MQ=255)
                       tGAATATATTCGCAATGCTGAACGTCCAGAAGGTCAGGGAGAAGCGCCTGCCCTTTCTGGTGATGCGCCGCTTGCAGTATTGCTGGTGACACTGGCGaaa                                                                       >  1:186157/1‑100 (MQ=255)
                        gAATATATTCGCAATGCTGAACGTCCAGAAGGTCAGGGAGAAGCGCCTGCCCTTTCTGGTGATGCGCCGCTTGCAGTATTGCTGGTGACACTGGCGAAag                                                                      <  2:983751/100‑1 (MQ=255)
                        gAATATATTCGCAATGCTGAACGTCCAGAAGGTCAGGGAGAAGCGCCAGCCCTTGCAGGTGATGCGCCGCTTGAAGTGTGGCTGGTGGCACTGGCGGAag                                                                      >  2:1241431/1‑100 (MQ=255)
                          atatatTCGCAATGCTGAACGTCCAGAAGGTCAGGGAGAAGCGCCTGCCCTTTCTGGTGATGCGCCGCTTGCAGTATTGCTGGTGACACTGGCGAAagag                                                                    <  2:943835/100‑1 (MQ=255)
                           tatatTCGCAATGCTGAACGTCCAGAAGGTCAGGGAGACGCGCCTGCCCTTTCTGGTGATGCGCCGCTTGCAGTATTGCTGGTGACACTGGCGAAAgaga                                                                   <  1:5323283/100‑1 (MQ=255)
                           tatatTCGCAATGCTGAACGTCCAGAAGGTCAGGGAGAAGCGCCTGCCCTTTCTGGTGATGCGCCGCTTGCAGTATTGCTGGTGACACTGGCGAAAgaga                                                                   <  1:4685146/100‑1 (MQ=255)
                            atatTCGCAATGCTGAACGTCCAGAAGGTCAGGGAGAAGCGCCTGCCCTTTCTGGTGATGCGCCGCTTGCAGTATTGCTGGTGACACTGGCGAAAGAGAt                                                                  <  1:3820053/100‑1 (MQ=255)
                             tatTCGCAATGCTGAACGTCCAGAAGGTCAGGGAGAAGCGCCTGCCCTTTCTGGTGATGCGCCGCTTGCAGTATTGCTGGTGACACTGGCGAAAGAGATg                                                                 <  2:2125272/100‑1 (MQ=255)
                             tatTCGCAATGCTGAACGTCCAGAAGGTCAGGGAGAAGCGCCTGCCCTTTCTGGTGATGCGCCGCTTGCAGTATTGCTGGTGACACTGGCGAAAGAGATg                                                                 <  1:2164092/100‑1 (MQ=255)
                             tatTCGCAATGCTGAACGTCCAGAAGGTCAGGGAGAAGCGCCTGCCCTTTCTGGTGATGCGCCGCTTGCAGTATTGCTGGTGACACTGGCGAAAGAGATg                                                                 >  1:2087632/1‑100 (MQ=255)
                              atTCGCAATGCTGAACGTCCAGAAGGTCAGGGAGAAGCGCCTGCCCTTTCTGGTGATGCGCCGCTTGCAGTATTGCTGGTGACACTGGCGAAAGAGATGa                                                                >  1:1641282/1‑100 (MQ=255)
                               ttCGCAATGCTGAACGTCCAGAAGGTCAGGGAGAAGCGCCTGCCCTTTCTGGTGATGCGCCGCTTGCAGTATTGCTGGTGACACTGGCGAAAGAGATGAc                                                               >  2:3550403/1‑100 (MQ=255)
                               ttCGCAATGCTGAACGTCCAGAAGGTCAGGGAGAAGCGCCTGCCCTTTCTGGTGATGCGCCGCTTGCAGTATTGCTGGTGACACTGGCGAAAGAGATGAc                                                               <  2:3843391/100‑1 (MQ=255)
                               ttCGCAATGCTGAACGTCCAGAAGGTCAGGGAGAAGCGCCTGCCCTTTCTGGTGATGCGCCGCTTGCAGTATTGCTGGTGACACTGGCGAAAGAGATGAc                                                               <  1:2516871/100‑1 (MQ=255)
                                 cGCAATGCTGAACGTCCAGAAGGTCAGGGAGAAGCGCCTGCCCTTTCTGGTGATGCGCCGCTTGCAGTATTGCTGGTGACACTGGCGAAAGAGATGACGc                                                             >  1:71928/1‑100 (MQ=255)
                                  gCAATGCTGAACGTCCAGAAGGTCAGGGAGAAGCGCCTGCCCTTTCTGGTGATGCGCCGCTTGCAGTATTGCTGGTGACACTGGCGAAAGAGATGACGcc                                                            >  2:3248978/1‑100 (MQ=255)
                                  gCAATGCTGAACGTCCAGAAGGTCAGGGAGAAGCGCCTGCCCTTTCTGGTGATGCGCCGCTTGCAGTATTGCTGGTGACACTGGCGAAAGAGATGACGcc                                                            <  1:5464196/100‑1 (MQ=255)
                                    aaTGCTGAACGTCCAGAAGGTCAGGGAGAAGCGCCTGCCCTTTCTGGTGATGCGCCGCTTGCAGTATTGCTGGTGACACTGGCGAAAGAGATGACGCCAg                                                          <  1:2546990/100‑1 (MQ=255)
                                    aaTGCTGAACGTCCAGAAGGTCAGGGAGAAGCGCCTGCCCTTTCTGGTGATGCGCCGCTTGCAGTATTGCTGGTGACACTGGCGAAAGAGATGACGCCAg                                                          <  1:3394203/100‑1 (MQ=255)
                                     aTGCTGAACGTCCAGAAGGTCAGGGAGAAGCGCCTGCCCTTTCTGGTGATGCGCCGCTTGCAGTATTGCTGGTGACACTGGCGAAAGAGATGACGCCAGt                                                         >  2:2276666/1‑100 (MQ=255)
                                     aTGCTGAACGTCCAGAAGGTCAGGGAGAAGCGCCTGCCCTTTCTGGTGATGCGCCGCTTGCAGTATTGCTGGTGACACTGGCGAAAGAGATGACGCCAGt                                                         <  2:4519259/100‑1 (MQ=255)
                                       gCTGAACGTCCAGAAGGTCAGGGAGAAGCGCCTGCCCTTTCTGGTGATGCGCCGCTTGCAGTATTGCTGGTGACACTGGCGAAAGAGATGACGCCAGTTg                                                       >  1:1740103/1‑100 (MQ=255)
                                        cTGAACGTCCAGAAGGTCAGGGAGAAGCGCCTGCCCTTTCTGGTGATGCGCCGCTTGCAGTATTGCTGGTGACACTGGCGAAAGAGATGACGCCAGTTGa                                                      >  2:1702987/1‑100 (MQ=255)
                                         tGAACGTCCAGAAGGTCAGGGAGAAGCGCCTGCCCTTTCTGGTGATGCGCCGCTTGCAGTATTGCTGGTGACACTGGCGAAAGAGATGACGCCAGTTGaa                                                     >  1:4240334/1‑100 (MQ=255)
                                         tGAACGTCCAGAAGGTCAGGGAGAAGCGCCTGCCCTTTCTGGTGATGCGCCGCTTGCAGTATTGCTGGTGACACTGGCGAAAGAGATGACGCCAGTTGaa                                                     <  1:3998145/100‑1 (MQ=255)
                                               tCCAGAAGGTCAGGGAGAAGCGCCTGCCCTTCCTGGTGATGCGCCGCTTGCAGTATTGCTGGTGACACTGGCGAAAGAGATGACGCCAGTTGAACaaaaa                                               >  1:4222229/1‑100 (MQ=255)
                                                    aaGGTCAGGGAGAAGCGCCTGCCCTTTCTGGTGATGCGCCGCTTGCAGTATTGCTGGTGACACTGGCGAAAGAGATGACGCCAGTTGAACAAAAACAGtt                                          <  2:4473967/100‑1 (MQ=255)
                                                    aaGGTCAGGGAGAAGCGCCTGCCCTTTCTGGTGATGCGCCGCTTGCAGTATTGCTGGTGACACTGGCGAAAGAGATGACGCCAGTTGAACAAAAACAGtt                                          <  2:2162723/100‑1 (MQ=255)
                                                     aGGTCAGGGAGAAGCGCCTGCCCTTTCTGGTGATGCGCCGCTTGCAGTATTGCTGGTGACACTGGCGAAAGAGATGACGCCAGTTGACCAAAAACAGttt                                         <  1:4797288/100‑1 (MQ=255)
                                                     aGGTCAGGGAGAAGCGCCTGCCCTTTCTGGTGATGCGCCGCTTGCAGTATTGCTGGTGACACTGGCGAAAGAGATGACGCCAGTTGAACAAAAACAGttt                                         <  2:4698757/100‑1 (MQ=255)
                                                       gTCAGGGAGAAGCGCCTGCCCTTTCTGGTGATGCGCCGCTTGCAGTATTGCTGGTGACACTGGCGAAAGAGATGACGCCAGTTGAACAAAAACAGTTTAc                                       >  1:3571817/1‑100 (MQ=255)
                                                        tCAGGGAGAAGCGCCTGCCCTTTCTGGTGATGCGCTGCTTGCAGTATTGCTGGTGACACTGGCGAAAGAGATGACGCCAGTTGAACAAAAACAGTTTACa                                      <  2:5416126/100‑1 (MQ=255)
                                                          aGGGAGAAGCGCCTGCCCTTTCTGGTGATGCGCCGCTTGCAGTATTGCTGGTGACACTGGCGAAAGAGATGACGCCAGTTGAACAAAAACAGTTTACATc                                    <  2:3349511/100‑1 (MQ=255)
                                                           gggAGAAGCGCCTGCCCTTTCTGGTGATGCGCCGCTTGCAGTATTGCTGGTGACACTGGCGAAAGAGATGACGCCAGTTGAACAAAAACAGTTTACagat                                   >  1:4670698/1‑97 (MQ=255)
                                                           gggAGAAGCGCCTGCCCTTTCTGGTGATGCGCCGCTTGCAGTATTGCTGGTGACACTGGCGAAAGAGATGACGCCAGTTGAACAAAAACAGTTTACATcc                                   <  1:443921/100‑1 (MQ=255)
                                                           gggAGAAGCGCCTGCCCTTTCTGGTGATGCGCCGCTTGCAGTATTGCTGGTGACACTGGCGAAAGAGATGACGCCAGTTGAACAAAAACAGTTTACATcc                                   >  1:2389486/1‑100 (MQ=255)
                                                            ggAGAAGCGCCTGCCCTTTCTGGCGATGCGCCGGTTGCAGTATTGCTGGTGACACTGGCGAAAGAGATGACGCCAGTTGAACAAAATCAGTTTACATCCt                                  <  2:3312650/100‑1 (MQ=255)
                                                             gagaAGCGCCTGCCCTTTCTGGTGATGCGCCGCTTGCAGTATTGCTGGTGACACTGGCGAAAGAGATGACGCCAGTTGAACAAAAACAGTTTACATCCtt                                 <  1:2628746/100‑1 (MQ=255)
                                                             gagaAGCGCCTGCCCTTTCTGGTGATGCGCCGCTTGCAGTATTGCTGGTGACACTGGCGAAAGAGATGACGCCAGTTGAACAAAAACAGTTTACATCCtt                                 >  1:1381611/1‑100 (MQ=255)
                                                             gagaAGCGCCTGCCCTTTCTGGTGATGCGCCGCTTGCAGTATTGCTGGTGACACTGGCGAAAGAGATGACGCCAGTTGAACAAAAACAGTTTACATCCtt                                 >  1:4499228/1‑100 (MQ=255)
                                                               gaAGCGCCTGCCCTTTCTGGTGATGCGCCGCTTGCAGTATTGCTGGTGACACTGGCGAAAGAGATGACGCCAGTTGAACAAAAACAGTTTACATCCTTGc                               >  2:590509/1‑100 (MQ=255)
                                                                  gcgcCTGCCCTTTCTGGTGATGCGCCGCTTGCAGTATTGCTGGTGACACTGGCGAAAGAGATGACGCCAGTTGAACAAAAACAGTTTACATCCTTGCTTc                            <  1:3405373/100‑1 (MQ=255)
                                                                  gcgcCTGCCCTTTCTGGTGATGCGCCGCTTGCAGTATTGCTGGTGACACTGGCGAAAGAGATGACGCCAGTTGAACAAAAACAGTTTACATCCTTGCTTc                            <  1:5248244/100‑1 (MQ=255)
                                                                   cgcCTGCCCTTTCTGGTGATGCGCCGCTTGCAGTATTGCTGGTGACACTGGCGAAAGAGATGACGCCAGTTGAACAAAAACAGTTTACATCCTTGCTTCt                           >  1:4265832/1‑100 (MQ=255)
                                                                   cgcCTGCCCTTTCTGGTGATGCGCCGCTTGCAGTATTGCTGGTGACACTGGCGAAAGAGATGACGCCAGTTGAACAAAAACAGTTTACATCCTTGCTTCt                           <  1:1072471/100‑1 (MQ=255)
                                                                     ccggCCCGTTCTGGTGATGCGCCGCCTGCAGTATTGCTGGTGACACTGGCGAAAGAGATGACGAAAAAAAAAAAAAAACAGTTTACATCCTTGCTTCTGc                         <  2:3000622/97‑1 (MQ=255)
                                                                       tGCCCTTTCTGGTGATGCGCCGCTTGCAGTATTGCTGGTGACACTGGCGAAAGAGATGACGCCAGTTGAACAAAAACAGTTTACATCCTTGCTTCTGCgg                       >  1:4303065/1‑100 (MQ=255)
                                                                       tGCCCTTTCTGGTGATGCGCCGCTTGCAGTATTGCTGGTGACACTGGCGAAAGAGATGACGCCAGTTGAACAAAAACAGTTTACATCCTTGCTTCTGCgg                       >  2:4310803/1‑100 (MQ=255)
                                                                           cTTTCTGGTGATGCGCCGCTTGCAGTATTGCTGGTGACACTGGCGAAAGAGATGACGCCAGTTGAACAAAAACAGTTTACATCCTTGCTTCTGCGGGAAg                   <  1:530740/100‑1 (MQ=255)
                                                                           cTTTCTGGTGATGCGCCGCTTGCAGTATTGCTGGTGACACTGGCGAAAGAGATGACGCCAGTTGAACAAAAACAGTTTACATCCTTGCTTCTGCGGGAAg                   >  2:2482833/1‑100 (MQ=255)
                                                                            tttCTGGTGATGCGCCGCTTGCAGTATTGCTGGTGACACTGGCGAAAGGGATGTCGCCAGTTGAACAAAAACAGTTTACATCCTTGCTTCTGCGGGaatg                  <  2:5101353/100‑3 (MQ=255)
                                                                            tttCTGGTGATGCGCCGCTTGCAGTATTGCTGGTGACACTGGCGAAAGAGATGACGCCAGTTGAACAAAAACAGTTTACATCCTTGCTTCTGCGGGAAgg                  >  1:4705151/1‑100 (MQ=255)
                                                                              tCTGGTGATGCGCCGCTTGCAGTATTGCTGGTGACACTGGCGAAAGAGATGACGAAAAAAAAAAAAAAACAGTTTACATCCTTGCTTCTGCGGGAAGGGa                <  2:2360435/100‑1 (MQ=255)
                                                                                 ggTGATGCGCCGCTTGCAGTATTGCTGGTGACACTGGCGAAAGAGATGACGCCAGTTGAACAAAAACAGTTTACATCCTTGCTTCTGCGGGAAGGGATTa             >  1:5380496/1‑100 (MQ=255)
                                                                                  gTGATGCGCCGCTTGCAGTATTGCTGGTGACACTGGCGAAAGAGATGACGCCAGTTGAACAAAAACAGTTTACATCCTTGCTTCTGCGGGAAGGGATTAt            >  1:2422203/1‑100 (MQ=255)
                                                                                    gttgggccgcTTGCAGTATTGCTGGTGACTCTGGCGAAAGAGATAACGCCAGTTGAACAAAAACAGTTTGCATTCTTGCTTCTGCGGGAAGGGATTATCg          <  2:482093/95‑1 (MQ=255)
                                                                                    gATGCGCCGCTTGCAGTATTGCTGGTGACACTGGCGAAAGAGATGGCGCCAGTTGAACAAAAACAGTTTACATCCTTGCTTCTGCGGGAAGGGATTATCg          >  2:4192229/1‑100 (MQ=255)
                                                                                    gATGCGCCGCTTGCAGTATTGCTGGTGACACTGGCGAAAGAGATGACGCCAGTTGAACAAAAACAGTTTACATCCTTGCTTCTGCGGGAAGGGATTATCg          >  2:1005256/1‑100 (MQ=255)
                                                                                    gATGCACCGCCTGCAGTATTGCTGGTGAAACTGGCGAAAGAGATGAAAAAAAAAAAAAAAAAACAGTTTACATCCTTGCTTCTGCGGGAAGGGATTATcn          <  2:4538340/100‑2 (MQ=255)
                                                                                     aTGCGCCGCTTGCAGTATTGCTGGTGACACTGGCGAAAGAGATGACGCCAGTTGAACAAAAACAGTTTACATCCTTGCTTCTGCGGGAAGGGATTATCgg         <  1:2836959/100‑1 (MQ=255)
                                                                                      tGCGCCGCTTGCAGTATTGCTGGTGACACTGGCGAAAGAGCGGACGCCAGTTGAACAAAAACAGTTTACATCCTTGCTTCTGCGGGAAGGGAGTATCGGa        >  2:3627951/1‑100 (MQ=255)
                                                                                      tGCGCCGCTTGCAGTATTGCTGGTGACACTGGCGAAAGAGATGACGCCAGTTGAACAAAAACAGTTTACATCCTTGCTTCTGCGGGAAGGGATTATCggc        >  1:5567943/1‑99 (MQ=255)
                                                                                      tGCGCCGCTTGCAGTATTGCTGGTGACACTGGCGAAAGAGATGACGCCAGTTGAACAAAAACAGTTTACATCCTTGCTTCTGCGGGAAGGGATTATCGGa        <  1:2827152/100‑1 (MQ=255)
                                                                                        cgccgcTTGCAGTATTGCTGGTGACACTGGCGAAAGAGATGACGCCAGTTGAACAAAAACAGTTTACATCCTTGCTTCTGCGGGAAGGGATTATCGGAtt      >  2:4260786/1‑100 (MQ=255)
                                                                                         gccgcTTGCAGTATTGCTGGTGACACTGGCGAAAGAGATGACGCCAGTTGAACAAAAACAGTTTACATCCTTGCTTCTGCGGGAAGGGATTATCGGATTg     >  1:18135/1‑100 (MQ=255)
                                                                                         gccgcTTGCAGTATTGCTGGTGACACTGGCGAAAGAGATGACGCCAGTTGAACAAAAACAGTTTACATCCTTGCTTCTGCGGGAAGGGATTATCGGATTg     <  2:4976899/100‑1 (MQ=255)
                                                                                         gccgcTTGCAGTATTGCTGGTGACACTGGCGAAAGAGATGACGCCAGTTGAACAAAAACAGTTTACATCCTTGCTTCTGCGGGAAGGGATTATCGGATTg     >  2:2027139/1‑100 (MQ=255)
                                                                                          ccgcTTGCAGTATTGCTGGTGACACTGGCGAAAGAGATGACGCCAGTTGAACAAAAACAGTTTACATCCTTGCTTCTGCGGGAAGGGATTATCGGATTGt    <  1:1507330/100‑1 (MQ=255)
                                                                                           cgcTTGCAGTATTGCTGGTGACACTGGCGAAAGAGATGACGCCAGTTGAACAAAAACAGTTTACATCCTTGCTTCTGCGGGAAGGGATTATCGGATTGtt   >  2:830655/1‑100 (MQ=255)
                                                                                            gcTTGCAGTATTGCTGGTGACACTGGCGAAAGAGATGACGCCAGTTGAACAAAAACAGTTTACATCCTTGCTTCTGCGGGAAGGGATTATCGGATTGTTa  >  2:1214445/1‑100 (MQ=255)
                                                                                                 |                                                                                              
GTTCACGTCAATGAACAGGTTCGTGAATATATTCGCAATGCTGAACGTCCAGAAGGTCAGGGAGAAGCGCCTGCCCTTTCTGGTGATGCGCCGCTTGAAGTATTGCTGGTGACACTGGCGAAAGAGATGACGCCAGTTGAACAAAAACAGTTTACATCCTTGCTTCTGCGGGAAGGGATTATCGGATTGTTA  >  NC_000913/2518590‑2518781

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