Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,503,125 G→A *223* (TAG→TAA ydcL → lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,503,1250GA97.6% 258.3 / ‑2.3 83*223* (TAG→TAAydcLlipoprotein
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base A (42/39);  total (43/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.04e-01

GGTTATCAAAGTCGTTCGTCAGGGCGAAGGTAAAGACCTGAATAACGAAAGTACGCCAATGGCTTTCGAAAATATTAAACAAGTTATTGATGACATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAGTTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAACCCCATCCTTCAGGGCGGGGATATAAGGCGCGGTTTTCCACCTGACCAGGTGTTTGCTGTTGTTCAATGTATTGGCGGATG  >  NC_000913/1502990‑1503257
                                                                                                                                       |                                                                                                                                    
ggTTATCAAAGTCGTTCGTCAGGGCGAAGGTAAAGACCTGAATAACGAAAGTACGCCAATGGCTTTCGAAAATATTAAACAAGTTATTGATGACATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTc                                                                                                                                   >  1:508331/1‑139 (MQ=255)
 gTTATCAAAGTCGTTCGTCAGGGCGAAGGTAAAGACCTGAATAACGAAAGTACGCCAATGGCTTTCGAAAATATTAAACAAGTTATTGATGACATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTcc                                                                                                                                  >  1:477648/1‑139 (MQ=255)
 gTTATCAAAGTCGTTCGTCAGGGCGAAGGTAAAGACCTGAATAACGAAAGTACGCCAATGGCTTTCGAAAATATTAAACAAGTTATTGATGACATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTcc                                                                                                                                  >  2:287568/1‑139 (MQ=255)
   tATCAAAGTCGTTCGTCAGGGCGAAGGTAAAGACCTGAATAACGAAAGTACGCCAATGGCTTTCGAAAATATTAAACAAGTTATTGATGACATGGCGACCGATGCCACCATGTTTGACGTTAACAAGAAGTAATTCCAg                                                                                                                                >  2:977470/1‑139 (MQ=255)
    aTCAAAGTCGTTCGTCAGGGCGAAGGTAAAGACCTGAATAACGAAAGTACGCCAATGGCTTTCGAAAATATTAAACAAGTTATTGATGACATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGa                                                                                                                               >  1:592215/1‑139 (MQ=255)
         aGTCGTTCGTCAGGGCGAAGGTAAAGACCTGAATAACGAAAGTACGCCAATGGCTTTCGAAAATATTAAACAAGTTATTGATGACATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACgcgc                                                                                                                          <  2:592215/139‑1 (MQ=255)
         aGTCGTTCGTCAGGGCGAAGGTAAAGACCTGAATAACGAAAGTACGCCAATGGCTTTCGAAAATATTAAACAAGTTATTGATGACATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACgcgc                                                                                                                          <  1:287568/139‑1 (MQ=255)
         aGTCGTTCGTCAGGGCGAAGGTAAAGACCTGAATAACGAAAGTACGCCAATGGCTTTCGAAAATATTAAACAAGTTATTGATGACATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACgcgc                                                                                                                          <  2:477648/139‑1 (MQ=255)
         aGTCGTTCGTCAGGGCGAAGGTAAAGACCTGAATAACGAAAGTACGCCAATGGCTTTCGAAAATATTAAACAAGTTATTGATGACATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACgcgc                                                                                                                          <  1:1444982/139‑1 (MQ=255)
         aGTCGTTCGTCAGGGCGAAGGTAAAGACCTGAATAACGAAAGTACGCCAATGGCTTTCGAAAATATTAAACAAGTTATTGATGACATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACgcgc                                                                                                                          <  2:1365102/139‑1 (MQ=255)
              ttcgtCAGGGCGAAGGTAAAGACCTGAATAACGAAAGTACGCCAATGGCTTTCGAAAATATTAAACAAGTTATTGATGACATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACgcgc                                                                                                                          <  1:1814516/134‑1 (MQ=255)
              ttcgtCAGGGCGAAGGTAAAGACCTGAATAACGAAAGTACGCCAATGGCTTTCGAAAATATTAAACAAGTTATTGATGACATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACgcgc                                                                                                                          >  2:1814516/1‑134 (MQ=255)
               tcgtCAGGGCGAAGGTAAAGACCTGAATAACGAAAGTACGCCAATGGCTTTCGAAAATATTAAACAAGTTATTGATGACATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGt                                                                                                                    >  1:1784124/1‑139 (MQ=255)
                      ggCGAAGGTAAAGACCTGAATAACGAAAGTACGCCAATGGCTTTCGAAAATATTAAACAAGTTATTGATGACATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATgg                                                                                                             <  2:294109/139‑1 (MQ=255)
                                       gAATAACGAAAGTACGCCAATGGCTTTCGAAAATATTAAACAAGTTATTGATGACATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATgg                                                                                                             <  2:1733847/122‑1 (MQ=255)
                                       gAATAACGAAAGTACGCCAATGGCTTTCGAAAATATTAAACAAGTTATTGATGACATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATgg                                                                                                             >  1:1733847/1‑122 (MQ=255)
                                       gAATAACGAAAGTACGCCAATGGCTTTCGAAAATATTAAACAAGTTATTGATGACATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTgc                                                                                            >  2:756770/1‑139 (MQ=255)
                                        aaTAACGAAAGTACGCCAATGGCTTTCGAAAATATTAAACAAGTTATTGATGACATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTCATCCGGTgcg                                                                                           >  2:1627316/1‑139 (MQ=255)
                                             cGAAAGTACGCCAATGGCTTTCGAAAATATTAAACAAGTTATTGATGACATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTa                                                                                      >  2:1023781/1‑139 (MQ=255)
                                                  gTACGCCAATGGCTTTCGAAAATATTAAACAAGTTATTGATGACATGGCGACCGATGCCACCATGTTTGACGTTAACAAGAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAAcc                                                                                 <  1:977470/139‑1 (MQ=255)
                                                      gCCAATGGCTTTCGAAAATATTAAACAAGTTATTGATGACATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACGCGcc                                                                                                                         <  1:978148/95‑1 (MQ=255)
                                                      gCCAATGGCTTTCGAAAATATTAAACAAGTTATTGATGACATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACGCGcc                                                                                                                         >  2:978148/1‑95 (MQ=255)
                                                          aTGGCTTTCGAAAATATTAAACAAGTTATTGATGACATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACgcgc                                                                                                                          >  2:1407740/1‑90 (MQ=255)
                                                          aTGGCTTTCGAAAATATTAAACAAGTTATTGATGACATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACgcgc                                                                                                                          <  1:1407740/90‑1 (MQ=255)
                                                            ggCTTTCGAAAATATTAAACAAGTTATTGATGACATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAACCCCATCCTTCa                                                                       <  2:508331/139‑1 (MQ=255)
                                                               tttCGAAAATATTAAACAAGTTATTGATGACATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAACCCCATCCTTCAggg                                                                    <  2:1784124/139‑1 (MQ=255)
                                                                ttCGAAAATATTAAACAAGTTATTGATGACATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAACCCCATCCTTCAGGGc                                                                   >  1:209956/1‑139 (MQ=255)
                                                                 tCGAAAATATTAAGCAAGTTATTGATGACATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGgc                                                                                                            <  2:999183/97‑1 (MQ=255)
                                                                 tCGAAAATATTAAACAAGTTATTGATGACATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGgc                                                                                                            >  1:999183/1‑97 (MQ=255)
                                                                 tCGAAAATATTAAACAAGTTATTGATGACATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAACCCCATCCTTCAGGGCg                                                                  <  2:970159/139‑1 (MQ=255)
                                                                       atatTAAACAAGTTATTGATGACATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGAt                                                                                                               <  1:1033261/88‑1 (MQ=255)
                                                                       atatTAAACAAGTTATTGATGACATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGAt                                                                                                               >  2:1033261/1‑88 (MQ=255)
                                                                                  gTTATTGATGACATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAACCCCATCCTTCAggg                                                                    <  1:632407/120‑1 (MQ=255)
                                                                                  gTTATTGATGACATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAACCCCATCCTTCAggg                                                                    >  2:632407/1‑120 (MQ=255)
                                                                                  gTTATTGATGACATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAACCCCATCCTTCAGGGCGGGGATATAAGGCGCGGt                                                 >  2:429429/1‑139 (MQ=255)
                                                                                    tATTGATGACATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAACCCCATCCTTCAGGGCGGGGATATAAGGCGCGGttt                                               <  1:756770/139‑1 (MQ=255)
                                                                                         atgaCATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACgcgc                                                                                                                          >  2:1170336/1‑59 (MQ=255)
                                                                                         atgaCATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACgcgc                                                                                                                          <  1:1170336/59‑1 (MQ=255)
                                                                                         atgaCATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAACCCCATCCTTCAGGGCGGGGATATAAGGCGCGGTTTTCCAc                                          <  1:6840/139‑1 (MQ=255)
                                                                                          tgaCATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAACCCCATCCTTCAGGGCGGGGATATAAGGCGCGGTTTTCCAcc                                         <  1:1023781/139‑1 (MQ=255)
                                                                                            aCATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAGTTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAACCCCATCCTTCAGGGCGGGGATATAAGGCGCGGTTTTCCACCTg                                       >  1:31374/1‑139 (MQ=255)
                                                                                            aCATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAACCCCATCCTTCAGGGCGGGGATATAAGGCGCGGTTTTCCACCTg                                       >  1:928491/1‑139 (MQ=255)
                                                                                            aCATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAACCCCATCCTTCAGGGCGGGGATATAAGGCGCGGTTTTCCACCTg                                       >  2:1480312/1‑139 (MQ=255)
                                                                                             cATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAACCCCATCCTTCAGGGCGGGGATATAAGgcgc                                                    <  1:105131/125‑1 (MQ=255)
                                                                                             cATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAACCCCATCCTTCAGGGCGGGGATATAAGgcgc                                                    >  2:105131/1‑125 (MQ=255)
                                                                                             cATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAACCCCATCCTTCAGGGCGGGGATATAAGGCGCGGTTTTCCACCTGa                                      >  1:343239/1‑139 (MQ=255)
                                                                                                ggCGTCCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGAcgcg                                                                                                                           >  1:1094948/1‑51 (MQ=38)
                                                                                                ggCGTCCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAAcccc                                                                               >  1:1093917/1‑95 (MQ=255)
                                                                                                ggCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAGTTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAACCCCATCCTTCAGGGCGGGGATATAAGGCGCGGTTTTCCACCTGACCa                                   <  2:31374/139‑1 (MQ=255)
                                                                                                ggCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGa                                                                                                                               <  2:1760481/47‑1 (MQ=255)
                                                                                                ggCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGa                                                                                                                               >  1:1760481/1‑47 (MQ=255)
                                                                                                ggCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGAcgcg                                                                                                                           <  2:1094948/51‑1 (MQ=255)
                                                                                                ggCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGgc                                                                                                            >  2:1408623/1‑66 (MQ=255)
                                                                                                ggCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGgc                                                                                                            <  1:1408623/66‑1 (MQ=255)
                                                                                                ggCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAAcccc                                                                               <  2:1093917/95‑1 (MQ=255)
                                                                                                ggCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAACCCCATCCTTCAggg                                                                    >  2:644889/1‑106 (MQ=255)
                                                                                                ggCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAACCCCATCCTTCAggg                                                                    <  1:644889/106‑1 (MQ=255)
                                                                                                ggCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAACCCCATCCTTCAGGGCGGGGATATAAGGCGCGGTTTTCCa                                           <  2:31852/131‑1 (MQ=255)
                                                                                                ggCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAACCCCATCCTTCAGGGCGGGGATATAAGGCGCGGTTTTCCa                                           >  1:31852/1‑131 (MQ=255)
                                                                                                    aCCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAACCCCATCCTTCAGGGCg                                                                  >  1:1691941/1‑104 (MQ=255)
                                                                                                    aCCGATGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAACCCCATCCTTCAGGGCg                                                                  <  2:1691941/104‑1 (MQ=255)
                                                                                                        aTGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAACCCCATCCTTCAggg                                                                    <  2:1619881/98‑1 (MQ=255)
                                                                                                        aTGCCACCATGTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAACCCCATCCTTCAggg                                                                    >  1:1619881/1‑98 (MQ=255)
                                                                                                           ccaccaTGTTTGACGTTAACAAAAAGTAATTCCAGACgcgc                                                                                                                          >  1:1252700/1‑41 (MQ=255)
                                                                                                           ccaccaTGTTTGACGTTAACAAAAAGTAATTCCAGACgcgc                                                                                                                          <  2:1252700/41‑1 (MQ=255)
                                                                                                           ccaccaTGTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAACCCCATCCTTCAGGTCGTGGATATAAGGGGCGGTTTTCCACCTGACCAGGTGTTTGCtg                        >  1:865242/1‑139 (MQ=255)
                                                                                                               caTGTTTGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAACCCCATCCTTCAGGGCGGGGATATAAGGCGCGGTTTTCCACCTGACCAGGTGTTTGCtgttgt                    >  1:369541/1‑139 (MQ=255)
                                                                                                                   tttGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATgg                                                                                                             >  1:136830/1‑46 (MQ=255)
                                                                                                                   tttGACGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATgg                                                                                                             <  2:136830/46‑1 (MQ=255)
                                                                                                                       aCGTTAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAACCCCATCCTTCAGGGCGGGGATATAAGGCGCGGTTTTCCACCTGACCAGGTGTTTGCTGTTGTTCAATGTa            <  1:578180/139‑1 (MQ=255)
                                                                                                                           tAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAACCCCATCCTTCAGGGCGGGGATATAAGGCGCGGTTTTCCACCTGAcc                                    <  2:1540593/111‑1 (MQ=255)
                                                                                                                           tAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAACCCCATCCTTCAGGGCGGGGATATAAGGCGCGGTTTTCCACCTGAcc                                    >  1:1540593/1‑111 (MQ=255)
                                                                                                                           tAACAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAACCCCATCCTTCAGGGCGGGGATATAAGGCGCGGTTTTCCACCTGACCAGGTGTTTGCTGTTGTTCAATGTATTgg        >  2:767493/1‑139 (MQ=255)
                                                                                                                            aaCAAAAAGTAATTCCAGACGCGCCATCGTTAGATgg                                                                                                             <  1:1226752/37‑1 (MQ=255)
                                                                                                                            aaCAAAAAGTAATTCCAGACGCGCCATCGTTAGATgg                                                                                                             >  2:1226752/1‑37 (MQ=255)
                                                                                                                            aaCAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAACCCCATCCtt                                                                         <  2:551478/73‑1 (MQ=255)
                                                                                                                            aaCAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAACCCCATCCtt                                                                         >  1:551478/1‑73 (MQ=255)
                                                                                                                            aaCAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAACCCCATCCTTCAGGGCGGGGATATAAGgcgc                                                    >  2:1694068/1‑94 (MQ=255)
                                                                                                                            aaCAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAACCCCATCCTTCAGGGCGGGGATATAAGgcgc                                                    <  1:1694068/94‑1 (MQ=255)
                                                                                                                            aaCAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAACCCCATCCTTCAGGGCGGGGATATAAGGCGCGGTTTTCCACCTGACCAGGTGTTTGCtg                        <  2:503711/122‑1 (MQ=255)
                                                                                                                            aaCAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAACCCCATCCTTCAGGGCGGGGATATAAGGCGCGGTTTTCCACCTGACCAGGTGTTTGCtg                        >  1:503711/1‑122 (MQ=255)
                                                                                                                            aaCAAAAAGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAACCCCATCCTTCAGGGCGGGGATATAAGGCGCGGTTTTCCACCTGACCAGGTGTTTGCTGTTGTTCAATGTATTGGc       >  2:1454386/1‑139 (MQ=255)
                                                                                                                                 aaaGTAATTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAACCCCATCCTTCAGGGCGGGGATATAAGGCGCGGTTTTCCACCTGACCAGGTGTTTGCTGTTGTTCAATGTATTGGCGGATg  <  1:1480312/139‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                    
GGTTATCAAAGTCGTTCGTCAGGGCGAAGGTAAAGACCTGAATAACGAAAGTACGCCAATGGCTTTCGAAAATATTAAACAAGTTATTGATGACATGGCGACCGATGCCACCATGTTTGACGTTAACAAAAAGTAGTTCCAGACGCGCCATCGTTAGATGGCGCTTTTTATCCGGTGCGCCGTAAAACCCCATCCTTCAGGGCGGGGATATAAGGCGCGGTTTTCCACCTGACCAGGTGTTTGCTGTTGTTCAATGTATTGGCGGATG  >  NC_000913/1502990‑1503257

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