Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,358,482 T→A 100% I46F (ATC→TTC)  cadA ← lysine decarboxylase, acid‑inducible

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,358,4820TA96.6% 290.0 / 0.5 90I46F (ATC→TTC) cadAlysine decarboxylase, acid‑inducible
Reads supporting (aligned to +/- strand):  ref base T (0/3);  new base A (42/44);  total (42/48)
Fisher's exact test for biased strand distribution p-value = 2.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.76e-01

AGGTTCTCGTTCATTTTGCTAATTTCTTCGCACAGCTCGAGATTATATTTATCCCAGTCAAAAATAACGCCGCACAGACGCGCATTGTTTTCGATCAGTTTTAATAAGTCGTCACGGTCGTTCGGGTAAACAATCTGGAAGTTCAGACGTTCAAGCGCGCGATGAAGTTCACGGATGGGTTCTTCTTTAAA  >  NC_000913/4358388‑4358578
                                                                                              |                                                                                                
aGGTTCTCGTTCATTTTGCTAATTTCTTCGCACAGCTCGAGATTATATTTATCCCAGTCAAAAATAACGCCGCACAGACGCGCATTGTTTTCGAACAGtt                                                                                             <  2:2156765/100‑1 (MQ=255)
aGGTTCTCGTTCATTTTGCTAATTTCTTCGCACAGCTCGAGATTATATTTATCCCAGTCAAAAATAACGCCGCACAGACGCGCATTGTTTTCGAACAGtt                                                                                             <  2:828395/100‑1 (MQ=255)
  gTTCTCGTTCATTTTGCTAATTTCTTCGCACAGCTCGAGATTATATTTATCCCAGTCAAAAATAACGCCGCACAGACGCGCATTGTTTTCGAACAGtttt                                                                                           <  1:489967/100‑1 (MQ=255)
  gTTCTCGTTCATTTTGCTAATTTCTTCGCACAGCTCGAGATTATATTTATCCCAGTCAAAAATAACGCCGCACAGACGCGCATTGTTTTCGAACAGtttt                                                                                           <  1:758268/100‑1 (MQ=255)
   ttCTCGTTCATTTTGCTAATTTCTTCGCACAGCTTGAGATTCTATTTATCCCAGTCAAAAATAACGCAGCAAAGACGCGCATTGTTTTCGAACAGTTTta                                                                                          >  2:97097/1‑100 (MQ=255)
      tcGTTCATTTTTCTAATTTCTTCGCACAGCTCTAGATTATATTTATCCCAGTCAAAAATAACGCCGCACAGACGCGCATTGTTTTCGAACAGTTTtaata                                                                                       <  2:1159593/100‑1 (MQ=255)
      tcGTTCATTTTGCTAATTTCTTCGCATAGCTCGAGATTATATTTATCCCAGTCAAAAATAACGCCGCACAGACGCGCATTGTTTTCGAACAGTTTtaata                                                                                       <  2:2189741/100‑1 (MQ=255)
      tcGTTCATTTTGCTAATTTCTTCGCACAGCTCGAGATTATATTTATCCCAGTCAAAAATAACGCCGCACAGACGCGCATTGTTTTCGAACAGTTTtaata                                                                                       >  1:2443109/1‑100 (MQ=255)
       cGTTCATTTTGCTAATTTCTTCGCACAGCTCGAGATTATATTTATCCCAGTCAAAAATAACGCCGCACAGACGCGCATTGTTTTCGAACAGTTTtaataa                                                                                      >  1:1353082/1‑100 (MQ=255)
         ttCATTTTGCTAATTTCTTCGCACAGCTCGAGATTATATTTATCCCAGTCAAAAATAACGCCGCACAGGCGCGCATTGTTTTCGAAGAGTTTTAATAAgt                                                                                    <  1:1724052/100‑1 (MQ=255)
         ttCATTTTAATAATTTCTTCGAACAGCTCTAGATTATAAAAAAAAAAAAAAAAAATAACGCCGCACAGACGCGCATTGTTTTCGAACAGTTTtaata                                                                                       <  2:416258/97‑1 (MQ=255)
          tCATTTTGCTAATTTCTTCGCACAGCTCGAGATTATATTTATCCTAGTCTAAAATAACGCCGCACACACTCGAATTGTTTTCGAAAAGATTtaat                                                                                        >  2:2336148/1‑95 (MQ=255)
               tttcTAATGTCTTCGCACAGGTCGAGATTATATTCATTCCGGTCAAAACTAACGCCGCACAGACGCGCATTGTTTTCGAACAGGTTTCATAAGTCGTCa                                                                               <  2:2356359/96‑1 (MQ=255)
                 gCTAATTTCTTCGCACAGCTCGAGATTATATTTATCCCAGTCAAAAATAACGCCGCACAGACGCGCATTGTTTTCGAACAGTTTTAATAAGTCGTCACgg                                                                            <  2:2244334/100‑1 (MQ=255)
                  cTAATTTCTTCGCACAGCTCGAGATTATATTTATCCCAGTCAAAAATAACGCCGCACAGACGCGCATTGTTTTCGAACAGTTTTAATAAGTCGTCACGGt                                                                           >  2:654341/1‑100 (MQ=255)
                  cTAATTTCTTCGCACAGCTCGAGATTATATTTATCCCAGTCAAAAATAACGCCGCACAGACGCGCATTGTTTTCGAACAGTTTTAATAAGTCGTCACGGt                                                                           <  2:947269/100‑1 (MQ=255)
                    aaTTTCTTCGCACAGCTCGAGATTATATTTATCCCAGTCAAAAATAACGCCGCACAGACGCGCATTGTTTTCGAACAGTTTTAATAAGTCGTCACGGTCg                                                                         >  1:2286980/1‑100 (MQ=255)
                     aTTTCTTCGCACAGCTCGAGATTATATTTATCCCAGTCAAAAATAACGCCGCACAGACGCGCATTGTTTTCGAACAGTTTTAATAAGTCGTCACGGTCGt                                                                        <  1:1793256/100‑1 (MQ=255)
                      tttCTTCGCACAGCTCGAGATTATATTTATCCCAGTCAAAAATAACGCCGCACAGACGCGCATTGTTTTCGATCAGTTTTAATAAGTCGTCACGGTCGtt                                                                       <  1:816116/100‑1 (MQ=255)
                      tttCTTCGCACAGCTCGAGATTATATTTATCCCAGTCAAAAATAACGCCGCACAGACGCGCATTGTTTTCGAACAGTTTTAATAAGTCGTCACGGTCGtt                                                                       >  1:1566450/1‑100 (MQ=255)
                       ttcttcGCACAGCTCGAGATTATATTTATCCCAGTCAAAAATAACGCCGCACAGACGCGCATTGTTTTCGAACAGTTTTAATAAGTCGTCACGGTCGTTc                                                                      >  2:2149077/1‑100 (MQ=255)
                         cttcGCACAGCTCGAGATTATATTTATCCCAGTCAAAAATAACGCCGCACAGACGCGCATTGTTTTCGAACAGTTTTAATAAGTCGTCACGGTCGTTCgg                                                                    >  1:235465/1‑100 (MQ=255)
                          ttcGCACAGCTCGAGATTATATTTATCCCAGTCAAAAATAACGCCGCACAGACGCGCATTGTTTTCGAACAGTTTTAATAAGTCGTCACGGTCGTTCggg                                                                   <  1:841982/100‑1 (MQ=255)
                          ttcGCACAGCTCGAGATTATATTTATCCCAGTCAAAAATAACGCCGCACAGACGCGCATTGTTTTCGAACAGTTTTAATAAGTCGTCACGGTCGTTCggg                                                                   <  2:1769478/100‑1 (MQ=255)
                          ttcGCACAGCTCGAGATTATATTTATCCCAGTCAAAAATAACGCCGCACAGACGCGCATTGTTTTCGAACAGTTTTAATAAGTCGTCACGGTCGTTCggg                                                                   <  2:960596/100‑1 (MQ=255)
                           tcGCACAGCTCGAGATTATATTTATCCCAGTCAAAAATAACGCCGCACAGACGCGCATTGTTTTCGAACAGTTTTAATAAGTCGTCACGGTCGTTCGGGt                                                                  <  2:1712538/100‑1 (MQ=255)
                           tcGCACAGCTCGAGATTATATTTATCCCAGTCAAAAATAACGCCGCACAGACGCGCATTGTTTTCGAACAGTTTTAATAAGTCGTCACGGTCGTTCGGGt                                                                  <  1:270902/100‑1 (MQ=255)
                           tcGCACAGCTCGAGATTATATTTATCCCAGTCAAAAATAACCCCGCACAGACGCGCATTGTTTTCGAACAGTTTTAATAAGTCGTCACGGTCGTTCGGGt                                                                  <  1:2294223/100‑1 (MQ=255)
                             gCACAGCTCGAGATTATATTTATCCCAGTCAAAAATAACGCCGCACAGACGCGCATTGTTTTCGAACAGTTTTAATAAGTCGTCACGGTCGTTCGGGTaa                                                                <  2:1195724/100‑1 (MQ=255)
                              cacaGCTCGAGATTATATTTATCCCAGTCAAAAATAACGCCGCACAGACGCGCATTGTTTTCGAACAGTTTTAATAAGTCGTCACGGTCGTTCGGGTaaa                                                               >  2:1769708/1‑100 (MQ=255)
                               acaGCTCGAGATTATATTTATCCCAGTCAAAAATAACGCCGCACAGACGCGCATTGTTTTCGAACAGTTTTAATAAGTCGTCACGGTCGTTCGGGTAAAc                                                              <  1:713803/100‑1 (MQ=255)
                                  gCTCGAGATTATATTTATCCCAGTCAAAAATAACGCCGCACAGACGCGCATTGTTTTCGAACAGTTTTAATAAGTCGTCACGGTCGTTCGGGTAAACAAt                                                           >  2:491231/1‑100 (MQ=255)
                                  gCTCGAGATTATATTTATCCCAGTCAAAAATAACGCCGCACAGACGCGCATTGTTTTCGAACAGTTTTAATAAGTCGTCACGGTCGTTCGGGTAAACAAt                                                           >  1:1249219/1‑100 (MQ=255)
                                  gCTCGAGATTATATTTATCCCAGTCAAAAATAACGCCGCACAGACGCGCATTGTTTTCGAACAGTTTTAATAAGTCGTCACGGTCGTTCGGGTAAACAAt                                                           >  1:1333235/1‑100 (MQ=255)
                                   cTCGAGATTATATTTATCCCAGTCAAAAATAACGCCGCACAGACGCGCATTGTTTTCGAACAGTTTTAATAAGTCGTCACGGTCGTTCGGGTAAACAATc                                                          <  2:1521892/100‑1 (MQ=255)
                                    tCGAGATTATATTTATCCCAGTCAAAAATAACGCCGCACAGACGCGCATTGTTTTCGAACAGTTTTAATAAGTCGTCACGGTCGTTCGGGTAAACAATCt                                                         >  1:313834/1‑100 (MQ=255)
                                     cGAGATTATATTTATCCCAGTCAAAAATAACGCCGCACAGACGCGCATTGTTTTCGAACAGTTTTAATAAGTCGTCACGGTCGTTCGGGTAAACAATCt                                                         <  1:2487409/99‑1 (MQ=255)
                                     cGAGATTATATTTATCCCAGTCAAAAATAACGCCGCACAGACGCGCATTGTTTTCGAACAGTTTTAATAAGTCGTCACGGTCGTTCGGGTAAACAATCTg                                                        >  2:149501/1‑100 (MQ=255)
                                     cGAGATTATATTTATCCCAGTCAAAAATAACGCCGCACAGACGCGCATTGTTTTCGAACAGTTTTAATAAGTCGTCACGGTCGTTCGGGTAAACAATCTg                                                        >  2:836968/1‑100 (MQ=255)
                                     cGAGATTATATTTATCCAAGTAAAAAATAACGCCGCACAGACGCGCATTGTTTTCGAACAGTTATAATAAGTCGTCACGGTCGTTCGGGTAAACAATCt                                                         >  2:1064629/1‑99 (MQ=255)
                                       agaTTATATTTATCCCAGTCAAAAATAACGCCGCACAGACGCGCATTGTTTTCGAACAGTTTTAATAAGTCGTCACGGTCGTTCGGGTAAACAATCTGGa                                                      <  1:768991/100‑1 (MQ=255)
                                        gaTTATATTTATCCCAGTCAAAAATAACGCCGCACAGACGCGCATTGTTTTCGAACAGTTTTAATAAGTCGTCACGGTCGTTCGGGTAAACAATCTGGaa                                                     <  2:2236139/100‑1 (MQ=255)
                                         aTTATATTTATCCCAGTCAAAAATAACGCCGCACAGACGCGCATTGTTTTCGATCAGTTTTAATAAGTCGTCACGGTCGTTCGGGTAAACAATCTGGAAg                                                    <  1:2502307/100‑1 (MQ=255)
                                         aTTATATTTATCCCAGTCAAAAATAACGCCGCACAGACGCGCATTGTTTTCGAACAGTTTTAATAAGTCGTCACGGTCGTTCGGGTAAACAATCTGGAAg                                                    <  1:404469/100‑1 (MQ=255)
                                          ttATATTTATCCCAGTCAAAAATAACGCCGCACAGACGCTCATTGTTTTCGAACAGTTTTAATAAGTCGTCACGGTCGTTCGGGTAAACAATCTGGAAg                                                    >  2:2099825/1‑99 (MQ=255)
                                            atatTTATCCCAGTCAAAAATAACGCCGCACAGACGCGCATTGTTTTCGAACAGTTTTAATAAGTCGTCACGGTCGTTCGGGTAAACAATCTGGAAGTTc                                                 >  1:2333556/1‑100 (MQ=255)
                                             tatTTATCCCAGTCAAAAATAACGCCGCACAGACGCGCATTGTTTTCGAACAGTTTTAATAAGTCGTCACGGTCGTTCGGGTAAACAATCTGGAAGTTCa                                                >  2:675989/1‑100 (MQ=255)
                                               tttATCCCAGTCAAAAATAACGCCGCACAGACGCGCATTGTTTTCGAACAGTTTTAATAAGTCGTCACGGTCGTTCGGGTAAACAATCTGGAAGTTCAGa                                              <  1:370085/100‑1 (MQ=255)
                                                 tATCCCAGTCAAAAATAACGCCGCACAGACGCGCATTGTTTTCGAACAGTTTTAATAAGTCGTCACGGTCGTTCGGGTAAACAATCTGGAAGTTCAGACg                                            >  1:1517021/1‑100 (MQ=255)
                                                 tATCCCAGTCAAAAATAACGCCGCACAGACGCGCATTGTTTTCGAACAGTTTTAATAAGTCGTCACGGTCGTTCGGGTAAACAATCTGGAAGTTCAGACg                                            >  2:1606337/1‑100 (MQ=255)
                                                   tCCCAGTCAAAAATAACGCCGCACAGACGCGCATTGTTTTCGAACAGTTTTAATAAGTCGTCACGGTCGTTCGGGTAAACAATCTGGAAGTTCAGACGtt                                          >  1:1261498/1‑100 (MQ=255)
                                                    cccAGTCAAAAATAACGCCGCACAGACGCGCATTGTTTTCGAACAGTTTTAATAAGTCGTCACGGTCGTTCGGGTAAACAATCTGGAAGTTCAGACGTTc                                         <  1:1016816/100‑1 (MQ=255)
                                                     ccAGTCAAAAATAACGCCGCACAGACGCGCATTGTTTTCGAACAGTTTTAATAAGTCGTCACGGTCGTTCGGGTAAACAATCTGGAAGTTCAGACGTTCa                                        <  1:1385309/100‑1 (MQ=255)
                                                     ccAGTCAAAAATAACGCCGCACAGACGCGCATTGTTTTCGAACAGTTTTAATAAGTCGTCACGGTCGTTCGGGTAAACAATCTGGAAGTTCAGACGTTCa                                        >  1:150808/1‑100 (MQ=255)
                                                     ccAGTCAAAAATAACGCCGCACAGACGCGCATTGTTTTCGAACAGTTTTAATAAGTCGTCACGGTCGTTCGGGTAAACAATCTGGAAGTTCAGACGTTCa                                        <  1:1510038/100‑1 (MQ=255)
                                                     ccAGTCAAAAATAACGCCGCACAGACGCGCATTGTTTTCGAACAGTTTTAATAAGTCGTCACGGTCGTTCGGGTAAACAATCTGGAAGTTCAGACGTTCa                                        <  2:1028172/100‑1 (MQ=255)
                                                     ccAGTCAAAAATAACGCCGCACAGACGCGCATTGTTTTCGAACAGTTTTAATAAGTCGTCACAGTCGTTCGGGTAAACAATCTGGAAGTTCAGACGTTCa                                        >  2:2236538/1‑100 (MQ=255)
                                                      cAGTCAAAAATAACGCCGCACAGACGCGCATTGTTTTCGAACAGTTTTAATAAGTCGTCACGGTCGTTCGGGTAAACAATCTGGAAGTTCAGACGTTCaa                                       >  1:248097/1‑100 (MQ=255)
                                                      cAGTCAAAAATAACGCCGCACAGACGCGCATTGTTTTCGAACAGTTTTAATAAGTCGTCACGGTCGTTCCGGTAAACAATCTGGAAGATCAGACGTTCaa                                       >  1:1535371/1‑100 (MQ=255)
                                                         tCAAGAATAACGGCGCACAGATGCGCATTGTTTTCGAGCAGTTTTAATAGGTCGTCACGGTCGTTCTGGTAAATAATCTGGAAGTTCAGACGTTCAAgcg                                    <  2:2091256/100‑1 (MQ=255)
                                                          cAAAAATAACGCCGCACAGACGCGCATTGTTTTCGAACAGTTTTAATAAGTCGTCACGGTCGTTCGGGTACACAATCTGGAAGTTCAGACGTTCAAgcgc                                   >  1:1264089/1‑100 (MQ=255)
                                                           aaaaaTAACGCCGCACAGACGCGCATTGTTTTCGAACAGTTTTAATAAGTCGTCACGGTCGTTCGGGTAAACAATCTGGAAGTTCAGACGTTCAAgcgcg                                  >  2:705229/1‑100 (MQ=255)
                                                           aaaaaTAACGCCGCACAGACGCGCATTGTTTTCGAACAGTTTTAATAAGTCGTCACGGTCGTTCGGGTAAACAATCTGGAAGTTCAGACGTTCAAgcgcg                                  >  2:2508671/1‑100 (MQ=255)
                                                            aaaaTAACGCCGCACAGACGCGCATTGTTTTCGAACAGTTTTAATAAGTCGTCACGGTCGTTCGGGTAAACAATCTGGAAGTTCAGACGTTCAAgcgcgc                                 <  2:1855333/100‑1 (MQ=255)
                                                               aTAACGCCGCACAGACGCGCATTGTTTTCGAACAGTTTTAATAAGTCGTCACGGTCGTTCGGGTAAACAATCTGGAAGTTCAGACGTTCAAGCGCGCGAt                              >  1:1506313/1‑100 (MQ=255)
                                                                 aaCGCCGCACAGACGCGCATTGTTTTCGAACAGTTTTAATAAGTCGTCACGGTCGTTCGGGTAAACAATCTGGAAGTTCAGACGTTCAAGCGCGCGATGa                            <  2:1513385/100‑1 (MQ=255)
                                                                  aCGCCGCACAGACGCGCATTGTTTTCGAACAGTTTTAATAAGTCGTCACGGTCGTTCGGGTAAACAATCTGGAAGTTCAGACGTTCAAGCGCGCGATGaa                           >  2:1284706/1‑100 (MQ=255)
                                                                       gcACAGACGCGCATTGTTTTCGAACAGTTTTAATAAGTCGTCACGGTCGTTCGGGTAAACAATCTGGAAGTTCAGACGTTCAAGCGCGCGATGAAGTTCa                      >  2:468378/1‑100 (MQ=255)
                                                                       gcACAGACGCGCATTGTTTTCGAACAGTTTTAATAAGTCGCCACGGTCGTTCGGGTAAACAATCTGGAAGTTCAGACGTTCAAGCGCGCGATGAAGTTCa                      >  2:1845100/1‑100 (MQ=255)
                                                                        cacaGACGCGCATTGTTTTCGAACAGTTTTAATAAGTCGTCACGGTCGTTCGGGTCAACCATCTGGAAGTTCAGACGTTCAAGCGCGCGATGAAGTTCAc                     >  2:66167/1‑100 (MQ=255)
                                                                        cacaGACGCGCATTGTTTTCGAACAGTTTTAATAAGTCGTCACGGTCGTTCGGGTAAACAATCTGGAAGTTCAGACGTTCAAGCGCGCGATGAAGTTCAc                     >  1:692678/1‑100 (MQ=255)
                                                                        cacaGACGCGCATTGTTTTCGAACAGTTTTAATAAGTCGTCACGGTCGTTCGGGTAAACAATCTGGAAGTTCAGACGTTCAAGCGCGCGATGAAGTTCAc                     <  1:2546370/100‑1 (MQ=255)
                                                                         acaGACGCGCATTGTTTTCGAACAGTTTTAATAAGTCGTCACGGTCGTTCGGGTAAACAATCTGGAAGTTCAGACGTTCAAGCGCGCGATGAAGTTCACg                    <  1:1994861/100‑1 (MQ=255)
                                                                         acaGACGCGCATTGTTTTCGAACAGTTTTAATAAGTCGTCACGGTCGTTCGGGTAAACAATCTGGAAGTTCAGACGTTCAAGCGCGCGATGAAGTTCACg                    <  2:2130703/100‑1 (MQ=255)
                                                                          caGACGCGCATTGTTTTCGATCAGTTTTAATTAGTCGTCACGGTCGTTCGGGTAAACAATCTGGAAGTTCAGACGTTCAAGCGCGCGATGAAGTTCACg                    <  2:1644420/99‑1 (MQ=255)
                                                                          caGACGCGCATTGTTTTCGAACAGTTTTAATAAGTCGTCACGGTCGTTCGGGTAAACAATCTGGAAGTTCAGACGTTCAAGCGCGCGATGAAGTTCACgg                   >  1:1715906/1‑100 (MQ=255)
                                                                            gACGCGCATTGTTTTCGAACAGTTTTAATAAGTCGTCACGGTCGTTCGGGTAAACAATCTGGAAGTTCAGACGTTCAAGCGCGCGATGAAGTTCACGGAt                 >  2:1064653/1‑100 (MQ=255)
                                                                             aCGCGCATTGTTTTCGAACAGTTTTAATAAGTCGTCACGGTCGTTCGGGTAAACAATCTGGAAGTTCAGACGTTCAAGCGCGCGATGAAGTTCACGGATg                <  1:2364597/100‑1 (MQ=255)
                                                                             aCGCGCATTGTTTTCGAACAGTTTTAATAAGTCGTCACGGTCGTTCGGGTAAACAATCTGGAAGTTCAGACGTTCAAGCGCGCGATGAAGTTCACGGATg                <  2:2023314/100‑1 (MQ=255)
                                                                             aCGCGCATTGTTTTCGAACAGTTTTAATAAGTCGTCACGGTCGTTCGGGTAAACAATCTGGAAGTTCAGACGTTCAAGCGCGCGATGAAGTTCACGGATg                <  2:2066304/100‑1 (MQ=255)
                                                                                    ttGTTTTCGAACAGTTTTAATAAGTCGTCACGGTCGTTCGGGTAAACAATCTGGAAGTTCAGACGTTCAAGCGCGCGATGAAGTTCACGGATGGGttcn          <  2:352925/99‑2 (MQ=255)
                                                                                      gTTTTCGAACAGTTTTAATAAGTCGTCACGGTCGTTCGGGTAAACAATCTGGAAGTTCAGACGTTCAAGCGCGCGATGAAGTTCACGGATGGGttcttc        >  2:1021748/1‑99 (MQ=255)
                                                                                      gTTTTCGAACAGTTTTAATAAGTCGTCACGGTCGTTCGGGTAAACAATCTGGAAGTTCAGACGTTCAAGCGCGCGATGAAGTTCACGGATGGGTtcttct       <  1:2458311/100‑1 (MQ=255)
                                                                                      gTTTTCGAACAGTTTTAATAAGTCGTCACGGTCGTTCGGGTAAACAATCTGGAAGTTCAGACGTTCAAGCGCGCGATGAAGTTCACGGATGGGTtcttct       <  1:800907/100‑1 (MQ=255)
                                                                                       ttttCGAACAGTTTTACTAAGTCGTCCAGGTCGTTTGGGTCAACAATCTGGAAGTTCAGACGGTCAAGCGTGCGTTGAACTTTACGCATGGGTTcttctt      >  2:2135946/1‑100 (MQ=255)
                                                                                         ttCGAACAGTTTTAATAAGTCGTCACGGTCGTTCGGGTAAACAATCTGGAAGTTCAGACGTTCAAGCGCGCGATGAAGTTCACGGATGGGTTCTTCTTTa    >  1:2345877/1‑100 (MQ=255)
                                                                                         ttCGAACAGTTTTAATAAGTCGTCACGGTCGTTCGGGTAAACAATCTGGAAGTTCAGACGTTCAAGCGCGCGATGAAGTTCACGGATGGGTTCTTCTTTa    <  1:296600/100‑1 (MQ=255)
                                                                                          tCGAACAGTTTTCATAAGTCGTAAGGGTCGTTCGGGTGAAACATCTAGAAGTTCAGATGTTAAAGCACGCCATCAAGTTCACGGATGGGTTCTTCGTTa  <  2:619162/100‑2 (MQ=255)
                                                                                           cGAACAGTTTTAATAAGTCGTCACGGTCGTTCGGGTAAACAATCTGGAAGTTCAGACGTTCAAGCGCGCGATGAAGTTCACGGATGGGTTCTTCTTTaaa  <  1:2552401/100‑1 (MQ=255)
                                                                                           cGAACAGTTTTAATAAGTCGTCACGGTCGTTCGGGTAAACAATCTGGAAGTTCAGACGTTCAAGCGCGCGATGAAGTTCACGGATGGGTTCTTCTTTaaa  >  2:2356178/1‑100 (MQ=255)
                                                                                              |                                                                                                
AGGTTCTCGTTCATTTTGCTAATTTCTTCGCACAGCTCGAGATTATATTTATCCCAGTCAAAAATAACGCCGCACAGACGCGCATTGTTTTCGATCAGTTTTAATAAGTCGTCACGGTCGTTCGGGTAAACAATCTGGAAGTTCAGACGTTCAAGCGCGCGATGAAGTTCACGGATGGGTTCTTCTTTAAA  >  NC_000913/4358388‑4358578

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: