Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,803,648 T→C S271S (AGT→AGC)  proV → ATP‑binding component of transport system for glycine, betaine and proline

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,803,6480TC99.0% 734.8 / ‑2.4 206S271S (AGT→AGC) proVATP‑binding component of transport system for glycine, betaine and proline
Reads supporting (aligned to +/- strand):  ref base T (0/2);  new base C (101/103);  total (101/105)
Fisher's exact test for biased strand distribution p-value = 4.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.85e-01

Only 100 of 208 total aligned reads displayed.

TGAAGTGGTACAGGTCGGCACACCGGATGAAATTCTCAATAATCCGGCGAATGATTATGTCCGTACCTTCTTCCGTGGCGTTGATATTAGTCAGGTATTCAGTGCGAAAGATATTGCCCGCCGGACACCGAATGGCTTAATTCGTAAAACCCCTGGCTTCGGCCCACGTTCGGCACTGAA  >  NC_000913/2803558‑2803737
                                                                                          |                                                                                         
tGAAGTGGTACAGGTCGGCACACCGGATGAAATTCTCAATAATCCGGCGAATGATTATGTCCGTACCTTCTTCCGTGGCGTTGATATTAGCCa                                                                                         >  2:3890889/1‑93 (MQ=255)
  aaGTGGTACAGGTCGGCACACCGGATGAAATTCTCAATAATCCGGCGAATGATTATGTCCGTACCTTCTTCCGTGGCGTTGATATTAGCCAgg                                                                                       >  2:2500405/1‑93 (MQ=255)
    gTGGTACAGGTCGGCACACCGGATGAAATTCTCAATAATCCGGCGAATGATTATGTCCGTACCTTCTTCCGTGGCGTTGATATTAGCCAGGTa                                                                                     >  2:5511196/1‑93 (MQ=255)
       gTACAGGTCGGCACACCGGATGAAATTCTCAATAATCCGGCGAATGATTATGTCCGTACCTTCTTCCGTGGCGTTGATATTAGCCAGGTATTc                                                                                  <  1:4266173/93‑1 (MQ=255)
        tACAGGTCGGCACACCGGATGAAATTCTCAATAATCCGGCGAATGATTATGTCCGTACCTTCTTCCGTGGCGTTGATATTAGCCAGGTATTCa                                                                                 >  1:2719715/1‑93 (MQ=255)
        tACAGGTCGGCACACCGGATGAAATTCTCAATAATCCGGCGAATGATTATGTCCGTACCTTCTTCCGTGGCGTTGATATTAGCCAGGTATTCa                                                                                 >  2:2316537/1‑93 (MQ=255)
        tACAGGTCGGCACACCGGATGAAATTCTCAATAATCCGGCGAATGATTATGTCCGTACCTTCTTCCGTGGCGTTGATATTAGCCAGGTATTCa                                                                                 >  2:4856640/1‑93 (MQ=255)
          cAGGTCGGCACACCGGATGAAATTCTCAATAATCCGGCGAATGATTATGTCCGTACCTTCTTCCGTGGCGTTGATATTAGCCAGGTATTCAGt                                                                               <  2:3412031/93‑1 (MQ=255)
           aGGTCGGCACACCGGATGAAATTCTCAATAATCCGGCGAATGATTATGTCCGTACCTTCTTCCGTGGCGTTGATATTAGCCAGGTATTCAGTg                                                                              <  1:5793470/93‑1 (MQ=255)
            ggTCGGCACACCGGATGAAATTCTCAATAATCCGGCGAATGATTATGTCCGTACCTTCTTCCGTGGCGTTGATATTAGCCAGGTATTCAGTGc                                                                             <  1:4136934/93‑1 (MQ=255)
             gTCGGCACACCGGATGAAATTCTCAATAATCCGGCGAATGATTATGTCCGTACCTTCTTCCGTGGCGTTGATATTAGCCAGGTATTCAGTGCg                                                                            >  1:4826944/1‑93 (MQ=255)
              tCGGCACACCGGATGAAATTCTCAATAATCCGGCGAATGATTATGTCCGTACCTTCTTCCGTGGCGTTGATATTAGCCAGGTATTCAGTGCGa                                                                           <  1:1257544/93‑1 (MQ=255)
               cGGCACACCGGATGAAATTCTCAATAATCCGGCGAATGATTATGTCCGTACCTTCTTCCGTGGCGTTGATATTAGCCAGGTATTCAGTGCGaa                                                                          <  2:2594028/93‑1 (MQ=255)
                ggCACACCGGATGAAATTCTCAATAATCCGGCGAATGATTATGTCCGTACCTTCTTCCGTGGCGTTGATATTAGCCAGGTATTCAGTGCGaaa                                                                         <  1:4918878/93‑1 (MQ=255)
                  cacaCCGGATGAAATTCTCAATAATCCGGCGAATGATTATGTCCGTACCTTCTTCCGTGGCGTTGATATTAGCCAGGTATTCAGTGCGAAAGa                                                                       <  1:3299472/93‑1 (MQ=255)
                   acacCGGATGAAATTCTCAATAATCCGGCGAATGATTATGTCCGTACCTTCTTCCGTGGCGTTGATATTAGCCAGGTATTCAGTGCGAAAGat                                                                      >  2:773350/1‑93 (MQ=255)
                     acCGGATGAAATTCTCAATAATCCGGCGAATGATTATGTCCGTACCTTCTTCCGTGGCGTTGATATTAGCCAGGTATTCAGTGCGAAAGatat                                                                    <  1:2893289/93‑1 (MQ=255)
                      ccGGATGAAATTCTCAATAATCCGGCGAATGATTATGTCCGTACCTTCTTCCGTGGCGTTGATATTAGCCAGGTATTCAGTGCGAAAGATAtt                                                                   <  1:784764/93‑1 (MQ=255)
                      ccGGATGAAATTCTCAATAATCCGGCGAATGATTATGTCCGTACCTTCTTCCGTGGCGTTGATATTAGCCAGGTATTCAGTGCGAAAGATAtt                                                                   <  1:3853358/93‑1 (MQ=255)
                       cGGATGAAATTCTCAATAATCCGGCGAATGATTATGTCCGTACCTTCTTCCGTGGCGTTGATATTAGCCAGGTATTCAGTGCGAAAGATATTg                                                                  >  2:5748114/1‑93 (MQ=255)
                        ggATGAAATTCTCAATAATCCGGCGAATGATTATGTCCGTACCTTCTTCCGTGGCGTTGATATTAGCCAGGTATTCAGTGCGAAAGATATTGc                                                                 <  1:1073430/93‑1 (MQ=255)
                         gATGAAATTCTCAATAATCCGGCGAATGATTATGTCCGTACCTTCTTCCGTGGCGTTGATATTAGCCAGGTATTCAGTGCGAAAGATATTGcc                                                                <  1:1876489/93‑1 (MQ=255)
                         gATGAAATTCTCAATAATCCGGCGAATGATTATGTCCGTACCTTCTTCCGTGGCGTTGATATTAGCCAGGTATTCAGTGCGAAAGATATTGcc                                                                >  2:5883988/1‑93 (MQ=255)
                          aTGAAATTCTCAATAATCCGGCGAATGATTATGTCCGTACCTTCTTCCGTGGCGTTGATATTAGCCAGGTATTCAGTGCGAAAGATATTGccc                                                               >  1:1825117/1‑93 (MQ=255)
                          aTGAAATTCTCAATAATCCGGCGAATGATTATGTCCGTACCTTCTTCCGTGGCGTTGATATTAGCCAGGTATTCAGTGCGAAAGATATTGccc                                                               >  2:1465590/1‑93 (MQ=255)
                          aTGAAATTCTCAATAATCCGGCGAATGATTATGTCCGTACCTTCTTCCGTGGCGTTGATATTAGCCAGGTATTCAGTGCGAAAGATATTGccc                                                               <  2:4347041/93‑1 (MQ=255)
                           tGAAATTCTCAATAATCCGGCGAATGATTATGTCCGTACCTTCTTCCGTGGCGTTGATATTAGCCAGGTATTCAGTGCGAAAGATATTGCccg                                                              >  2:5096561/1‑93 (MQ=255)
                           tGAAATTCTCAATAATCCGGCGAATGATTATGTCCGTACCTTCTTCCGTGGCGTTGATATTAGCCAGGTATTCAGTGCGAAAGATATTGCccg                                                              >  1:3726493/1‑93 (MQ=255)
                           tGAAATTCTCAATAATCCGGCGAATGATTATGTCCGTACCTTCTTCCGTGGCGTTGATATTAGCCAGGTATTCAGTGCGAAAGATATTGCccg                                                              >  1:3770635/1‑93 (MQ=255)
                                ttCTCAATAATCCGGCGAATGATTATGTCCGTACCTTCTTCCGTGGCGTTGATATTAGCCAGGTATTCAGTGCGAAAGATATTGCCCGCCGGa                                                         <  2:4213605/93‑1 (MQ=255)
                                ttCTCAATAATCCGGCGAATGATTATGTCCGTACCTTCTTCCGTGGCGTTGATATTAGCCAGGTATTCAGTGCGAAAGATATTGCCCGCCGGa                                                         >  1:5568680/1‑93 (MQ=255)
                                ttCTCAATAATCCGGCGAATGATTATGTCCGTACCTTCTTCCGTGGCGTTGATATTAGCCAGGTATTCAGTGCGAAAGATATTGCCCGCCGGa                                                         <  1:3535649/93‑1 (MQ=255)
                                   tcAATACTCCGGCGACTGATTATGTCCGTACCTTCTTCCGTGGCGTTGATATTAGCCAGGTATTCAGTGCGAAAGTTATTGCCCGCCGGacac                                                      <  1:5478626/93‑1 (MQ=255)
                                   tcAATAATCCGGCGAATGATTATGTCCGTACCTTCTTCCGTGGCGTTGATATTAGCCAGGTATTCAGTGCGAAAGATATTGCCCGCCGGacac                                                      <  2:671761/93‑1 (MQ=255)
                                   tcAATAATCCGGCGAATGATTATGTCCGTACCTTCTTCCGTGGCGTTGATATTAGCCAGGTATTCAGTGCGAAAGATATTGCCCGCCGGacac                                                      <  1:5581110/93‑1 (MQ=255)
                                   tcAATAATCCGGCGAATGATTATGTCCGTACCTTCTTCCGTGGCGTTGATATTAGCCAGGTATTCAGTGCGAAAGATATTGCCCGCCGGacac                                                      <  1:2935449/93‑1 (MQ=255)
                                   tcAATAATCCGGCGAATGATTATGTCCGTACCTTCTTCCGTGGCGTTGATATTAGCCAGGTATTCAGTGCGAAAGATATTGCCCGCCGGacac                                                      <  1:1718752/93‑1 (MQ=255)
                                   tcAATAATCCGGCGAATGATTATGTCCGTACCTTCTTCCGTGGCGTTGATATTAGCCAGGTATTCAGTGCGAAAGATATTGCCCGCCGGacac                                                      <  1:5198606/93‑1 (MQ=255)
                                       taatCCGGCGAATGATTATGTCCGTACCTTCTTCCGTGGCGTTGATATTAGCCAGGTATTCAGTGCGAAAGATATTGCCCGCCGGACACCGaa                                                  >  2:548201/1‑93 (MQ=255)
                                        aatCCGGCGAATGATTATGTCCGTACCTTCTTCCGTGGCGTTGATATTAGCCAGGTATTCAGTGCGAAAGATATTGCCCGCCGGACACCGAAt                                                 <  1:454400/93‑1 (MQ=255)
                                         atCCGGCGAATGATTATGTCCGTACCTTCTTCCGTGGCGTTGATATTAGCCAGGTATTCAGTGCGAAAGATATTGCCCGCCGGACACCGAATg                                                <  1:1848373/93‑1 (MQ=255)
                                          tCCGGCGAATGATTATGTCCGTACCTTCTTCCGTGGCGTTGATATTAGCCAGGTATTCAGTGCGAAAGATATTGCCCGCCGGACACCGAATgg                                               >  1:5593927/1‑93 (MQ=255)
                                           ccGGCGAATGATTATGTCCGTACCTTCTTCCGTGGCGTTGATATTAGCCAGGTATTCAGTGCGAAAGATATTGCCCGCCGGACACCGAATGGc                                              <  2:2392109/93‑1 (MQ=255)
                                              gCGAATGATTATGTCCGTACCTTCTTCCGTGGCGTTGATATTAGCCAGGTATTCAGTGCGAAAGATATTGCCCGCCGGACACCGAATGGCTTa                                           >  1:335898/1‑93 (MQ=255)
                                                  aTGATTATGTCCGTACCTTCTTCCGTGGCGTTGATATTAGCCAGGTATTCAGTGCGAAAGATATTGCCCGCCGGACACCGAATGGCTTAATTc                                       >  2:4330482/1‑93 (MQ=255)
                                                    gATTATGTCCGTACCTTCTTCCGTGGCGTTGATATTAGCCAGGTATTCAGTGCGAAAGATATTGCCCGCCGGACACCGAATGGCTTAATTCGt                                     >  2:3460755/1‑93 (MQ=255)
                                                    gATTATGTCCGTACCTTCTTCCGTGGCGTTGATATTAGCCAGGTATTCAGTGCGAAAGATATTGCCCGCCGGACACCGAATGGCTTAATTCGt                                     <  1:5501955/93‑1 (MQ=255)
                                                    gATTATGTCCGTACCTTCTTCCGTGGCGTTGATATTAGCCAGGTATTCAGTGCGAAAGATATTGCCCGCCGGACACCGAATGGCTTAATTCGt                                     >  2:1982944/1‑93 (MQ=255)
                                                       tATGTCCGTACCTTCTTCCGTGGCGTTGATATTAGCCAGGTATTCAGTGCGAAAGATATTGCCCGCCGGACACCGAATGGCTTAATTCGTaaa                                  >  2:4933736/1‑93 (MQ=255)
                                                        aTGTCCGTACCTTCTTCCGTGGCGTTGATATTAGCCAGGTATTCAGTGCGAAAGATATTGCCCGCCGGACACCGAATGGCTTAATTCGTaaaa                                 <  2:5357176/93‑1 (MQ=255)
                                                        aTGTCCGTACCTTCTTCCGTGGCGTTGATATTAGCCAGGTATTCAGTGCGAAAGATATTGCCCGCCGGACACCGAATGGCTTAATTCGTaaaa                                 <  1:1653989/93‑1 (MQ=255)
                                                        aTGTCCGTACCTTCTTCCGTGGCGTTGATATTAGCCAGGTATTCAGTGCGAAAGATATTGCCCGCCGGACACCGAATGGCTTAATTCGTaaaa                                 <  2:1837011/93‑1 (MQ=255)
                                                        aTGTCCGTACCTTCTTCCGTGGCGTTGATATTAGCCAGGTATTCAGTGCGAAAGATATTGCCCGCCGGACACCGAATGGCTTAATTCGTaaaa                                 >  2:2344324/1‑93 (MQ=255)
                                                          gTCCGTACCTTCTTCCGTGGCGTTGATATTAGCCAGGTATTCAGTGCGAAAGATATTGCCCGCCGGACACCGAATGGCTTAATTCGTAAAAcc                               <  2:1157763/93‑1 (MQ=255)
                                                          gTCCGTACCTTCTTCCGTGGCGTTGATATTAGCCAGGTATTCAGTGCGAAAGATATTGCCCGCCGGACACCGAATGGCTTAATTCGTAAAAcc                               >  1:5096728/1‑93 (MQ=255)
                                                            ccGTACCTTCTTCCGTGGCGTTGATATTAGCCAGGTATTCAGTGCGAAAGATATTGCCCGCCGGACACCGAATGGCTTAATTCGTAAAAcccc                             <  2:3210271/93‑1 (MQ=255)
                                                            ccGTACCTTCTTCCGTGGCGTTGATATTAGCCAGGTATTCAGTGCGAAAGATATTGCCCGCCGGACACCGAATGGCTTAATTCGTAAAAcccc                             <  1:3578271/93‑1 (MQ=255)
                                                            ccGTACCTTCTTCCGTGGCGTTGATATTAGCCAGGTATTCAGTGCGAAAGATATTGCCCGCCGGACACCGAATGGCTTAATTCGTAAAAcccc                             <  1:308196/93‑1 (MQ=255)
                                                            ccGTACCTTCTTCCGTGGCGTTAATATTAGCCAGGTATTCAGTGCGAAAGATATTGCCCGCCGGACACCGAATGGCTTAATTCGTAAAAcccc                             <  2:1478208/93‑1 (MQ=255)
                                                             cGTACCTTTTTTCGTGGCGTTCATATTAGCCAGGTATTCAGTGCGAAAGATATTGCCCGCCGGACACCGAATGGCTTAATTCGTAAAACCCCt                            <  1:1501387/93‑1 (MQ=255)
                                                             cGTACCTTCTTCCGTGGCGTTGATATTAGCCAGGTATTCAGTGCGAAAGATATTGCCCGCCGGACACCGAATGGCTTAATTCGTAAAACCCCt                            >  1:4490785/1‑93 (MQ=255)
                                                             cGTACCTTCTTCCGTGGCGTTGATATTAGCCAGGTATTCAGTGCGAAAGATATTGCCCGCCGGACACCGAATGGCTTAATTCGTAAAACCCCt                            >  2:2179926/1‑93 (MQ=255)
                                                             cGTACCTTCTTCCGTGGCGTTGATATTAGCCAGGTATTCAGTGCGAAAGATATTGCCCGCCGGACACCGAATGGCTTAATTCGTAAAACCCCt                            <  1:6027584/93‑1 (MQ=255)
                                                              gTACCTTCTTCCGTGGCGTTGATATTAGCCAGGTATTCAGTGCGAAAGATATTGCCCGCCGGACACCGAATGGCTTAATTCGTAAAACCCCTg                           <  1:4874742/93‑1 (MQ=255)
                                                              gTACCTTCTTCCGTGGCGTTGATATTAGCCAGGTATTCAGTGCGAAAGATATTGCCCGCCGGACACCGAATGGCTTAATTCGTAAAACCCCTg                           <  1:4094425/93‑1 (MQ=255)
                                                              gTACCTTCTTCCGTGGCGTTGATATTAGCCAGGTATTCAGTGCGAAAGATATTGCCCGCCGGACACCGAATGGCTTAATTCGTAAAACCCCTg                           >  1:5939252/1‑93 (MQ=255)
                                                              gTACCTTCTTCCGTGGCGTTGATATTAGCCAGGTATTCAGTGCGAAAGATATTGCCCGCCGGACACCGAATGGCTTAATTCGTAAAACCCCTg                           >  2:143222/1‑93 (MQ=255)
                                                               tACCTTCTTCCGTGGCGTTGATATTAGCCAGGTATTCAGTGCGAAAGATATTGCCCGCCGGACACCGAATGGCTTAATTCGTAAAACCCCTgg                          >  2:6008683/1‑93 (MQ=255)
                                                                aCCTTCTTCCGTGGCGTTGATATTAGCCAGGTATTCAGTGCGAAAGATATTGCCCGCCGGACACCGAATGGCTTAATTCGTAAAACCCCTGGc                         >  2:5619485/1‑93 (MQ=255)
                                                                aCCTTCTTCCGTGGCGTTGATATTAGCCAGGTATTCAGTGCGAAAGATATTGCCCGCCGGACACCGAATGGCTTAATTCGTAAAACCCCTGGc                         >  2:857683/1‑93 (MQ=255)
                                                                aCCTTCTTCCGTGGCGTTGATATTAGCCAGGTATTCAGTGCGAAAGATATTGCCCGCCGGACACCGAATGGCTTAATTCGTAAAACCCCTGGc                         <  1:2205850/93‑1 (MQ=255)
                                                                 ccTTCTTCCGTGGCGTTGATATTAGCCAGGTATTCAGTGCGAAAGATATTGCCCGCCGGACACCGAATGGCTTAATTCGTAAAACCCCTGGCt                        <  2:1915973/93‑1 (MQ=255)
                                                                  cttcttCCGTGGCGTTGATATTAGCCAGGTATTCAGTGCGAAAGATATTGCCCGCCGGACACCGAATGGCTTAATTCGTAAAACCCCTGGCtt                       <  1:3153147/93‑1 (MQ=255)
                                                                   ttcttcCGTGGCGTTGATATTAGCCAGGTATTCAGTGCGAAAGATATTGCCCGCCGGACACCGAATGGCTTAATTCGTAAAACCCCTGGCTTc                      >  2:336346/1‑93 (MQ=255)
                                                                   ttcttcCGTGGCGTTGATATTAGCCAGGTATTCAGTGCGAAAGATATTGCCCGCCGGACACCGAATGGCTTAATTCGTAAAACCCCTGGCTTc                      >  1:2601923/1‑93 (MQ=255)
                                                                   ttcttcCGTGGCGTTGATATTAGCCAGGTATTCAGTGCGAAAGATATTGCCCGCCGGACACCGAATGGCTTAATTCGTAAAACCCCTGGCTTc                      >  2:6041211/1‑93 (MQ=255)
                                                                    tcttcCGTGGCGTTGATATTAGCCAGGTATTCAGTGCGAAAGATATTGCCCGCCGGACACCGAATGGCTTAATTCGTAACACCCCTGGCTTCg                     >  2:4019390/1‑93 (MQ=255)
                                                                    tcttcCGTGGCGTTGATATTAGCCAGGTATTCAGTGCGAAAGATATTGCCCGCCGGACACCGAATGGCTTAATTCGTAAAACCCCTGGCTTCg                     >  2:2718650/1‑93 (MQ=255)
                                                                     cttcCGTGGCGTTGATATTAGCCAGGTATTCAGTGCGAAAGATATTGCCCGCCGGACACCGAATGGCTTAATTCGTAAAACCCCTGGCTTCgg                    >  1:2005587/1‑93 (MQ=255)
                                                                     cttcCGTGGCGTTGATATTAGCCAGGTATTCAGTGCGAAAGATATTGCCCGCCGGACACCGAATGGCTTAATTCGTAAAACCCCTGGCTTCgg                    >  2:1341950/1‑93 (MQ=255)
                                                                         cGTGGCGTTGATATTAGCCAGGTATTCAGTGCGAAAGATATTGCCCGCCGGACACCGAATGGCTTAATTCGTAAAACCCCTGGCTTCGGCCCa                <  2:3483993/93‑1 (MQ=255)
                                                                         cGTGGCGTTGATATTAGCCAGGTATTCAGTGCGAAAGATATTGCCCGCCGGACACCGAATGGCTTAATTCGTAAAACCCCTGGCTTCGGCCCa                <  1:1742121/93‑1 (MQ=255)
                                                                         cGTGGCGTTGATATTAGCCAGGTATTCAGTGCGAAAGATATTGCCCGCCGGACACCGAATGGCTTAATTCGTAAAACCCCTGGCTTCGGCCCa                >  2:3331777/1‑93 (MQ=255)
                                                                           tGGCGTTGATATTAGCCAGGTATTCAGTGCGAAAGATATTGCCCGCCGGACACCGAATGGCTTAATTCGTAAAACCCCTGGCTTCGGCCCACg              >  2:2972602/1‑93 (MQ=255)
                                                                           tGGCGTTGATATTAGCCAGGTATTCAGTGCGAAAGATATTGCCCGCCGGACACCGAATGGCTTAATTCGTAAAACCCCTGGCTTCGGCCCACg              >  2:5429129/1‑93 (MQ=255)
                                                                            ggCGTTGATATTAGCCAGGTATTCAGTGCGAAAGATATTGCCCGCCGGACACCGAATGGCTTAATTCGTAAAACCCCTGGCTTCGGCCCACGt             <  2:5899119/93‑1 (MQ=255)
                                                                             gCGTTGATATTAGCCAGGTATTCAGTGCGAAAGATATTGCCCGCCGGACACCGAATGGCTTAATTCGTAAAACCCCTGGCTTCGGCCCACGtt            >  1:3806770/1‑93 (MQ=255)
                                                                             gCGTTGATATTAGCCAGGTATTCAGTGCGAAAGATATTGCCCGCCGGACACCGAATGGCTTAATTCGTAAAACCCCTGGCTTCGGCCCACGtt            >  2:4735865/1‑93 (MQ=255)
                                                                             gCGTTGATATTAGCCAGGTATTCAGTGCGAAAGATATTGCCCGCCGGACACCGAATGGCTTAATTCGTAAAACCCCTGGCTTCGGCCCACGtt            <  2:1642027/93‑1 (MQ=255)
                                                                             gCGTTGATATTAGCCAGGTATTCAGTGCGAAAGATATTGCCCGCCGGACACCGAATGGCTTAATTCGTAAAACCCCTGGCTTCGGCCCACGtt            <  2:4117524/93‑1 (MQ=255)
                                                                              cGTTGATATTAGCCAGGTATTCAGTGCGAAAGATATTGCCCGCCGGACACCGAATGGCTTAATTCGTAAAACCCCTGGCTTCGGCCCACGTTc           <  2:2784758/93‑1 (MQ=255)
                                                                               gTTGATATTAGCCAGGTATTCAGTGCGAAAGATATTGCCCGCCGGACACCGAATGGCTTAATTCGTAAAACCCCTGGCTTCGGCCCACGTTCg          >  1:2656204/1‑93 (MQ=255)
                                                                               gTTGATATTAGCCAGGTATTCAGTGCGAAAGATATTGCCCGCCGGACACCGAATGGCTTAATTCGTAAAACCCCTGGCTTCGGCCCACGTTCg          >  1:139547/1‑93 (MQ=255)
                                                                                ttGATATTAGCCAGGTATTCAGTGCGAAAGATATTGCCCGCCGGACACCGAATGGCTTAATTCGTAAAACCCCTGGCTTCGGCCCACGTTCgg         >  1:4205918/1‑93 (MQ=255)
                                                                                ttGATATTAGCCAGGTATTCAGTGCGAAAGATATTGCCCGCCGGACACCGAATGGCTTAATTCGTAAAACCCCTGGCTTCGGCCCACGTTCgg         >  2:4559442/1‑93 (MQ=255)
                                                                                 tGATATTAGCCAGGTATTCAGTGCGAAAGATATTGCCCGCCGGACACCGAATGGCTTAATTCGTAAAACCCCTGGCTTCGGCCCACGTTCGGc        <  1:1140357/93‑1 (MQ=255)
                                                                                    tatTAGCCTGGTATTCAGTGCGAAAGATATTGCCCGCCGGACACCGAATGGCTTAATTCGTAAAACCCCTGGCTTCGGCCCACGTTCGGCACt     >  1:5635520/1‑93 (MQ=255)
                                                                                    tatTAGCCAGGTATTCAGTGCGAAAGATATTGCCCGCCGGACACCGAATGGCTTAATTCGTAAAACCCCTGGCTTCGGCCCACGTTCGGCACt     >  2:728791/1‑93 (MQ=255)
                                                                                    tatTAGCCAGGTATTCAGTGCGAAAGATATTGCCCGCCGGACACCGAATGGCTTAATTCGTAAAACCCCTGGCTTCGGCCCACGTTCGGCACt     <  1:1200181/93‑1 (MQ=255)
                                                                                       tAGCCAGGTATTCAGTGCGAAAGATATTGCCCGCCGGACACCGAATGGCTTAATTCGTAAAACCCCTGGCTTCGGCCCACGTTCGGCACTGaa  <  2:2907991/93‑1 (MQ=255)
                                                                                          |                                                                                         
TGAAGTGGTACAGGTCGGCACACCGGATGAAATTCTCAATAATCCGGCGAATGATTATGTCCGTACCTTCTTCCGTGGCGTTGATATTAGTCAGGTATTCAGTGCGAAAGATATTGCCCGCCGGACACCGAATGGCTTAATTCGTAAAACCCCTGGCTTCGGCCCACGTTCGGCACTGAA  >  NC_000913/2803558‑2803737

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