Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,525,586 A→G intergenic (‑38/+93) sgcA ← / ← sgcQ putative PTS system enzyme IIA component/putative nucleoside triphosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,525,5860AG99.7% 1021.5 / ‑5.4 323intergenic (‑38/+93)sgcA/sgcQputative PTS system enzyme IIA component/putative nucleoside triphosphatase
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (151/171);  total (152/171)
Fisher's exact test for biased strand distribution p-value = 4.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.70e-01

Only 100 of 335 total aligned reads displayed.

CGCCAGTCGGTGGCTTTACGTTGTGCCTGTACCCATTTAATATCGTTAATCATCCGTCGTTGACTCCATGCCGATTCGGGTTAATCTGGTAGCGATCCCCGTCGATACTTTTGACGAAGGCGGCAGGGATCGCAGAAGGGTTATTGCTCTTTGCCGGGGATAAGCGCTCTGCTTATCCCCA  >  NC_000913/4525496‑4525676
                                                                                          |                                                                                          
cGCCAGTCGGTGGCTTTACGTTGTGCCTGTACCCATTTAATATCGTTAATCATCCGTCGTTGACTCCATGCCGATTCGGGTTAATCTGGTGGc                                                                                          <  2:631227/93‑1 (MQ=255)
 gCCAGTCGGTGGCTTTACGTTGTGCCTGTACCCATTTAATATCGTTAATCATCCGTCGTTGACTCCATGCCGATTCGGGTTAATCTGGTGGCg                                                                                         >  2:2149520/1‑93 (MQ=255)
 gCCAGGCGGTGGCTTTACGTTGTGCCTGTACCCATTTAATATCGTTAATCATCCGTCGTTGACTCCATGCCGATTCGGGTTAATCTGGTGGCg                                                                                         <  1:1020140/93‑1 (MQ=255)
   cAGTCGGTGGCTTTACGTTGTGCCTGTACCCATTTAATATCGTTAATCATCCGTCGTTGACTCCATGCCGATTCGGGTTAATCTGGTGGCGAt                                                                                       >  2:5818286/1‑93 (MQ=255)
   cAGTCGGTGGCTTTACGTTGTGCCTGTACCCATTTAATATCGTTAATCATCCGTCGTTGACTCCATGCCGATTCGGGTTAATCTGGTGGCGAt                                                                                       >  2:2305452/1‑93 (MQ=255)
   cAGTCGGTGGCTTTACGTTGTGCCTGTACCCATTTAATATCGTTAATCATCCGTCGTTGACTCCATGCCGATTCGGGTTAATCTGGTGGCGAt                                                                                       >  2:3262227/1‑93 (MQ=255)
   cAGTCGGTGGCTTTACGTTGTGCCTGTACCCATTTAATATCGTTAATCATCCGTCGTTGACTCCATGCCGATTCGGGTTAATCTGGTGGCGAt                                                                                       >  2:3128756/1‑93 (MQ=255)
       cGGTGGCTTTACGTTGTGCCTGTACCCATTTAATATCGTTAATCATCCGTCGTTGACTCCATGCCGATTCGGGTTAATCTGGTGGCGATcccc                                                                                   <  1:2143811/93‑1 (MQ=255)
         gTGGCTTTACGTTGTGCCTGTACCCATTTAATATCGTTAATCATCCGTCGTTGACTCCATGCCGATTCGGGTTAATCTGGTGGCGATCCCCGt                                                                                 <  2:3188684/93‑1 (MQ=255)
          tGGCTTTACGTTGTGCCTGTACCCATTTAATATCGTTAATCATCCATCGTTGACTCCATGCCGATTCGGGTTAATCTGGTGGCGATCCCCGTc                                                                                >  2:6413287/1‑93 (MQ=255)
           ggCTTTACGTTGTGCCTGTACCCATTTAATATCGTTAATCATCCGTCGTTGACTCCATGCCGATTCGGGTTAATCTGGTGGCGATCCCCGTCg                                                                               <  1:3149329/93‑1 (MQ=255)
            gCTTTACGTTGTGCCTGTACCCATTTAATATCGTTAATCATCCGTCGTTGACTCCATGCCGATTCGGGTTAATCTGGTGGCGATCCCCGTCGa                                                                              <  1:4063858/93‑1 (MQ=255)
             cTTTACGTTGTGCCTGTACCCATTTAATATCGTTAATCATCCGTCGTTGACTCCATGCCGATTCGGGTTAATCTGGTGGCGATCCCCGTCGAt                                                                             <  2:5606329/93‑1 (MQ=255)
             cTTTACGTTGTGCCTGTACCCATTTAATATCGTTAATCATCCGTCGTTGACTCCATGCCGATTCGGGTTAATCTGGTGGCGATCCCCGTCGAt                                                                             >  2:1172101/1‑93 (MQ=255)
              tttACGTTGTGCCTGTACCCATTTAATATCGTTAATCATCCGTCGTTGACTCCATGCCGATTCGGGTTAATCTGGTGGCGATCCCCGTCGATa                                                                            >  1:6812558/1‑93 (MQ=255)
               ttACGTTGTGCCTGTACCCATTTAATATCGTTAATCATCCGTCGTTGACTCCATGCCGATTCGGGTTAATCTGGTGGCGATCCCCGTCGATAc                                                                           <  2:4751284/93‑1 (MQ=255)
               ttACGTTGTGCCTGTACCCATTTAATATCGTTAATCATCCGTCGTTGACTCCATGCCGATTCGGGTTAATCTGGTGGAGATCCCCGTCGATAc                                                                           >  1:5791726/1‑93 (MQ=255)
                 aCGTTGTGCCTGTACCCATTTAATATCGTTAATCATCCGTCGTTGACTCCATGCCGATTCGGGTTAATCTGGTGGCGATCCCCGTCGATACtt                                                                         >  1:4361148/1‑93 (MQ=255)
                   gTTGCGCCTGTACCCATTTAATATCGTTAATCATCCGTCGTTGACTCCATGCCGATTCGGGTTAATCTGGTGGCGATCCCCGTCGATACtttt                                                                       <  1:6398241/93‑1 (MQ=255)
                     tgtgCCTGTACCCATTTAATATCGTTAATCATCCGTCGTTGACTCCATGCCGATTCGGGTTAATCTGGTGGCGATCCCCGTCGATACTTTTGa                                                                     <  1:5695915/93‑1 (MQ=255)
                       tgCCTGTACCCATTTAATATCGTTAATCATCCGTCGTTGACTCCATGCCGATTCGGGTTAATCTGGTGGCGATCCCCGTCGATACTTTTGACg                                                                   <  2:2978360/93‑1 (MQ=255)
                          cTGTACCCATTTAATATCGTTAATCATCCGTCGTTGACTCCATGCCGATTCGGGTTAATCTGGTGGCGATCCCCGTCGATACTTTTGACGAAg                                                                <  2:4575164/93‑1 (MQ=255)
                          cTGTACCCATTTAATATCGTTAATCATCCGTCGTTGACTCCATGCCGATTCGGGTTAATCTGGTGGCGATCCCCGTCGATACTTTTGACGAAg                                                                >  2:2646890/1‑93 (MQ=255)
                           tGTACCCATTTAATATCGTTAATCATCCGTCGTTGACTCCATGCCGATTCGGGTTAATCTGGTGGCGATCCCCGTCGATACTTTTGACGAAgg                                                               >  1:3860522/1‑93 (MQ=255)
                           tGTACCCATTTAATATCGTTAATCATCCGTCGTTGACTCCATGCCGATTCGGGTTAATCTGGTGGCGATCCCCGTCGATACTTTTGACGAAgg                                                               >  2:6616/1‑93 (MQ=255)
                           tGTACCCATTTAATATCGTTAATCATCCGTCGTTGACTCCATGCCGATTCGGGTTAATCTGGTGGCGATCCCCGTCGATACTTTTGACGAAgg                                                               >  2:1413609/1‑93 (MQ=255)
                           tGTACCCATTTAATATCGTTAATCATCCGTCGTTGACTCCATGCCGATTCGGGTTAATCTGGTGGCGATCCCCGTCGATACTTTTGACGAAgg                                                               >  1:6925654/1‑93 (MQ=255)
                               cccATTTAATATCGTTAATCATCCGTCGTTGACTCCATGCCGATTCGGGTTAATCTGGTGGCGATCCCCGTCGATACTTTTGACGAAggcggc                                                           >  2:7406840/1‑93 (MQ=255)
                                ccATTTAATATCGTTAATCATCCGTCGTTGACTCCATGCCGATTCGGGTTAATCTGGTGGCGATCCCCGTCGATACTTTTGACGAAGGCGGCa                                                          >  1:1653935/1‑93 (MQ=255)
                                ccATTTAATATCGTTAATCATCCGTCGTTGACTCCATGCCGATTCGGGTTAATCTGGTGGCGATCCCCGTCGATACTTTTGACGAAGGCGGCa                                                          <  2:3890132/93‑1 (MQ=255)
                                 cATTTAATATCTTTAATCATCCGTCGTTGCCTCCATGCCGATTCGGGTTAATCTGGTGGCGATCCCCGTCGATACTTTTGACGAAGGAGGCAg                                                         <  2:4863810/93‑1 (MQ=255)
                                 cATTTAATATCGTTAATCATCCGTCGTTGACTCCATGCCGATTCGGGTTAATCTGGTGGCGATCCCCGTCGATACTTTTGACGAAGGCGGCAg                                                         >  2:5518654/1‑93 (MQ=255)
                                 cATTTAATATCGTTAATCATCCGTCGTTGACTCCATGCCGATTCGGGTTAATCTGGTGGCGATCCCCGTCGATACTTTTGACGAAGGCGGCAg                                                         >  2:2523331/1‑93 (MQ=255)
                                 cATTTAATATCGTTAATCATCCGTCGTTGACTCCATGCCGATTCGGGTTAATCTGGTGGCGATCCCCGTCGATACTTTTGACGAAGGCGGCAg                                                         >  2:4522301/1‑93 (MQ=255)
                                 cATTTAATATCGTTAATCATCCGTCGTTGACTCCATGCCGATTCGGGTTAATCTGGTGGCGATCCCCGTCGATACTTTTGACGAAGGCGGCAg                                                         <  1:4577247/93‑1 (MQ=255)
                                    ttAATATCGTTAATCATCCGTCGTTGACTCCATGCCGATTCGGGTTAATCTGGTGGCGATCCCCGTCGATACTTTTGACGAAGGCGGCAGGGa                                                      >  1:6523460/1‑93 (MQ=255)
                                    ttAATATCGTTAATCATCCGTCGTTGACTCCATGCCGATTCGGGTTAATCTGGTGGCGATCCCCGTCGATACTTTTGACGAAGGCGGCAGGGa                                                      <  1:1818202/93‑1 (MQ=255)
                                     tAATATCGTTAATCATCCGTCGTTGACTCCATGCCGATTCGGGTTAATCTGGTGGCGATCCCCGTCGATACTTTTGACGAAGGCGGCAGGGAt                                                     <  2:5097014/93‑1 (MQ=255)
                                     tAATATCGTTAATCATCCGTCGTTGACTCCATGACGATTCGGGTTAATCTGGTGGCGATCCCCGTCGATACTTTTTACGAAGGCGGCAGGGAt                                                     >  2:7342480/1‑93 (MQ=255)
                                      aaTATCGTTAATCATCCGTCGTTGACTCCATGCCGATTCGGGTTAATCTGGTGGCGATCCCCGTCGATACTTTTGACGAAGGCGGCAGGGATc                                                    <  1:4203013/93‑1 (MQ=255)
                                      aaTATCGTTAATCATCCGTCGTTGACTCCATGCCGATTCGGGTTAATCTGGTGGCGATCCCCGTCGATACTTTTGACGAAGGCGGCAGGGATc                                                    >  1:3051991/1‑93 (MQ=255)
                                      aaTATCGTTAATCATCCGTCGTTGACTCCATGCCGATTCGGGTTAATCTGGTGGCGATCCCCGTCGATACTTTTGACGAAGGCGGCAGGGATc                                                    >  1:4826368/1‑93 (MQ=255)
                                       atatCGTTAATCATCCGTCGTTGACTCCATGCCGATTCGGGTTAATCTGGTGGCGATCCCCGTCGATACTTTTGACGAAGGCGGCAGGGATCg                                                   <  2:4171665/93‑1 (MQ=255)
                                        tatCGTTAATCATCCGTCGTTGACTCCATGCCGATTCGGGTTAATCTGGTGGCGATCCCCGTCGATACTTTTGACGAAGGCGGCAGGGATCGc                                                  <  2:3456688/93‑1 (MQ=255)
                                        tatCGTTAATCATCCGTCGTTGACTCCATGCCGATTCGGGTTAATCTGGTGGCGATCCCCGTCGATACTCTTGAGGAAGGCGGCAGGGATCGc                                                  >  2:2692686/1‑93 (MQ=255)
                                         atCGTTAATCATCCGTCGTTGACTCCATGCCGATTCGGGTTAATCTGGTGGCGATCCCCGTCGATACTTTTGACGAAGGCGGCAGGGATCGCa                                                 <  1:7495189/93‑1 (MQ=255)
                                         atCGTTAATCATCCGTCGTTGACTCCATGCCGATTCGGGTTAATCTGGTGGCGATCCCCGTCGATACTTTTGACGAAGGCGGCAGGGATCGCa                                                 >  2:7920219/1‑93 (MQ=255)
                                          tCGTTAATCATCCGTCGTTGACTCCATGCCGATTCGGGTTAATCTGGTGGCGATCCCCGTCGATACTTTTGACGAAGGCGGCAGGGATCGCAg                                                >  1:3785635/1‑93 (MQ=255)
                                          tCGTTAATCATCCGTCGTTGACTCCATGCCGATTCGGGTTAATCTGGTGGCGATCCCCGTCGATACTTTTGACGAAGGCGGCAGGGATCGCAg                                                <  1:6628163/93‑1 (MQ=255)
                                          tCGTTAATCATCCGTCGTTAACTCCATGCCGATTCGGGTTAATCTGGTGGCGATCCCCGTCGATACTTTTGACGAAGGCGGCAGGGATCACAg                                                <  2:7028097/93‑1 (MQ=255)
                                            gTTAATCATCCGTCGTTGACTCCATGCCGATTCGGGTTAATCTGGTGGCGATCCCCGTCGATACTTTTGACGAAGGCGGCAGGGATCGCAGaa                                              <  1:3657745/93‑1 (MQ=255)
                                             ttAATCATCCGTCGTTGACTCCATGCCGATTCGGGTTAATCTGGTGGCGATCCCCGTCGATACTTTTGACGAAGGCGGCAGGGATCGCAGAAg                                             >  1:3947494/1‑93 (MQ=255)
                                             ttAATCATCCGTCGTTGACTCCATGCCGATTCGGGTTAATCTGGTGGCGATCCCCGTCGATACTTTTGACGAAGGCGGCAGGGATCGCAGAAg                                             >  2:6258356/1‑93 (MQ=255)
                                              tactcatcCGTCGTTGACTCCATGCCGATTCGGGTTAATCTGGTGGCGATCCCCGTCGATACTTTTGACGAAGGCGGCAGGGATCGCAGAAgg                                            <  1:1431002/90‑1 (MQ=255)
                                              tAATCATCCGTCGTTGACTCCATGCCGATTCGGGTTAATCTGGTGGCGATCCCCGTCGATACTTTTGACGAAGGCGGCAGGGATCGCAGAAgg                                            >  1:4110055/1‑93 (MQ=255)
                                              tAATCATCCGTCGTTGACTCCATGCCGATTCGGGTTAATCTGGTGGCGATCCCCGTCGATACTTTTGACGAAGGCGGCAGGGATCGCAGAAgg                                            <  1:3590992/93‑1 (MQ=255)
                                              tAATCATCCGTCGTTGACTCCATGCCGATTCGGGTTAATCTGGTGGCGATCCCCGTCGATACTTTTGACGAAGGCGGCAGGGATCGCAGAAgg                                            <  1:7775919/93‑1 (MQ=255)
                                               aaTCATCCGTCGTTGACTCCATGCCGATTCGGGTTAATCTGGTGGCGATCCCCGTCGATACTTTTGACGAAGGCGGCAGGGATCGCAGAAggg                                           <  1:5188285/93‑1 (MQ=255)
                                                atcatcCGTCGTCGACTACATGCCGATTCGGGTTAATCTGGTGGCGATCCCTGTCGATACTTTTGACGAAGGCGGCAGGGATCGCAGAAGGGt                                          <  2:1359097/93‑1 (MQ=255)
                                                 tcatcCGTCGTTGACTCCATGCCGATTCGGGTTAATCTGGTGGCGATCCCCGTCGATACTTTTGACGAAGGCGGCAGGGATCGCAGAAGGGtt                                         >  2:719270/1‑93 (MQ=255)
                                                 tcatcCGTCGTTGACTCCATGCCGATTCGGGTTAATCTGGTGGCGATCCCCGTCGATACTTTTGACGAAGGCGGCAGGGATCGCAGAAGGGtt                                         >  2:4008046/1‑93 (MQ=255)
                                                  catcCGTCGTTGACTCCATGCCGATTCGGGTTAATCTGGTGGCGATCCCCGTCGATACTTTTGACGAAGGCGGCAGGGATCGCAGAAGGGTTa                                        <  2:5182960/93‑1 (MQ=255)
                                                  catcCGTCGTTGACTCCATGCCGATTCGGGTTAATCTGGTGGCGATCCCCGTCGATACTTTTGACGAAGGCGGCAGGGATCGCAGAAGGGTTa                                        <  2:6103137/93‑1 (MQ=255)
                                                   atcCGTCGTTGACTCCATGCCGATTCGGGTTAATCTGGTGGCGATCCCCGTCGATACTTTTGACGAAGGCGGCAGGGATCGCAGAAGGGTTAt                                       <  2:1613327/93‑1 (MQ=255)
                                                   atcCGTCGTTGACTCCATGCCGATTCGGGTTAATCTGGTGGCGATCCCCGTCGATACTTTTGACGAAGGCGGCAGGGATCGCAGAAGGGTTAt                                       <  1:5123750/93‑1 (MQ=255)
                                                   atcCGTCGTTGACTCCATGCCGATTCGGGTTAATCTGGTGGCGATCCCCGTCGATACTTTTGACGAAGGCGGCAGGGATCGCAGAAGGGTTAt                                       >  1:2368803/1‑93 (MQ=255)
                                                   atcCGTCGTTGACTCCATGCCGATTCGGGTTAATCTGGTGGCGATCCCCGTCGATACTTTTGACGAAGGCGGCAGGGATCGCAGAAGGGTTAt                                       >  1:3315814/1‑93 (MQ=255)
                                                       gtcgtTGACTCCATGCCGATTCGGGTTAATCTGGTGGCGATCCCCGTCGATACTTTTGACGAAGGCGGCAGGGATCGCAGAAGGGTTATTGct                                   >  1:1951175/1‑93 (MQ=255)
                                                       gtcgtTGACTCCATGCCGATTCGGGTTAATCTGGTGGCGATCCCCGTCGATACTTTTGACGAAGGCGGCAGGGATCGCAGAAGGGTTATTGct                                   <  1:491568/93‑1 (MQ=255)
                                                       gtcgtTGACTCCATGCCGATTCGGGTTAATCTGGTGGCGATCCCCGTCGATACTTTTGACGAAGGCGGCAGGGATCGCAGAAGGGTTATTGct                                   >  2:5284302/1‑93 (MQ=255)
                                                          gtTGACTCCATGCCGATTCGGGTTAATCTGGTGGCGATCCCCGTCGATACTTTTGACGAAGGCGGCAGGGATCGCAGAAGGGTTATTGCTCtt                                <  2:1259238/93‑1 (MQ=255)
                                                           ttGACTCCATGCCGATTCGGGTTAATCTGGTGGCGATCCCCGTCGATACTTTTGACGAAGGCGGCAGGGATCGCAGAAGGGTTATTGCTCttt                               <  1:3554063/93‑1 (MQ=255)
                                                            tGACTCCATGCCGATTCGGGTTAATCTGGTGGCGATCCCCGTCGATACTTTTGACGAAGGCGGCAGGGATCGCAGAAGGGTTATTGCTCTTTg                              <  2:1494493/93‑1 (MQ=255)
                                                              aCTCCATGCCGATTCGGGTTAATCTGGTGGCGATCCCCGTCGATACTTTTGACGAAGGCGGCAGGGATCGCAGAAGGGTTATTGCTCTTTGcc                            >  2:4384890/1‑93 (MQ=255)
                                                                tCCATGCCGATTCGGGTTAATCTGGTGGCGATCCCCGTCGATACTTTTGACGAAGGCGGCAGGGATCGCAGAAGGGTTATTGCTCTTTGCCgg                          >  1:7354179/1‑93 (MQ=255)
                                                                 ccATGCCGATTCGGGTTAATCTGGTGGCGATCCCCGTCGATACTTTTGACGAAGGCGGCAGGGATCGCAGAAGGGTTATTGCTCTTTGCCggg                         >  2:5873304/1‑93 (MQ=255)
                                                                   aTGCCGATTCGGGTTAATCTGGTGGCGATCCCCGTCGATACTTTTGACGAAGGCGGCAGGGATCGCAGAAGGGCTATTGCTCTTTGCCggggg                       >  2:4255644/1‑92 (MQ=255)
                                                                    tGCCGATTCGGGTTAATCTGGTGGCGATCCCCGTCGATACTTTTGACGAAGGCGGCAGGGATCGCAGAAGGGTTATTGCTCTTTGCCGGGGAt                      >  2:431773/1‑93 (MQ=255)
                                                                    tGCCGATTCGGGTTAATCTGGTGGCGATCCCCGTCGATACTTTTGACGAAGGCGGCAGGGATCGCAGAAGGGTTATTGCTCTTTGCCGGGGAt                      >  2:2730030/1‑93 (MQ=255)
                                                                     gCCGATTCGGGTTAATCTGGTGGCGATCCCCGTCGATACTTTTGACGAAGGCGGCAGGGATCGCAGAAGGGTTATTGCTCTTTGCCGGGGATa                     <  1:6238253/93‑1 (MQ=255)
                                                                      ccGATTCGGGTTAATCTGGTGGCGATCCCCGTCGATACTTTTGACGAAGGCGGCAGGGATCGCAGAAGGGTTATTGCTCTTTGCCGGGGATaa                    <  2:3654213/93‑1 (MQ=255)
                                                                      ccGATTCGGGTTAATCTGGTGGCGATCCCCGTCGATACTTTTGACGAAGGCGGCAGGGATCGCAGAAGGGTTATTGCTCTTTGCCGGGGATaa                    <  2:782958/93‑1 (MQ=255)
                                                                      ccGATTCGGGTTAATCTGGTGGCGATCCCCGTCGATACTTTTGACGAAGGCGGCAGGGATCGCAGAAGGGTTATTGCTCTTTGCCGGGGATaa                    <  2:7135772/93‑1 (MQ=255)
                                                                      ccGATTCGGGTTAATCTGGTGGCGATCCCCGTCGATACTTTTGACGAAGGCGGCAGGGATCGCAGAAGGGTTATTGCTCTTTGCCGGGGATaa                    <  2:6759243/93‑1 (MQ=255)
                                                                       cGATTCGGGTTAATCTGGTGGCGATCCCCGTCGATACTTTTGACGAAGGCGGCAGGGATCGCAGAAGGGTTATTGCTCTTTGCCGGGGATAAg                   <  1:7290920/93‑1 (MQ=255)
                                                                         aTTCGGGTTAATCTGGTGGCGATCCCCGTCGATACTTTTGACGAAGGCGGCAGGGATCGCAGAAGGGTTATTGCTCTTTGCCGGGGATAGgcg                 >  1:6541377/1‑93 (MQ=255)
                                                                          ttCGGGTTAATCTGGTGGCGATCCCCGTCGATACTTTTGACGAAGGCGGCAGGGATCGCAGAAGGGTTATTGCTCTTTGCCGGGGATAAgcgc                >  2:7541468/1‑93 (MQ=255)
                                                                          ttCGGGTTAATCTGGTGGCGATCCCCGTCGATACTTGTGACGAAGGCGGCAGGGATCGCAGAAGGGTTATTGCTCTTTGCCGGGGATAAgcgc                <  1:5259148/93‑1 (MQ=255)
                                                                            cGGGTTAATCTGGTGGCGATCCCCGTCGATACTTTTGACGAAGGCGGCAGGGATCGCAGAAGGGTTATTGCTCTTTGCCGGGGATAAGCGctc              >  1:911030/1‑93 (MQ=255)
                                                                            cGGGTTAATCTGGTGGCGATCCCCGTCGATACTTTTGACGAAGGCGGCAGGGATCGCAGAAGGGTTATTGCTCTTTGCCGGGGATAAGCGctc              <  1:445577/93‑1 (MQ=255)
                                                                              ggTTAATCTGGTGGCGATCCCCGTCGATACTTTTGACGAAGGCGGCAGGGATCGCAGAAGGGTTATTGCTCTTTGCCGGGGATAAGCGCTCTg            >  2:4122641/1‑93 (MQ=255)
                                                                               gTTAATCTGGTGGCGATCCCCGTCGATACTTTTGACGAAGGCGGCAGGGATCGCAGAAGGGTTATTGCTCTTTGCCGGGGATAAGCGCTCTGc           >  2:3340180/1‑93 (MQ=255)
                                                                                  aaTCTGGTGGCGATCCCCGTCGATACTTTTGACGAAGGCGGCAGGGATCGCAGAAGGGTTATTGCTCTTTGCCGGGGATAAGCGCTCTGCTTa        >  1:1200155/1‑93 (MQ=255)
                                                                                  aaTCTGGTGGCGATCCCCGTCGATACTTTTGACGAAGGCGGCAGGGATCGCAGAAGGGTTATTGCTCTTTGCCGGGGATAAGCGCTCTGCTTa        >  2:6478872/1‑93 (MQ=255)
                                                                                  aaTCTGGTGGCGATCCCCGTCGATACTTTTGACGAAGGCGGCAGGGATCGCAGAAGGGTTATTGCTCTTTGCCGGGGATAAGCGCTCTGCTTa        <  1:7211220/93‑1 (MQ=255)
                                                                                   aTCTGGTGGCGATCCCCGTCGATACTTTTGACGAAGGCGGCAGGGATCGCAGAAGGGTTATTGCTCTTTGCCGGGGATAAGCGCTCTGCTTAt       <  1:2557928/93‑1 (MQ=255)
                                                                                   aTCTGGTGGCGATCCCCGTCGATACTTTTGACGAAGGCGGCAGGGATCGCAGAAGGGTTATTGCTCTTTGCCGGGGATAAGCGCTCTGCTTAt       <  2:1742677/93‑1 (MQ=255)
                                                                                    tCTGGTGGCGATCCCCGTCGATACTTTTGACGAAGGCGGCAGGGATCGCAGAAGGGTTATTGCTCTTTGCCGGGGATAAGCGCTCTGCTTATc      >  1:6066681/1‑93 (MQ=255)
                                                                                    tCTGGTGGCGATCCCCGTCGATACTTTTGACGAAGGCGGCAGGGATCGCAGAAGGGTTATTGCTCTTTGCCGGGGATAAGCGCTCTGCTTATc      <  2:5912339/93‑1 (MQ=255)
                                                                                      tggtggCGATCCCCGTCGATACTTTTGACGAAGGCGGCAGGGATCGCAGAAGGGTTATTGCTCTTTGCCGGGGATAAGCGCTCTGCTTATccc    >  1:5054928/1‑93 (MQ=255)
                                                                                          |                                                                                          
CGCCAGTCGGTGGCTTTACGTTGTGCCTGTACCCATTTAATATCGTTAATCATCCGTCGTTGACTCCATGCCGATTCGGGTTAATCTGGTAGCGATCCCCGTCGATACTTTTGACGAAGGCGGCAGGGATCGCAGAAGGGTTATTGCTCTTTGCCGGGGATAAGCGCTCTGCTTATCCCCA  >  NC_000913/4525496‑4525676

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: