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.1% 351.5 / ‑4.6 108intergenic (+39/‑63)ycdY/ycdZredox enzyme maturation protein (REMP) chaperone for YcdX/DUF1097 family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (48/58);  total (50/58)
Fisher's exact test for biased strand distribution p-value = 4.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.35e-01

Only 100 of 119 total aligned reads displayed.

TAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTGTTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGAACATACTTCTCTCTATTGCAATCACAACAG  >  NC_000913/1100136‑1100329
                                                                                                 |                                                                                                
tAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATttatt                                                                                                <  1:1323329/100‑1 (MQ=255)
tAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATttatt                                                                                                >  2:2727521/1‑100 (MQ=255)
 aaCCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATttattg                                                                                               <  2:2902614/100‑1 (MQ=255)
  aCCCGAGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATttattgt                                                                                              <  2:2557699/100‑1 (MQ=255)
    ccGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATttattgttg                                                                                            <  2:3653594/100‑1 (MQ=255)
    ccGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATttattgttg                                                                                            <  2:3296938/100‑1 (MQ=255)
    ccGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATttattgttg                                                                                            <  1:2421388/100‑1 (MQ=255)
    ccGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATttattgttg                                                                                            <  1:1909915/100‑1 (MQ=255)
     cgcgATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATttattgtt                                                                                             <  2:3968971/98‑1 (MQ=255)
      gagATGCCATTAGTGCAATGTGGGACGAGCTCGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCATTCATAACTGCAACGGATAATttattgttg                                                                                            <  2:934578/96‑1 (MQ=255)
         aTGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATaa                                                                                       <  2:2163892/100‑1 (MQ=255)
         aTGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATaa                                                                                       >  1:3293108/1‑100 (MQ=255)
          tGCCATTAGTGCAATGTGGGACGAGCTCGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATa                                                                                        >  2:1079644/1‑98 (MQ=255)
           cgatctAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATaaaa                                                                                     <  1:748369/95‑1 (MQ=255)
            ccATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAAt                                                                                    >  1:136064/1‑100 (MQ=255)
             cATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATATAAtg                                                                                   >  2:3796041/1‑100 (MQ=255)
             cATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAAtg                                                                                   <  2:2980652/100‑1 (MQ=255)
                 aatGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAAtgtgtg                                                                               <  1:2496611/98‑1 (MQ=255)
                 aGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATaaaaagatcg                                                                               >  2:748365/1‑94 (MQ=255)
                   tGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCt                                                                             <  1:3127670/100‑1 (MQ=255)
                    gCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTc                                                                            >  1:3578732/1‑100 (MQ=255)
                    gCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTc                                                                            <  2:2675825/100‑1 (MQ=255)
                    gCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTc                                                                            >  2:3350703/1‑100 (MQ=255)
                     cAATTTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGc                                                                              <  2:2058125/97‑1 (MQ=255)
                     cAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTGTTGTTGCATAAAATGTGTGCTCg                                                                           >  1:3915026/1‑100 (MQ=255)
                     cAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCg                                                                           >  1:3524997/1‑100 (MQ=255)
                        tgtgGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGGAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGa                                                                          >  2:932272/1‑98 (MQ=255)
                         gtgGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACAGATAATTTATTGTTGCATAAAATGTGTGCTCGAtct                                                                       <  2:2686762/100‑1 (MQ=255)
                            ggACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCtcat                                                                    >  2:1056992/1‑100 (MQ=255)
                                 aGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATg                                                               >  2:630464/1‑100 (MQ=255)
                                 aGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATg                                                               >  1:3489912/1‑100 (MQ=255)
                                 aGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATg                                                               >  2:2143850/1‑100 (MQ=255)
                                  ggttgaagaaGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATgg                                                              <  2:3847169/96‑1 (MQ=255)
                                   cTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGc                                                             <  2:1630763/100‑1 (MQ=255)
                                   cTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGc                                                             <  2:1421863/100‑1 (MQ=255)
                                    tGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGc                                                             <  2:2539929/99‑1 (MQ=255)
                                      gaagaaGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCcgc                                                          <  1:3482667/100‑1 (MQ=255)
                                      gaagaaGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCcgc                                                          >  1:2124070/1‑100 (MQ=255)
                                       aagaagATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCcgcg                                                         <  1:3231912/100‑1 (MQ=255)
                                       aagaagATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCcgcg                                                         >  2:1109269/1‑100 (MQ=255)
                                       aagaagATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCcgcg                                                         <  1:426530/100‑1 (MQ=255)
                                       aagaagATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCcgcg                                                         <  2:3329330/100‑1 (MQ=255)
                                       aagaagATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCcgcg                                                         <  2:940229/100‑1 (MQ=255)
                                       aagaagATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCcgcg                                                         <  1:1159155/100‑1 (MQ=255)
                                        agaagaTTCTGAAGAATAAGTGTGATCTGCGGCACTCATGACTGCGACGGCTAATTTATTGTTGCATAAAATGTGTGCTGGATCTCATTCATGGGcgcg                                                         <  2:2912660/99‑1 (MQ=255)
                                        agaagaTTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGt                                                        >  1:995670/1‑100 (MQ=255)
                                          aagaTTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGttt                                                      >  1:1185831/1‑100 (MQ=255)
                                          aagaTTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCCCGATCTCATTCATGGCCGCGttt                                                      >  2:3646822/1‑100 (MQ=255)
                                           agaTTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGtttt                                                     >  1:3615090/1‑100 (MQ=255)
                                           agaTTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGtttt                                                     >  2:2100481/1‑100 (MQ=255)
                                            gaTTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTc                                                    >  1:1611026/1‑100 (MQ=255)
                                            gaTTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTc                                                    <  2:2273123/100‑1 (MQ=255)
                                              ttCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTTCATAAAATGTGTGCTCGTTCTCGTTCATGGCCGCGTTTTCTg                                                  >  1:3323197/1‑100 (MQ=255)
                                              ttCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTg                                                  <  2:4126023/100‑1 (MQ=255)
                                               tCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGc                                                 >  2:1991368/1‑100 (MQ=255)
                                                cTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCt                                                >  2:592699/1‑100 (MQ=255)
                                                 taaagaaTAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGc                                                 <  2:3368060/96‑1 (MQ=255)
                                                 tGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTa                                               >  2:2883653/1‑100 (MQ=255)
                                                 tGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTa                                               <  2:2005416/100‑1 (MQ=255)
                                                 tGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAAGTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTa                                               >  1:2483455/1‑100 (MQ=255)
                                                  gaagaaTAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTAt                                              >  2:254933/1‑100 (MQ=255)
                                                   aagaaTAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATc                                             <  1:884570/100‑1 (MQ=255)
                                                    agaaTAAGTGTGATCTATGTCACTCATAACTGCAAATGATAATTTATTGTTTCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCa                                            <  1:3348253/100‑1 (MQ=255)
                                                    agaaTAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCa                                            <  2:3743205/100‑1 (MQ=255)
                                                    agaaTAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCa                                            <  2:729774/100‑1 (MQ=255)
                                                      aaTAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTGATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATg                                          <  2:1523945/100‑1 (MQ=255)
                                                       aTAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATgc                                         <  2:3323172/100‑1 (MQ=255)
                                                       aTAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATgc                                         <  1:1827651/100‑1 (MQ=255)
                                                       aTAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATgc                                         <  1:3071079/100‑1 (MQ=255)
                                                       aTAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATgc                                         >  2:1977923/1‑100 (MQ=255)
                                                           gtgtGATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGcgcg                                      <  1:1227208/99‑1 (MQ=255)
                                                           gtgtGATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCgg                                     >  1:3413830/1‑100 (MQ=255)
                                                              tgATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCAt                                  <  2:478619/100‑1 (MQ=255)
                                                               gATCTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATg                                 >  1:2495015/1‑100 (MQ=255)
                                                                  cTACGTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGAAc                              <  2:2122038/100‑1 (MQ=255)
                                                                      gTCACTCATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATAATGCGCGGCATGAACATAc                          >  2:1611037/1‑100 (MQ=255)
                                                                         aCTCATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGAACATACTtc                       >  1:863060/1‑100 (MQ=255)
                                                                          cTCATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGAACATACTtct                      >  2:428426/1‑100 (MQ=255)
                                                                            cATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGAACATACTtctct                    >  1:515613/1‑100 (MQ=255)
                                                                            cATAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGAACATACTtctct                    >  1:3103024/1‑100 (MQ=255)
                                                                             aTAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGAACATACTtctctc                   <  2:3481754/100‑1 (MQ=255)
                                                                             aTAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGAACATACTtctctc                   >  1:1803987/1‑100 (MQ=255)
                                                                             aTAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGAACATACTtctctc                   <  1:2667180/100‑1 (MQ=255)
                                                                              tAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCTCGGCATGAACATACTTctctct                  <  2:2421432/100‑1 (MQ=255)
                                                                              tAACTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGAACATACTTctctct                  <  1:2894826/100‑1 (MQ=255)
                                                                                aCTGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGAACATACTTCTCTCTAt                <  1:3608895/100‑1 (MQ=255)
                                                                                 cTGCAACGGATAATTTATTGTTGCATAAAAGGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGAACATACTTCTCTCTa                 <  1:3943561/98‑1 (MQ=255)
                                                                                  tGCAACGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGAACATACTTCTCTCTATTg              <  2:4046293/100‑1 (MQ=255)
                                                                                   ncAACGGATAATCTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGAACATACTTCTCTCTATTGc             >  1:2737483/2‑100 (MQ=255)
                                                                                      aCGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGAACATACTTCTCTCTATTGCAAt          <  1:3486764/100‑1 (MQ=255)
                                                                                      aCGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGAACATACTTCTCTCTATTGCAAt          >  1:1964374/1‑100 (MQ=255)
                                                                                       cGGATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGAACATACTTCTCTCTATTGCAATc         >  2:628140/1‑100 (MQ=255)
                                                                                        ggATAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGAACATACTTCTCTCTATTGCAATca        <  2:2106933/100‑1 (MQ=255)
                                                                                          aTAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGAACATACTTCTCTCTATTGCAATcaca      >  1:795428/1‑100 (MQ=255)
                                                                                           tAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGAACATACTTCTCTCTATTGCAATcac       <  2:123100/98‑1 (MQ=255)
                                                                                           tAATTTATTGTTGCATAAAATGTGTGCTCGATCTCATTAATGGCCGCGTTTTCTGCTATCATGCGCGGCATGAACATACTTCTCTCTCTTGCACTCacaa     >  2:3628515/1‑100 (MQ=255)
                                                                                             aTTTATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGAACATACTTCTCTCTATTGCAATCacaaca   <  1:2959845/100‑1 (MQ=255)
                                                                                              tttATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGAACATACTTCTCTCTATTGCagatcggaaga  >  1:1806233/1‑90 (MQ=255)
                                                                                              tttATTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGAACATACTTCTCTCTATTGCAATCACAACAg  >  2:3648067/1‑100 (MQ=255)
                                                                                              tttATTATTGCATAAAATGTGTGCTCGATCTCATTCATTGCCGCGTATTCTGCTATCATGACCGACATGAACATACTTCTCTCGATTGCAATGacaaca   >  2:1492254/1‑99 (MQ=255)
                                                                                                 |                                                                                                
TAACCCGCGATGCCATTAGTGCAATGTGGGACGAGCTGGAAGAAGATTCTGAAGAATAAGTGTGATCTACGTCACTCATAACTGCAACGGATAATTTGTTGTTGCATAAAATGTGTGCTCGATCTCATTCATGGCCGCGTTTTCTGCTATCATGCGCGGCATGAACATACTTCTCTCTATTGCAATCACAACAG  >  NC_000913/1100136‑1100329

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: