Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,612,329 A→G 8.8% intergenic (+40/‑82) ytjA → / → yjjU uncharacterized protein/putative patatin‑like family phospholipase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,612,3290AG 8.8% 182.4 / 8.4 68intergenic (+40/‑82)ytjA/yjjUuncharacterized protein/putative patatin‑like family phospholipase
Reads supporting (aligned to +/- strand):  ref base A (24/38);  new base G (3/3);  total (27/41)
Fisher's exact test for biased strand distribution p-value = 6.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.33e-01

TGGTGGTCTGGCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTGCGCCGTTGTCGCTAAGGCCGTTCCAGCCTG  >  NC_000913/4612187‑4612477
                                                                                                                                              |                                                                                                                                                    
tggtggTCTGGCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGAGGAc                                                                                                                                                >  2:429335/1‑149 (MQ=255)
  gtggtCTGGCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAATCCAGTACATCGGAATg                                                                                                                                              >  2:61272/1‑149 (MQ=255)
       cTGGCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGGCGGGATTCTTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCGCTGGACTGGCttt                                                                                                                                         <  2:482557/149‑1 (MQ=255)
       cTGGCCGGTACCGCTGCAGGCGCAGCTAAAATTGGCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGGCCGCTGGACTGGCttt                                                                                                                                         <  1:163139/149‑1 (MQ=255)
        tGGCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTCTCGTCGGGATTCTTCTGTTCCTGGTTAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCtttt                                                                                                                                        <  2:496889/149‑1 (MQ=255)
         ggCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCGCTGGACTGGCtttgg                                                                                                                                       <  1:868986/149‑3 (MQ=255)
          gCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGc                                                                                                                                      <  1:681738/149‑1 (MQ=255)
          gCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGc                                                                                                                                      >  2:1018490/1‑149 (MQ=255)
          gCCGGGACCGCTGCAGGCGCAGCTAACATTGTCCTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGCCCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGc                                                                                                                                      <  1:879172/149‑1 (MQ=255)
            cGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAACCGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCgg                                                                                                                                    <  2:728754/149‑1 (MQ=255)
            cGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCgg                                                                                                                                    >  1:815983/1‑149 (MQ=255)
            cGGTACCGCTGCAGGCGCAGCTAAAATTGTCTGTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCGAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCgg                                                                                                                                    >  1:190460/1‑149 (MQ=255)
             gggaCCGCTGCTGGCGCAGCTAAAATTGTCTTGGGAGTCGCGATTATTCTGTTCCTGGTGAGTTTGAGCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGt                                                                                                                                   <  2:874335/146‑1 (MQ=255)
             ggTACCGCTGCAGGCGCATCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGt                                                                                                                                   <  2:484068/149‑1 (MQ=255)
                 ccGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTa                                                                                                                               <  2:815983/149‑1 (MQ=255)
                  cgcgGCAGGCGCAGCTAAAAGTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAg                                                                                                                              <  1:796595/149‑1 (MQ=255)
                   gCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGGGAGTTTGAGCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGc                                                                                                                             <  2:351875/149‑1 (MQ=255)
                     ggCATGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTTATCTGTCCCTGGTGAGTTTGTGCATGGGCCGAAAACGACCCTTGATTTCAGACCATTACGAAAAATTTATCCAACGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGa                                                                                                                           <  2:589433/148‑1 (MQ=255)
                                cTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGCCCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTg                                                                                                                >  2:969036/1‑149 (MQ=255)
                                 tAAAATTGTCTTTGTCGTCGGGATTATTCTGTTACTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTAAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTATGCGGTTTTAGCGAATAATAAATTGc                                                                                                               >  1:278410/1‑149 (MQ=255)
                                  aaaaTTGGCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTAAGCGAATAATAAATTGCg                                                                                                              <  1:16254/149‑1 (MQ=255)
                                      ttGTCTTGGTCGTCGGCATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTAc                                                                                                          <  1:850830/149‑1 (MQ=255)
                                          cTTTGTCGGCGGGATTATTCTGTGCCTGGGGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTg                                                                                                      <  1:511181/149‑1 (MQ=255)
                                           tttGTCGGCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGt                                                                                                     <  1:280989/149‑1 (MQ=255)
                                              gtcgtcGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCATTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCAt                                                                                                  >  2:701882/1‑149 (MQ=255)
                                              gtcgtcGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCAt                                                                                                  >  2:971003/1‑149 (MQ=255)
                                                 gtcgGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATTATAAATTGCGTTACTTTGTCATTCt                                                                                               >  2:236320/1‑149 (MQ=255)
                                                  gccggATTATTCTGATCCTGGGGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCtt                                                                                              <  2:289660/146‑1 (MQ=255)
                                                   cgGGATTATTCTGTTCCTGGTGATTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCTGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCttt                                                                                             <  1:256820/149‑1 (MQ=255)
                                                     ggATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTca                                                                                           <  2:1018440/149‑1 (MQ=255)
                                                               gTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAATAAAATTTTATCAAAAGCCAGTACAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAATAAc                                                                                 <  1:157619/149‑1 (MQ=255)
                                                                ttCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACa                                                                                <  1:764286/149‑1 (MQ=255)
                                                                 tCCTGGTGAGTTTGTTCATGGGCCGTTAACGACGCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAg                                                                               <  1:447617/149‑1 (MQ=255)
                                                                 tCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAg                                                                               <  1:1018490/149‑1 (MQ=255)
                                                                 tCCTGGTGAGTGTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAACTTGCGTTACTTTGTCATTCTTTCACAACAAAAACAg                                                                               <  2:749426/149‑1 (MQ=255)
                                                                      gTGAGTTTGTTCATGGGCCGAAAACTACCCTAGATTTCATACCATTAAGAAAATTTTATCCAAATCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAgg                                                                          <  1:583926/149‑1 (MQ=255)
                                                                         aGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCag                                                                       <  1:191970/149‑1 (MQ=255)
                                                                         aGTTTGATCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCag                                                                       <  1:741852/149‑1 (MQ=255)
                                                                         aGTTGGTTCATGGGCCGAAAACGACCCTAGATATCAGACCATTAAGAAAATTTTATCCAAAGCCCGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCag                                                                       <  2:680233/149‑1 (MQ=255)
                                                                           tttGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCagag                                                                     <  1:967630/149‑1 (MQ=255)
                                                                              gTTCATGGTCCGAAAACGACCCTACATTTAAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGATTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGtt                                                                  <  2:362454/149‑2 (MQ=255)
                                                                                 cATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTgggg                                                               <  1:280179/149‑1 (MQ=255)
                                                                                      gCCGAAAACGCCCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCTGACTGGCTTTTGAGGTTTTAGCGATTTATAAATTGCGTTCCTTTGTGATTCTTTCAGAACAAAATCAGGAAGGCTGAGGTGTGGAAGCg                                                          >  2:707551/1‑149 (MQ=255)
                                                                                          aaaaCGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATa                                                      >  2:421504/1‑149 (MQ=255)
                                                                                             aCGAGCCTAGATTTCCTACCATTAAGAAAATTTTCTCCAACTCCAGTCCAGCTGACTGGCTTTTGCGGTTTTAGCGAATAATAATTTGCGATACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCt                                                   <  1:436093/149‑1 (MQ=255)
                                                                                              cGACCCTAGATTTCACACCATTAAGAAAATTTTATCCAAAGCCCGGCCAGCGGACGGGCTTTTGCGGTGTTAGCGAATAATAAATTGCCTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGGGGGGCAGCGAATACCTg                                                  <  2:914189/149‑1 (MQ=255)
                                                                                               gACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGTCAGCGAATACCTGt                                                 >  1:278024/1‑149 (MQ=255)
                                                                                                 cccTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTaa                                               >  1:103562/1‑149 (MQ=255)
                                                                                                  ccTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTTCGTTACTTTGTCATTCTTTCACAACAAAAAAAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTaaa                                              >  1:834560/1‑148 (MQ=255)
                                                                                                   cTAGATTACAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATAGCGTTACTTTGTCAATCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCTAATCCATGTAACg                                             >  2:434202/1‑149 (MQ=255)
                                                                                                     aGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCt                                           >  2:389515/1‑149 (MQ=255)
                                                                                                       aTTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTg                                         >  2:820956/1‑149 (MQ=255)
                                                                                                          tCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTa                                      >  2:799585/1‑149 (MQ=255)
                                                                                                                cATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTg                                >  1:858405/1‑149 (MQ=255)
                                                                                                                cATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGGCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTg                                <  2:639487/149‑1 (MQ=255)
                                                                                                                 aTTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTgc                               >  2:1026174/1‑149 (MQ=255)
                                                                                                                   tAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTgcgc                             <  1:869252/149‑1 (MQ=255)
                                                                                                                      gAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTATGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTGAGCAGt                          >  1:368114/1‑149 (MQ=255)
                                                                                                                           ttttATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTGCGCCGTTGTCg                     <  1:308512/149‑1 (MQ=255)
                                                                                                                           ttttATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTGCGCCGTTGTCg                     >  2:406033/1‑149 (MQ=255)
                                                                                                                                 ccAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTGTTAGCGAATAATAAATTTCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTGCGCCGTTGTCGCTAAgg               <  2:190460/149‑1 (MQ=255)
                                                                                                                                 caaaaGCCAGTCCGCTGGACTGGCTTTTGCGGGTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTGCGCCGTTGTCGCTTAgg               >  1:321097/3‑149 (MQ=255)
                                                                                                                                    aaGCCAGTCCGCTGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTGCGCCGTTGTCGCTAAGGCCg            >  1:833014/1‑149 (MQ=255)
                                                                                                                                    aaGCCAGTCCGCTGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTGCGCCGTTGTCGCTAAGGCCg            >  1:409303/1‑149 (MQ=255)
                                                                                                                                      gCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGATTTATTAATTGCGTTAGTTTGTCATTCTTTCACAACAAAATCAGGAAGGGAGTGGTGGGGGTTCGAATACCTGTAACGCTTTGTATTATTGCGCCGTTGTCGCTAAGGCCGtt          <  2:322981/149‑1 (MQ=255)
                                                                                                                                           tCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTGCGCCGTTGTCGCTAAGGCCGTTCCAGc     >  1:957551/1‑149 (MQ=255)
                                                                                                                                           gccAGCGTTCTGGCTTTAGCGCTGTTAGCGAAAAATAAAATGCGTTACTTTGTCAGTCTTTCAGAACAAAAACAGGTAGGCAGAGGTGGGTCAGGGAATACCTGTAACGCTTGGTAATATTGCGCCGTTGTCGCTAAGGGCGTTCCacc     <  2:16107/148‑3 (MQ=255)
                                                                                                                                             cAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTGCGCCGTTGTCGCTAAGGCCGTTCCAGCCt   <  1:432773/149‑1 (MQ=255)
                                                                                                                                              aGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAACCAGGCAGGCAGAGGGGGGGCAGCGAATACCTGTAACGCTTGGTAATATTGCGCCGTTGTCGCTAAGGCCGTTCCAGCCTg  <  1:971003/149‑1 (MQ=255)
                                                                                                                                              aGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTGCGCCGTTGTCGCTAAGGCCGTTCCAGCCTg  <  1:701882/149‑1 (MQ=255)
                                                                                                                                              aGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACACGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTGCGCCGTTGTCGCTAAGGCCGTTCCAGCCTg  <  1:540638/149‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                    
TGGTGGTCTGGCCGGTACCGCTGCAGGCGCAGCTAAAATTGTCTTTGTCGTCGGGATTATTCTGTTCCTGGTGAGTTTGTTCATGGGCCGAAAACGACCCTAGATTTCAGACCATTAAGAAAATTTTATCCAAAGCCAGTCCAGCGGACTGGCTTTTGCGGTTTTAGCGAATAATAAATTGCGTTACTTTGTCATTCTTTCACAACAAAAACAGGAAGGCAGAGGTGGGGCAGCGAATACCTGTAACGCTTGGTAATATTGCGCCGTTGTCGCTAAGGCCGTTCCAGCCTG  >  NC_000913/4612187‑4612477

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: