Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,407,293 T→A 18.1% intergenic (+27/‑82) 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,2930TA18.1% 456.0 / 85.9 192intergenic (+27/‑82)accC/yhdTacetyl‑CoA carboxylase, biotin carboxylase subunit/DUF997 family putative inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (93/63);  new base A (23/12);  total (117/75)
Fisher's exact test for biased strand distribution p-value = 5.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01

Only 100 of 202 total aligned reads displayed.

TGACGAGAACTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGGACACTCGTTTTGTT  >  NC_000913/3407194‑3407392
                                                                                                   |                                                                                                   
tgaCGAGAACTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGAt                                                                                                     <  1:2077478/100‑1 (MQ=255)
   cGAGAACTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGAtttt                                                                                                  >  2:72911/1‑100 (MQ=255)
   cGAGAACTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGAtttt                                                                                                  >  1:237303/1‑100 (MQ=255)
   cGAGAACTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGAtttt                                                                                                  >  1:1377347/1‑100 (MQ=255)
   cGAGAACTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGAtttt                                                                                                  >  2:508308/1‑100 (MQ=255)
    gagaACTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTc                                                                                                 >  2:284250/1‑100 (MQ=255)
     agaACTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTcc                                                                                                <  2:2784117/100‑1 (MQ=255)
         cTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGcc                                                                                            <  2:4156433/100‑1 (MQ=255)
         cTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGcc                                                                                            >  1:2150803/1‑100 (MQ=255)
          ttCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCt                                                                                           >  1:1356666/1‑100 (MQ=255)
           tCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCtt                                                                                          >  1:322870/1‑100 (MQ=255)
            ccAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCttt                                                                                         >  1:2708232/1‑100 (MQ=255)
            ccAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGAAAATCCGGCCttt                                                                                         <  1:1771358/100‑1 (MQ=255)
            ccAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGAAAATCCGGCCttt                                                                                         <  2:575893/100‑1 (MQ=255)
            ccAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCAGATTTTCCGGCCttt                                                                                         <  2:3916111/100‑1 (MQ=255)
            ccAGCATGGTGGCACTAACATCCACTATCTCGATAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTATCCGGCCttt                                                                                         >  2:2765208/1‑100 (MQ=255)
             cAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCATGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCtttt                                                                                        <  1:3634560/100‑1 (MQ=255)
              aGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCttttt                                                                                       <  1:980079/100‑1 (MQ=255)
              aGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGACCttttt                                                                                       <  1:2200916/100‑1 (MQ=255)
               gCATGTTGGCACTAACATCTACTATCTAGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTTCTAAAGCGTCAAAAGGCCGGAAAATCCGGCCttttga                                                                                      <  2:905431/100‑3 (MQ=37)
               gCATGGTGGCACTAATATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCtttttt                                                                                      >  2:2069328/1‑100 (MQ=255)
                cATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTtta                                                                                     <  1:1108407/100‑1 (MQ=255)
                cATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTtta                                                                                     <  1:3520814/100‑1 (MQ=255)
                 aTGGTGGCACTAACATCCACTCTCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTtta                                                                                     <  2:2703403/99‑1 (MQ=255)
                 aTGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGAAAATCCGGCCttttgac                                                                                     <  2:2507678/99‑4 (MQ=38)
                  tggtggCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTttatt                                                                                   >  2:1243189/1‑100 (MQ=255)
                     tggCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACt                                                                                >  1:2470935/1‑100 (MQ=255)
                      ggCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTg                                                                               >  1:72897/1‑100 (MQ=255)
                      ggCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTg                                                                               <  1:4019016/100‑1 (MQ=255)
                        cACTAACATCCACTATCTTGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTggg                                                                             <  1:3736317/100‑1 (MQ=255)
                          cTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGa                                                                           >  1:2276390/1‑100 (MQ=255)
                           tAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGAt                                                                          >  1:3366886/1‑100 (MQ=255)
                               aTCCACTATCTGGAGAAAAAACTCGGTCTTCAGGCAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTAATGGTGATCGAc                                                                      <  2:2197151/100‑1 (MQ=255)
                               aTCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATcn                                                                        <  2:1427922/98‑2 (MQ=255)
                               aTCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGAc                                                                      >  2:2354680/1‑100 (MQ=255)
                                  cACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAAc                                                                   >  1:1760047/1‑100 (MQ=255)
                                     tATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAAcccc                                                                >  1:658642/1‑100 (MQ=255)
                                      aTCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAAcccc                                                                >  1:1596265/1‑99 (MQ=255)
                                        cTGGAGAAAAAACTCGGTCTTCAGGAAAACTAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCAt                                                             >  2:301080/1‑100 (MQ=255)
                                          ggAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCTTTTTTTATTACTAGGGATCGACAcccccc                                                               >  2:3096468/1‑96 (MQ=255)
                                          ggAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATaa                                                           <  2:3171647/100‑1 (MQ=255)
                                            agaAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAgg                                                         >  2:2436186/1‑100 (MQ=255)
                                             gaAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGt                                                        >  1:2259006/1‑100 (MQ=255)
                                               aaaaaCTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTAc                                                      >  1:339787/1‑100 (MQ=255)
                                               aaaaaCTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTAc                                                      >  1:499145/1‑100 (MQ=255)
                                                      cGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATcccc                                               <  1:1684536/100‑1 (MQ=255)
                                                       ggTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCg                                              <  1:1333916/100‑1 (MQ=255)
                                                        gTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGc                                             >  2:1615561/1‑100 (MQ=255)
                                                        gTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGc                                             >  2:3069589/1‑100 (MQ=255)
                                                          cTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCtt                                           >  1:2166542/1‑100 (MQ=255)
                                                             cAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTttct                                        >  1:2584291/1‑100 (MQ=255)
                                                             cAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTttct                                        <  2:2163831/100‑1 (MQ=255)
                                                               ggAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTACGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTttctt                                       >  2:3262157/1‑99 (MQ=255)
                                                                gAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCa                                     >  2:3796111/1‑100 (MQ=255)
                                                                 aaaaaTAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCAc                                    >  2:265669/1‑100 (MQ=255)
                                                                  aaaaTAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCAcc                                   <  1:1632092/100‑1 (MQ=255)
                                                                  aaaaTAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCAcc                                   >  1:2735160/1‑100 (MQ=255)
                                                                    aaTAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCa                                 >  2:450722/1‑100 (MQ=255)
                                                                    aaTAAGACTGCTAAAGCGTCAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCAt                                >  1:3344983/1‑100 (MQ=255)
                                                                       aaGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCa                              <  2:3577024/100‑1 (MQ=255)
                                                                         gACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAgg                            <  1:1639686/100‑1 (MQ=255)
                                                                         gACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAgg                            <  2:1517426/100‑1 (MQ=255)
                                                                          aCTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAggg                           <  2:1175130/100‑1 (MQ=255)
                                                                           cTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAATGTACAATCCCCGCTTTCTTCACCCATCAGGGa                          <  1:3982078/100‑1 (MQ=255)
                                                                            tGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGAc                         <  1:4182567/100‑1 (MQ=255)
                                                                             gCTAAAGCGTCAAACGGCCGGATTTTCCGGCCTTTTTTATTACCGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACa                        <  2:3498980/100‑1 (MQ=255)
                                                                             gCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGTTACAATCCCCGCTTTCTTCACCCATCAgtg                           >  2:1898829/1‑95 (MQ=255)
                                                                             gCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACa                        >  1:2964960/1‑100 (MQ=255)
                                                                             gCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACa                        >  2:3404098/1‑100 (MQ=255)
                                                                             gCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACa                        <  2:2714822/100‑1 (MQ=255)
                                                                             gCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGAACa                        >  2:2980786/1‑100 (MQ=255)
                                                                               tAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACaaa                      >  1:2291723/1‑100 (MQ=255)
                                                                               tAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACaaa                      >  2:1317878/1‑100 (MQ=255)
                                                                                aacgCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACaaa                      <  1:1806063/96‑1 (MQ=255)
                                                                                aaaGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACaaaa                     >  2:2361281/1‑100 (MQ=255)
                                                                                aaaGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACaaaa                     >  2:657640/1‑100 (MQ=255)
                                                                                aaaGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACaaaa                     >  2:3679333/1‑100 (MQ=255)
                                                                                    cGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATCAGGTATAATCCCCGCTTTCTTCACCCATCAGGGTCCAAAAATg                 >  2:3874032/1‑100 (MQ=255)
                                                                                    cGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATg                 <  1:564703/100‑1 (MQ=255)
                                                                                    cGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATg                 <  2:1526951/100‑1 (MQ=255)
                                                                                    cGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACACCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAAt                  >  2:1243936/1‑99 (MQ=255)
                                                                                    ataaaaaaGGCCGGAAAATCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATg                 >  1:1279959/5‑100 (MQ=38)
                                                                                     gTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATgg                <  1:3893284/100‑1 (MQ=255)
                                                                                      aaaaaaGGCCGGAAAATCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATCAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGGa               >  2:1662678/3‑100 (MQ=38)
                                                                                      aaaaaaGGCCGGAAAATCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGGa               >  1:1011995/3‑100 (MQ=39)
                                                                                      aaaaaaGGCCGGAAAATCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGGa               >  1:2227590/3‑100 (MQ=38)
                                                                                      aaaaaaGGCCGGAAAATCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGGa               >  1:4152265/3‑100 (MQ=39)
                                                                                      aaaaaaGGCCGGAAAATCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTCCTTAACCCATCAGGGACAAAAAATGGa               >  2:3534124/3‑100 (MQ=38)
                                                                                       cAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGGac              >  1:348406/1‑100 (MQ=255)
                                                                                       cAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGCGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGGac              <  1:2942584/100‑1 (MQ=255)
                                                                                       aaaaaGGCCGGAAAATCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGGac              >  2:2518593/2‑100 (MQ=39)
                                                                                       aaaaaGGCCGGAAAATCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGGac              >  1:3828054/2‑100 (MQ=39)
                                                                                       aaaaaGGCCGGAAAATCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGGac              >  2:2320061/2‑100 (MQ=39)
                                                                                       aaaaaGGCCGGAAAATCCGGCATTTGTTATTACTGGGGATCGACAACCACCATAAGGTACAATCTCCGCTTTCTTCACCGATCAGGGGCAAAAAATGGac              >  1:454445/2‑100 (MQ=255)
                                                                                        aaaaGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGGaca             >  2:229279/1‑100 (MQ=255)
                                                                                          aaGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGGACACt           >  1:1989394/1‑100 (MQ=255)
                                                                                              ccGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTGACCCATCAGGGACAAAAAATGGACACTCGtt       >  1:1114350/1‑100 (MQ=255)
                                                                                              ccGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGGACACTCGtt       >  2:196818/1‑100 (MQ=255)
                                                                                              ccGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGGACACTCGtt       >  2:4070654/1‑100 (MQ=255)
                                                                                                 gATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGGACACTCGTTTTg    <  1:2646685/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
TGACGAGAACTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGGACACTCGTTTTGTT  >  NC_000913/3407194‑3407392

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