Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,107,129 A→G E239G (GAG→GGG)  sbp → periplasmic sulfate‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,107,1290AG99.6% 815.4 / NA 258E239G (GAG→GGG) sbpperiplasmic sulfate‑binding protein
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (143/114);  total (144/114)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

Only 100 of 266 total aligned reads displayed.

TTGCCTGGGAAAACGAAGCTCTGCTGGCAGCGAATGAACTGGGGAAAGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGAGCCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTACTAAAGAGGTGGCGGAAGCCTACCTGAAATATCTCTACTCGCCA  >  NC_000913/4107039‑4107217
                                                                                          |                                                                                        
ttGCCTGGGAAAACGAAGCTCTGCTGGCAGCGAATGAACTGGGGAAAGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGGGc                                                                                        >  2:4587671/1‑93 (MQ=255)
   ccTGGGAAAACGAAGCTCTGCTGGCAGCGAATGAACTGGGGAAAGATAAATTCGAAATCGTCACGCCGAGCGAGTCTATCCTCGCAGGGCCaa                                                                                     >  2:3131364/1‑93 (MQ=255)
    cTGGGAAAACGAAGCTCTGCTGGCAGCGAATGAACTGGGGAAAGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGGGCCAAc                                                                                    <  1:2564657/93‑1 (MQ=255)
    cTGGGAAAACGAAGCTCTGCGGGCAGCGAATGAACTGGGGAAAGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGGGCCAAc                                                                                    <  1:1390028/93‑1 (MQ=255)
     tGGGAAAACGAAGCTCTGCTGGCAGCGAATGGACTGGGGAAAGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGGGCCAAcc                                                                                   >  2:155902/1‑93 (MQ=255)
     tGGGAAAACGAAGCTCTGCTGGCAGCGAATGAACTGGGGAAAGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGGGCCAAcc                                                                                   <  1:2822128/93‑1 (MQ=255)
     tGGGAAAACGAAGCTCTGCTGGCAGCGAATGAACTGGGGAAAGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGGGCCAAcc                                                                                   >  2:3794884/1‑93 (MQ=255)
      gggAAAACGAAGCTCTGCTGGCAGCGAATGAACTGGGGAAAGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGGGCCAACCg                                                                                  >  1:2912446/1‑93 (MQ=255)
       ggAAAACGAAGCTCTGCTGGCAGCGAATGAACTGGGGAAAGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGGGCCAACCgt                                                                                 <  1:1183075/93‑1 (MQ=255)
       ggAAAACGAAGCTCTGCTGGCAGCGAATGAACTGGGGAAAGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGGGCCAACCgt                                                                                 >  1:3865216/1‑93 (MQ=255)
         aaaaCGAAGCTCTGCTGGCAGCGAATGAACTGGGGAAAGATAAATTCGAAAGCGTCACGCCGAGTGAGTCTATCCTCGCAGGGCCAACCgtgt                                                                               <  1:1049863/93‑1 (MQ=255)
          aaaCGAAGCTCTGCTGGCAGCGAATGAACTGGGGAAAGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGGGCCAACCGTGTc                                                                              <  1:4223598/93‑1 (MQ=255)
             cGAAGCTCTGCTGGCAGCGAATGAACTGGGGAAAGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGGGCCAACCGTGTCggt                                                                           <  1:1703887/93‑1 (MQ=255)
               aaGCTCTGCTGGCAGCGAATGAACTGGGGAAAGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGGGCCAACCGTGTCggtgg                                                                         <  2:3577662/93‑1 (MQ=255)
                aGCTCTGCTGGCAGCGAATGAACTGGGGAAAGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGGGCCAACCGTGTCggtggt                                                                        <  2:2070676/93‑1 (MQ=255)
                 gCTCTGCTGGCAGCGAATGAACTGGGGAAAGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGGGCCAACCGTGTCGGTGGTc                                                                       <  2:797356/93‑1 (MQ=255)
                 gCTCTGCTGGCAGCGAATGAACTGGGGAAAGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGGGCCAACCGTGTCGGTGGTc                                                                       <  1:339419/93‑1 (MQ=255)
                   tctGCTGGCAGCGAATGAACTGGGGAAAGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGGGCCAACCGTGTCGGTGGTCGa                                                                     >  1:5071110/1‑93 (MQ=255)
                   tctGCTGGCAGCGAATGAACTGGGGAAAGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGGGCCAACCGTGTCGGTGGTCGa                                                                     <  2:5887000/93‑1 (MQ=255)
                    ctgctgGCAGCGAATGAACTGGGGAAAGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGGGCCAACCGTGTCGGTGGTCGAt                                                                    <  1:5933852/93‑1 (MQ=255)
                       ctgGCAGCGAATGAACTGGGGAAAGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGGGCCAACCGTGTCGGTGGTCGATaaa                                                                 >  1:1565662/1‑93 (MQ=255)
                       ctgGCAGCGAATGAACTGGGGAAAGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGGGCCAACCGTGTCGGTGGTCGATaaa                                                                 <  2:6182181/93‑1 (MQ=255)
                       ctgGCAGCGAATGAACTGGGGAAAGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGGGCCAACCGGGTCGGTGGTCGATacg                                                                 >  2:3409165/1‑91 (MQ=255)
                        tgGCAGCGAATGAACTGGGGAAAGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGGGCCAACCGTGTCGGTGGTCGATAAAg                                                                >  2:2905746/1‑93 (MQ=255)
                        tgGCAGCGAATGAACTGGGGAAAGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGGGCCAACCGTGTCGGTGGTCGATAAAg                                                                <  1:5732868/93‑1 (MQ=255)
                         ggCAGCGAATGAACTGGGGAAAGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGGGCCAACCGTGTCGGTGGTCGATAAAGt                                                               <  2:3346413/93‑1 (MQ=255)
                         ggCAGCGAATGAACTGGGGAAAGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGGGCCAACCGTGTCGGTGGTCGATAAAGt                                                               >  1:2092472/1‑93 (MQ=255)
                          gCAGCGAATGAACTGGGGAAAGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGGGCCAACCGTGTCGGTGGTCGATAAAGTg                                                              <  1:568322/93‑1 (MQ=255)
                          gCAGCGAATGAACTGGGGAAAGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGGGCCAACCGTGTCGGTGGTCGATAAAGTg                                                              >  1:5268175/1‑93 (MQ=255)
                           cAGCGAATGAACTGGGGAAAGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGGGCCAACCGTGTCGGTGGTCGATAAAGTgg                                                             <  2:3008360/93‑1 (MQ=255)
                            aGCGAATGAACTGGGGAAAGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGGGCCAACCGTGTCGGTGGTCGATAAAGTGGt                                                            <  1:6034361/93‑1 (MQ=255)
                            aGCGAATGAACTGGGGAAAGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGGGCCAACCGTGTCGGTGGTCGATAAAGTGGt                                                            <  2:4238553/93‑1 (MQ=255)
                             gCGAATGAACTGGGGAAAGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGGGCCAACCGTGTCGGTGGTCGATAAAGTGGTc                                                           >  1:2396927/1‑93 (MQ=255)
                             gCGAATGAACTGGGGAAAGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGGGCCAACCGTGTCGGTGGTCGATAAAGTGGTc                                                           <  1:4016830/93‑1 (MQ=255)
                               gAATGAACTGGGGAAAGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGGGCCAACCGTGTCGGTGGTCGATAAAGTGGTCga                                                         >  2:887209/1‑93 (MQ=255)
                               gAATGAACTGGGGAAAGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGGGCCAACCGTGTCGGTGGTCGATAAAGTGGTCga                                                         >  2:3844884/1‑93 (MQ=255)
                                aaTGAACTGGGGAAAGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGGGCCAACCGTGTCGGTGGTCGATAAAGTGGTCgag                                                        >  1:3532378/1‑93 (MQ=255)
                                aaTGAACTGGGCAAAGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGGGCCAACCGTGTCGGTGGTCGATAAAGTGGTCgag                                                        <  2:1648916/93‑1 (MQ=255)
                                 aTGAACTGGGGAAAGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGGGCCAACCGTGTCGGTGGTCGATAAAGTGGTCgaga                                                       >  1:3113660/1‑93 (MQ=255)
                                 aTGAACTGGGGAAAGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGGGCCAACCGTGTCGGTGGTCGATAAAGTGGTCgaga                                                       >  1:840113/1‑93 (MQ=255)
                                  tGAACTGGGGAAAGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGGGCCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGaa                                                      <  2:4461157/93‑1 (MQ=255)
                                    aaCTGGGGAAAGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGGGCCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGaaaa                                                    >  1:5856348/1‑93 (MQ=255)
                                     aCTGGGGAAAGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGGGCCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGaaaaa                                                   >  2:4608122/1‑93 (MQ=255)
                                      cTGGGGAAAGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGGGCCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGaaaaaa                                                  <  1:502792/93‑1 (MQ=255)
                                       tGGGGAAAGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGGGCCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAg                                                 >  2:4772602/1‑93 (MQ=255)
                                       tGGGGAAAGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGGGCCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAg                                                 >  2:4933955/1‑93 (MQ=255)
                                        ggggAAAGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGGGCCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAgg                                                <  1:4464869/93‑1 (MQ=255)
                                        ggggAAAGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGGGCCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAgg                                                >  1:1611543/1‑93 (MQ=255)
                                         gggAAAGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGGGCCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGt                                               <  2:1563673/93‑1 (MQ=255)
                                          ggAAAGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGGGCCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTa                                              >  1:615164/1‑93 (MQ=255)
                                          ggAAAGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGGGCCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTa                                              >  2:3600493/1‑93 (MQ=255)
                                           gAAAGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGGGCCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTAc                                             >  2:832941/1‑93 (MQ=255)
                                           gAAAGATAAATGCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGGGCCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTAc                                             <  1:887514/93‑1 (MQ=255)
                                               gATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGGGCCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTACTaaa                                         >  2:5771157/1‑93 (MQ=255)
                                                aTAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGGGCCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTACTAAag                                        <  1:1822664/93‑1 (MQ=255)
                                                 tAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGGGCCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTACTAAaga                                       >  2:2856882/1‑93 (MQ=255)
                                                 tAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGGGCCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTACTAAaga                                       <  1:3342179/93‑1 (MQ=255)
                                                 tAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGGGCCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTACTAAaga                                       >  2:387702/1‑93 (MQ=255)
                                                 tAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGGGCCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTACTAAaga                                       >  2:4825428/1‑93 (MQ=255)
                                                  aaaTTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGGGCCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTACTAAagag                                      >  1:2308929/1‑93 (MQ=255)
                                                   aaTTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGGGCCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTACTAAAGAgg                                     <  2:5244847/93‑1 (MQ=255)
                                                   aaTTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGGGCCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTACTAAAGAgg                                     <  1:3831381/93‑1 (MQ=255)
                                                     ttCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGGGCCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTACTAAAGAGGTg                                   <  1:2201260/93‑1 (MQ=255)
                                                     ttCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGGGCCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTACTAAAGAGGTg                                   >  2:177895/1‑93 (MQ=255)
                                                      tCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGGGCCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTACTAAAGAGGTgg                                  <  2:1535133/93‑1 (MQ=255)
                                                       cGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGGGCCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTACTAAAGAGGTGGc                                 <  1:1548036/93‑1 (MQ=255)
                                                        gAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGGGCCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTACTAAAGAGGTGGCg                                <  2:2715104/93‑1 (MQ=255)
                                                        gAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGGGCCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTACTAAAGAGGTGGCg                                <  2:1317474/93‑1 (MQ=255)
                                                        gAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGGGCCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTACTAAAGAGGTGGCg                                >  1:2652816/1‑93 (MQ=255)
                                                           aTCGGCACGCCGAGTGAGTCTATCCTCGCAGGGCCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTACTAAAGAGGTGGCGGaa                             <  1:4949765/93‑1 (MQ=255)
                                                            tCGTCACGCCGAGTGAGTCTATCCTCGCAGGGCCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTACTAAAGAGGTGGCGGAAg                            >  2:1273395/1‑93 (MQ=255)
                                                             cGTCACGCCGAGTGAGTCTATCCTCGCAGGGCCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTACTAAAGAGGTGGCGGAAGc                           >  2:5428536/1‑93 (MQ=255)
                                                              gTCACGCCGAGTGAGTCTATCCTCGCAGGGCCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTACTAAAGAGGTGGCGGAAGcc                          >  1:4909949/1‑93 (MQ=255)
                                                              gTCACGCCGAGTGAGTCTATCCTCGCAGGGCCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTACTAAAGAGGTGGCGGAAGcc                          <  2:5746534/93‑1 (MQ=255)
                                                               tCACGCCGAGTGAGTCTATCCTCGCAGGGCCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTACTAAAGAGGTGGCGGAAGCCt                         >  2:3679823/1‑93 (MQ=255)
                                                               tCACGCCGAGTGAGTCTATCCTCGCAGGGCCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTACTAAAGAGGTGGCGGAAGCCt                         >  1:483603/1‑93 (MQ=255)
                                                               tCACGCCGAGTGAGTCAATCCTCGCAGGGCCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTACTAAAGAGGTGGCGGAAGCCt                         <  2:5025696/93‑1 (MQ=255)
                                                                 aCGCCGAGTGAGTCTATCCTCGCAGGGCCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTACTAAAGAGGTGGCGGAAGCCTAc                       >  2:4511752/1‑93 (MQ=255)
                                                                 aCGCCGAGTGAGTCTATCCTCGCAGGGCCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTACTAAAGAGGTGGCGGAAGCCTAc                       >  2:1967486/1‑93 (MQ=255)
                                                                  cGCCGAGTGAGTCTATCCTCGCAGGGCCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTACTAAAGAGGTGGCGGAAGCCTAcc                      <  2:4102956/93‑1 (MQ=255)
                                                                  cGCCGAGTGAGTCTATCCTCGCAGGGCCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTACTAAAGAGGTGGCGGAAGCCTAcc                      <  1:1260553/93‑1 (MQ=255)
                                                                  cGCCGAGTGAGTCTATCCTCGCAGGGCCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTACTAAAGAGGTGGCGGAAGCCTAcc                      >  2:1213064/1‑93 (MQ=255)
                                                                  cGCCGAGTGAGTCTATCCTCGCAGGGCCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTACAAAAGAGGTGGCGGAAGCCTAcc                      <  1:3997274/93‑1 (MQ=255)
                                                                     cGAGTGAGTCTATCCTCGCAGGGCCAACCGTGTCGGTGGTCGATAAGGTGGTCGAGAAAAAAGGTACTAAAGAGGTGGCGGAAGCCTACCTGa                   <  2:5073582/93‑1 (MQ=255)
                                                                     cGAGTGAGTCTATCCTCGCAGGGCCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTACTAAAGAGGTGGCGGAAGCCTACCTGa                   >  1:4059983/1‑93 (MQ=255)
                                                                      gagtgagtCTATCCTCGCAGGGCCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTACTAAAGAGGTGGCGGAAGCCTACCTGaa                  >  2:5599874/1‑93 (MQ=255)
                                                                      gagtgagtCTATCCTCGCAGGGCCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTACTAAAGAGGTGGCGGAAGCCTACCTGaa                  >  2:1025998/1‑93 (MQ=255)
                                                                      gagtgagtCTATCCTCGCAGGGCCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTACTAAAGAGGTGGCGGAAGCCTACCTGaa                  >  2:321645/1‑93 (MQ=255)
                                                                       agtgagtCTATCCTCGCAGGGCCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTACTAAAGAGGTGGCGGAAGCCTACCTGaaa                 <  1:3279392/93‑1 (MQ=255)
                                                                       agtgagtCTATCCTCGCAGGGCCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTACTAAAGAGGTGGCGGAAGCCTACCTGaaa                 >  2:5050198/1‑93 (MQ=255)
                                                                       agtgagtCTATCCTCGCAGGGCCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTACTAAAGAGGTGGCGGAAGCCTACCTGaaa                 >  1:3674166/1‑93 (MQ=255)
                                                                        gtgagtCTATCCTCGCAGGGCCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTACTAAAGAGGTGGCGGAAGCCTACCTGAAat                >  2:618535/1‑93 (MQ=255)
                                                                         tgagtCTATCCTCGCAGGGCCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTACTAAAGAGGTGGCGGAAGCCTACCTGAAata               <  2:3125199/93‑1 (MQ=255)
                                                                            gtCTATCCTCGCAGGGCCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTACTAAAGAGGTGGCGGAAGCCTACCTGAAATAtct            <  1:436609/93‑1 (MQ=255)
                                                                             tCTATCCTCGCAGGGCCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTACTAAAGAGGTGGCGGAAGCCTACCTGAAATAtctc           <  1:1016445/93‑1 (MQ=255)
                                                                               tATCCTCGCAGGGCCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTACTAAAGAGGTGGCGGAAGCCTACCTGAAATATCTCTa         >  2:220920/1‑93 (MQ=255)
                                                                               tATCCTCGCAGGGCCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTACTAAAGAGGTGGCGGAAGCCTACCTGAAATATCTCTa         >  2:67456/1‑93 (MQ=255)
                                                                                aTCCTCGCAGGGCCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTACTAAAGAGGTGGCGGAAGCCTACCTGAAATATCTCTAc        >  2:4655742/1‑93 (MQ=255)
                                                                                    tCGCAGGGCCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTACTAAAGAGGTGGCGGAAGCCTACCTGAAATATCTCTACTCGc    <  2:1848528/93‑1 (MQ=255)
                                                                                      gCAGGGCAAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTACTAAAGAGGTGGCGGAAGCCTACCTGAAATATCTCTACTCGCCa  <  2:1519079/93‑1 (MQ=255)
                                                                                          |                                                                                        
TTGCCTGGGAAAACGAAGCTCTGCTGGCAGCGAATGAACTGGGGAAAGATAAATTCGAAATCGTCACGCCGAGTGAGTCTATCCTCGCAGAGCCAACCGTGTCGGTGGTCGATAAAGTGGTCGAGAAAAAAGGTACTAAAGAGGTGGCGGAAGCCTACCTGAAATATCTCTACTCGCCA  >  NC_000913/4107039‑4107217

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