Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,519,367 C→T intergenic (‑303/+98) aroK ← / ← hofQ shikimate kinase I/DNA catabolic putative fimbrial transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,519,3670CT99.3% 483.0 / ‑5.0 140intergenic (‑303/+98)aroK/hofQshikimate kinase I/DNA catabolic putative fimbrial transporter
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (58/81);  total (58/82)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

Only 100 of 140 total aligned reads displayed.

CTGGCGATAAATCACCGGGTAAGATTAGCGTAAAAAAGACAGCAAAATGCCGCCTGAATGATAAATCATCATCATGGAACGCAGGCAGCGAAGAGCACCAACTCTCAAATCGGTACTCCTTGTATGCTAAATACCTGCGCGCGTCAAATAGATGAAACACGTTCAGCGTAAAAACATTACCGCTTTTACGGC  >  NC_000913/3519270‑3519461
                                                                                                 |                                                                                              
ctGGCGATAAATCACCGGGTAAGATTAGCGTAAAAAAGACAGCAAAATGCCGCCTGAATGATAAATCATCATCATGGAACGCAGGCAGCGAAGAGCATCa                                                                                              <  1:1959812/100‑1 (MQ=255)
ctGGCGATAAATCACCGGGTAAGATTAGCGTAAAAAAGACAGCAAAATGCCGCCTGAATGATAAATCATCATCATGGAACGCAGGCAGCGAAGAGCATCa                                                                                              <  1:1864052/100‑1 (MQ=255)
ctGGCGATAAATCACCGGGTAAGATTAGCGTAAAAAAGACAGCAAAATGCCGCCTGAATGATAAATCATCATCATGGAACGCAGGCAGCGAAGAGCATCa                                                                                              <  1:2808831/100‑1 (MQ=255)
ctGGCGATAAATCACCGGGTAAGATTAGCGTAAAAAAGACAGCAAAATGCCGCCTGAATGATAAATCATCATCATGGAACGCAGGCAGCGAAGAGCATCa                                                                                              <  1:910485/100‑1 (MQ=255)
ctGGCGATAAATCACCGGGTAAGATTAGCGTAAAAAAGACAGCAAAATGCCGCCTGAATGATAAATCATCATCATGGAACGCAGGCAGCGAAGAGCATCa                                                                                              <  1:3155237/100‑1 (MQ=255)
  ggCGATAAATCACCGGGTAAGATTAGCGTAAAAAAGACAGCAAAATGCCGCCTGAATGATAAATCATCATCATGGAACGCAGGCAGCGAAGAGCATCAAc                                                                                            >  1:2366892/1‑100 (MQ=255)
  ggCGATAAATCACCGGGTAAGATTAGCGTAAAAAAGACAGCAAAATGCCGCCTGAATGATAAATCATCATCATGGAACGCAGGCAGCGAAGAGCATCAAc                                                                                            <  1:1911336/100‑1 (MQ=255)
  ggCGATAAATCACCGGGTAAGATTAGCGTAAAAAAGACAGCAAAATGCCGCCTGAATGATAAATCATCATCATGGAACGCAGGCAGCGAAGAGCATCAAc                                                                                            >  2:1849745/1‑100 (MQ=255)
      aTAAATCACCGGGTAAGATTAGCGTAAAAAAGACAGCAAAATGCCGCCTGAATGATAAATCATCATCATGGAACGCAGGCAGCGAAGAGCATCAACtctc                                                                                        <  1:1724721/100‑1 (MQ=255)
         aaTCACCGGGTAAGATTAGCGTAAAAAAGACAGCAAAATGCCGCCTGAATGATAAATCATCATCATGGAACGCAGGCAGCGAAGAGCATCAACTCTCaaa                                                                                     <  1:382934/100‑1 (MQ=255)
         aaTCACCGGGTAAGATTAGCGTAAAAAAGACAGCAAAATGCCGCCTGAATGATAAATCATCATCATGGAACGCAGGCAGCGAAGAGCATCAACTCTCaaa                                                                                     >  2:2222418/1‑100 (MQ=255)
          aTCACCGGGTAAGATTAGCGTAAAAAAGACAGCAAAATGCCGCCTGAATGATAAATCATCATCATGGAACGCAGGCAGCGAAGAGCATCAACTCTCAAAt                                                                                    >  1:3087492/1‑100 (MQ=255)
          aTCACCGGGTAAGATTAGCGTAAAAAAGACAGCAAAATGCCGCCTGAATGATAAATCATCATCATGGAACGCAGGCAGCGAAGAGCATCAACTCTCAAAt                                                                                    >  2:1175056/1‑100 (MQ=255)
          aTCACCGGGTAAGATTAGCGTAAAAAAGACAGCAAAATGCCGCCTGAATGATAAATCATCATCATGGAACGCAGGCAGCGAAGAGCATCAACTCTCAAAt                                                                                    >  2:2491363/1‑100 (MQ=255)
          aTCACCGGGTAAGATTAGCGTAAAAAAGACAGCAAAATGCCGCCTGAATGATAAATCATCATCATGGAACGCAGGCAGCGAAGAGCATCAACTCTCAAAt                                                                                    >  1:2644860/1‑100 (MQ=255)
            cACCGGGTAAGATTAGCGTAAAAAAGACAGCAAAATGCCGCCTGAATGATAAATCATCATCATGGAACGCAGGCAGCGAAGAGCATCAACTCTCAAATCg                                                                                  >  1:1598468/1‑100 (MQ=255)
             aCCGGGTAAGATTAGCGTAAAAAAGACAGCAAAATGCCGCCTGAATGATAAATCATCATCATGGAACGCAGGCAGCGAAGAGCATCAACTCTCAAATCgg                                                                                 <  1:805794/100‑1 (MQ=255)
             aCCGGGTAAGATTAGCGTAAAAAAGACAGCAAAATGCCGCCTGAATGATAAATCATCATCATGGAACGCAGGCAGCGAAGAGCATCAACTCTCAAATCgg                                                                                 <  1:995885/100‑1 (MQ=255)
              ccGGGTAAGATTAGCGTAAAAAAGACAGCAAAATGCCGCCTGAATGATAAATCATCATCATGGAACGCAGGCAGCGAAGAGCATCAACTCTCAAATCGGt                                                                                <  2:890805/100‑1 (MQ=255)
              ccGGGTAAGATTAGCGTAAAAAAGACAGCAAAATGCCGCCTGAATGATAAATCATCATCATGGAACGCAGGCAGCGAAGAGCATCAACTCTCAAATCGGt                                                                                >  1:2542128/1‑100 (MQ=255)
              ccGGGTAAGATTAGCGTAAAAAAGACAGCAAAATGCCGCCTGAATGATAAATCATCATCATGGAACGCAGGCAGCGAAGAGCATCAACTCTCAAATCGGt                                                                                <  1:1065704/100‑1 (MQ=255)
              ccGGGTAAGATTAGCGTAAAAAAGACAGCAAAATGCCGCCTGAATGATAAATCATCATCATGGAACGCAGGCAGCGAAGAGCATCAACTCTCAAATCGGt                                                                                <  1:2019099/100‑1 (MQ=255)
               cGGGTAAGATTAGCGTAAAAAAGACAGCAAAATGCCGCCTGAATGATAAATCATCATCATGGAACGCAGGCAGCGAAGAGCATCAACTCTCAAATCGGTa                                                                               <  2:258734/100‑1 (MQ=255)
                gggTAAGATTAGCGTAAAAAAGACAGCAAAATGCCGCCTGAATGATAAATCATCATCATGGAACGCAGGCAGCGAAGAGCATCAACTCTCAAATCGGTAc                                                                              <  1:2932752/100‑1 (MQ=255)
                    aaGATTAGCGTAAAAAAGACAGCAAAATGCCGCCTGAATGATAAATCATCATCATGGAACGCAGGCAGCGAAGAGCATCAACTCTCAAATCGGTACTCCt                                                                          <  2:140168/100‑1 (MQ=255)
                     aGATTAGCGTAAAAAAGACAGCAAAATGCCGCCTGAATGATAAATCATCATCATGGAACGCAGGCAGCGAAGAGCATCAACTCTCAAATCGGTACTCCtt                                                                         <  1:2911931/100‑1 (MQ=255)
                       aTTAGCGTAAAAAAGACAGCAAAATGCCGCCTGAATGATAAATCATCATCATGGAACGCAGGCAGCGAAGAGCATCAACTCTCAAATCGGTACTCCTTGt                                                                       <  2:2846727/100‑1 (MQ=255)
                        ttAGCGTAAAAAAGACAGCAAAATGCCGCCTGAATGATAAATCATCATCATGGAACGCAGGCAGCGAAGAGCATCAACTCTCAAATCGGTACTCCTTGTa                                                                      <  1:2289828/100‑1 (MQ=255)
                          aGCGTAAAAAAGACAGCAAAATGCCGCCTGAATGATAAATCATCATCATGGAACGCAGGCAGCGAAGAGCATCAACTCTCAAATCGGTACTCCTTGTATg                                                                    <  2:2399327/100‑1 (MQ=255)
                          aGCGTAAAAAAGACAGCAAAATGCCGCCTGAATGATAAATCATCATCATGGAACGCAGGCAGCGAAGAGCATCAACTCTCAAATCGGTACTCCTTGTATg                                                                    <  2:3045852/100‑1 (MQ=255)
                           gCGTAAAAAAGACAGCAAAATGCCGCCTGAATGATAAATCATCATCATGGAACGCAGGCAGCGAAGAGCATCAACTCTCAAATCGGTACTCCTTGTATGc                                                                   >  2:2636090/1‑100 (MQ=255)
                           gCGTAAAAAAGACAGCAAAATGCCGCCTGAATGATAAATCATCATCATGGAACGCAGGCAGCGAAGAGCATCAACTCTCAAATCGGTACTCCTTGTATGc                                                                   >  2:711810/1‑100 (MQ=255)
                           gCGTAAAAAAGACAGCAAAATGCCGCCTGAATGATAAATCATCATCATGGAACGCAGGCAGCGAAGAGCATCAACTCTCAAATCGGTACTCCTTGTATGc                                                                   <  2:3180571/100‑1 (MQ=255)
                            cGTAAAAAAGACAGCAAAATGCCGCCTGAATGATAAATCATCATCATGGAACGCAGGCAGCGAAGAGCATCAACTCTCAAATCGGTACTCCTTGTATGCt                                                                  >  1:3066646/1‑100 (MQ=255)
                               aaaaaaGACAGCAAAATGCCGCCTGAATGATAAATCATCATCATGGAACGCAGGCAGCGAAGAGCATCAACTCTCAAATCGGTACTCCTTGTATGCTaaa                                                               >  1:76379/1‑100 (MQ=255)
                               aaaaaaGACAGCAAAATGCCGCCTGAATGATAAATCATCATCATGGAACGCAGGCAGCGAAGAGCATCAACTCTCAAATCGGTACTCCTTGTATGCTaaa                                                               >  1:1784551/1‑100 (MQ=255)
                               aaaaaaGACAGCAAAATGCCGCCTGAATGATAAATCATCATCATGGAACGCAGGCAGCGAAGAGCATCAACTCTCAAATCGGTACTCCTTGTATGCTaaa                                                               >  2:565299/1‑100 (MQ=255)
                               aaaaaaGACAGCAAAATGCCGCCTGAATGATAAATCATCATCATGGAACGCAGGCAGCGAAGAGCATCAACTCTCAAATCGGTACTCCTTGTATGCTaaa                                                               >  2:2450614/1‑100 (MQ=255)
                                 aaaaGACAGCAAAATGCCGCCTGAATGATAAATCATCATCATGGAACGCAGGCAGCGAAGAGCATCAACTCTCAAATCGGTACTCCTTGTATGCTAAATa                                                             <  2:2608961/100‑1 (MQ=255)
                                  aaaGACAGCAAAATGCCGCCTGAATGATAAATCATCATCATGGAACGCAGGCAGCGAAGAGCATCAACTCTCAAATCGGTACTCCTTGTATGCTAAATAc                                                            <  1:2874629/100‑1 (MQ=255)
                                  aaaGACAGCAAAATGCCGCCTGAATGATAAATCATCATCATGGAACGCAGGCAGCGAAGAGCATCAACTCTCAAATCGGTACTCCTTGTATGCTAAATAc                                                            <  1:573780/100‑1 (MQ=255)
                                   aaGACAGCAAAATGCCGCCTGAATGATAAATCATCATCATGGAACGCAGGCAGCGAAGAGCATCAACTCTCAAATCGGTACTCCTTGTATGCTAAATAcc                                                           <  2:560263/100‑1 (MQ=255)
                                    aGACAGCAAAATGCCGCCTGAATGATAAATCATCATCATGGAACGCAGGCAGCGAAGAGCATCAACTCTCAAATCGGTACTCCTTGTATGCTAAATACCt                                                          >  2:213382/1‑100 (MQ=255)
                                     gACAGCAAAATGCCGCCTGAATGATAAATCATCATCATGGACCGCAGGCAGCGAAGAGCATCAACTCTCAAATCGGTACTCCTTGTATGCTAAATACCTg                                                         <  2:2053582/100‑1 (MQ=255)
                                       cAGCAAAATGCCGCCTGAATGATAAATCATCATCATGGAACGCAGGCAGCGAAGAGCATCAACTCTCAAATCGGTACTCCTTGTATGCTAAATACCTgcg                                                       <  2:1963891/100‑1 (MQ=255)
                                       cAGCAAAATGCCGCCTGAATGATAAATCATCATCATGGAACGCAGGCAGCGAAGAGCATCAACTCTCAAATCGGTACTCCTTGTATGCTAAATACCTgcg                                                       >  1:1367502/1‑100 (MQ=255)
                                       cAGCAAAATGCCGCCTGAATGATAAATCATCATCATGGAACGCAGGCAGCGAAGAGCATCAACTCTCAAATCGGTACTCCTTGTATGCTAAATACCTgcg                                                       <  2:328524/100‑1 (MQ=255)
                                       cAGCAAAATGCCGCCTGAATGATAAATCATCATCATGGAACGCAGGCAGCGAAGAGCATCAACTCTCAAATCGGTACTCCTTGTATGCTAAATACCTgcg                                                       >  2:1479999/1‑100 (MQ=255)
                                       cAGCAAAATGCCGCCTGAATGATAAATCATCATCATGGAACGCAGGCAGCGAAGAGCATCAACTCTCAAATCGGTACTCCTTGTATGCTAAATACCTgc                                                        >  2:427062/1‑99 (MQ=255)
                                        aGCAAACTGCCGCCTGAATGATAAATCATCATCATGGAACGCAGGCAGCGAAGAGCATCAACTCTCAAATCGGTACTCCTTGTATGCTAAATACCTgcgc                                                      <  1:142655/100‑1 (MQ=255)
                                         gCAAAATGCCGCCTGAATGATAAATCATCGTCATGGAACGCAGGCAGCGAAGAGCATCAACTCTCAAATCGGTACTCCTTGTATGCTAAATACCTgcgcg                                                     >  2:338722/1‑100 (MQ=255)
                                         gCAAAATGCCGCCTGAATGATAAATCATCATCATGGAACGCAGGCAGCGAAGAGCATCAACTCTCAAATCGGTACTCCTTGTATGCTAAATACCTgcgcg                                                     >  2:2158358/1‑100 (MQ=255)
                                         gCAAAATGCCGCCTGAATGATAAATCATCATCATGGAACGCAGGCAGCGAAGAGCATCAACTCTCAAATCGGTACTCCTTGTATGCTAAATACCTgcgcg                                                     >  2:1881590/1‑100 (MQ=255)
                                         gCAAAATGCCGCCTGAATGATAAATCATCATCATGGAACGCAGGCAGCGAAGAGCATCAACTCTCAAATCGGTACTCCTTGTATGCTAAATACCTgcgcg                                                     <  2:2905515/100‑1 (MQ=255)
                                         gCAAAATGCCGCCTGAATGATAAATCATCATCATGGAACGCAGGCAGCGAAGAGCATCAACTCTCAAATCGGTACTCCTTGTATGCTAAATACCTgcgcg                                                     >  1:2628340/1‑100 (MQ=255)
                                         gCAAAATGCCGCCTGAATGATAAATCATCATCATGGAACGCAGGCAGCGAAGAGCATCAACTCTCAAATCGGTACTCCTTGTATGCTAAATACCTgcgcg                                                     >  1:1524624/1‑100 (MQ=255)
                                         gCAAAATGCCGCCTGAATGATAAATCATCATCATGGAACGCAGGCAGCGAAGAGCATCAACTCTCAAATCGGTACTCCTTGTATGCTAAATACCTgcgcg                                                     >  1:2438193/1‑100 (MQ=255)
                                           aaaaTGCCGCCTGAATGATAAATCATCATCATGGAACGCAGGCAGCGAAGAGCATCAACTCTCAAATCGGTACTCCTTGTATGCTAAATACCTGcgcgcg                                                   >  1:2060827/1‑100 (MQ=255)
                                           aaaaTGCCGCCTGAATGATAAATCATCATCATGGAACGCAGGCAGCGAAGAGCATCAACTCTCAAATCGGTACTCCTTGTATGCTAAATACCTGcgcgcg                                                   <  2:1625357/100‑1 (MQ=255)
                                           aaaaTGCCGCCTGAATGATAAATCATCATCATGGAACGCAGGCAGCGAAGAGCATCAACTCTCAAATCGGTACTCCTTGTATGCTAAATACCTGcgcgcg                                                   >  1:2814817/1‑100 (MQ=255)
                                             aaTGCCGCCTGAATGATAAATCATCATCATGGAACGCAGGCAGCGAAGAGCATCAACTCTCAAATCGGTACTCCTTGTATGCTAAATACCTGCGCGCGTc                                                 <  2:2343371/100‑1 (MQ=255)
                                              aTGCCGCCTGAATGATAAATCATCATCATGGAACGCAGGCAGCGAAGAGCATCAACTCTCAAATCGGTACTCCTTGTATGCTAAATACCTGCGCGCGTCa                                                <  1:3139426/100‑1 (MQ=255)
                                                nccgccTGAATGATAAATCATCATCATGGAACGCAGGCAGCGAAGAGCATCAACTCTCAAATCGGTACTCCTTGTATGCTAAATACCTGCGCGCGTCaaa                                              >  2:3238038/2‑100 (MQ=255)
                                                gccgccTGAATGATAAATCATCATCATGGAACGCAGGCAGCGAAGAGCATCAACTCTCAAATCGGTACTCCTTGTATGCTAAATACCTGCGCGCGTCaaa                                              <  1:958006/100‑1 (MQ=255)
                                                gccgccTGAATGATAAATCATCATCATGGAACGCAGGCAGCGAAGAGCATCAACTCTCAAATCGGTACTCCTTGTATGCTAAATACCTGCGCGCGTCaaa                                              <  1:2478805/100‑1 (MQ=255)
                                                 ccgccTGAATGATAAATCATCATCATGGAACGCAGGCAGCGAAGAGCATCAACTCTCAAAGCGGGACTCCTTGTATGCTAAATACCTGCGCGCGTc                                                 <  2:46672/96‑1 (MQ=255)
                                                  cgccTGAATGATAAATCATCATCATGGAACGCAGGCAGCGAAGAGCATCAACTCTCAAATCGGTACTCCTTGTATGCTACATACCTGCGCGCGTCAAATa                                            <  2:1988190/100‑1 (MQ=255)
                                                   gccTGAATGATAAATCATCATCATGGAACGCAGGCAGCGAAGAGCATCAACTCTCAAATCGGTACTCCTTGTATGCTAAATACCTGCGCGCGTCAAATAg                                           <  2:71244/100‑1 (MQ=255)
                                                     cTGAATGATAAATCATCATCATGGAACGCAGGCAGCGAAGAGCATCAACTCTCAAATCGGTACTCCTTGTATGCTAAATACCTGCGCGCGTCAAATAGAt                                         >  1:912109/1‑100 (MQ=255)
                                                      tGAATGATAAATCATCATCATGGAACGCAGGCAGCGAAGAGCATCAACTCTCAAATCGGTACTCCTTGTATGCTAAATACCTGCGCGCGTCAAATAGATg                                        <  2:2814782/100‑1 (MQ=255)
                                                            aTAAATCATCATCATGGAACGCAGGCAGCGAAGAGCATCAACTCTCAAATCGGTACTCCTTGTATGCTAAATACCTGCGCGCGTCAAATAGATGAAacac                                  <  1:426337/100‑1 (MQ=255)
                                                            aTAAATCATCATCATGGAACGCAGGCAGCGAAGAGCATCAACTCTCAAATCGGTACTCCTTGTATGCTAAATACCTGCGCGCGTCAAATAGATGAAacac                                  >  2:849640/1‑100 (MQ=255)
                                                             tAAATCATCATCATGGAACGCAGGCAGCGAAGAGCATCAACTCTCAAATCGGTACTCCTTGTATGCTAAATACCTGCGCGCGTCAAATAGATGAAACACg                                 <  2:3190023/100‑1 (MQ=255)
                                                             tAAATCATCATCATGGAACGCAGGCAGCGAAGAGCATCAACTCTCAAATCGGTACTCCTTGTATGCTAAATACCTGCGCGCGTCAAATAGATGAAACACg                                 >  1:1359851/1‑100 (MQ=255)
                                                              aaaTCATCATCATGGAACGCAGGCAGCGAAGAGCATCAGCTCTCAAATCGGTACTCCTTGTATGCTAAATACCTGCGCGCGTCAAATAGATGAAACACGt                                <  2:1784529/100‑1 (MQ=255)
                                                               aaTCATCATCATGGAACGCAGGCAGCGAAGAGCATCAACTCTCAAATCGGTACTCCTTGTATGCTAAATACCTGCGCGCGTCAAATAGATGAAACCCGtn                               <  2:812205/100‑2 (MQ=255)
                                                               aaTCATCATCATGGAACGCAGGCAGCGAAGAGCATCAACTCTCAAATCGGTACTCCTTGTATGCTAAATACCTGCGCGCGTCAAATAGATGAAACACGtt                               >  1:548078/1‑100 (MQ=255)
                                                                  catcatcatGGAACGCAGGCAGCGAAGAGCATCAACTCTCAAATCGGTACTCCTTGTATGCTAAATACCTGCGCGCGTCAAATAGATGAAACACGTTCAg                            <  1:1810511/100‑1 (MQ=255)
                                                                  catcatcatGGAACGCAGGCAGCGAAGAGCATCAACTCGCAAAGCGGTACTCCTTGTATGCTAAATACCTGCGCGCGTCAAATAGATGAAACACGTTCAg                            <  1:1089220/100‑1 (MQ=255)
                                                                    tcattatGGAACGCAGGCAGCGAAGAGCATCAACTCTCAAATCGGTACTCCTTGTATGCTAAATACCTGCGCGCGTCAAATAAATGAAACACGTTCAGc                           <  2:1344892/99‑1 (MQ=255)
                                                                      atcatGGAACGCAGGCAGCGAAGAGCATCAACTCTCAAATCGGTACTCCTTGTATGCTAAATACCTGCGCGCGTCAAATAGATGAAACACGTTCAGCGTa                        >  1:95347/1‑100 (MQ=255)
                                                                      atcatGGAACGCAGGCAGCGAAGAGCATCAACTCTCAAATCGGTACTCCTTGTATGCTAAATACCTGCGCGCGTCAAATAGATGAAACACGTTCAGCGTa                        >  1:1603572/1‑100 (MQ=255)
                                                                      atcatGGAACGCAGGCAGCGAAGAGCATCAACTCTCAAATCGGTACTCCTTGTATGCTAAATACCTGCGCGCGTCAAATAGATGAAACACGTTCAGCGTa                        >  1:2458546/1‑100 (MQ=255)
                                                                        catGGAACGCAGGCAGCGAAGAGCATCAACTCTCAAATCGGTACTCCTTGTATGCTAAATACCTGCGCGCGTCAAATAGATGAAACACGTTCAGCGTaaa                      <  1:2491392/100‑1 (MQ=255)
                                                                        catGGAACGCAGGCAGCGAAGAGCATCAACTCTCAAATCGGTACTCCTTGTATGCTAAATACCTGCGCGCGTCAAATAGATGAAACACGTTCAGCGTaaa                      <  1:3111384/100‑1 (MQ=255)
                                                                         atGGAACGCAGGCAGCGAAGAGCATCAACTCTCAAATCGGTACTCCTTGTATGCTAAATACCTGCGCGCGTCAAATAGATGAAACACGTTCAGCGTaaaa                     >  1:2315631/1‑100 (MQ=255)
                                                                         atGGAACGCAGGCAGCGAAGAGCATCAACTCTCAAATCGGTACTCCTTGTATGCTAAATACCTGCGCGCGTCAAATAGATGAAACACGTTCAGCGTaaaa                     <  2:178551/100‑1 (MQ=255)
                                                                         atGGAACGCAGGCAGCGAAGAGCATCAACTCTCAAATCGGTACTCCTTGTATGCTAAATACCTGCGCGCGTCAAATAGATGAAACACGTTCAGCGTaaaa                     >  1:1407477/1‑100 (MQ=255)
                                                                            gAACGCAGGCAGCGAAGAGCATCAACTCTCAAATCGGTACTCCTTGTATGCTAAATACCTGCGCGCGTCAAATAGATGAAACACGTTCAGCGTAAAAACa                  <  1:2274732/100‑1 (MQ=255)
                                                                              aCGCAGGCAGCGAAGAGCATCAACTCTCAAATCGGTACTCCTTGTATGCTAAATACCTGCGCGCGTCAAATAGATGAAACACGTTCAGCGTAAAAACAtt                <  1:2348078/100‑1 (MQ=255)
                                                                              aCGCAGGCAGCGAAGAGCATCAACTCTCAAATCGGTACTCCTTGTATGCTAAATACCTGCGCGCGTCAAATAGATGAAACACGTTCAGCGTAAAAACAtt                <  1:2611229/100‑1 (MQ=255)
                                                                                  aggcagCGAAGAGCATCAACTCTCAAATCGGTACTCCTTGTATGCTAAATACCTGCGCGCGTCAAATAGATGAAACACGTTCAGCGTAAAAACATTACCg            >  1:69871/1‑100 (MQ=255)
                                                                                    gcagCGAAGAGCATCAACTCTCAAATCGGTACTCCTTGTATGCTAAATACCTGCGCGCGTCAAATAGATGAAACACGTTCAGCGTAAAAACATTACCGc           >  2:2276140/1‑99 (MQ=255)
                                                                                    gcagCGAAGAGCATCAACTCTCAAATCGGTACTCCTTGTATGCTAAATACCTGCGCGCGTCAAATAGATGAAACACGTTCAGCGTAAAAACATTACCGCt          <  1:519908/100‑1 (MQ=255)
                                                                                       gCGAAGAGCATCAACTCTCAAATCGGTACTCCTTGTATGCTAAATACCTGCGCGCGTCAAATAGATGAAACACGTTCAGCGTAAAAACATTACCGCtttt       <  1:1849767/100‑1 (MQ=255)
                                                                                        cGAAGAGCATCAACTCTCAAATCGGTACTCCTTGTATGCTAAATACCTGCGCGCGTCAAATAGATGAAACACGTTCAGCGTAAAAACATTACCGCtttt       <  2:3095655/99‑1 (MQ=255)
                                                                                        cGAAGAGCATCAACTCTCAAATCGGTACTCCTTGTATGCTAAATACCTGCGCGCGTCAAATAGATGAAACACGTTCAGCGTAAAAACATTACCGCTTTTa      <  1:1141546/100‑1 (MQ=255)
                                                                                         gAAGAGCACCAACTCTCAAATCGGTACTCCTTGTATGCTAAATACCTGCGCGCGTCAAATAGATGAAACACGTTCAGCGTAAAAACATTACCGCTTTTAc     <  2:4133/100‑1 (MQ=255)
                                                                                           agagCATCAACTCTCAAATCGGTACTCCTTGTATGCTAAATACCTGCGCGCGTCAAATAGATGAAACACGTTCAGCGTAAAAACATTACCGCTTTTACgg   <  2:2542097/100‑1 (MQ=255)
                                                                                            gagCATCAACTCTCAAATCGGTACTCCTTGTATGCTAAATACCTGCGCGCGTCAAATAGATGAAACACGTTCAGCGTAAAAACATTACCGCTTTTACGGc  <  1:1698656/100‑1 (MQ=255)
                                                                                                 |                                                                                              
CTGGCGATAAATCACCGGGTAAGATTAGCGTAAAAAAGACAGCAAAATGCCGCCTGAATGATAAATCATCATCATGGAACGCAGGCAGCGAAGAGCACCAACTCTCAAATCGGTACTCCTTGTATGCTAAATACCTGCGCGCGTCAAATAGATGAAACACGTTCAGCGTAAAAACATTACCGCTTTTACGGC  >  NC_000913/3519270‑3519461

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: