Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,370,616 A→G 27.9% intergenic (‑204/‑72) yjeH ← / → groS putative transporter/Cpn10 chaperonin GroES, small subunit of GroESL

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,370,6160AG27.9% 131.5 / 62.0 96intergenic (‑204/‑72)yjeH/groSputative transporter/Cpn10 chaperonin GroES, small subunit of GroESL
Reads supporting (aligned to +/- strand):  ref base A (27/42);  new base G (10/17);  total (37/59)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGAGGTGAATCAGGGTTTTCACCCGATTTTGTGCTGATCAGAATTTTTTTTCTTTTTCCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGTCACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGGAGAGTTATCAATGAATATTCGTCCATTGCATGATC  >  NC_000913/4370518‑4370712
                                                                                                  |                                                                                                
aatgtgaggtgaatcagggttttcacccgattttgtgctgatcagaattttttttctttttcccccttgaaggggcgaagcctcatccccatttGGTCAc                                                                                                 <  2:821579‑M2/6‑1 (MQ=255)
 atgtgaggtgaatcagggttttcacccgattttgtgctgatcagaattttttttctttttcccccttgaaggggcgaagcctcatccccatttGGTCAcc                                                                                                <  2:878975‑M2/7‑1 (MQ=255)
      gAATCAGGGTTTTCACCCGATTTTGTGCTGATCAGAATTTTTTTTCTTTTTCCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGTCGCCAGCCg                                                                                           <  1:1289447/100‑1 (MQ=255)
      gAATCAGGGTTTTCACCCGATTTTGTGCTGATCAGAATTTTTTTTCTTTTTCCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGTCACCAGCCg                                                                                           <  2:302615/100‑1 (MQ=255)
       ggtgaatcagggttttcacccgattttgtgctgatcagaattttttttctttttcccccttgaaggggcgaagcctcatccccatttGGTCACCAGCCgg                                                                                          >  1:1549403‑M2/88‑100 (MQ=255)
       ggtgaatcagggttttcacccgattttgtgctgatcagaattttttttctttttcccccttgaaggggcgaagcctcatccccatttGGTCACCAGCCgg                                                                                          >  2:1307543‑M2/88‑100 (MQ=255)
       aaTCAGGGTTTTCAACCGAGTTTGTGCTGATCAGAATTTTTTTTCTTTTTCCCCCTTGAAGGGGCGAAGCCTCTTCCCCATTTCTCTGGTCGCCAGCCg                                                                                           <  2:5183/99‑1 (MQ=255)
        aTCAGGGTTTTCACCCGATTTTGTGCTGATCAGAATTTTTTTTCTTTTTCCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGTCACCAGCCggg                                                                                         <  2:1091528/100‑1 (MQ=255)
        aTCAGGGTTTTCACCCGATTTTGTGCTGATCAGAATTTTTTTTCTTTTTCCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGTCACCAGCCggg                                                                                         <  2:326562/100‑1 (MQ=255)
         tgaatcagggttttcacccgattttgtgctgatcagaattttttttctttttcccccgtgaaggggcgaagcctcatccccatttGGTCACCAGCCggg                                                                                         <  1:1768952‑M2/14‑1 (MQ=255)
           aGGGTTTTCACCCGATTTTGTGCTGATCAGAATTTTTTTTCTTTCTCCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGTCACCAGCCGGGaaa                                                                                      >  2:893669/1‑100 (MQ=255)
            gggTTTTCACCCGATTTTGTGCTGATCAGAATTTTTTTTCTTTTTCCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGTCACCAGCCGGGAAAc                                                                                     <  2:2189321/100‑1 (MQ=255)
            atcagggttttcacccgattttgtgctgatcagaattttttttctttttcccccttgaaggggcgaagcctcatccccatttGGTCACCAGCCGGGAAAc                                                                                     <  2:194226‑M2/18‑1 (MQ=255)
             ggTTTTCACCCGATTTTGTGCTGATCAGAATTTTTTTTCTTTTTCCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGTCGCCAGCCGGGAAAcc                                                                                    <  1:477005/100‑1 (MQ=255)
               agggttttcacccgattttgtgctgatcagaattttttttctttttcccccttgaaggggcgaagcctcatccccatttGGTCACCAGCCGGGAAACCAc                                                                                  >  2:338802‑M2/80‑100 (MQ=255)
                  gttttctcccgattttgtgctgatcagaattttttttctttttcccccttgaaggggcgaagcctcatccccatttGGTCACCAGCCGGGAAACCACGTa                                                                               <  2:1218128‑M2/24‑1 (MQ=255)
                  gttttcacccgattttgtgctgatcagaattttttttctttttcccccttgaaggggcgaagcctcatccccatttGGTCACCAGCCGGGAAACCACGTa                                                                               <  2:1684998‑M2/24‑1 (MQ=255)
                    aCCCGATTTTGTGCTGATCAGAATTTTTTTTCTTTTTCCCCCTTGAAGGGTATACGCCTCATCCCCATTTC‑‑‑‑GTCAAAAGCCGGAAAACAACGTAATCTcc                                                                         >  1:1161336/1‑100 (MQ=255)
                     cccGATTTTGTGCTGATCAGAATTTTTTTTCTTTTTCCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGTCACCAGCCGGGAAACCACGTAAGc                                                                            <  2:852658/100‑1 (MQ=255)
                        gATTTTGTGCTGATCAGAATTTTTTTTCTTTTTCCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGTCACCAGCCGGGAAACCACGTAAGCTc                                                                          <  1:1954354/99‑1 (MQ=255)
                        acccgattttgtgctgatcagaattttttttctttttcccccttgaaggggcgaagcctcatccccatttGGTCACCAGCCGGGAAACCACGTAAGc                                                                            >  1:2231396‑M2/71‑97 (MQ=255)
                         cccgattttgtgctgatcagaattttttttctttttcccccttgaaggggcgaagcctcatccccatttGGTCACCAGCCGGGAAACCACGTAAGCTCCg                                                                        >  1:1415267‑M2/70‑100 (MQ=255)
                          ccgattttgtgctgatcagaattttttttctttttcccccttgaaggggcgaagcctcatccccatttGGTCACCAGCCGGGAAACCACGTAAGCTCCgg                                                                       <  1:1605138‑M2/32‑1 (MQ=255)
                           cgattttgtgctgatcagaattttttttctttttcccccttgaaggggcgaagcctcatccccatttGGTCACCAGCCGGGAAACCACGTAAGCTCCGGc                                                                      <  2:1689963‑M2/33‑1 (MQ=255)
                            ttGTGCTGATCAGAATTTTTTTTCTTTTTCCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGTCGCCAGCCGGGAAACCACGTAAGCTCCGGCg                                                                     <  1:1864954/100‑1 (MQ=255)
                            ttGTGCTGATCAGAATTTTTTTTCTTTTTCCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGTCACCAGCCGGGAAACCACGTAAGCTCCGGCg                                                                     <  2:2085094/100‑1 (MQ=255)
                            ttGTGCTGATCAGAATTTTTTTTCTTTTTCCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGTCACCAGCCGGGAAACCACGTAAGCTCCGGCg                                                                     <  1:942159/100‑1 (MQ=255)
                            ttGTGCTGATCAGAATTTTTTTTCTTTTTCCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGTCACCAGCCGGGAAACCACGTAAGCTCCGGCg                                                                     <  2:1846226/100‑1 (MQ=255)
                            gattttgtgctgatcagaattttttttctttttcccccttgaaggggcgaagcctcatccccatttGGTCACCAGCCGGGAAACCACGTAAGCTCCGGCg                                                                     <  1:2072886‑M2/34‑1 (MQ=255)
                             tgtgCTGATCAGAATTTTTTTTCTTTTTCCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGTCGCCAGCCGGGAAACCACGTAAGCTCCGGCGt                                                                    >  1:142805/1‑100 (MQ=255)
                               tgCTGATCAGAATTTTTTTTCTTTTTCCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGTCGCCAGCCGGGAAACCACGTAAGCTCCGGCGTCa                                                                  <  1:842723/100‑1 (MQ=255)
                                ttgtgctgatcagaattttttttctttttcccccttgaaggggcgaagcctcatccccatttGGTCACCAGCCGGGAAACCACGTAAGCTCCGGCGTCAc                                                                 >  2:1898305‑M2/63‑100 (MQ=255)
                                gCTGATCAGAATTTTTTTTTTTTTTTTTTCTTTAATAGGCGAAGCATCATCTCCATTTCTCTGTTCACCAGCCGGGGAACCACGTAAGCTCCGGCGTCAc                                                                 >  2:1338936/1‑100 (MQ=255)
                                 tgtgctgatcagaattttttttctttttcccccttgaaggggcgaagcctcatccccatttGGTCACCAGCCGGGAAACCACGTAAGCTCCGGCGTCAc                                                                 <  2:2224295‑M2/38‑1 (MQ=255)
                                 cTGATCAGAATTTTTTTTCTTTTTCCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGTCACCAGCCGGGAAACCACGTAAGCTCCGGCGTCAcc                                                                <  1:492073/100‑1 (MQ=255)
                                 cTGATCAGAATTTTTTTTCTTTTTCCACCTTAAAAGGGCAAAGCCTCATCTCCATTTCTCTGCTCGCCAGCCGGGAAACCACGTAACCTCCGGCTTCAc                                                                 >  2:557825/1‑99 (MQ=255)
                                  tGATCAGAATTTTTTTTCTTTTTCCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGTCGCCAGCCGGGAAACCACGTAAGCTCCGGCGTCAccc                                                               >  2:1287135/1‑100 (MQ=255)
                                   gATCAGAATTTTTTTTCTTTTTCCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGTCGCCAGCCGGGAAACCACGTACGCTCCGGCGCCACCCa                                                              >  2:1305991/1‑100 (MQ=255)
                                   gATCAGAATTTTTTTTCTTTTTCCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGTCGCCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCa                                                              <  2:2201149/100‑1 (MQ=255)
                                    aTCAGAATTTTTTTTTTTTTTTTTTCTTGAAGGGGAGAAGCCTCATCCCCATTTCTCTGGTCACCAGCCGGGAAACCACGTAAGCTCCTGCGTCACCCAt                                                             >  2:1452792/1‑100 (MQ=255)
                                    aTCAGAATTTTTTTTCTTTTTCCTCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGTCACCAGCCGGGAAACCACGTAAGCTCCTGCGTCACCCAt                                                             >  2:1452917/1‑100 (MQ=255)
                                      cAGAATTTTTTTTCTTTTTCCCCCTTGGAGGGGCGAAGCCTCATCCCCATTTCTCTGGTCGCCAGCCGGGAAACCACGTAAGCTCCGGCGTCAccc                                                               <  2:1071479/96‑1 (MQ=255)
                                        gAATTTTTTTTCTTTTTTTTTTTTTAAGTGGCGAAGCCTCATCCCCATTTCTCTGGTCACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATa                                                            >  2:151684/1‑97 (MQ=255)
                                        gAATTTTTTTTCTTTTTCCCCCTTGAAGGGGCGAAGCCTCATCCCTATTTCTCTGGTCACCAGCCGGGAAACCTCGTAAGCTCCGGCGTCACCCATAACa                                                         <  2:68876/100‑1 (MQ=255)
                                        gAATTTTTTTTCTTTTTCCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGTCACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACa                                                         <  1:227400/100‑1 (MQ=255)
                                        aatttttttttCTTTTTCCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGTCGCCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACa                                                         <  1:876719/97‑1 (MQ=255)
                                          cagaattttttttctttttcccccttgaaggggcgaagcctcatccccatttGGTCACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGa                                                       <  1:338801‑M2/48‑1 (MQ=255)
                                           agaagtttttgtattgttcccccttgaaggggcgaagcctcatccccatgtGGTCACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAGCa                                                         <  1:486136‑M2/46‑1 (MQ=255)
                                            tttttttCTTTTTCCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGTCGCCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATa                                                     <  1:1741366/100‑1 (MQ=255)
                                              attttttttctttttcccccttgaaggggcgaagcctcatccccatttGGTCACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACg                                                   >  1:2189170‑M2/49‑100 (MQ=255)
                                                tttctttttcCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGTCGCCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGa                                                 <  1:1715295/100‑1 (MQ=255)
                                                tttctttttcCCACTTGAATGGAAAAACCCACATATATATTAATCTGGTCGCCAGCCGGTGAACCACGTAATCTCAGGCGTCACCCATAACAGAt                                                      >  2:531532/1‑95 (MQ=255)
                                                     ttttcCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGTCACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTtc                                            <  1:678148/100‑1 (MQ=255)
                                                     ttctttttctcccttgaaggggcgaagcctcatccccatttGGTCACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTtc                                            <  1:103051‑M2/59‑1 (MQ=255)
                                                     ttctttttcccccttgaaggggcgaagcctcatccccatttGGTCACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTtc                                            >  1:211252‑M2/42‑100 (MQ=255)
                                                     ttctttttcccccttgaaggggcgaagcctcatccccatttGGTCACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTtc                                            <  1:453028‑M2/59‑1 (MQ=255)
                                                       ttcCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGTCACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTtctc                                          <  1:515360/100‑1 (MQ=255)
                                                        tttttcccccttgaaggggcgaagcctcatccccatttGGTCACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCa                                         <  1:2027822‑M2/62‑1 (MQ=255)
                                                        tcCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGTCGCCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCa                                         <  1:474669/100‑1 (MQ=255)
                                                           ttcccccttgaaggggcgaagcctcatccccatttGGTCACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAg                                      <  1:8911‑M2/65‑1 (MQ=255)
                                                            tcccccttgaaggggcgaagcctcatccccatttGGTCACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAgg                                     >  2:1244835‑M2/35‑100 (MQ=255)
                                                            ccTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGTCACCAGCCGGGAAACCACGTAAGCTCCGGCGCCACCCATAACAGATACGGACTTTCTCAAAgg                                     >  1:860783/1‑100 (MQ=255)
                                                             cccccttgaaggggcgaagcctcctccccatttGGTCACCAGCCGGGAAACCACGTAAGCTCCGGCGGCACCCATAACACATACGGACTTACTCAAAGga                                    <  2:2231078‑M2/67‑1 (MQ=255)
                                                             cTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGTCGCCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAgg                                     <  2:132997/99‑1 (MQ=255)
                                                              ttGAAGGGGCGAAGCCTCATCCCCATTTCTCTGTTCGCCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGgag                                   >  2:1482655/1‑100 (MQ=255)
                                                               cccttgaaggggcgaagcatcatccccatttGGTCACCAGCCGGGAAACCACGTAAGCTCCGTCATAACCCATAACAGAAACGGACTTCGTCAAAGgag                                   >  2:447790‑M2/32‑99 (MQ=255)
                                                                 aaGGGGCGAAGCCTCATCCCCATTTCTCTGGTCGCCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGGAGAGt                                >  2:1927965/1‑100 (MQ=255)
                                                                     ngCGAAGCAACATCCCCATTTCTCTAGTCACCAGCCGGGAAACAACATAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTAAAAGGAGAGtt                               >  2:2097131/2‑97 (MQ=255)
                                                                     ggCGAAGCCTCATCCCCATTTCTCTGGTCGCCAGCCGGGAAACCACGTAAGCTCCGGCGTCACTCATAACAGATACGGACTTTCTCAAAGGAGAGTTATc                            <  2:1541629/100‑1 (MQ=255)
                                                                      gCGAAGCCTCATCCCCATTTCTCTGGTCGCCAGCCGGGAAACCACGTAAGCACCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGGAGAGTTATCa                           <  1:759490/100‑1 (MQ=255)
                                                                         aaGCCTCATCCCCATTTCTCTGGTCACCAGCCTGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGGAGAGTTATCAATg                        <  2:819268/100‑1 (MQ=255)
                                                                         aaGCCTCATCCCCATTTCTCTGGTCACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGGAGAGTTATCAATg                        >  2:341297/1‑100 (MQ=255)
                                                                          gcgaagcctcatccccatttGGTCACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGGAGAGTTATCAATGa                       <  1:1653019‑M2/80‑1 (MQ=255)
                                                                           gCCTCATCCCCATTTCTCTGGTCGCCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGGAGAGTTATCAATGaa                      >  1:613096/1‑100 (MQ=255)
                                                                           gCCTCATCCCCATTTCTCTGGTCGCCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGGAGAGTTATCAATGa                       <  1:640596/99‑1 (MQ=255)
                                                                           gCCTCATCCCCATTTCTCTGGTCGCAAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGGAGAGTTATCAATGaa                      >  2:588388/1‑100 (MQ=255)
                                                                             cTCATCCCCATTTCTCTGGTCACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGGAGAGTTATCAATGAata                    >  1:1614436/1‑100 (MQ=255)
                                                                               cATCCCCATTTCTCTGGTCACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGGAGAGTTATCAATGAATAtt                  <  1:1634114/100‑1 (MQ=255)
                                                                                 tCCCCATTTCTCTGGTCGCCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGGAGAGTTATCAATGAATATTCg                >  1:1657820/1‑100 (MQ=255)
                                                                                 tCCCCATTTCTCTGGTCACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGGAGAGTTATCAATGAATATTCg                <  1:893670/100‑1 (MQ=255)
                                                                                 ctcatccccatttGGTCACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCAATAACAGATACACACTTTCTCAAAGTAGAGTTATCAATAAatat                   >  1:852079‑M2/14‑97 (MQ=255)
                                                                                   cccATTTCTCTGGTCACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACTGACTTTCTCAAAGGAGAGTTATCAATGAATATTCGTc              >  1:2173689/1‑100 (MQ=255)
                                                                                   catccccatttGGTCACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGGAGAGTTATCAATGAATATTCGTc              <  2:211252‑M2/89‑1 (MQ=255)
                                                                                   catccccatttGGTCACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGGAGAGTTATCAATGAATATTCGTc              >  1:2128082‑M2/12‑100 (MQ=255)
                                                                                   catccccatttGGTCACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGGAGAGTTATCAATGAATATTCGTc              <  1:1244840‑M2/89‑1 (MQ=255)
                                                                                    ccATTTCTCTGGTCACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGGAGAGTTATCAATGAATATTCGTcc             <  1:1248029/100‑1 (MQ=255)
                                                                                     tccccatttGGTCACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGGAGAGTTATCAATGAATATTCGTCCa            <  2:1225759‑M2/91‑1 (MQ=255)
                                                                                       cccatttGGTCACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGGAGAGTTATCAATGAATATTCGTCCAtt          >  2:513776‑M2/8‑100 (MQ=255)
                                                                                         catttGGTCACCAGCCGGGAAACCACGTAAGCTCTGGCGTCACCCATAACAGATACGGACTTTCTCAAAGGAGAGTTATCAATGAATATTCGTCCATTGc        >  2:587330‑M2/6‑100 (MQ=255)
                                                                                         catttGGTCACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGGAGAGTTATCAATGAATATTCGTCCATTGc        <  1:578260‑M2/95‑1 (MQ=255)
                                                                                          ctctGGTCACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGGAGAGTTATCAATGAATATTCGTCCATTGCa       <  1:411092/100‑1 (MQ=255)
                                                                                           tctGGTCGCCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGGAGAGTTATCAATGAATATTCGTCCATTGCAt      <  1:1891385/100‑1 (MQ=255)
                                                                                             tGGTCGCCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAATAGATACGGACTTTCTCAAAGGAGAGTTATCAATGAATATTCGTCCATTGCATg     <  2:312909‑M2/98‑1 (MQ=255)
                                                                                             tGGTCACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGGAGAGTTATCAATGAATATTCGTCCATTGCATGa    >  1:41771‑M2/2‑100 (MQ=255)
                                                                                             nggTCACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGGAGAGTTATCAATGAATATTCGTCCATTGCATGa    >  2:244596/2‑100 (MQ=255)
                                                                                               gTCACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGGAGAGTTATCAATGAATATTCGTCCATTGCATGATc  <  2:609725/100‑1 (MQ=255)
                                                                                                  |                                                                                                
TGAGGTGAATCAGGGTTTTCACCCGATTTTGTGCTGATCAGAATTTTTTTTCTTTTTCCCCCTTGAAGGGGCGAAGCCTCATCCCCATTTCTCTGGTCACCAGCCGGGAAACCACGTAAGCTCCGGCGTCACCCATAACAGATACGGACTTTCTCAAAGGAGAGTTATCAATGAATATTCGTCCATTGCATGATC  >  NC_000913/4370518‑4370712

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