Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,038,298 A→G intergenic (+82/+188) rrfA → / ← mobB 5S ribosomal RNA/molybdopterin‑guanine dinucleotide biosynthesis protein B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,038,2980AG82.7% 193.7 / 62.7 139intergenic (+82/+188)rrfA/mobB5S ribosomal RNA/molybdopterin‑guanine dinucleotide biosynthesis protein B
Reads supporting (aligned to +/- strand):  ref base A (13/11);  new base G (66/49);  total (79/60)
Fisher's exact test for biased strand distribution p-value = 8.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.17e-02

Only 100 of 142 total aligned reads displayed.

GCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAG  >  NC_000913/4038208‑4038333
                                                                                          |                                   
gCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGGCa                                   >  1:2285175/1‑93 (MQ=21)
  cAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGGCaaa                                 >  1:3573282/1‑93 (MQ=18)
  cAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGGCaaa                                 >  1:5929766/1‑93 (MQ=21)
  cAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACaaa                                 <  2:3628866/93‑1 (MQ=33)
    ggCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGGCACATc                               >  2:1810433/1‑93 (MQ=255)
    ggCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGGCAAATc                               <  1:958420/93‑1 (MQ=18)
    ggCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGGCAAATc                               >  1:608001/1‑93 (MQ=21)
    ggCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATc                               <  2:1496723/93‑1 (MQ=33)
    ggCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATc                               >  1:3659495/1‑93 (MQ=34)
     gCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATcc                              <  1:2586176/93‑1 (MQ=33)
      cATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGGCAAATccg                             >  1:266513/1‑93 (MQ=21)
      cATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGGCAAATccg                             <  2:4840854/93‑1 (MQ=18)
      cATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGGCAAATccg                             >  1:2809639/1‑93 (MQ=21)
      cATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGGCAAATccg                             <  2:5731111/93‑1 (MQ=18)
       aTCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGGCAAATccgc                            >  1:4154984/1‑93 (MQ=21)
       aTCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGGCAAATccgc                            >  2:232421/1‑93 (MQ=21)
       aTCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGGCAAATccgc                            <  2:845379/93‑1 (MQ=18)
       aTCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTGCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGGCAAATccgc                            <  2:3912491/93‑1 (MQ=16)
       aTCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTATCGTTTTAGCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGGCAAATccgc                            <  1:3103118/93‑1 (MQ=11)
        tCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGGCAAATccgcc                           <  2:4084059/93‑1 (MQ=18)
        tCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGGCAAATccgcc                           >  1:3711611/1‑93 (MQ=21)
        acAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATccgcc                           <  1:514143/92‑1 (MQ=18)
         cAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGGCAAATccgccg                          <  1:4312733/93‑1 (MQ=18)
         cAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGGCAAATccgccg                          <  1:4748825/93‑1 (MQ=18)
         cAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGGCAAATccgccg                          >  2:869719/1‑93 (MQ=18)
         cAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGGCAAATccgccg                          <  2:1041525/93‑1 (MQ=18)
         cAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGGCAAATccgccg                          <  1:1927498/93‑1 (MQ=18)
         cAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGGCAAATccgccg                          <  2:6109413/93‑1 (MQ=18)
          aaaTAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGGCAAATCCGCCgg                         <  2:2236324/93‑1 (MQ=18)
          aaaTAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGGCAAATCCGCCgg                         <  2:3022434/93‑1 (MQ=18)
          aaaTAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGGCAAATCCGCCgg                         <  2:4072229/93‑1 (MQ=18)
          aaaTAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGGCAAATCCGCCgg                         >  1:6541938/1‑93 (MQ=21)
          aaaTAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTAAGTAGGGCAAATCCGCCgg                         <  1:2379358/93‑1 (MQ=16)
           aaTAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGGCAAATCCGCCggg                        <  1:1094424/93‑1 (MQ=18)
            aTAAAACGAAAGGCTCAGTCGGAATACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGGCAAATCCGCCGGGa                       <  1:2785496/93‑1 (MQ=16)
            aTAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGGCAAATCCGCCGGGa                       <  1:5373198/93‑1 (MQ=18)
             tAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGGCAAATCCGCCGGGAg                      <  1:2450783/93‑1 (MQ=18)
             tAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGGCAAATCCGCCGGGAg                      >  2:6066094/1‑93 (MQ=21)
              aaaaCGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGGCAAATCCGCCGGGAGc                     <  1:4508470/93‑1 (MQ=18)
               aaaCGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGGCAAATCCGCCGGGAGCg                    >  2:1754558/1‑93 (MQ=18)
               aaaCGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGGCAAATCCGCCGGGAGCg                    <  1:6636878/93‑1 (MQ=18)
                aaCGAAAGGCTCAGTCGGAAGACTGGGCGTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGGCAAATCCGCCGGGAGCgg                   <  2:3284515/93‑1 (MQ=255)
                aaCGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGGCAAATCCGCCGGGAGCgg                   >  2:5380831/1‑93 (MQ=21)
                aaCGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCgg                   >  2:6468441/1‑93 (MQ=34)
                 aCGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGGCAAATCCGCCGGGAGCGGa                  <  1:2710979/93‑1 (MQ=18)
                  cGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGGCAAATCCGCCGGGAGCGGAt                 >  1:1117461/1‑93 (MQ=21)
                   gAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGGCAAATCCGCCGGGAGCGGAtt                <  2:2764613/93‑1 (MQ=18)
                   gAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGGCAAATCCGCCGGGAGCGGAtt                <  2:1197190/93‑1 (MQ=18)
                    aaaGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGGCAAATCCGCCGGGAGCGGAttt               <  1:3357949/93‑1 (MQ=18)
                    aaaGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGGCAAATCCGCCGGGAGCGGAttt               >  1:4140309/1‑93 (MQ=21)
                    aaaGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGGCAAATCCGCCGGGAGCGGAttt               >  1:1043693/1‑93 (MQ=21)
                     aaGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGGCAAATCCGCCGGGAGCGGATTTg              >  1:4363874/1‑93 (MQ=18)
                     aaGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGGCAAATCCGCCGGGAGCGGATTTg              >  2:2318904/1‑93 (MQ=21)
                     aaGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGGCAAATCCGCCGGGAGCGGATTTg              >  1:3470354/1‑93 (MQ=21)
                     aaGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTg              >  1:6625952/1‑93 (MQ=34)
                      aGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGGCAAATCCGCCGGGAGCGGATTTGa             >  1:5529512/1‑93 (MQ=21)
                      aGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGGCAAATCCGCCGGGAGCGGATTTGa             >  1:3894326/1‑93 (MQ=21)
                      aGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGa             >  1:4793191/1‑93 (MQ=34)
                       ggCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGGCAAATCCGCCGGGAGCGGATTTGaa            >  1:5019372/1‑93 (MQ=21)
                       ggCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGGCAAATCCGCCGGGAGCGGATTTGaa            <  1:6235426/93‑1 (MQ=18)
                       ggCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGaa            >  1:783228/1‑93 (MQ=34)
                        gCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGGCAAATCCGCCGGGAGCGGATTTGAAc           <  2:4327701/93‑1 (MQ=18)
                        gCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGGCAAATCCGCCGGGAGCGGATTTGAAc           >  2:3851526/1‑93 (MQ=18)
                        gCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAAc           >  1:3145098/1‑93 (MQ=34)
                         cTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGGCAAATCCGCCGGGAGCGGATTTGAACg          <  2:4376676/93‑1 (MQ=18)
                         cTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGGCAAATCCGCCGGGAGCGGATTTGAACg          >  2:1906076/1‑93 (MQ=18)
                         cTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACg          >  1:793542/1‑93 (MQ=34)
                          tCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGGCAAATCCGCCGGGAGCGGATTTGAACGt         >  1:1499769/1‑93 (MQ=21)
                          tCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGGCAAATCCGCCGGGAGCGGATTTGAACGt         >  2:5038762/1‑93 (MQ=21)
                          tCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGGCAAATCCGCCGGGAGCGGATTTGAACGt         <  2:3097682/93‑1 (MQ=18)
                          tCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGGCAAATCCGCCGGGAGCGGATTTGAACGt         >  1:5039715/1‑93 (MQ=18)
                          tCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGGCAAATCCGCCGGGAGCGGATTTGAACGt         >  2:1604143/1‑93 (MQ=21)
                          tCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGt         <  2:6098222/93‑1 (MQ=33)
                          tCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGt         <  2:5910275/93‑1 (MQ=33)
                          tCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTAGGTGAACGCTCTCCTGAGTAGGGCAAATCCGCCGGGAGCGGATTTGAACGt         >  1:6389866/1‑93 (MQ=16)
                           cAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGGCAAATCCGCCGGGAGCGGATTTGAACGtt        >  1:6219389/1‑93 (MQ=18)
                           cAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGGCAAATCCGCCGGGAGCGGATTTGAACGtt        >  1:1294590/1‑93 (MQ=21)
                           cAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATGCGCCGGGAGAGGATTTGAACGtt        >  2:2846391/1‑93 (MQ=16)
                            aGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGGCAAATCCGCCGGGAGCGGATTTGAACGTTg       <  1:3110802/93‑1 (MQ=18)
                             gTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGGCAAATCCGCCGGGAGCGGATTTGAACGTTGc      <  1:5495370/93‑1 (MQ=18)
                             gTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGGCAAATCCGCCGGGAGCGGATTTGAACGTTGc      >  1:4684888/1‑93 (MQ=18)
                             gTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGGCAAATCCGCCGGGAGCGGATTTGAACGTTGc      >  2:2243350/1‑93 (MQ=18)
                             gTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGGCAAATCCGCCGGGAGCGGATTTGAACGTTGc      <  1:647057/93‑1 (MQ=18)
                             gTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGGCAAATCCGCCGGGAGCGGATTTGAACGTTGc      >  2:5304363/1‑93 (MQ=18)
                             gTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGc      <  2:3385653/93‑1 (MQ=33)
                              tCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCg     <  2:1453518/93‑1 (MQ=33)
                               cGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGGCAAATCCGCCGGGAGCGGATTTGAACGTTGCGa    >  1:4074382/1‑93 (MQ=18)
                                ggAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGGCAAATCCGCCGGGAGCGGATTTGAACGTTGCGaa   <  2:64501/93‑1 (MQ=17)
                                ggAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGGCAAATCCGCCGGGAGCGGATTTGAACATTGCGaa   >  2:2309065/1‑93 (MQ=12)
                                 gAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGGCAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAg  >  1:5107734/1‑93 (MQ=14)
                                 gAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGGCAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAg  >  2:577678/1‑93 (MQ=14)
                                 gAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGGCAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAg  >  1:1339317/1‑93 (MQ=14)
                                 gAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGGCAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAg  >  2:5689043/1‑93 (MQ=14)
                                 gAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGGCAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAg  <  1:885359/93‑1 (MQ=14)
                                 gAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGGCAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAg  >  2:1308722/1‑93 (MQ=14)
                                 gAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGGCAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAg  >  2:4270509/1‑93 (MQ=14)
                                 gAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGGCAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAg  >  2:1971462/1‑93 (MQ=14)
                                 gAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGGCAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAg  >  2:1987874/1‑93 (MQ=14)
                                 gAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGGCAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAg  >  1:3272171/1‑93 (MQ=14)
                                 gAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAg  >  2:5370749/1‑93 (MQ=31)
                                                                                          |                                   
GCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAG  >  NC_000913/4038208‑4038333

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