Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 383,477 G→A 23.1% G153G (GGC→GGT yaiP ← putative family 2 glycosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913383,4770GA23.1% 174.3 / 47.9 91G153G (GGC→GGTyaiPputative family 2 glycosyltransferase
Reads supporting (aligned to +/- strand):  ref base G (33/37);  new base A (13/8);  total (46/45)
Fisher's exact test for biased strand distribution p-value = 3.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.42e-01

TAATAAACGGTGCGCCACCCAGGAGCTGCTGTAGCGTGCGCTTCATAACAATCATCGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCGCCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAATTTCTGCCAGCATATAGCCCATTCCGTCTTGATCGGGCGGAACATAG  >  NC_000913/383380‑383576
                                                                                                 |                                                                                                   
taataaACGGTGCGCCACCCAGGAGCTGCTGTAGCGTGCGCTTCATAACAATCATCGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCAGCGcc                                                                                                   >  1:990037/1‑100 (MQ=255)
  ataaACGGTGCGCCACCCAGGAGCTGCTGTAGCGTGCGCTTCATAACAATCATCGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCGCCtt                                                                                                 <  2:1979982/100‑1 (MQ=255)
  ataaACGGTGCGCCACCCAGGAGCTGCTGTAGCGTGCGCTTCATAACAATCATCGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCACCtt                                                                                                 >  2:1754387/1‑100 (MQ=255)
    aaaCGGTGCGCCACCCAGGAGCTGCTGTAGCGTGCGCTTCATAACAATCATCGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCGCCTTTc                                                                                               <  1:1554948/100‑1 (MQ=255)
     aaCGGTGCGCCACCCAGGAGCTGCTGTAGCGTGCGCTTCATAACAATCATCGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCGCCTTTCa                                                                                              <  2:605434/100‑1 (MQ=255)
      aCGGTGCGCCACCCAGGAGCTGCTGTAGCGTGCGCTTCATAACAATCATCGGCAACTGTACGGTCGCGCGGATGTGCGGTAACAGACCCGCACCTTTCaa                                                                                             <  2:1153189/100‑1 (MQ=255)
       cGGTGCGCCACCCAGGAGCTGCTGTAGCGTGCGCTTCATAACAATCATCGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCGCCTTTCaa                                                                                             <  2:582451/99‑1 (MQ=255)
       cGGTGCGCCACCCAGGAGCTGCTGTAGCGTGCGCTTCATAACAATCATCGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCGCCTTTCAAc                                                                                            <  1:510488/100‑1 (MQ=255)
       cGGTGCGCCACCCAGGAGCTGCTGTAGCGTGCGCTTCATAACAATCATCGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCACCTTTCaa                                                                                             <  1:1332448/99‑1 (MQ=255)
         gTGCGCCACCCAGGAGCTGCTGTAGCGTGCGCTTCATAACAATCATCGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCGCCTTTCAACGc                                                                                          <  1:129847/100‑1 (MQ=255)
         gTGCGCCACCCAGGAGCTGCTGTAGCGTGCGCTTCATAACAATCATCGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCGCCTTTCAACGc                                                                                          <  1:1522158/100‑1 (MQ=255)
          tGCGCCACCCAGGAGCTGCTGTAGCGTGCGCTTCATAACAATCATCGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCGCCTTTCAACGCa                                                                                         >  2:1959171/1‑100 (MQ=255)
           gcgcCACCCAGGAGCTGCTGTAGCGTGCGCTTCATAACAATCATCGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCGCCTTTCAACGCAg                                                                                        <  1:188920/100‑1 (MQ=255)
              ccACCCAGGAGCTGCTGTAGCGTGCGCTTCATAACAATCATCGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCGCCTTTCAACGCAGTag                                                                                     >  1:2079013/1‑100 (MQ=255)
              ccACCCAGGAGCTGCTGTAGCGTGCGCTTCATAACAATCATCGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCGCCTTTCAACGCAGTag                                                                                     >  1:1318084/1‑100 (MQ=255)
               naCCCAGGAGCTGCTGTAGCGTGCGCTTCATAACAATCATCGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCGCCTTTCAACGCAGTag                                                                                     >  2:1712301/2‑99 (MQ=255)
                aCCCAGGAGCTGCTGTAGCGTGCGCTTCATAACAATCATCGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCGCCTTTCAACGCAGTaga                                                                                    >  2:279739/1‑99 (MQ=255)
                 cccAGGAGCTGCTGTAGCGTGCGCTTCATAACAATCATCGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCGCCTTTCAACGCAGTAGAgg                                                                                  <  1:1391300/100‑1 (MQ=255)
                 cccAGGAGCTGCTGTAGCGTGCGCTTCATAACAATCATCGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCACCTTTCAACGCAGTAGAgg                                                                                  <  1:1534564/100‑1 (MQ=255)
                    aGGAGCTGCTGTAGCGTGCGCTTCATAACAATCATCGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCACCTTTCAACGCAGTAGAGGGaa                                                                               >  1:254623/1‑100 (MQ=255)
                    aGGAGCTGCTGTAGCGTGCGCTTAATAACAATCATCGGCAACTTTACGGTAGCGCGGATGTGCGGTAACAGACCCGCACCTTTCAACGCAGTAGAGGGaa                                                                               >  2:997431/1‑100 (MQ=255)
                       aGCTGCTGTAGCGTGCGCTTCATAACAATCATCGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCGCCTTTCAACGCAGTAGAGGGAATgc                                                                            <  1:148075/100‑1 (MQ=255)
                         ctgctgTAGCGTGCGCTTCATAACAATCATCGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCGCCTTTCAACGCAGTAGAGGGAATgccg                                                                          <  2:1126975/100‑1 (MQ=255)
                           gctgTAGCGTGCGCTTCATAACAATCATCGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCGCCTTTCAACGCAGTAGAGGGAATgccgcc                                                                        >  1:1817543/1‑100 (MQ=255)
                           gctgTAGCGTGCGCTTCATAACAATCATCGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCGCCTTTCAACGCAGTAGAGGGAATgccg                                                                          >  2:901737/1‑98 (MQ=255)
                             tgTAGCGTGCGCTTCATAACAATCATCGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCGCCTTTCAACGCAGTAGAGGGAATGCCGCCTa                                                                      <  2:140507/100‑1 (MQ=255)
                                aGCGTTCTCTTCATAACAATCATCGGCAACTTTACTGTCGCGCGGATGTGCGGTAACAGACCCGCGCCTTTCAACGCAGTAGAGGGAATGCCGCCTACgg                                                                   <  2:1039387/100‑1 (MQ=255)
                                 gCGTGCGCTTCATAACAATCATCGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCGCCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGc                                                                  >  1:1769590/1‑100 (MQ=255)
                                  cGTGCGCTTCATAACAATCATCGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCGCCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCa                                                                 <  1:1840616/100‑1 (MQ=255)
                                    tGCGCTTCATAACAATCATCGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCACCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATc                                                               >  2:1519634/1‑100 (MQ=255)
                                       gcTTCATAACAATCATCGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCGCCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGc                                                            <  2:1862605/100‑1 (MQ=255)
                                       gcTTCATAACAATCATCGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACACGCGCCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGc                                                            >  2:2187560/1‑100 (MQ=255)
                                        cTTCATAACAATCATCGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGTGCCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCa                                                           >  2:1527281/1‑100 (MQ=255)
                                           cATAACAATCATCGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCGCCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACcg                                                        <  1:2378252/100‑1 (MQ=255)
                                           cATAACAATCATCGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCGCCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACcg                                                        <  1:1390245/100‑1 (MQ=255)
                                                cAATCATCGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCGCCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTc                                                   >  1:2113483/1‑100 (MQ=255)
                                                cAATCATCGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCACCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCt                                                    <  1:638722/99‑1 (MQ=255)
                                                cAATCATCGGCAACTTTACCGTAACGCAGATATGCCGTAACAGACCTACGCCTTCCAAAGCAGTAGAGCTAATGCCGCCTACGGCATAGGCACcgc                                                       >  2:40500/1‑96 (MQ=255)
                                                  atcatcGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCGCCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCaa                                                 <  1:2390058/100‑1 (MQ=255)
                                                  atcatcGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCGCCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCaa                                                 >  1:368736/1‑100 (MQ=255)
                                                    catcGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCACCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAAtt                                               >  2:1066080/1‑100 (MQ=255)
                                                       cGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCGCCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAATTTCt                                            >  1:498263/1‑100 (MQ=255)
                                                        ggCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCGCCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAATTTCTg                                           >  1:1321207/1‑100 (MQ=255)
                                                         gCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCGCCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAATTTCTGc                                          >  2:1324305/1‑100 (MQ=255)
                                                         gCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCACCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAATTTCTGc                                          <  1:1577275/100‑1 (MQ=255)
                                                          cAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCGCCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAATTTCTGcc                                         <  1:2046681/100‑1 (MQ=255)
                                                          cAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCGCCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAATTTCTGc                                          >  1:681889/1‑99 (MQ=255)
                                                          cAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCACCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAATTTCTGcc                                         >  1:1281367/1‑100 (MQ=255)
                                                            aCTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCGCCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAATTTCTGCCAg                                       <  2:314413/100‑1 (MQ=255)
                                                             cTTTACGGTCGCGGGGATGTGCGGTAACAGACCCGCGCCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAATTTCTGcc                                         <  1:1204211/97‑1 (MQ=255)
                                                             cTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCGCCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAATTTCTGCCAg                                       >  1:857820/1‑99 (MQ=255)
                                                             cTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCGCCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAATTTCTGCCAGc                                      <  2:2020924/100‑1 (MQ=255)
                                                             cTTTACGGTCGCGCGGATGTGCGGTAACAGAACCGCGCCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCAGCGGCACCGCGCTCAATTTCTGCCAGc                                      >  2:1306915/1‑100 (MQ=255)
                                                              tttACGGTCGCGCGGATGTGCGGTAACAGACCCGCGCCTTGCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAATTTCTGCCAg                                       <  1:1959193/98‑1 (MQ=255)
                                                               ttACGGTCGCGCGGATGTGCGGTAACAGACCCGCGCCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAATTTCTGCCAGCat                                    <  1:90942/100‑1 (MQ=255)
                                                               ttACGGTCGCGCGGATGTGCGGTAACAGACCCGCGCCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAATTTCTGCCAGCat                                    <  1:959989/100‑1 (MQ=255)
                                                                  cGGTCGCGCGGATGTGCGGTAACAGACCCGCGCCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAATTTCTGCCAGCAtata                                 >  1:1637182/1‑100 (MQ=255)
                                                                   ggTCGCGCGGATGTGCGGTAACAGACCCGCGCCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAATTTCTGCCAGCATATAg                                >  1:1682249/1‑100 (MQ=255)
                                                                   ggTCGCGCGGATGTGCGGTAACAGACCCGCGCCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAATTTCTGCCAGCATATAg                                <  1:1332042/100‑1 (MQ=255)
                                                                   ggTCGCGCGGATGTGCGGTAACAGACCCGCGCCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAATTTCTGCCAGCATATAg                                >  1:1565055/1‑100 (MQ=255)
                                                                    gTCGCGCGGATGTGCGGTAACAGACCCGCGCCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAATTTCTGCCAGCATATAGc                               <  1:691996/100‑1 (MQ=255)
                                                                    gTCGCGCGGATGTGCGGTAACAGACCCGCGCCTTTCAACGCAGTAGAGGGAATGCCGACTACGGCATCGGCACCGCGCTCAATTTCTGCCAGCATATAGc                               >  2:870825/1‑100 (MQ=255)
                                                                    gTCGCGCGGATGTGCGGTAACAGACCCGCACCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCGCAATTTCTGCCAGCATATAGc                               <  1:1165560/100‑1 (MQ=255)
                                                                            gATGTGCGGTAACAGACCCGCGCCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAATTTCTGCCAGCATATAGCCCATTCCg                       >  1:874027/1‑100 (MQ=255)
                                                                             aTGTGCGGTAACAGACCCGCGCCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAATTTCTGCCAGCATATAGCCCATTCCg                       <  1:1562315/99‑1 (MQ=255)
                                                                             aTGTGCGGTAACAGACCCGCGCCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAATTTCTGCCAGCATATAGCCCATTCCGt                      <  1:2051614/100‑1 (MQ=255)
                                                                              tgtgCGGTAACAGACCCGCACCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAATTTCTGCCAGCATATAGCCCATTCCGTc                     <  1:1924716/100‑1 (MQ=255)
                                                                                ttcGGTAGCAGACCCGCGCCTTTCAATGCAGTCGAGGGAATTCCGCCTACGGCATCGGCACGGCGCTCAATTTCTGCCGTCATATAGCCCAGTCCGTCtt                   <  2:2255033/98‑1 (MQ=255)
                                                                                tgCGGTAACAGACCCGCGCCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAATTTCTGCCAGCATATAGCCCATTCCGTCtt                   >  2:297634/1‑100 (MQ=255)
                                                                                tgCGGTAACAGACCCGCGCCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAATTTCTGCCAGCATATAGCCCATTCCGTCtt                   <  1:1703282/100‑1 (MQ=255)
                                                                                tgCGGTAACAGACCCGCGCCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAATTTCTGCCAGCATATAGCCCATTCCGTCtt                   >  1:1973742/1‑100 (MQ=255)
                                                                                 gCGGTAACAGACCCGCACCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAATTTCTGCCAGCATATAGCCCATTCCGTCTTg                  >  1:591736/1‑100 (MQ=255)
                                                                                 gCGGTAACAGACCCGCACCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAATTTCTGCCAGCATATAGCCCATTCCGTCTTg                  >  1:1580913/1‑100 (MQ=255)
                                                                                  cGGTAACAGACCCGCACCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAATTTCTGCCAGCATATAGCCCATTCCGTCTTGa                 >  2:787156/1‑100 (MQ=255)
                                                                                   ggTAACAGACCCGCACCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAATTTCTGCCAGCATATAGCCCATTCCGTCTTGAt                >  2:238438/1‑100 (MQ=255)
                                                                                    gTAACAGACCCGCGCCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAATTTCTGCCAGCATATAGCCCATTCCGTCTTGATc               <  1:1364484/100‑1 (MQ=255)
                                                                                     tAACAGACCCGCGCCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAATTTCTGCCAGCATATAGCCCATTCCGTCTTGATCg              <  2:1663501/100‑1 (MQ=255)
                                                                                      aaCAGACCCGCGCCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAATTTCTGCCAGCATATAGCCCATTCCGTCTTGATCgg             >  1:524262/1‑100 (MQ=255)
                                                                                      aaCAGACCCGCGCCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAATTTCTGCCAGCATATAGCCCATTCCGTCTTGATCgg             >  2:1950267/1‑100 (MQ=255)
                                                                                        cAGACCCGCGCCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAATTTCTGCCAGCATATAGCCCATTCCGTCTTGATCGGGc           >  1:2371792/1‑100 (MQ=255)
                                                                                        cAGACCCGCACCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAATTTCTGCCAGCATATAGCCCATTCCGTCTTGATCGGGc           >  1:1338494/1‑100 (MQ=255)
                                                                                        cAGACCCGCACCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAATTTCTGCCAGCATATAGCCCATTCCGTCTTGATCGGGc           <  1:1556966/100‑1 (MQ=255)
                                                                                         aGACCCGCGCCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAATTTCTGCCAGCATATAGCCCATTCCGTCTTGATCGGGc           <  1:990812/99‑1 (MQ=255)
                                                                                         aGACCCGCGCCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAATTTCTGCCAGCATATAGCCCATTCCGTCTTGATCGGGCg          <  2:1135329/100‑1 (MQ=255)
                                                                                         aGACCCGCGCCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAATTTCTGCCAGCATATAGCCCATTCCGTCTTGATCGGGCg          <  2:666970/100‑1 (MQ=255)
                                                                                          gACCCGCGCCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAATTTCTGCCAGCATATAGCCCATTCCGTCTTGATCGGGCgg         <  1:1452315/100‑1 (MQ=255)
                                                                                          gACCCGCACCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAATTTCTGCCAGCATATCGCCCATTCCGTCTTGATCGGGCgg         >  1:622233/1‑100 (MQ=255)
                                                                                            cccGCGCCTTTCAACGCAGTAGCGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAATTTCTGCCAGCATATAGCCCATTCCGTCTTGATCGGGCGGaa       >  2:1390273/1‑100 (MQ=255)
                                                                                            cccGCGCCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAATTTCTGCCAGCATATAGCCCATTCCGTCTTGATCGGGCGGaa       >  2:1721792/1‑100 (MQ=255)
                                                                                            cccGCGCCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAATTTCTGCCAGCATATAGCCCATTCCGTCTTGATCGGGCGGaa       >  1:2338223/1‑100 (MQ=255)
                                                                                            cccGCACCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAATTTCTGCCAGCATATAGCCCATTCCGTCTTGATCGGGCGGaa       >  2:190492/1‑100 (MQ=255)
                                                                                              cgcgCCGTTCAACGCAGTAGAGGGAATGCGGCCTACGGCATCGGCACCGCGCTCAAGTTCTGCCAGCATATAGCCCATTCCGTCTTGATCGGGCGGAACa     <  2:34503/100‑1 (MQ=255)
                                                                                               gcgcCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAATTTCTGCCAGCATATAGCCCATTCCGTCTTGATCGGGCGGAACAt    <  2:528831/100‑1 (MQ=255)
                                                                                                 gcCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAATTTCTGCCAGCATATAGCCCATTCCGTCTTGATCGGGCGGAACATAg  >  1:2064174/1‑100 (MQ=255)
                                                                                                 |                                                                                                   
TAATAAACGGTGCGCCACCCAGGAGCTGCTGTAGCGTGCGCTTCATAACAATCATCGGCAACTTTACGGTCGCGCGGATGTGCGGTAACAGACCCGCGCCTTTCAACGCAGTAGAGGGAATGCCGCCTACGGCATCGGCACCGCGCTCAATTTCTGCCAGCATATAGCCCATTCCGTCTTGATCGGGCGGAACATAG  >  NC_000913/383380‑383576

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: