Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,705,971 C→A 21.0% intergenic (‑181/+127) dppB ← / ← dppA dipeptide/heme ABC transporter permease/dipeptide/heme ABC transporter periplasmic binding protein; dipeptide chemotaxis receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,705,9710CA21.0% 487.3 / 113.8 187intergenic (‑181/+127)dppB/dppAdipeptide/heme ABC transporter permease/dipeptide/heme ABC transporter periplasmic binding protein; dipeptide chemotaxis receptor
Reads supporting (aligned to +/- strand):  ref base C (63/84);  new base A (34/5);  total (98/89)
Fisher's exact test for biased strand distribution p-value = 4.85e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.47e-01

Only 100 of 197 total aligned reads displayed.

GGAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCT  >  NC_000913/3705874‑3706070
                                                                                                 |                                                                                                    
ggAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATgcgg                                                                                                    <  1:572417/100‑1 (MQ=255)
 gAATGAAGTATTGCTGCCGGATGGTGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATgcggc                                                                                                   >  2:3320276/1‑100 (MQ=255)
 gAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATgcggc                                                                                                   >  2:556278/1‑100 (MQ=255)
  aaTGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATgcggcg                                                                                                  <  1:121065/100‑1 (MQ=255)
  aaTGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATgcggcg                                                                                                  <  2:2232082/100‑1 (MQ=255)
    tGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCCCCGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGtt                                                                                                >  1:881979/1‑100 (MQ=255)
     gAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTa                                                                                               <  2:1807929/100‑1 (MQ=255)
     gAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTa                                                                                               >  2:2120707/1‑100 (MQ=255)
     gAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTa                                                                                               <  2:2447375/100‑1 (MQ=255)
        gTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACg                                                                                            <  1:2390087/100‑1 (MQ=255)
        gTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACg                                                                                            <  2:401802/100‑1 (MQ=255)
        gTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACg                                                                                            <  2:1847870/100‑1 (MQ=255)
         tATTGCTGCCGGATGGCGCTGCACGCAGCGGCACAGCGCCCCAGGGTATCCCGATTGTTGAAAAAGTGCGCTTTCTTTATGCCGGATGCGGCGATAAcac                                                                                           >  1:1681203/1‑98 (MQ=25)
         tATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGc                                                                                           <  2:1913299/100‑1 (MQ=255)
          aTTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGcc                                                                                          <  1:3556715/100‑1 (MQ=255)
             gctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTa                                                                                       >  2:75278/1‑100 (MQ=255)
               tgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATc                                                                                     >  2:2914380/1‑100 (MQ=255)
               tgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATc                                                                                     >  1:1907426/1‑100 (MQ=255)
                gcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATcc                                                                                    <  2:1130074/100‑1 (MQ=255)
                gcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATcc                                                                                    >  2:3550257/1‑100 (MQ=255)
                 ccGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGTTATACAGATTGTTGACAAAGTGCGCTTTGTTTATTCCGGATGCGGCGTTAACGCCTTATCCg                                                                                   <  1:577404/100‑1 (MQ=255)
                  cGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTTCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCgg                                                                                  <  1:2074348/100‑1 (MQ=255)
                  cGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATccag                                                                                  <  1:1556217/100‑3 (MQ=255)
                   ggATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGc                                                                                 >  2:2198290/1‑100 (MQ=255)
                    gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGcc                                                                                <  1:678891/100‑1 (MQ=255)
                    gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGcc                                                                                >  2:3297105/1‑100 (MQ=255)
                     aTGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTACCGCCTTATCCGGcc                                                                                >  2:815701/1‑99 (MQ=255)
                     aTGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCt                                                                               >  2:1610332/1‑100 (MQ=255)
                      tGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTa                                                                              >  1:2546238/1‑100 (MQ=255)
                      tGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTa                                                                              >  1:1480729/1‑100 (MQ=255)
                      tGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTa                                                                              >  1:1835129/1‑100 (MQ=255)
                       ggCGCGGCACGCAGGGGCACAGCGCCACCGTGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGcc                                                                                <  1:3196189/97‑1 (MQ=255)
                         cgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACCAAGCGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAcca                                                                           >  1:2996871/1‑98 (MQ=255)
                         cgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGcataaac                                                                           <  2:2960953/100‑7 (MQ=25)
                         cgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGAGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑aa                                                                          >  2:198668/1‑100 (MQ=255)
                          gcTGCACGCAGAGGCACAGCGCCACCGTGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGcataaac                                                                           <  1:2370949/99‑7 (MQ=25)
                          gcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑aaa                                                                         >  2:2656896/1‑100 (MQ=255)
                          gcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGcataaaca                                                                          <  1:764793/100‑8 (MQ=25)
                             gCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGGCTTATCCGGCCTAC‑AAAATc                                                                      >  1:528286/1‑100 (MQ=255)
                             gCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATc                                                                      <  2:615006/100‑1 (MQ=255)
                              cACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCg                                                                     <  2:1395981/100‑1 (MQ=255)
                                cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTg                                                                   >  1:311113/1‑100 (MQ=255)
                                cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTg                                                                   <  1:1709244/100‑1 (MQ=255)
                                 gCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGc                                                                  <  1:2704921/100‑1 (MQ=255)
                                  cAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGGCTTATCCGGCCTAC‑AAAATCGTGCa                                                                 >  1:452354/1‑100 (MQ=255)
                                      ggCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAAtt                                                             >  2:3724192/1‑100 (MQ=255)
                                       gCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTc                                                            <  1:3262211/100‑1 (MQ=255)
                                        cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCa                                                           <  1:3183636/100‑1 (MQ=255)
                                        cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGATAACGACTTATCCGCCCTAC‑AAAATCGTGCAAATTCa                                                           >  1:2563651/1‑100 (MQ=255)
                                          caGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAAt                                                         >  2:1048401/1‑100 (MQ=255)
                                          caGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAAt                                                         <  2:2331287/100‑1 (MQ=255)
                                          caGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAAt                                                         <  1:2869746/100‑1 (MQ=255)
                                           aGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAAtn                                                        <  1:1120755/100‑2 (MQ=255)
                                                   cGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACAtt                                                    <  1:2555084/96‑1 (MQ=255)
                                                   cGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGAAAAC‑AAACTCGTGCAAATTCCATACCTTg                                                   <  1:3023817/97‑1 (MQ=17)
                                                    gggTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGa                                               >  2:3093577/1‑100 (MQ=255)
                                                         taCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAAtttt                                          <  1:1610325/100‑1 (MQ=255)
                                                           cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAATTTTCg                                        >  2:3323367/1‑100 (MQ=255)
                                                           cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAATTTTCg                                        <  1:3212783/100‑1 (MQ=255)
                                                           cAGATTGTTGACAAAGTGCGCTTTGTTTATGCAGGATGCAGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAATTTTCg                                        <  1:2664375/100‑1 (MQ=255)
                                                           cAGATTGTTGACAAAGTGCACTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAATTTTCg                                        <  1:2145348/100‑1 (MQ=255)
                                                                 gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGcc                                  <  2:2313728/100‑1 (MQ=255)
                                                                 gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGcc                                  >  1:1370161/1‑100 (MQ=255)
                                                                  ttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCt                                 >  1:892696/1‑100 (MQ=255)
                                                                  ttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCt                                 >  1:1765383/1‑100 (MQ=255)
                                                                   tgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTg                                <  2:1226554/100‑1 (MQ=255)
                                                                    gACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGa                               >  1:1163939/1‑100 (MQ=255)
                                                                         aGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGc                          >  2:989331/1‑100 (MQ=255)
                                                                         aGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGc                          <  1:1048396/100‑1 (MQ=255)
                                                                              gcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAgc                     <  1:1224607/100‑1 (MQ=255)
                                                                              gcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAgc                     >  1:1953889/1‑100 (MQ=255)
                                                                              gcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAgc                     <  2:1180077/100‑1 (MQ=255)
                                                                               cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAgcg                    >  2:3353863/1‑100 (MQ=255)
                                                                               cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAgcg                    <  2:774409/100‑1 (MQ=255)
                                                                                 ttttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAAACCTGCAAATTCAATATATTGCAGGATACGTGTAGGCCTGATAAGCGTAGCGCa                  >  2:882405/9‑100 (MQ=17)
                                                                                  tGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCAt                 <  1:3753949/100‑1 (MQ=255)
                                                                                  cttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCATAC‑AAGATCGTGCAAATTCAATATATTGCATGAGCCACGTAGGCCTGATAAGCGTAGCGCAt                 >  1:1307495/8‑100 (MQ=17)
                                                                                   ttgtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATc                >  1:1135662/7‑100 (MQ=18)
                                                                                   gTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATc                >  1:2274810/1‑100 (MQ=255)
                                                                                    tttATGCCGGACGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCa               <  2:3401971/100‑1 (MQ=255)
                                                                                     ttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAg              <  1:136937/100‑1 (MQ=255)
                                                                                     gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAGATCGTGCAAATTCAATATATTGCATGAGCCACGTAGGCCTGATAAGCGTAGCGCATCAg              >  2:1549921/5‑100 (MQ=17)
                                                                                     gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAg              >  1:2117020/5‑100 (MQ=18)
                                                                                     gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑‑AAATCATGCAAATCCAAAATATTGCAGGAGTTATGTAGGCCTGATAAGCGTAGCGCATCAgg             >  2:3120682/5‑100 (MQ=21)
                                                                                       aggCCGGATAAGGCGTTTACGCCGCATCTGGCCTAC‑AAATTCTTCCAAATTCAATATATTGCAGAAATCATGTAGGCCTGATAAGCGTAGCGCATCAGGc            >  1:3698658/3‑100 (MQ=18)
                                                                                       aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAATCGCTTAAATTCAATATATTGCACCATCCTCGTAGGCCTGATAAGCGTAGCGCATCAGGc            >  2:1019220/3‑100 (MQ=21)
                                                                                       aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACTTTGCAAATTCAATATATTGCAATCTCCGTGTAGGCCTGATAAGCGTAGCGCATCAGGc            >  1:722526/3‑100 (MQ=18)
                                                                                       aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGc            >  1:1848525/3‑100 (MQ=18)
                                                                                       aggCCGGATAAGGCGTTAACGCCGCATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGc            >  2:618701/3‑100 (MQ=18)
                                                                                       aggCCGGATAAGGCGTTAACGCCGCATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGc            >  2:2151747/3‑100 (MQ=18)
                                                                                        tGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCATGATAAGCGGAGCGCATCAGGc            >  2:2834524/1‑99 (MQ=255)
                                                                                        ggCCGGATAAGGCGTTTACGCCGCATCTGGCCTAC‑AAATTCTTCCAAATTCAATATATTGCAGAAATCATGTAGGCCTGATAAGCGTAGCGCATCAGGCa           >  2:1260548/2‑100 (MQ=17)
                                                                                        ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAATCGCTTAAATTCAATATATTGCACCATCCTCGTAGGCCTGATAAGCGTAGCGCATCAGGCa           >  2:1784483/2‑100 (MQ=21)
                                                                                        ggCCGGATAAGGCGTTAACGCCGCATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCa           >  2:3661315/2‑100 (MQ=18)
                                                                                         gCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAATCGCTTAAATTCAATATATTGCACCATCCTCGTAGGCCTGATAAGCGTAGCGCATCAGGCaa          >  1:3347711/1‑99 (MQ=21)
                                                                                          ccGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCAGGACCAT‑GTAGGCCTGATAAGCGTAGCGCATCAGGCagt         >  1:929404/1‑98 (MQ=18)
                                                                                            ggATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCAtttt       >  1:1090170/1‑100 (MQ=255)
                                                                                              aTGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg     <  1:1032340/100‑1 (MQ=255)
                                                                                                gcggcgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcc   >  2:1877741/1‑100 (MQ=255)
                                                                                                 cggcgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcct  >  2:151445/1‑100 (MQ=255)
                                                                                                 |                                                                                                    
GGAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCT  >  NC_000913/3705874‑3706070

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