Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,705,970 G→A 24.6% intergenic (‑180/+128) 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,9700GA24.6% 523.3 / 162.2 262intergenic (‑180/+128)dppB/dppAdipeptide/heme ABC transporter permease/dipeptide/heme ABC transporter periplasmic binding protein; dipeptide chemotaxis receptor
Reads supporting (aligned to +/- strand):  ref base G (82/114);  new base A (56/8);  total (140/122)
Fisher's exact test for biased strand distribution p-value = 5.14e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.06e-03

Only 100 of 277 total aligned reads displayed.

CAGGGAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCC  >  NC_000913/3705871‑3706069
                                                                                                   |                                                                                                       
cAGGGAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATg                                                                                                         >  2:1925872/1‑100 (MQ=255)
     aaTGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATgcggcg                                                                                                    <  2:1468712/100‑1 (MQ=255)
     aaTGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATgcggcg                                                                                                    <  2:471075/100‑1 (MQ=255)
      aTGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATgcggcg                                                                                                    <  2:2700613/99‑1 (MQ=255)
        gAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTa                                                                                                 >  1:1009198/1‑100 (MQ=255)
         aaGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTaa                                                                                                >  1:4040551/1‑100 (MQ=255)
           gTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACg                                                                                              >  1:3405010/1‑100 (MQ=255)
             aTTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGcc                                                                                            >  1:2534256/1‑100 (MQ=255)
             aTTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGcc                                                                                            <  1:825032/100‑1 (MQ=255)
              ttGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCt                                                                                           >  1:1601513/1‑100 (MQ=255)
                gctgcCGGATGGCGTTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTa                                                                                         <  2:3423004/100‑1 (MQ=255)
                 ctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTAt                                                                                        <  2:1636553/100‑1 (MQ=255)
                    ccGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCCTATcc                                                                                      <  2:3943944/99‑1 (MQ=255)
                      ggATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGc                                                                                   <  1:3481210/100‑1 (MQ=37)
                       gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGcc                                                                                  <  2:1259998/100‑1 (MQ=255)
                        aTGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCt                                                                                 <  2:3340254/100‑1 (MQ=255)
                        aTGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCt                                                                                 <  2:2824479/100‑1 (MQ=255)
                        aTGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCt                                                                                 <  2:2038150/100‑1 (MQ=255)
                        aTGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTAGACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGcat                                                                                 <  1:797060/100‑3 (MQ=25)
                         tGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGCATa                                                                                <  2:307625/100‑1 (MQ=37)
                          ggCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAc                                                                               <  2:1716758/100‑1 (MQ=255)
                           gcgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACa                                                                             >  1:2450178/1‑100 (MQ=255)
                            cgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACaa                                                                            <  1:656862/100‑1 (MQ=255)
                                 cACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCg                                                                       <  1:1055079/100‑1 (MQ=255)
                                 cACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCg                                                                       >  1:2077268/1‑100 (MQ=255)
                                  aCGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGt                                                                      >  1:389396/1‑100 (MQ=255)
                                   cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTg                                                                     <  1:964967/100‑1 (MQ=255)
                                   cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTg                                                                     <  1:3290523/100‑1 (MQ=255)
                                    gCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGc                                                                    <  2:2258673/100‑1 (MQ=255)
                                     cAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCa                                                                   <  1:1124989/100‑1 (MQ=255)
                                       gagGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCaaa                                                                 <  2:3053229/100‑1 (MQ=255)
                                       gagGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCaaa                                                                 >  1:476592/1‑100 (MQ=255)
                                        agGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAAt                                                                <  2:924853/100‑1 (MQ=255)
                                         ggCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCCTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACCAAATCGTGCAAAtt                                                               >  2:2773549/1‑100 (MQ=255)
                                           cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCa                                                             >  1:769454/1‑100 (MQ=255)
                                           cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCa                                                             <  2:351687/100‑1 (MQ=255)
                                           cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCa                                                             >  1:435164/1‑100 (MQ=255)
                                             caGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAAt                                                           <  1:2090335/100‑1 (MQ=255)
                                              aGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATa                                                          >  2:3576802/1‑100 (MQ=255)
                                              aGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATa                                                          <  1:2835950/100‑1 (MQ=255)
                                                  ccACCGGGTCTACACATTGTTGCCAACGTGCGCTTTGTTTATGCCAGATGCGGCGTTAACGCCTTATCCGGCCTACAAACTCGTCCAAATTCAATACAtt                                                      >  2:3822937/1‑100 (MQ=17)
                                                  ccACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACAtt                                                      <  2:3299980/100‑1 (MQ=255)
                                                    aCCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGc                                                    <  1:3576825/100‑1 (MQ=255)
                                                     ccGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCa                                                   >  2:1830516/1‑100 (MQ=255)
                                                             aCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTc                                        <  1:2229846/100‑1 (MQ=255)
                                                             aCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTc                                        <  2:2077562/100‑1 (MQ=255)
                                                             aCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTc                                        <  1:140372/100‑1 (MQ=255)
                                                              cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCg                                       <  1:2975676/100‑1 (MQ=255)
                                                                gATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTa                                     <  2:699352/100‑1 (MQ=255)
                                                                    gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGcc                                 >  2:301001/1‑100 (MQ=255)
                                                                    gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGAcc                                 <  1:2827052/100‑1 (MQ=255)
                                                                    gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAAGcc                                 >  2:682579/1‑100 (MQ=255)
                                                                     ttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCAGGCCTACAAAATCGTGCAAATTCAAGACATTGAAGAATTGTCGtgtg                                   >  2:3701934/1‑94 (MQ=14)
                                                                      tgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCTGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTg                               <  1:274349/100‑1 (MQ=255)
                                                                      tgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGcc                                 >  2:2669298/1‑98 (MQ=255)
                                                                        aCAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGAt                             <  2:2601955/100‑1 (MQ=255)
                                                                        aCAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGAt                             <  1:1830529/100‑1 (MQ=255)
                                                                         cAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATa                            <  1:334005/100‑1 (MQ=255)
                                                                           aaGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAg                          >  2:2634767/1‑100 (MQ=255)
                                                                           aaGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAg                          <  1:2299190/100‑1 (MQ=255)
                                                                            aGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGc                         >  1:3826771/1‑100 (MQ=255)
                                                                            aGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGc                         >  1:3009585/1‑100 (MQ=255)
                                                                             gTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCg                        <  2:955020/100‑1 (MQ=255)
                                                                              tGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCgg                       <  1:3476098/100‑1 (MQ=255)
                                                                                 gcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAgc                    >  2:1788527/1‑100 (MQ=255)
                                                                                 gcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAgc                    <  1:3691026/100‑1 (MQ=255)
                                                                                   tttttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCATACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAgcgc                  >  1:2555997/10‑100 (MQ=17)
                                                                                   attttgtaggCCGGATAAGGCGTTAACGCCGCATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAgcgc                  >  2:2437288/10‑100 (MQ=17)
                                                                                     tGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCAt                >  1:989729/1‑100 (MQ=255)
                                                                                      gTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATc               >  2:768741/1‑100 (MQ=255)
                                                                                       tgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCAGGACCATGTAGGCCTGATAAGCGTAGCGCATCa              >  2:586815/6‑100 (MQ=18)
                                                                                        ttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAg             <  1:4025592/100‑1 (MQ=255)
                                                                                        gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAg             >  1:619150/5‑100 (MQ=17)
                                                                                         taggCCGGATGCGGCGTTAACGCCTTATCCGGTCTACGAAATCGGGCAAAGTCAATACATTGCAGAAGTTTCGTAGGCCTGATAAGCGGAGCGCATc               <  2:3633767/94‑1 (MQ=17)
                                                                                         taggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCAGGACCATGTAGGCCTGATAAGCGTAGCGCATCAgg            >  1:1288580/4‑100 (MQ=18)
                                                                                         tATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTACGTAGGCCTGATAAGCGGAGCGCATCAg             >  2:2564800/1‑99 (MQ=255)
                                                                                          aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCAGGACCATGTAGGCCTGATAAGCGTAGCGCATCAGGc           >  1:3043500/3‑100 (MQ=18)
                                                                                          aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAACCTGCAAATTCAATATATTGCAGGATACGTGTAGGCCTGATAAGCGTAGCGCATCAGGc           >  1:3217352/3‑100 (MQ=18)
                                                                                          aggCCGGATAAGGCATTCACGCCGCATCCGGCCTACAAAACCTTGCTAATTCAATATATTGCAGGGACTATGTAGGCCTGATAAGCATAGCGCATCAGGc           >  1:1875500/3‑100 (MQ=18)
                                                                                          aTGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGc           <  2:4041211/100‑1 (MQ=255)
                                                                                           ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAATTCTGCAAATTCAACAAATTGCGAATCCCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAt         >  1:1985800/2‑100 (MQ=18)
                                                                                           ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCAGGACCATGTAGGCCTGATAAGCGTAGCGCATCAGGCa          >  1:1277254/2‑100 (MQ=18)
                                                                                           ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAGTTCGTGCAAATTCAATAAATTGCAATATGACGTAGGCCTGATAAGCGTAGCGCATCAGGcga         >  1:3767428/2‑98 (MQ=21)
                                                                                           ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCa          >  1:1462990/2‑100 (MQ=21)
                                                                                           ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCa          >  2:1411145/2‑100 (MQ=25)
                                                                                           ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAATCGCTTAAATTCAATATATTGCACCATCCTCGTAGGCCTGATAAGCGTAGCGCATCAGGCa          >  2:3176860/2‑100 (MQ=21)
                                                                                           ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAACCTGCAAATTCAATATATTGCAGGATACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCa          >  2:4063685/2‑100 (MQ=18)
                                                                                           ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCa          >  2:2372723/2‑100 (MQ=17)
                                                                                           ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAATCATGCAAATCCAAAATATTGCAGGAGTTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAt         >  2:2929945/2‑100 (MQ=21)
                                                                                           ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAATCATGCAAATCCAAAATATTGCAGGAGTTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAt         >  1:1210486/2‑100 (MQ=21)
                                                                                            gCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAAGCATGCAAATTCAATATATTGCAGAGTTCACGTAGGCCTGATAAGCGTAGCGCATCAGGCag         >  2:1055938/1‑99 (MQ=17)
                                                                                            gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCAGGACCATGTAGGCCTGATAAGCGTAGCGCATCAGGCag         >  1:2663455/1‑99 (MQ=18)
                                                                                            gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCAGGACCATGTAGGCCTGATAAGCGTAGCGCATCAGGCag         >  1:1511559/1‑99 (MQ=18)
                                                                                            gCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGCaa         >  2:3351670/1‑99 (MQ=21)
                                                                                            gCCGGATAAGGCGTTAACGCCGCATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCAt         >  1:2356037/1‑100 (MQ=18)
                                                                                             ccGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCGTGCAAATTCAATACATTGGAGAATATTCGTAGGCCTGATAAGCGGAGCGCATCAGGCAtt        >  2:1140797/1‑100 (MQ=255)
                                                                                               ggATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCAtttt      >  1:878592/1‑100 (MQ=255)
                                                                                                 acgcggcgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg    <  2:2505781/98‑1 (MQ=255)
                                                                                                 aTGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg    <  1:2634783/100‑1 (MQ=255)
                                                                                                   gcggcgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcc  <  1:2381726/100‑1 (MQ=255)
                                                                                                   |                                                                                                       
CAGGGAATGAAGTATTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCC  >  NC_000913/3705871‑3706069

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: