Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,705,971 C→A 29.6% 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,9710CA29.6% 405.4 / 182.4 200intergenic (‑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 (66/75);  new base A (49/10);  total (115/85)
Fisher's exact test for biased strand distribution p-value = 1.73e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.27e-01

Only 100 of 212 total aligned reads displayed.

AGGGAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCT  >  NC_000913/3705872‑3706070
                                                                                                   |                                                                                                           
  ggAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATgcgg                                                                                                           <  2:1704215/100‑1 (MQ=255)
     aTGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGt                                                                                                        >  1:101470/1‑100 (MQ=255)
     aTGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGt                                                                                                        <  2:2643177/100‑1 (MQ=255)
     aTGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAATGCGCTTTGTTTATGCCGGATGCGGCGt                                                                                                        >  2:2279716/1‑100 (MQ=255)
      tGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGtt                                                                                                       <  2:1470437/100‑1 (MQ=255)
      tGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGtt                                                                                                       <  2:684624/100‑1 (MQ=255)
       gAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTa                                                                                                      >  1:1489292/1‑100 (MQ=255)
        aaGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTACGCCGCATccggcataaa                                                                                                     <  1:865209/100‑6 (MQ=38)
            aTTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGcc                                                                                                 <  1:2448528/100‑1 (MQ=255)
               gctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTa                                                                                              <  1:1709789/100‑1 (MQ=255)
                ctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTAt                                                                                             <  2:494654/100‑1 (MQ=255)
                     ggATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGc                                                                                        >  2:1824644/1‑100 (MQ=255)
                     ggATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGc                                                                                        <  1:3293593/100‑1 (MQ=37)
                       aTGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCt                                                                                      <  2:1683930/100‑1 (MQ=255)
                       aTGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGATTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCt                                                                                      <  2:1992241/100‑1 (MQ=255)
                        tGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGCAT‑‑a                                                                                   <  1:2517435/100‑1 (MQ=25)
                         ggCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑Ac                                                                                  >  1:161693/1‑100 (MQ=255)
                          gcgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑Ac                                                                                  >  2:3460089/1‑99 (MQ=255)
                          gcgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑a                                                                                <  2:1529982/100‑1 (MQ=255)
                           cgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGcataaac                                                                                  <  2:3415184/100‑7 (MQ=25)
                            gcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑aaa                                                                              >  1:2248387/1‑100 (MQ=255)
                               gCACGCAGAGGCACAGCGCCATCGGGTAGACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGcataaacaaag                                                                              <  2:2831839/100‑11 (MQ=18)
                                cACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCg                                                                          <  1:2295274/100‑1 (MQ=255)
                                  cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTg                                                                        <  1:962469/100‑1 (MQ=255)
                                  cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTg                                                                        >  2:2162169/1‑100 (MQ=255)
                                  cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTg                                                                        <  1:1370036/100‑1 (MQ=255)
                                   gCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGc                                                                       >  1:1874661/1‑100 (MQ=255)
                                   gCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGc                                                                       <  1:1974731/100‑1 (MQ=255)
                                   gCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGc                                                                       <  2:3582598/100‑1 (MQ=255)
                                      gagGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCaaa                                                                    <  2:1811987/100‑1 (MQ=255)
                                       agGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAAt                                                                   <  1:628829/100‑1 (MQ=255)
                                          cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTc                                                                 <  1:3035695/99‑1 (MQ=255)
                                           acaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCaa                                                               >  1:1794866/1‑100 (MQ=255)
                                           acaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCaa                                                               <  1:938356/100‑1 (MQ=255)
                                            caGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAAt                                                              <  2:310388/100‑1 (MQ=255)
                                            caGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAAt                                                              <  2:671520/100‑1 (MQ=255)
                                             aGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCTGCCT‑‑AC‑AAAATCGTGCAAATTCAATa                                                             <  1:2219017/100‑1 (MQ=255)
                                             aGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATa                                                             >  1:740880/1‑100 (MQ=255)
                                             aGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATa                                                             <  1:338575/100‑1 (MQ=255)
                                              gcgcCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATAc                                                            >  1:2019660/1‑100 (MQ=255)
                                               cgcCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACa                                                           >  2:3128333/1‑100 (MQ=255)
                                               cgcCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACa                                                           <  1:213036/100‑1 (MQ=255)
                                                gcCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACAt                                                          >  1:134892/1‑100 (MQ=255)
                                                    ccGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGc                                                       <  2:269972/99‑1 (MQ=255)
                                                         tataCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑At‑t                                              >  1:2340773/1‑100 (MQ=255)
                                                          ataCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑At‑tt                                             >  1:2943577/1‑100 (MQ=255)
                                                            aCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTc                                           <  1:317356/100‑1 (MQ=255)
                                                             cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑g                                        <  1:3068839/100‑1 (MQ=255)
                                                             cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑g                                        <  1:1469346/100‑1 (MQ=255)
                                                             cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑g                                        >  1:1918879/1‑100 (MQ=255)
                                                                 ttgttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAgg                                    <  1:3513273/100‑1 (MQ=255)
                                                                 gcataaaCAAAGCGCACTTCGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAgg                                    >  2:1160602/7‑100 (MQ=18)
                                                                   gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGcc                                  <  2:740877/100‑1 (MQ=255)
                                                                   gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGcc                                  <  1:3311824/100‑1 (MQ=255)
                                                                    ttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCt                                 <  2:1056185/100‑1 (MQ=255)
                                                                     tgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTg                                >  1:355915/1‑100 (MQ=255)
                                                                      gACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGa                               >  2:1263608/1‑100 (MQ=255)
                                                                           aGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGAA‑GA‑AT‑TTTC‑‑GTGGGACTGATatg                           >  2:3579623/1‑97 (MQ=17)
                                                                            gTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCg                         >  1:3020061/1‑100 (MQ=255)
                                                                            gTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCg                         >  1:115850/1‑100 (MQ=255)
                                                                              gcgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGa                       <  1:805363/100‑1 (MQ=255)
                                                                                gcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAGTCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAg                      <  2:1587848/99‑1 (MQ=255)
                                                                                gcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAgc                     <  2:617577/100‑1 (MQ=255)
                                                                                 tccgatctaggCCGGATAAGGCGTTTACGCCGCATCCGGCCT‑‑AC‑AAATTC‑TGCAAATTCAACAAATTGC‑‑GA‑ATCCCTT‑‑GTAGGCCTGATAAGCGTAgcgc                   <  1:1140327/90‑1 (MQ=18)
                                                                                 gaatttgtaggCCGGATAAGGCGTTTACGCCGCATCTGGCCT‑‑AC‑AAATTCTTCCAAATTCAGTATATTGCA‑GA‑AA‑TCAT‑‑GTAGGCCTGATAAGCGTAgcg                    >  1:1752813/11‑100 (MQ=12)
                                                                                 cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTTCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAgcg                    <  1:2451639/100‑1 (MQ=255)
                                                                                 cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAgcg                    >  1:446404/1‑100 (MQ=255)
                                                                                 cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAgcg                    >  2:113461/1‑100 (MQ=255)
                                                                                 agttttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCATAAAC‑AAAACTGTGCAAATTCAATATATTGCA‑GG‑AA‑ACAC‑‑GTAGGCCTGATAAGCGAAg                      >  2:754061/11‑100 (MQ=2)
                                                                                  tttGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAgcgc                   >  2:328768/1‑100 (MQ=255)
                                                                                  tttGGTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCTT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAgcg                    <  2:3058806/99‑1 (MQ=255)
                                                                                   ttGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCa                  >  1:1578042/1‑100 (MQ=255)
                                                                                    tttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑ACAAAAATCTTGCCAATTCAATATATTGCA‑GG‑AC‑CAT‑‑‑GTAGGCCTGATAAGTGTAGCGCAt                 >  1:1271158/8‑100 (MQ=18)
                                                                                      tttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCa               >  1:1692890/1‑100 (MQ=255)
                                                                                       gtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCT‑‑AC‑AAATTC‑TGCAAATTCAACAAATTGC‑‑GA‑ATCCCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAgg             >  1:2192925/5‑100 (MQ=18)
                                                                                       gtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCT‑‑AC‑AAAACAGTGCAAATTCAATATATTGCA‑AT‑CT‑TCTG‑‑GTAGGCCTGATAAGCATCGCGCATCAg              >  2:2763245/5‑100 (MQ=18)
                                                                                        taggCCGGATAAGGCGTTTACGCCGCATCCGGCCT‑‑AC‑AAAATCGTGCTAATTCAATATATTGCA‑GA‑AA‑CCTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAgg             >  1:3458623/4‑100 (MQ=18)
                                                                                        taggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑AC‑AAAGTCGTATAAATTCAATGTATTGCA‑GT‑TT‑ATTG‑‑GTAGGCCTGATAAGCGCAGCGCATCAgg             >  1:2786673/4‑100 (MQ=21)
                                                                                         aggCCGGATAAGGCGTTTACGCCGCATCCGGCCT‑‑AC‑AAAACAGTGCAAATTCAATATATTGCA‑AT‑CT‑TCTG‑‑GTAGGCCTGATAAGCATCGCGCATCAGGc            >  1:1833242/3‑100 (MQ=18)
                                                                                         aggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑AC‑AAAACCTTGCTAATTCAATATATTGCAGGG‑AC‑TAT‑‑‑GTAGGCCTGATAAGCATAGCGCATCAGGc            >  1:549715/3‑100 (MQ=18)
                                                                                         aggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑AC‑AAAAACATGCAAATTCAATAGATTGCA‑GAGAT‑TAT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc            >  1:2900952/3‑100 (MQ=17)
                                                                                         aTGCCGGATGCGGCGTGAACGCCTTATCCGGCCT‑‑AC‑AAAGTCGTATAAATTCAATGTATTGCA‑GT‑TT‑ATTG‑‑GTAGGCCTGATGCGCTGCGCTTATCAGGc            <  2:653849/100‑1 (MQ=21)
                                                                                          ggCCGGATAAGGCGTTTACGCCGCATCCGGCCT‑‑AC‑AAAAGCATGCAAATTCAATATATTGCA‑GA‑GT‑TCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa           >  1:2813960/2‑100 (MQ=17)
                                                                                          ggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑ACAAAAATCTTGCCAATTCAATATATTGCA‑GG‑AC‑CAT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa           >  2:2065962/2‑100 (MQ=18)
                                                                                          ggCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑AC‑AAAACTTTGCAAATTCAATATATTGC‑‑‑A‑AT‑CTCCGTGTAGGCCTGATAAGCGTAGCGCATCAGGc            >  2:2463804/2‑99 (MQ=18)
                                                                                           gCCGGATCAGGCGTTAACGCCGCATCCGGCCT‑‑AC‑AAAGTCGAGCAAATTCACTACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGATCGCATCAGGCa           >  2:2439315/1‑99 (MQ=17)
                                                                                           gCCGGATAAGGCGTTTACGCCGCATCCGGCCT‑‑AC‑AAAAGCATGCAAATTCAATATATTGCA‑GA‑GT‑TCAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCag          >  1:1986198/1‑99 (MQ=17)
                                                                                           gCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑AC‑AAGTTCGTGCAAATTCAATAAATTGCA‑‑A‑TA‑TGAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGcgat         >  2:582636/1‑97 (MQ=21)
                                                                                           gCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑AC‑AAAACTTTGCAAATTCAATATATTGCA‑AT‑CT‑CCGT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaa          >  2:2319068/1‑99 (MQ=18)
                                                                                           gCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑AC‑AAAACTGTGCAAATTCAATATATTGCA‑GG‑AA‑ACAC‑‑GTAGGCCTGATAAGCGAAGCGCATCAGGCag          >  1:2153613/1‑99 (MQ=18)
                                                                                           gCCGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑AC‑‑AAATCATGCAAATCCAAAATATTGCA‑GG‑AG‑TTAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAtt         >  1:1666181/1‑100 (MQ=21)
                                                                                           gCCGGATAAGGCGTTAACGCCGCATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACACTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAt          >  1:2869837/1‑100 (MQ=18)
                                                                                            ccGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAtt         >  1:2078441/1‑100 (MQ=255)
                                                                                            ccGGATAAGGCGTTCACGCCGCATCCGGCCT‑‑AC‑AAAAGCATGAAAATTCAATATATTGCA‑GG‑AG‑CTGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGACagt         >  1:3240046/1‑98 (MQ=17)
                                                                                                aTGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg     <  1:1065965/100‑1 (MQ=255)
                                                                                                aTGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg     <  1:2571908/100‑1 (MQ=255)
                                                                                                  gcggcgTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAttttt      >  2:3405976/1‑97 (MQ=255)
                                                                                                  gcggcgTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcn   <  2:135915/100‑2 (MQ=255)
                                                                                                  gcggcgTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcc   <  2:1505022/100‑1 (MQ=255)
                                                                                                   cggcgTTAACAGCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCACTACATTGCA‑GA‑AT‑TTTC‑‑GTAGACCTGATAAGCGAAGCGCATCAGGCATTtttacc   >  2:1882599/1‑99 (MQ=17)
                                                                                                   |                                                                                                           
AGGGAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCT‑‑AC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑TTTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCT  >  NC_000913/3705872‑3706070

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: ‑