Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 981,566 C→T R433C (CGC→TGC)  ycbB → putative amidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913981,5660CT99.7% 1047.7 / NA 295R433C (CGC→TGC) ycbBputative amidase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (141/153);  minor base G (1/0);  total (142/153)
Fisher's exact test for biased strand distribution p-value = 4.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 300 total aligned reads displayed.

GGTCGCCCCGATCGCAAAACGCCGATGATGAGCAGTGCCCTTAACAACGTAGTGGTAAACCCGCCGTGGAACGTACCTCCAACTCTGGCACGCAAAGATATTCTGCCAAAAGTGCGCAACGATCCGGGATATCTCGAAAGCCATGGCTATACGGTGATGCGCGGCTGGAACAGCAGAGAA  >  NC_000913/981476‑981655
                                                                                           |                                                                                         
ggtcgCCCCGATCGCAAAACGCCGATGATGAGCAGTGCCCTTAACAACGTAGTGGTAAACCCGCCGTGGAACGTACCTCCAACTCTGGCATGc                                                                                         >  1:7040480/1‑93 (MQ=255)
 gtcgCCCCGATCGCAAAACGCCGATGATGAGCAGTGCCCTTAACAACGTAGTGGTAAACCCGCCGTGGAACGTACCTCCAACTCTGGCATGCa                                                                                        <  1:5350563/93‑1 (MQ=255)
  tcgCCCCGATCGCAAAACGCCGATGATGAGCAGTGCCCTTAACAACGTAGTGGTAAACCCGCCGTGGAACGTACCTCCAACTCTGGCATGCaa                                                                                       >  2:876400/1‑93 (MQ=255)
      cccGATCGCAAAACGCCGATGATGAGCAGTGCCCTTAACAACGTAGTGGTAAACCCGCCGTGGAACGTACCTCCAACTCTGGCATGCAAAGat                                                                                   >  1:1293212/1‑93 (MQ=255)
       ccGATCGCAAAACGCCGATGATGAGCAGTGGCCTTAACAACGTAGTGGTAAACCCGCCGTGGAACGTACCTCCAACTCTGGCATGCAAAGata                                                                                  >  1:2691545/1‑93 (MQ=255)
       ccGATCGCAAAACGCCGATGATGAGCAGTGCCCTTAACAACGTAGTGGTAAACCCGCCGTGGAACGTACCTCCAACTCTGGCATGCAAAGata                                                                                  <  1:4286981/93‑1 (MQ=255)
       ccGATCGCAAAACGCCGATGATGAGCAGTGCCCTTAACAACGTAGTGGTAAACCCGCCGTGGAACGTACCTCCAACTCTGGCATGCAAAGata                                                                                  >  1:670958/1‑93 (MQ=255)
       ccGATCGCAAAACGCCGATGATGAGCAGTGCCCTTAACAACGTAGTGGTAAACCCGCCGTGGAACGTACCTCCAACTCTGGCATGCAAAGata                                                                                  >  1:6646426/1‑93 (MQ=255)
        cGATCGCAAAACGCCGATGATGAGCAGTGCCCTTAACAACGTAGTGGTAAACCCGCCGTGGAACGTACCTCCAACTCTGGCATGCAAAGatat                                                                                 <  2:604284/93‑1 (MQ=255)
         gATCGCAAAACGCCGATGATGAGCAGTGCCCTTAACAACGTAGTGGTAAACCCGCCGTGGAACGTACCTCCAACTCTGGCATGCAAAGATAtt                                                                                >  2:6543928/1‑93 (MQ=255)
          aTCGCAAAACGCCGATGATGAGCAGTGCCCTTAACAACGTAGTGGTAAACCCGCCGTGGAACGTACCTCCAACTCTGCCATGCAAAGATATTc                                                                               >  1:4538704/1‑93 (MQ=255)
          aTCGCAAAACGCCGATGATGAGCAGGGCCCTTAACAACGTAGTGGTAAACCCGCCGTGGAACGTACCTCCAACTCTCGCATGCAAAGatatac                                                                               >  2:1575262/1‑91 (MQ=255)
            cGCAAAACGCCGATGATGAGCAGTGCCCTTAACAACGTAGTGGTAAACCCGCCGTGGAACGTCCCTCCAACTCTGGCATGCAAAGATATTCTg                                                                             <  2:6714399/93‑1 (MQ=255)
             gCAAAACGCCGATGATTAGCAGTGCCCTTAACAACGTAGTGGTAATCCCTTCGTGGAACGTACCTCCAACTCTGGCATGCTAAGATAATCTGc                                                                            >  2:7836448/1‑93 (MQ=255)
               aaaaCGCCGATGATGAGCAGTGCCCTTAACAACGTAGTGGTAAACCCGCCGTGGAACGTACCTCCAACTCTGGCATGCAAAGATATTCTGCCa                                                                          >  1:41207/1‑93 (MQ=255)
                aaaCGCCGATGATGAGCAGTGCCCTTAACAAGGTAGGGGTAACCCCGCCGTGGAACGTACCTCCAACTCTGGCATGCAAAGATATTCTGCCaa                                                                         <  1:566964/93‑1 (MQ=255)
                aaaCGCCGATGATGAGCAGTGCCCTTAACAACGTAGTGGTAAACCCGCCGTGGAACGTACCTCCAACTCTGGCATGCAAAGATATTCTGCCaa                                                                         >  2:6265299/1‑93 (MQ=255)
                aaaCGCCGATGATGAGCAGTGCCCTTAACAACGTAGTGGTAAACCCGCCGTGGAACGTACCTCCAACTCTGGCATGCAAAGATATTCTGCCaa                                                                         >  1:7440692/1‑93 (MQ=255)
                  aCGCCGATGATGAGCAGTGCCCTTAACAACGTAGTGGTAAACCCGCCGTGGAACGTACCTCCAACTCTGGCATGCAAAGATATTCTGCCaaaa                                                                       >  1:8777346/1‑93 (MQ=255)
                  aCGCCGATGAGGGAGAGTGCCCTTAACAACGTAGTGGTAAACCCCGCCGGGAACGTACCTCCAACTCTGGCATGCAAAGATATTCTGCCaaaa                                                                       <  1:5755052/93‑1 (MQ=255)
                   cGCCGATGATGAGCAGTGCCCTTAACAACGTAGTGGTAAACCCGCCGTGGAACGTACCTCCAACTCTGGCATGCAAAGATATTCTGCCAAAAg                                                                      <  2:7952243/93‑1 (MQ=255)
                    gCCGATGATGAGCAGTGCCCTTAACAACGTAGTGGTAACCCCGCCGTGGAACGTACCTCCAACTCTGGCATGCAAAGATATTCTGCCAAAAGt                                                                     <  1:8212275/93‑1 (MQ=255)
                    gCCGATGATGAGCAGTGCCCTTAACAACGTAGTGGTAAACCCGCCGTGGAACGTACCTCCAACTCTGGCATGCAAAGATATTCTGCCAAAAGt                                                                     >  1:5202091/1‑93 (MQ=255)
                    gCCGATAATGGGCAGTGCCCTAAACAACGAAGTGGTAACCCCGCCGTGGAACGTACCTCCAACTCTGGCATGCAAAGATATTCTGCCAAAAGt                                                                     <  2:7008686/93‑1 (MQ=255)
                      cGATGATGAGCAGTGCCCTTAACAACGTAGTGGTAAACCCGCCGTGGAACGTACCTCCAACTCTGGCATGCAAAGATATTCTGCCAAAAGTgc                                                                   >  1:2250415/1‑93 (MQ=255)
                        atgatgAGCAGTGCCCTTAACAACGTAGTGGTAACCCCGCCGTGGAACGTACCTCCAACTCTGGCATGCAAAGATATTCTGCCAAAAGTgcgc                                                                 <  2:260704/93‑1 (MQ=255)
                        atgatgAGCAGTGCCCTTAACAACGTAGTGGTAAACCCGCCGTGGAACGTACCTCCAACTCTGGCATGCAAAGATATTCTGCCAAAAGTgcgc                                                                 <  2:4965911/93‑1 (MQ=255)
                          gatgaGCAGTGCCCTTAACAACGTAGTGGTAAACCCGCCGTGGAACGTACCTCCAACTCTGGCATGCAAAGATATTCTGCCAAAAGTGCGCaa                                                               >  2:2729172/1‑93 (MQ=255)
                           atgaGCAGTGCCCTTAACAACGTAGTGGTAACCCCGCCGTGGAACGTACCTCCAACTCTGGCATGCAAAGATATTCTGCCAAAAGTGCGCAAc                                                              <  1:5885964/93‑1 (MQ=255)
                             gaGCAGTGCCCTTAACAACGTAGTGGTAAACCCGCCGTGGAACGTACCTCCAACTCTGGCATGCAAAGATATTCTGCCAAAAGTGCGCAACGa                                                            >  1:537266/1‑93 (MQ=255)
                             gaGCAGTGCCCTTAACAACGTAGTGGTAAACCCGCCGTGGAACGTACCTCCAACTCTGGCATGCAAAGATATTCTGCCAAAAGTGCGCAACGa                                                            >  1:1859746/1‑93 (MQ=255)
                             gaGCAGTCCCCTTAACAACGTAGTGGTAAACCCGCCGTGGAACGTACCTCCAACTCTGGCATGCAAAGATATTCTGCCAAAAGTGCGCAACGa                                                            <  2:4015705/93‑1 (MQ=255)
                              atcggTGCCCTTAACAACGTAGTGGAAAGCCCGCCGTGGAACGTACCTCCAACTCTGGCATGCAAAGATATTCTGCCAAAAGTGCGCAACGAt                                                           <  1:7595343/89‑1 (MQ=255)
                              aGCAGTGCCCTTAACAACGTAGTGGTAACCCCGCCGTGGAACGTACCTCCAACTCTGGCATGCAAAGATATTCTGCCAAAAGTGCGCAACGAt                                                           <  1:6844837/93‑1 (MQ=255)
                               gCAGTGCCCTTAACAACGTAGTGGTAACCCCGCCGTGGAACGTACCTCCAACTCTGGCATGCAAAGATATTCTGCCAAAAGTGCGCAACGATc                                                          <  2:1909708/93‑1 (MQ=255)
                               gCAGTGCCCTTAACAACGTAGTGGTAAACCCGCCGTGGAACGTACCTCCAACTCTGGCATGCAAAGATATTCTGCCAAAAGTGCGCAACGATc                                                          <  2:7036442/93‑1 (MQ=255)
                                cggTGCCCTTAACATCTTGGTGGAAACCCCGCCGTGGAACGTACCTCCAACCCTGGCATGCAAAGATATTCTGCCAAAAGTGCGCAACGATcc                                                         <  2:4779264/91‑1 (MQ=255)
                                 aGTGCCCTTAACAACGTAGTGGTAAACCCGCCGTGGAACGTACCTCCAACTCTGGCATGCAAAGATATTCTGCCAAAAGTGCGCAACGATACg                                                        <  1:8365424/93‑1 (MQ=255)
                                  gTGCCCTTAACAACGTAGTGGTAAACCCGCCGTGGAACGTACCTCCAACTCTGGCATGCAAAGATATTCTGCCAAAAGTGCGCAACGATCCgg                                                       <  1:4346193/93‑1 (MQ=255)
                                  gTGCCCTTAACAACGTAGTGGTAAACCCGCCGTGGAACGTACCTCCAACTCTGGCATGCAAAGATATTCTGCCAAAAGTGCGCAACGATCCgg                                                       >  1:1680689/1‑93 (MQ=255)
                                   tGCCCTTAACAACGTAGTGGTAAACCCGCCGTGGAACGTACCTCCAACTCTGGCATGCAAAGATATTCTGCCAAAAGTGCGCAACGATCCggg                                                      <  1:4470952/93‑1 (MQ=255)
                                    ccccTTAACAACGTAGTGGTAACCCCGCCGTGGAACGTACCTCCAACTCTGGCATGCAAAGATATTCTGCCAAAAGTGCGCAACGATCCGGGa                                                     <  2:2143006/92‑1 (MQ=255)
                                     cccTTAACAACGTAGTGGTAAACCCGCCGTGGAACGTACCTCCAACTCTGGCATGCAAAGATATTCTGCCAAAAGTGCGCAACGATCCGGGat                                                    >  2:2306753/1‑93 (MQ=255)
                                     cccTTAACAACGTAGTGGTAAACCCGCCGTGGAACGTACCTCCAACTCTGGCATGCAAAGATATTCTGCCAAAAGTGCGCAACGATCCGGGat                                                    <  2:1320430/93‑1 (MQ=255)
                                       cTTAACAACGTAGTGGTAACCCCGCCGTGGAACGTACCTCCAACTCTGGCATGCAAAGATATTCTGCCAAAAGTGCGCAACGATCCGGGatat                                                  <  2:4318734/93‑1 (MQ=255)
                                        ttAACAACGTAGTGGTAAACCCGCCGTGGAACGTACCTCCAACTCTGGCATGCAAAGATATTCTGCCAAAAGTGCGCAACGATCCGGGATAtc                                                 >  1:2367097/1‑93 (MQ=255)
                                        ttAACAACGTAGTGGTAAACCCGCCGTGGAACGTACCTCCAACTCTGGCATGCAAAGATATTCTGCCAAAAGTGCGCAACGATCCGGGATAtc                                                 <  2:7358100/93‑1 (MQ=255)
                                         tAACAACGTAGTGGTAAACCCGCCGTGGAACGTACCTCCAACTCTGGCATGCAAAGATATTCTGCCAAAAGTGCGCAACGATCCGGGATAtct                                                <  1:3848832/93‑1 (MQ=255)
                                           acaccGTAGTGGTAACCCCGCCGTGGAACGTACCTCCAACTCTGGCATGCAAAGATATTCTGCCAAAAGTGCGCAACGATCCGGGATATCTCg                                              <  2:5070571/93‑1 (MQ=255)
                                           acaacGTAGTGGTAAACCCGCCGTGGAACGTACCTCCAACTCTGGCATGCAAAGATATTCTGCCAAAAGTGCGCAACGATCCGGGATATCTCg                                              <  1:714984/93‑1 (MQ=255)
                                             agcGTAGTGGTAACCCCGCCGTGGAACGTACCTCCAACTCTGGCATGCAAAGATATTCTGCCAAAAGTGCGCAACGATCCGGGATATCTCGaa                                            <  2:6843800/91‑1 (MQ=255)
                                             aacTTAGTGGTAACCCCGCCGTGGAACGTACCTCCAACTCTGGCATGCAAAGATATTCTGCCAAAAGTGCGCAACGATCCGGGATATCTCGaa                                            <  2:412698/93‑1 (MQ=255)
                                              acGTAGTGGTAAACCCGCCGTGGAACGTACCTCCAACTCTGGCATGCAAAGATATTCTGCCAAAAGTGCGCAACGATCCGGGATATCTCGaaa                                           >  1:2580781/1‑93 (MQ=255)
                                                gTAGTGGTAAACCCGCCGTGGAACGTACCTCCAACTCTGGCATGCAAAGATATTCTGCCAAAAGTGCGCAACGATCCGGGATATCTCGAAAGc                                         >  2:1280742/1‑93 (MQ=255)
                                                gTAGTGGTAAACCCGCCGTGGAACGTACCTCCAACTCTGGCATGCAAAGATATTCTGCCAAAAGTGCGCAACGATCCGGGATATCTCGAAAGc                                         <  1:3234335/93‑1 (MQ=255)
                                                gTAGTGGAAAACCCGCCGTGGAACGTACCTCCAACTCTGGCATGCAAAGATATTCTGCCAAAAGTGCGCAACGATCCGGGATATCTCGAAAGc                                         <  1:7551149/93‑1 (MQ=255)
                                                 tAGTGGTAAACCCGCCGTGGAACGTACCTCCAACTCTGGCATGCAAAGATATTCTGCCAAAAGTGCGCAACGATCCGGGATATCTCGAAAGcc                                        >  1:2112854/1‑93 (MQ=255)
                                                 tAGTGGTAAACCCGCCGTGGAACGTACCTCCAACTCTGGCATGCAAAGATATTCTGCCAAAAGTGCGCAACGATCCGGGATATCTCGAAAGcc                                        <  2:8709635/93‑1 (MQ=255)
                                                  aGTGGTAAACCCGCCGTGGAACGTGCCTCCAACTCTGGCATGCAAAGATATTCTGCCAAAAGTGCGCAACGATCCGGGATATCTCGAAAGCCa                                       <  1:7780253/93‑1 (MQ=255)
                                                  aGTGGTAAACCCGCCGTGGAACGTACCTGCAACTCTGGCATGCAAAGATATTCTGCCAAAAGGGCGCAACGATCCGGGATATCTCGAAAGCCa                                       >  1:858627/1‑93 (MQ=255)
                                                  aGTGGTAAACCCGCCGTGGAACGTACCTCCAACTCTGGCATGCAAAGATATTCTGCCAAAAGTGCGCAACGATCCGGGATATCTCGAAAGcaa                                       >  1:2853413/1‑91 (MQ=255)
                                                   gTGGTAAACCCGCCGTGGAACGTACCTCCAACTCTGGCATGCAAAGATATTCTGCCAAAAGTGCGCAACGATCCGGGATATCTCGAAAGCCAt                                      >  1:7883100/1‑93 (MQ=255)
                                                     ggTAAACCCGCCGTGGAACGTACCTCCAACTCTGGCATGCAAAGATATTCTGCCAAAAGTGCGCAACGATCCGGGATATCTCGAAAGCCATgg                                    >  2:7574996/1‑93 (MQ=255)
                                                     ggTAAACCCGCCGTGGAACGTACCTCCAACTCTGGCATGCAAAGATATTCTGCCAAAAGTGCGCAACGATCCGGGATATCTCGAAAGCCATgg                                    <  2:3637455/93‑1 (MQ=255)
                                                      gTAAACCCGCCGTGGAACGTACCTCCAACTCTGGCATGCAAAGATATTCTGCCAAAAGTGCGCAACGATCCGGGATATCTCGAAAGCCATGGc                                   >  1:1071227/1‑93 (MQ=255)
                                                      gTAAACCCGCCGTGGAACGTACCTCCAACTCTGGCATGCAAAGATATTCTGCCAAAAGTGCGCAACGATCCGGGATATCTCGAAAGCCATGGc                                   <  2:7102692/93‑1 (MQ=255)
                                                       gaaaCCCGCCGTGGAACGTACCTCCAACTCTGGCATGCAAAGATATTCTGCCAAAAGTGCGCAACGATCCGGGATATCTCGAAAGCCATGGCt                                  <  1:1822716/92‑1 (MQ=255)
                                                         aaCCCGCCGTGGAACGTACCTCCAACTCTGGCATGCAAAGATATTCTGCCAAAAGTGCGCAACGATCCGGGATATCTCGAAAGCCATGGCtat                                >  1:5092763/1‑93 (MQ=255)
                                                          aCCCGCCGTGGAACGTACCTCCAACTCTGGCATGCAAAGATATTCTGCCAAAAGTGCGCAACGATCCGGGATATCTCGAAAGCCATGGCtata                               <  1:6159783/93‑1 (MQ=255)
                                                               gccgTGGAGCGTACCTCCAACTCTGGCATGCAAAGATATTCTGCCAAAAGTGCGCAACGATCCGGGATATATCGAAAGCCATGGCTATACGGt                           <  2:8323803/93‑1 (MQ=255)
                                                               gccgTGGAACGTACCTCCAACTCTGGCATGCAAAGATATTCTGCCAAAAGTGCGCAACGATCCGGGATATCTCGAAAGCCATGGCTATACGGt                           >  2:6167981/1‑93 (MQ=255)
                                                                ccgTGGAACGTACCTCCAACTCTGGCATGCAAAGATATTCTGCCAAAAGTGCGCAACGATCCGGGATATCTCGAAAGCCATGGCTATACGGTg                          >  2:4721605/1‑93 (MQ=255)
                                                                 cgTGGAACGTACCTCCAACTCTGGCATGCAAAGATATTCTTCCAAAAGTGCGCAACGATCCGGGATATCTCGAAAGCCATGGCTATACGGTGa                         <  1:7371958/93‑1 (MQ=255)
                                                                 cgTGGAACGTACCTCCAACTCTGGCATGCAAAGATATTCTGCCAAAAGTGCGCAACGATCCGGGATATCTCGAAAGCCATGGCTATACGGTGa                         <  1:978622/93‑1 (MQ=255)
                                                                 cgTGAAACGTACCTCCAACTCTGGCATGCAAAGATATTCTGCCAAAAGTGCGCAACGATCCGGGATATCTCGAAAGCCATGGCTATACGGTGa                         <  1:119610/93‑1 (MQ=255)
                                                                  gTGGAACGTACCTCCAACTCTGGCATGCAAAGATATTCTGCCAAAAGTGCGCAACGATCCGGGATATCTCGAAAGCCATGGCTATACGGTGAt                        >  1:472543/1‑93 (MQ=255)
                                                                  gTGGAACGTACCTCCAACTCTGGCATGCAAAGATATTCTGCCAAAAGTGCGCAACGATCCGGGATATCTCGAAAGCCATGGCTATACGGTGAt                        >  2:3853380/1‑93 (MQ=255)
                                                                    ggAACGTACCTCCAACTCTGGCATGCAAAGATATTCTGCCAAAAGTGCGCAACGATCCGGGATATCTCGAAAGCCATGGCTATACGGTGATgc                      <  2:5420300/93‑1 (MQ=255)
                                                                    ggAACGTACCTCCAACTCTGGCATGCAAAGATATTCTGCCAAAAGTGCGCAACGATCCGGGATATCTCGAAAGCCATGGCTATACGGTGATgc                      <  2:3760926/93‑1 (MQ=255)
                                                                     gAACGTACCTCCAACTCTGGCATGCAAAGATATTCTGCCAAAAGTGCGCAACGATCCGGGATATCTCGAAAGCCATGGCTATACGGTGATgcg                     >  2:4820368/1‑93 (MQ=255)
                                                                     gAACGTACCTCCAACTCTGGCATGCAAAGATATTCTGCCAAAAGTGCGCAACGATCCGGGATATCTCGAAAGCCATGGCTATACGGTGATgcg                     >  2:1240756/1‑93 (MQ=255)
                                                                      aaCGTACCTCCAACTCTGGCATGCAAAGATATTCTGCCAAAAGTGCGCAACGATCCGGGATATCTCGAAAGCCATGGCTATACGGTGATgcgc                    >  2:240928/1‑93 (MQ=255)
                                                                       aCGTACCTCCAACTCTGGCATGCAAAGATATTCTGCCAAAAGTGCGCAACGATCCGGGATATCTCGAAAGCCATGGCTATACGGTGATGcgcg                   >  2:2900724/1‑93 (MQ=255)
                                                                       aCGTACCTCCAACTCTGGCATGCAAAGATATTCTGCCAAAAGTGCGCAACGATCCGGGATATCTCGAAAGCCATGGCTATACGGTGATGcgcg                   >  2:1127566/1‑93 (MQ=255)
                                                                        cGTACCTCCAACTCTGGCATGCAAAGATATTCTGCCAAAAGTGCGCAACGATCCGGGATATCTCGAAAGCCATGGCTATACGGTGATGCGCgg                  >  2:8498959/1‑93 (MQ=255)
                                                                         gTACCTCCAACTCTGGCATGCAAAGATATTCTGCCAAAAGTGCGCAACGATCCGGGATATCTCGAAAGCCATGGCTATACGGTGATGCGCGGc                 <  2:3063501/93‑1 (MQ=255)
                                                                          tACCTCCAACTCTGGCATGCAAAGATATTCTGCCAAAAGTGCGCAACGATCCGGGATATCTCGAAAGCCATGGCTATACGGTGATGCGCGGCt                <  2:5833895/93‑1 (MQ=255)
                                                                          tACCTCCAACTCTGGCATGCAAAGATATTCTGCCAAAAGTGCGCAACGATCCGGGATATCTCGAAAGCCATGGCTATACGGTGATGCGCGGCt                >  1:4058093/1‑93 (MQ=255)
                                                                          tACCTCCAACTCTGGCATGCAAAGATATTCTGCCAAAAGTGCGCAACGATCCGGGATATCTCGAAAGCCATGGCTATACGGTGATGCGCGGCt                >  2:3124782/1‑93 (MQ=255)
                                                                          tACCTCCAACTCTGGCATGCAAAGATATTCTGCCAAAAGTGCGCAACGATCCGGGATATCTCGAAAGCCATGGCTATACGGTGATGCGCGGCt                >  2:6996336/1‑93 (MQ=255)
                                                                          tACCTCCAACTCTGGCATGCAAAGATATTCTGCCAAAAGTGCGCAACGATCCGGGATATCTCGAAAGCCATGGCTATACGGTGATGCGCGGCt                >  1:3059836/1‑93 (MQ=255)
                                                                                 aaCTCTGGCATGCAAAGATATTCTGCCAAAAGTGCGCAACGATCCGGGATATCTCGAAAGCCATGGCTATACGGTGATGCGCGGCTGGAAcag         <  2:5299148/93‑1 (MQ=255)
                                                                                  aCTCTGGCATGCAAAGATATTCCGCCAAAAGTGCGCAACGATCCGGGATATCTCGAAAGCCATGGCTATACGGTGATGCGCGGCTGGAAcagc        >  2:7784294/1‑93 (MQ=255)
                                                                                   ctctGGCATGCAAAGATATTCTGCCAAAAGTGCGCAACGATCCGGGATATCTCGAAAGCCATGGCTATACGGTGATGCGCGGCTGGAAcagca       <  1:6321701/93‑1 (MQ=255)
                                                                                   ctctGGCATGCAAAGATATTCTGCCAAAAGTGCGCAACGATCCGGGATATCTCGAAAGCCATGGCTATACGGTGATGCGCGGCTGGAAcagca       <  2:4515610/93‑1 (MQ=255)
                                                                                    tctGGCATGCAAAGATATTCTGCCAAAAGTGCGCAACGATCCGGGATATCTCGAAAGCCATGGCTATACGGTGATGCGCGGCTGGAAcagcag      >  2:5583480/1‑93 (MQ=255)
                                                                                     ctGGCATGCAAAGATATTCTGCCAAAAGTGCGCAACGATCCGGGATATCTCGAAAGCCATGGCTATACGGTGATGCGCGGCTGGAACAGCaga     <  1:1468822/93‑1 (MQ=255)
                                                                                     caggCATGCAAAGATATTCTGCCAAAAGTCCGCATCGATCCGGGATATCACGAAAGCCATGGCTATACTGGGATGCGCGGCTGGAACAGCaga     <  2:7686420/91‑1 (MQ=255)
                                                                                      tGGCATGCAAAGATATTCTGCCAAAAGTGCGCAACGATCCGGGATATCTCGAAAGCCATGGCTATACGGTGATGCGCGGCTGGAACAGCagag    >  1:65646/1‑93 (MQ=255)
                                                                                       ggCATGCAAAGATATTCTGCCAAAAGTGCGCAACGATCCGGGATATCTCGAAAGCCATGGCTATACGGTGATGCGCGGCTGGAACAGCAgaga   <  2:334319/93‑1 (MQ=255)
                                                                                           |                                                                                         
GGTCGCCCCGATCGCAAAACGCCGATGATGAGCAGTGCCCTTAACAACGTAGTGGTAAACCCGCCGTGGAACGTACCTCCAACTCTGGCACGCAAAGATATTCTGCCAAAAGTGCGCAACGATCCGGGATATCTCGAAAGCCATGGCTATACGGTGATGCGCGGCTGGAACAGCAGAGAA  >  NC_000913/981476‑981655

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 18 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: