Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,407,295 T→A 19.2% intergenic (+29/‑80) accC → / → yhdT acetyl‑CoA carboxylase, biotin carboxylase subunit/DUF997 family putative inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,407,2950TA19.2% 427.1 / 90.9 192intergenic (+29/‑80)accC/yhdTacetyl‑CoA carboxylase, biotin carboxylase subunit/DUF997 family putative inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (92/62);  new base A (25/12);  total (117/75)
Fisher's exact test for biased strand distribution p-value = 4.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.13e-01

Only 100 of 200 total aligned reads displayed.

CGAGAACTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGGACACTCGTTTTGTT  >  NC_000913/3407197‑3407392
                                                                                                  |                                                                                                 
cGAGAACTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGAtttt                                                                                                  >  2:3584807/1‑100 (MQ=255)
cGAGAACTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGAtttt                                                                                                  >  1:1377347/1‑100 (MQ=255)
cGAGAACTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGAtttt                                                                                                  >  2:3158137/1‑100 (MQ=255)
 gagaACTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAGCTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTc                                                                                                 <  2:3568490/100‑1 (MQ=255)
  agaACTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAACACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTcc                                                                                                <  2:3070583/100‑1 (MQ=255)
    aaCTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCgg                                                                                              <  2:1482401/100‑1 (MQ=255)
       ttCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCt                                                                                           >  1:1356666/1‑100 (MQ=255)
        tCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCtt                                                                                          >  1:322870/1‑100 (MQ=255)
        tCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGAAAATCCGGCCtt                                                                                          <  1:2232053/100‑1 (MQ=255)
         ccAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCttt                                                                                         <  2:1424832/100‑1 (MQ=255)
         ccAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGAAAATCCGGCCttt                                                                                         <  1:1771358/100‑1 (MQ=255)
          cAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCATGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCtttt                                                                                        <  1:3634560/100‑1 (MQ=255)
             cATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTtta                                                                                     <  1:1108407/100‑1 (MQ=255)
             cATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTtta                                                                                     <  1:3520814/100‑1 (MQ=255)
             cATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTtta                                                                                     <  2:3305393/100‑1 (MQ=255)
             cATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGAAAATCCGGCCttttgac                                                                                     <  1:2223897/100‑4 (MQ=38)
             cATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGAAAATCCGGCCttttgac                                                                                     <  1:919004/100‑4 (MQ=38)
              aTGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTttat                                                                                    <  2:1243538/100‑1 (MQ=255)
              aTGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTtta                                                                                     <  2:4194683/99‑1 (MQ=255)
              aTGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGAAAATCCGGCCttttgac                                                                                     <  2:3047565/99‑4 (MQ=38)
                ggtggCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTttatta                                                                                  >  2:1519565/1‑100 (MQ=255)
                ggtggCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTttatta                                                                                  <  2:1877440/100‑1 (MQ=255)
                 gtagCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTAc                                                                                 <  2:2712102/97‑1 (MQ=255)
                   ggCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTg                                                                               <  1:4019016/100‑1 (MQ=255)
                   ggCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTg                                                                               >  1:72897/1‑100 (MQ=255)
                    gCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTgg                                                                              >  1:2114694/1‑100 (MQ=255)
                     cACTAACATCCACTATCTTGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTggg                                                                             <  1:3736317/100‑1 (MQ=255)
                     cACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTggg                                                                             <  1:2169964/100‑1 (MQ=255)
                       cTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTggg                                                                             >  2:636079/1‑98 (MQ=255)
                       cTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGa                                                                           <  2:1907133/100‑1 (MQ=255)
                       cTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGa                                                                           >  2:2737449/1‑100 (MQ=255)
                        tAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGAt                                                                          <  2:485177/100‑1 (MQ=255)
                        tAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGAt                                                                          >  1:3366886/1‑100 (MQ=255)
                            aTCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGAc                                                                      <  1:2159278/100‑1 (MQ=255)
                            aTCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGAc                                                                      >  1:2426645/1‑100 (MQ=255)
                             tCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAACACTGCTAAAGCGTCAAAAGGCCGGATCTTCCGGCCTTTTTTATTCCTGGGGa                                                                           >  2:2977600/1‑94 (MQ=255)
                               cACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAAc                                                                   >  1:1760047/1‑100 (MQ=255)
                                  tATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAAcccc                                                                >  1:658642/1‑100 (MQ=255)
                                   aTCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAAcccc                                                                >  1:1596265/1‑99 (MQ=255)
                                    tCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCa                                                              <  1:2733640/100‑1 (MQ=255)
                                     cTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCAt                                                             >  1:2278960/1‑100 (MQ=255)
                                       ggAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATaa                                                           >  2:1254185/1‑100 (MQ=255)
                                         agaaagttcTCGGTCTTCAGGAAAAATAAGACTGCTAAGGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAgg                                                         <  2:2375280/92‑1 (MQ=255)
                                         agaAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAgg                                                         >  2:2097304/1‑100 (MQ=255)
                                          gcaaaaaCTCGGTCCTCAGGAAAAATAAGACTGCTAAAGCGTCACAACGCCCGAGTATCCGCACTTTATTCCTGCTGCAGATCGCCAACCCCCCTACgg                                                         >  2:3888240/3‑99 (MQ=255)
                                          gaAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGt                                                        <  1:2563886/100‑1 (MQ=255)
                                            aaaaaCTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTAc                                                      >  1:339787/1‑100 (MQ=255)
                                            aaaaaCTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTAc                                                      >  2:2359702/1‑100 (MQ=255)
                                            aaaaaCTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTAc                                                      >  1:499145/1‑100 (MQ=255)
                                              aaaCTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACaa                                                    <  2:3258415/100‑1 (MQ=255)
                                               aaCTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAAt                                                   >  2:2458613/1‑100 (MQ=255)
                                                   cGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATcccc                                               <  1:1684536/100‑1 (MQ=255)
                                                    ggTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCg                                              <  1:1333916/100‑1 (MQ=255)
                                                     gTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGCTCGACAACCCCCATAAGGTACAATCCCCGc                                             >  2:3855330/1‑100 (MQ=255)
                                                     gTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCg                                              <  2:2040745/99‑1 (MQ=255)
                                                     gTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGc                                             >  2:1615561/1‑100 (MQ=255)
                                                            ggAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTttcttc                                      >  2:540211/1‑100 (MQ=255)
                                                            ggAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTttcttc                                      >  2:855764/1‑100 (MQ=255)
                                                            ggAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTttcttc                                      >  2:3701687/1‑100 (MQ=255)
                                                            ggAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATCTTCCGGCCTTTCTTATCACTGAGGATCGGCAACCCCCATAAGATACAAACCCCTCTttcttc                                      >  2:1215972/1‑100 (MQ=255)
                                                               aaaaTAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCAcc                                   <  1:1632092/100‑1 (MQ=255)
                                                                aaaTAAGACTGCTAAAGCGTCAAAAGGCCGGATATTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCAccc                                  >  2:2508781/1‑100 (MQ=255)
                                                                 aaTAAGACTGCTAAAGCGTCAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCAt                                >  1:3344983/1‑100 (MQ=255)
                                                                      gACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAgg                            <  1:1639686/100‑1 (MQ=255)
                                                                        cTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAATGTACAATCCCCGCTTTCTTCACCCATCAGGGa                          <  1:3982078/100‑1 (MQ=255)
                                                                         tGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGAc                         <  1:4182567/100‑1 (MQ=255)
                                                                           cTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACaa                       <  2:3981350/100‑1 (MQ=255)
                                                                             aacgCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACaaa                      <  1:1806063/96‑1 (MQ=255)
                                                                             aaaGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACaaaa                     <  2:2176639/100‑1 (MQ=255)
                                                                                gCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAAt                  >  1:2264128/1‑100 (MQ=255)
                                                                                aataaaaaaGGCCGGAAAATCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACaaaaaa                   >  2:2306630/6‑99 (MQ=38)
                                                                                aataaaaaaGGCCGGAAAATCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACaaaaa                    >  2:2542020/6‑98 (MQ=38)
                                                                                aataaaaaaGGCCGGAAAATCCGGCCTTTTTAATTACTGGAGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACaaaaagt                  >  1:2669052/6‑98 (MQ=37)
                                                                                 cGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATg                 <  1:564703/100‑1 (MQ=255)
                                                                                 cGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATg                 <  2:1526951/100‑1 (MQ=255)
                                                                                 ataaaaaaGGCCGGAAAATCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATg                 >  1:1279959/5‑100 (MQ=38)
                                                                                  taaaaaaGGCCGGAAAATCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATg                 >  2:907627/4‑99 (MQ=38)
                                                                                  gTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATgg                <  1:3893284/100‑1 (MQ=255)
                                                                                   aaaaaaGGCCGGAAAATCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGGa               >  1:1011995/3‑100 (MQ=39)
                                                                                   aaaaaaGGCCGGAAAATCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGGa               >  1:2348108/3‑100 (MQ=38)
                                                                                   aaaaaaGGCCGGAAAATCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGGa               >  1:4152265/3‑100 (MQ=39)
                                                                                    cAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGGac              >  1:348406/1‑100 (MQ=255)
                                                                                    cAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGGac              <  2:2775190/100‑1 (MQ=255)
                                                                                    cAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGCGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGGac              <  1:2942584/100‑1 (MQ=255)
                                                                                    cAAAAGGCCAGATTTTCCGGCGTTTTTTATTACTGTTGGTAGACAACCCCCAGAAGGTTAAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGGan              <  2:2289137/100‑2 (MQ=255)
                                                                                    aaaaaGGCCGGAAAATCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGGac              >  2:2659557/2‑100 (MQ=39)
                                                                                    aaaaaGGCCGGAAAATCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGGac              >  1:3828054/2‑100 (MQ=39)
                                                                                    aaaaaGGCCGGAAAATCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGGac              >  2:2320127/2‑100 (MQ=39)
                                                                                    aaaaaGGCCGGAAAATCCGGCATTTGTTATTACTGGGGATCGACAACCACCATAAGGTACAATCTCCGCTTTCTTCACCGATCAGGGGCAAAAAATGGac              >  1:454445/2‑100 (MQ=255)
                                                                                    aaaaaGGCAGGAAAATCCGGCCTTTTTTATTACTGGAGAGCGACAACCACCATACGGTACAATCCCTGCCTTCTTGACCCATCACGGCCAAAAATTGCac              >  2:2985614/2‑100 (MQ=255)
                                                                                     aaaaGGCCGGATTTTCCGGCCTTTTTTATTACTGGGTATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGGaca             >  1:2964056/1‑100 (MQ=255)
                                                                                     aaaaGGCCGGAAAATCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGGaca             >  2:723035/1‑100 (MQ=255)
                                                                                      aaaGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATg                 >  2:3515467/1‑95 (MQ=255)
                                                                                       aaGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGGACACt           >  1:1989394/1‑100 (MQ=255)
                                                                                         ggCCGGAAAATCCGGCCTTTTTCATTACCGGGGATCGACAACCCCCATAAGGTCCAATCCCCGCTTTCCTCACCCATCAGGGACAACAAATGGacaccc          >  2:1116863/1‑97 (MQ=255)
                                                                                          gCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGAACAAAAAATGGACACTCGt        >  2:4126529/1‑100 (MQ=255)
                                                                                           ccGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTGACCCATCAGGGACAAAAAATGGACACTCGtt       >  1:1114350/1‑100 (MQ=255)
                                                                                           ccGGATTTTCCGGCCTTTTTTAATACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAACAATGGACACTc          >  2:3483551/1‑97 (MQ=255)
                                                                                           ccGAATTCTCCCGGCTGTTCTATTACTCGGGATCGACACCCCCAATAAGGTACAATCGCCGCTTTCTTCACCCATCAGGGACAAAAAATGGACACTCg         <  2:2804523/98‑1 (MQ=255)
                                                                                                ttttCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGGACACTCGTTTTGtt  <  1:2590486/100‑1 (MQ=255)
                                                                                                  |                                                                                                 
CGAGAACTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGGACACTCGTTTTGTT  >  NC_000913/3407197‑3407392

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: