Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,706,017 C→T 18.6% intergenic (‑227/+81) 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,706,0170CT18.6% 524.3 / 91.5 213intergenic (‑227/+81)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/100);  new base T (35/3);  total (110/103)
Fisher's exact test for biased strand distribution p-value = 1.25e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.18e-01

Only 100 of 216 total aligned reads displayed.

GCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTT  >  NC_000913/3705918‑3706116
                                                                                                    |                                                                                                      
gcgcCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATAc                                                                                                        <  2:756339/100‑1 (MQ=255)
gcgcCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATAc                                                                                                        <  2:1717408/100‑1 (MQ=255)
      ccGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCa                                                                                                  <  1:1722107/100‑1 (MQ=255)
        gggTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGc                                                                                                   <  1:2540662/97‑1 (MQ=255)
        gggTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGa                                                                                               <  1:2726085/100‑1 (MQ=255)
          gTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGACCTACAAAATCGTGCAAATTCAATACATTGCAGaa                                                                                              >  2:1000886/1‑99 (MQ=255)
             taCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAAtttt                                                                                        >  1:96292/1‑100 (MQ=255)
              aCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTc                                                                                       <  2:2553549/100‑1 (MQ=255)
              aCAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTc                                                                                       >  2:445149/1‑100 (MQ=255)
               cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCg                                                                                      <  1:554843/100‑1 (MQ=255)
                 gATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTa                                                                                    >  1:71821/1‑100 (MQ=255)
                  ggcataaaCAAAGCGCACTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAg                                                                                   >  1:2269458/8‑100 (MQ=18)
                  aTTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAg                                                                                   <  2:56786/100‑1 (MQ=255)
                  aTTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAg                                                                                   <  1:3241177/100‑1 (MQ=255)
                   ttgttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGG                                                                                  <  1:120160‑M1/100‑1 (MQ=255)
                     gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCC                                                                                <  2:3057609‑M1/100‑1 (MQ=255)
                     gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCC                                                                                <  1:2453805‑M1/100‑1 (MQ=255)
                     gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCC                                                                                >  1:2298712‑M1/1‑100 (MQ=255)
                        gACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGA                                                                             >  2:1112880‑M1/1‑100 (MQ=255)
                        gACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTG                                                                              <  1:1419994‑M1/99‑1 (MQ=255)
                          cAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATa                                                                           <  2:496698/100‑1 (MQ=255)
                           aaaGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCtaataa                                                                          >  2:2980698/1‑100 (MQ=255)
                                 cgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAg                                                                    <  1:3028648/100‑1 (MQ=255)
                                 cgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAg                                                                    <  1:2768602/100‑1 (MQ=255)
                                   cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAgcg                                                                  >  2:2536623/1‑100 (MQ=255)
                                    gttttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAACTTTGCAAATTCAATATATTGCAATCTCCGTGTAGGCCTGATAAGCGTAgcgc                                                                 >  1:1381591/10‑100 (MQ=18)
                                     ttGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCa                                                                <  1:2319709‑M2/100‑2 (MQ=255)
                                     atttgtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAATTCTGCAAATTCAACAAATTGCGAATCCCTTGTAGGCCTGATAAGCGTAGCGCAt                                                               >  1:459037/9‑100 (MQ=18)
                                       gTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCatc                                                              <  1:1243170‑M2/100‑4 (MQ=255)
                                        tttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCatca                                                             >  1:2578396‑M2/1‑96 (MQ=255)
                                         ttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCatcag                                                            <  1:2027936‑M2/100‑6 (MQ=255)
                                          taggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAATCATGCAAATCCAAAATATTGCAGGAGTTATGTAGGCCTGATAAGCGTAGCGCATCAGGc                                                          >  2:2064445/4‑100 (MQ=21)
                                           ntgCCGGAAGCGGCGTTAACGCCTTATCCGGCCTACAAAATAGTGCAAATACAATACATTGCAGAATCTCCGTAGGCCGGATAAGCGGAGCGCatcaggc                                                          <  2:1296808‑M2/99‑8 (MQ=255)
                                           aggccGGATAAGGCGTTAACGACGCATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCatcaggc                                                          <  2:1491592‑M2/98‑8 (MQ=255)
                                           aggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTACCGCATCAGGc                                                          >  2:2132410/3‑100 (MQ=18)
                                           aggCCGGATAAGGCGTTCACGCCGCATCCGGCGTACGAAATCGTGCTAATTCAATATATTGTGGGCTCCACGTAGGCCTGATAAGCGTAGCGCATCAGGc                                                          >  1:1827863/3‑100 (MQ=21)
                                           aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAgg                                                           >  1:1066213/3‑99 (MQ=21)
                                           aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAGTCGTATAAATTCAATGTATTGCAGTTTATTGGTAGGCCTGATAAGCGCAGCGCATCAGGc                                                          >  2:1433463/3‑100 (MQ=25)
                                           aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAGCACGCAAACTCAATATATTGCAGAGATCATGTAGGCCTGATAAGCGTAGCGCATCAGGc                                                          >  1:2929041/3‑100 (MQ=14)
                                           aTGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATgcgctccgcttatcaggc                                                          <  1:739240‑M1/100‑19 (MQ=255)
                                           aTGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATgcgctccgcttatcaggc                                                          <  2:96290‑M1/100‑19 (MQ=255)
                                           aTGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATgcgctccgcttatcaggc                                                          <  1:1186655‑M1/100‑19 (MQ=255)
                                           aTGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATgcgctccgcttatcaggc                                                          <  1:1992686‑M1/100‑19 (MQ=255)
                                            tGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCa                                                         >  1:651762/1‑100 (MQ=255)
                                            tGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCa                                                         >  1:346974/1‑100 (MQ=255)
                                            tGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCa                                                         >  2:334523/1‑100 (MQ=255)
                                            tGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCa                                                         >  1:1574772/1‑100 (MQ=255)
                                            tGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCa                                                         >  1:2514116/1‑100 (MQ=255)
                                            ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAATTCTGCAAATTCAACAAATTGCGAATCCCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAt                                                        >  1:758786/2‑100 (MQ=18)
                                            ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                         >  2:628195/2‑100 (MQ=18)
                                            ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAACCTGCAAATTCAATATATTGCAGGATACGTGTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                         >  1:1097111/2‑100 (MQ=17)
                                            ggCCGGATAAGGCGTTAACGCCGCATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCa                                                         >  2:3142817/2‑100 (MQ=18)
                                             gCCGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCag                                                        >  2:2926657/1‑99 (MQ=18)
                                              ccGGATAAGGCGTTTACGCCGCATCCGGCCTACAAAATCGTGCTAATTCAATATATTGCAGAAACCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCagt                                                       >  2:2161661/1‑98 (MQ=18)
                                               cGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATgcgctccgcttatcaggcctac                                                      <  1:2065379‑M1/100‑23 (MQ=255)
                                                ggataaggcgTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCAGGACCATGTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                                                     >  2:3245867/7‑100 (MQ=18)
                                                ggataaggcgTTCACGCCGCATCCGGCCTACAAAACTGTGCAAATTCAATATATTGCAGGAAACACGTAGGCCTGATAAGCGAAGCGCATCAGGCAgttt                                                     >  2:3173169/7‑100 (MQ=18)
                                                ggataaggcgTTCACGCCGCATCCGGCCTACAAATCATGCAAATCCAAAATATTGCAGGAGTTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAtttat                                                    >  1:150825/7‑98 (MQ=21)
                                                ggATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATgcgctccgcttatcaggcctacg                                                     <  1:1738871‑M1/100‑24 (MQ=255)
                                                 gataaggcgTTTACGCCGCATCCGGCCTACAAATTCTGCAAATTCAACAAATTGCGAATCCCTTGTAGGCCTGATAAGCGTAGCGCATCAGGCAtttgcg                                                   >  2:1654229/6‑97 (MQ=21)
                                                 gataaggcgTTCACGCCGCATCCGGCCTACGAAATCGTGCTAATTCAATATATTGTGGGCTCCACGTAGGCCTGATAAGCGTAGCGCATCAGGCAgtttt                                                    <  1:478085/95‑1 (MQ=21)
                                                 gATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGACgcgctccgcttatcaggcctacga                                                    <  1:41579‑M1/100‑25 (MQ=255)
                                                  ataaggcgTTTACGCCTTATCCGGCCTACAAAAGCGTGCAAATTCAATACATTGCATGGGCCATGTAGGCCTGATAAGCGTAGCGCATCAGGCAatttta                                                   >  2:852679/5‑99 (MQ=18)
                                                  aTGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg                                                   <  1:493807/100‑1 (MQ=255)
                                                  aTGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg                                                   >  2:898528/1‑100 (MQ=255)
                                                  aTGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg                                                   <  2:1203661/100‑1 (MQ=255)
                                                  aTGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttg                                                   <  1:2862721/100‑1 (MQ=255)
                                                   tGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgc                                                  <  2:592659/100‑1 (MQ=255)
                                                   tGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgc                                                  >  1:1655618/1‑100 (MQ=255)
                                                   tGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgc                                                  >  1:1544192/1‑100 (MQ=255)
                                                       gcgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATgcgctccgcttatcaggcctacgaaaattc                                              <  1:3235040‑M1/100‑31 (MQ=255)
                                                        cgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttg                                             <  2:540356/100‑1 (MQ=255)
                                                        cgTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttg                                             <  2:138552/100‑1 (MQ=255)
                                                         gTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgc                                            >  2:1570183/1‑100 (MQ=255)
                                                           tAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCa                                          >  2:1455542/1‑100 (MQ=255)
                                                            aaCGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCATAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCa                                          >  1:317301/1‑99 (MQ=255)
                                                              cGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCa                                       >  2:244358/1‑100 (MQ=255)
                                                               gCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAg                                      <  2:1595774/100‑1 (MQ=255)
                                                                ccTTATCCGGCCTACAAAATCGATCAAATACAATAAATTGCAGAATATTCGTAGGCCTGATAAGCGGAGCGCATCAAGCATTTTTGCCTTTGCTAtta                                       >  2:2597756/1‑96 (MQ=18)
                                                                    aTCCGGCCTATAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTg                                 <  1:42761/100‑1 (MQ=255)
                                                                       cGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTAt                              <  2:2269442/100‑1 (MQ=255)
                                                                       cGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTAt                              <  2:1188646/100‑1 (MQ=255)
                                                                        ggCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATg                             <  2:284799/100‑1 (MQ=255)
                                                                         gCCTATAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATgg                            <  2:3126902/100‑1 (MQ=255)
                                                                           cTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCt                          >  1:1621840/1‑100 (MQ=255)
                                                                              cAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTTGGCCTGATAAGCGGAGCGCATCAGGCATTTTTTCTTTTTCTATCTGTCTTGTATCGCtt                         >  2:2004695/1‑98 (MQ=21)
                                                                               ttccgatctGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCtttt                       <  1:2420684/92‑1 (MQ=255)
                                                                                aaaaTCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTa                      <  2:1651247/100‑1 (MQ=255)
                                                                                  aaTCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAat                    <  2:684105/100‑1 (MQ=255)
                                                                                     cGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAattat                 <  2:176973/100‑1 (MQ=255)
                                                                                       tGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTc               >  2:1214875/1‑100 (MQ=255)
                                                                                         cAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTTTTTTTTTTTTCAATCTTGTAAGGCTTTTAATTATTCGa             >  2:798048/1‑100 (MQ=21)
                                                                                          aaaTTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGa             <  1:877622/99‑1 (MQ=255)
                                                                                          aaaTTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGAt            <  2:107808/100‑1 (MQ=255)
                                                                                           aaTTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATa           <  2:877360/100‑1 (MQ=255)
                                                                                            aTTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATag          <  1:1132154/100‑1 (MQ=255)
                                                                                             ttCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATaga         >  1:892122/1‑100 (MQ=255)
                                                                                               cAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAgaga       <  1:234274/100‑1 (MQ=255)
                                                                                                aaTACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCTATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGAc      >  2:1836600/1‑100 (MQ=255)
                                                                                                 aTACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACg     >  1:3094236/1‑100 (MQ=255)
                                                                                                    |                                                                                                      
GCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAATCGTGCAAATTCAATACATTGCAGAATTTTCGTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGTTT  >  NC_000913/3705918‑3706116

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: