Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,100,233 G→A 100% intergenic (+39/‑63) ycdY → / → ycdZ redox enzyme maturation protein (REMP) chaperone for YcdX/DUF1097 family inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,100,2330GA99.4% 533.9 / ‑4.7 159intergenic (+39/‑63)ycdY/ycdZredox enzyme maturation protein (REMP) chaperone for YcdX/DUF1097 family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (79/79);  total (79/80)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.20e-01

Only 100 of 162 total aligned reads displayed.

TAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTGTTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGAACATACTTCTCTCTATTGCAATCACAACA  >  NC_000913/1100136‑1100328
                                                                                                 |                                                                                               
   cccGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATttattgtt                                                                                            >  2:357614/1‑100 (MQ=255)
      gcgATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCa                                                                                         >  2:603421/1‑100 (MQ=255)
      gcgATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCa                                                                                         <  1:1466165/100‑1 (MQ=255)
         aTGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATaa                                                                                      <  1:806181/100‑1 (MQ=255)
           gCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATACTTTATTGGTGCATaac                                                                                     >  1:3888030/1‑98 (MQ=255)
           gCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATaaaa                                                                                    >  1:2076268/1‑100 (MQ=255)
            ccATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAAt                                                                                   >  1:2170983/1‑100 (MQ=255)
             cATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAAtg                                                                                  >  2:3840449/1‑100 (MQ=255)
             cATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAAtg                                                                                  >  1:2542721/1‑100 (MQ=255)
              aTTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAAtgt                                                                                 >  2:2473456/1‑100 (MQ=255)
               ttAGTGCAATGTGGGACGAGCTTGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCATTCATAACTGCAACGGATAATTTATTGTTGCATAAAAtgtg                                                                                <  1:4259148/100‑1 (MQ=255)
               ttAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAAtgtg                                                                                <  2:2917696/100‑1 (MQ=255)
                 aGTACAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATaaaa                                                                                    >  1:3037690/1‑94 (MQ=255)
                  gTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGc                                                                             >  1:3224712/1‑100 (MQ=255)
                   tGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCt                                                                            <  1:1403765/100‑1 (MQ=255)
                    gCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTc                                                                           >  2:2684653/1‑100 (MQ=255)
                    gCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTc                                                                           >  2:487983/1‑100 (MQ=255)
                        tgtgGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGAtc                                                                       <  2:3325919/100‑1 (MQ=255)
                         gtgGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGAtct                                                                      <  2:99265/100‑1 (MQ=255)
                          tgGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGAtctc                                                                     >  1:1053285/1‑100 (MQ=255)
                           gggACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTCCGTCACTCATAACTGCAAGGGATAATTTATTGTTGCATAAAATGTGTGCTCGATAtca                                                                    >  2:1249689/1‑100 (MQ=255)
                           gggACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCtca                                                                    >  2:1342931/1‑100 (MQ=255)
                            ggACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCtcat                                                                   <  1:605986/100‑1 (MQ=255)
                            ggACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCtcat                                                                   <  2:326059/100‑1 (MQ=255)
                            ggACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCtcat                                                                   <  2:2978569/100‑1 (MQ=255)
                             gACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCtcatt                                                                  >  2:644700/1‑100 (MQ=255)
                                gAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCtcattcat                                                               >  2:747975/1‑100 (MQ=255)
                                 aGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATg                                                              <  1:618393/100‑1 (MQ=255)
                                  gCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCtcattcat                                                               >  1:1856748/1‑98 (MQ=255)
                                   cTGGAAGTAGATTCTGAAGAATAAGTGTGATCGACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAATGCGTGCTCGATCTCATTCATGGc                                                            <  1:3630728/100‑1 (MQ=255)
                                   cTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGc                                                            <  2:4271235/100‑1 (MQ=255)
                                    tGGAAGAAGCTTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGCTCTCATTCATGGcc                                                           >  2:751095/1‑100 (MQ=255)
                                     ggAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGcc                                                           <  1:1912138/99‑1 (MQ=255)
                                     ggAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCcg                                                          >  1:327874/1‑100 (MQ=255)
                                     ggAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCcg                                                          >  1:3885452/1‑100 (MQ=255)
                                      gaataaGATTCTTAAGTATTAGTGTGTTCTACGTTACTCATAACTGCAACGCTTAATTTATTGTTTTATAAATTGTTTGTTCGAGCTCATTCATGGCcgc                                                         <  1:2686517/100‑1 (MQ=255)
                                      gaagaaGATTCTGAAGAATAAGTGTGATCTACGTCACCCATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCcgc                                                         >  1:3660584/1‑100 (MQ=255)
                                       aagaagATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCcgcg                                                        >  1:1768328/1‑100 (MQ=255)
                                       aagaagATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCcgcg                                                        >  2:595496/1‑100 (MQ=255)
                                        agaagaTTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGt                                                       >  1:1293914/1‑100 (MQ=255)
                                          aagaTTCTGAAGAATAAGTGTTATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAATGCGTGCTCGATCTCATTCATGGCCGCGtt                                                      <  1:558124/99‑1 (MQ=255)
                                           agaTTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGtttt                                                    >  1:1237559/1‑100 (MQ=255)
                                            naTTCTGAAGAATAAGTGAGATCTACCTCACTCATAACCGCAACGGATAATCTAGTGTCGCAGAAAATGTGTGCTCGAACGCATTCATGGCCGAGtttt                                                    >  1:1544750/2‑99 (MQ=255)
                                               tctgcggaaTAAGTGTGATCTACGTCTCTCATATGTGTATCGGATAAGTTATTGTTGCACTAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGc                                                <  1:1161910/94‑1 (MQ=255)
                                                cTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCt                                               >  1:1114977/1‑100 (MQ=255)
                                                 tGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTa                                              >  2:604769/1‑100 (MQ=255)
                                                  gaagaaTAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTAt                                             <  2:2316100/100‑1 (MQ=255)
                                                    agaaTAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTGTTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCa                                           <  2:2735438/100‑1 (MQ=255)
                                                    agaaTAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCa                                           <  2:3742209/100‑1 (MQ=255)
                                                     gaaTAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATc                                            <  1:336022/98‑1 (MQ=255)
                                                       aTAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATgc                                        >  1:1749142/1‑100 (MQ=255)
                                                          aGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGcgcg                                     <  2:1333442/100‑1 (MQ=255)
                                                           gtgtGATCTACGTCACTCATAACTGCACCGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCACGTTTTCTGCTATCATGcgcg                                     <  1:417571/99‑1 (MQ=255)
                                                           gtgtGATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCgg                                    <  1:3455296/100‑1 (MQ=255)
                                                             gtgATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCa                                  <  2:588363/100‑1 (MQ=255)
                                                             gtgATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCa                                  <  2:2274778/100‑1 (MQ=255)
                                                               gATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATg                                >  1:1298730/1‑100 (MQ=255)
                                                                 tCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGaa                              <  2:1016917/100‑1 (MQ=255)
                                                                 tCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGaa                              >  2:1703439/1‑100 (MQ=255)
                                                                  cTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGAAc                             <  1:2362203/100‑1 (MQ=255)
                                                                   tACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGAACa                            <  1:1215051/100‑1 (MQ=255)
                                                                    aCGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGAACAt                           <  2:116577/100‑1 (MQ=255)
                                                                    aCGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGAACAt                           >  1:305400/1‑100 (MQ=255)
                                                                    aCGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGAACAt                           <  1:357612/100‑1 (MQ=255)
                                                                      gTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAATGGGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGAACATAc                         >  2:3730500/1‑100 (MQ=255)
                                                                       tCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGAACATACt                        <  1:1015391/100‑1 (MQ=255)
                                                                       tCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGAACATACt                        <  2:2262531/100‑1 (MQ=255)
                                                                        cACTCATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGAACATACtt                       >  2:2063970/1‑100 (MQ=255)
                                                                         aCTCATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGAACATACTtc                      <  2:2435661/100‑1 (MQ=255)
                                                                         aCTCATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGAACATACTtc                      >  2:3967087/1‑100 (MQ=255)
                                                                         aCTCATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGAACATACTtc                      <  2:3088098/100‑1 (MQ=255)
                                                                          cTCATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGAACATACTtct                     >  2:1784195/1‑100 (MQ=255)
                                                                            cATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGAACATACTtctct                   >  1:821515/1‑100 (MQ=255)
                                                                            cATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGAACATACTtctct                   >  2:3378912/1‑100 (MQ=255)
                                                                            cATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGAACATACTtctc                    <  2:1298734/99‑1 (MQ=255)
                                                                             aTAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGAACATACTtctctc                  >  1:4227025/1‑100 (MQ=255)
                                                                              tAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGAACATACTTctctct                 <  1:1920463/100‑1 (MQ=255)
                                                                              tAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGAACATACTTctctct                 >  2:465243/1‑100 (MQ=255)
                                                                                aCTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGAACATACTTCTCTCTAt               <  2:3727852/100‑1 (MQ=255)
                                                                                aCTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGAACATACTTCTCTCTAt               >  1:1672941/1‑100 (MQ=255)
                                                                                aCTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGAACATACTTCTCTCTAt               >  1:4259400/1‑100 (MQ=255)
                                                                                aCTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGAACATACTTCTCTCTAt               <  2:1635276/100‑1 (MQ=255)
                                                                                 cTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGAACATACTTCTCTCTAtt              <  1:4013499/100‑1 (MQ=255)
                                                                                  tGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGAACATACTTCTCTCTATTg             <  2:112804/100‑1 (MQ=255)
                                                                                  tGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGAACATACTTCTCTCTATTg             <  1:1384919/100‑1 (MQ=255)
                                                                                  tGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGAACATACTTCTCTCTATTg             <  1:1340599/100‑1 (MQ=255)
                                                                                   gCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGAACATACTTCTCTCTATTGc            >  1:3194196/1‑100 (MQ=255)
                                                                                   gCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGAACATACTTCTCTCTATTGc            <  1:200326/100‑1 (MQ=255)
                                                                                       cGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGAACATACTTCTCTCTATTGCAATc        <  2:1376958/100‑1 (MQ=255)
                                                                                        ggATAATTTATTGTTGCATAAAATGTGTGCTCGGTCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGAACATACTTCTCTCTATTGCAATca       <  1:2518436/100‑1 (MQ=255)
                                                                                        ggATAATTTATTGTTGCATAAAATGTGTGCTCGATCTGATTCATGGCAGCGTTTTCTGCTATCATGCGCGGCATGAACATACTTCTCTCTATTGCaa          <  1:1100704/97‑1 (MQ=255)
                                                                                        ggATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATTATGCGCGGCATGAACATACTTCTCTCTATTGCAATca       <  1:2283127/100‑1 (MQ=255)
                                                                                         gATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGAACATACTTCTCTCTATTGCAATcac      <  1:2126702/100‑1 (MQ=255)
                                                                                          aTAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGAACATACTTCTCTCTATTGCAATcaca     <  1:3130647/100‑1 (MQ=255)
                                                                                            aaTTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGAACATACTTCTCTCTATTGCAATCacaac   >  1:3814233/1‑100 (MQ=255)
                                                                                            aaTTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGAACATACTTCTCTCTATTGCAATCacaac   <  2:1435209/100‑1 (MQ=255)
                                                                                             aTTTATTGTTGCATAAAATGTGTGCTCGATGTCCTTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGAACATACTTCTCTCTATTGAAATCacaaca  >  1:2801730/1‑100 (MQ=255)
                                                                                             aTTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGAACATACTTCTCTCTATTGCAATCacaaca  <  2:162157/100‑1 (MQ=255)
                                                                                             aTTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGAACATACTTCTCTCTATTGCAATCacaaca  >  2:508931/1‑100 (MQ=255)
                                                                                             aTTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTCTCTGCTATCATGCGCGGCATGACCATACTTCTCTCTATTGCAATCacaaca  >  2:3646431/1‑100 (MQ=255)
                                                                                                 |                                                                                               
TAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTGTTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGAACATACTTCTCTCTATTGCAATCACAACA  >  NC_000913/1100136‑1100328

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: