Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,905,748 A→G S83P (TCT→CCT)  pstB ← ATP‑binding component of high‑affinity phosphate‑specific transport system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,905,7480AG100.0% 809.8 / NA 252S83P (TCT→CCT) pstBATP‑binding component of high‑affinity phosphate‑specific transport system
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (125/127);  total (125/127)

Only 100 of 257 total aligned reads displayed.

AAGCGATGTTGTCGTAGATGGACATCGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGAGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCAAACATTTTGTTGAAGGTACGCA  >  NC_000913/3905658‑3905837
                                                                                          |                                                                                         
aaGCGATGTTGTCGTAGATGGACATCGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGt                                                                                         >  2:3876484/1‑93 (MQ=255)
  gCGATGTTGTCGTAGATGGACATCGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTg                                                                                       >  1:4123183/1‑93 (MQ=255)
   cGATGTTGTCGTAGATGGACATCGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTgg                                                                                      <  2:5355906/93‑1 (MQ=255)
      tgttgtCGTAGATGGACATCGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCa                                                                                   >  2:5489465/1‑93 (MQ=255)
       gttgtCGTAGATGGACATCGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAg                                                                                  >  1:6175836/1‑93 (MQ=255)
        ttgtCGTAGATGGACATCGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAgg                                                                                 <  2:4189353/93‑1 (MQ=255)
        ttgtCGTAGATGGACATCGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAgg                                                                                 >  1:1240850/1‑93 (MQ=255)
         tgtCGTAGATGGACATCGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGa                                                                                <  1:5046739/93‑1 (MQ=255)
         tgtCGTAGATGGACATCGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGa                                                                                <  2:185688/93‑1 (MQ=255)
            cGTAGATGGACATCGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGAtgt                                                                             <  1:1395743/93‑1 (MQ=255)
            cGTAGATGGACATCGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGAtgt                                                                             <  2:4472871/93‑1 (MQ=255)
             gTAGATGGACATCGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGGTATCCTGAGGGTTGGTCAGGAtgtt                                                                            >  1:6056746/1‑93 (MQ=255)
             gTAGATGGACATCGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGAtgtt                                                                            <  2:3281445/93‑1 (MQ=255)
              tAGATGGACATCGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGAtgttg                                                                           <  1:899838/93‑1 (MQ=255)
               aGATGGACATCGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGAtgttgt                                                                          <  2:2880855/93‑1 (MQ=255)
               aGATGGACAGCGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGAtgttgt                                                                          <  2:3565660/93‑1 (MQ=255)
                  tGGACATCGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGc                                                                       >  2:3817720/1‑93 (MQ=255)
                   ggACATCGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGcc                                                                      <  2:2746134/93‑1 (MQ=255)
                   ggACATCGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGcc                                                                      <  1:473093/93‑1 (MQ=255)
                    gACATCGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCa                                                                     >  2:5194833/1‑93 (MQ=255)
                     aCATCGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCAt                                                                    >  2:3982194/1‑93 (MQ=255)
                        tCGGAAACGGCGTCGGTTTCTGGAACTCCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCaa                                                                 <  1:4975665/93‑1 (MQ=255)
                        tCGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCaa                                                                 <  2:2765835/93‑1 (MQ=255)
                          ggAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGc                                                               <  1:6045021/93‑1 (MQ=255)
                          ggAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGc                                                               <  2:5889207/93‑1 (MQ=255)
                          ggAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGc                                                               <  2:3345401/93‑1 (MQ=255)
                           gAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCa                                                              <  2:298957/93‑1 (MQ=255)
                           gAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCa                                                              >  1:926257/1‑93 (MQ=255)
                            aaaCGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAg                                                             >  1:2761867/1‑93 (MQ=255)
                             aaCGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGa                                                            >  1:531476/1‑93 (MQ=255)
                              aCGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGaa                                                           <  2:5424911/93‑1 (MQ=255)
                               cGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAAt                                                          >  2:151652/1‑93 (MQ=255)
                                ggCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAAtt                                                         >  1:3800031/1‑93 (MQ=255)
                                 gCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAAttt                                                        <  1:2016294/93‑1 (MQ=255)
                                 gCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAAttt                                                        >  1:2120760/1‑93 (MQ=255)
                                    tCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCAc                                                     >  1:4776008/1‑93 (MQ=255)
                                    tCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCAc                                                     >  1:3980702/1‑93 (MQ=255)
                                    tCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCAc                                                     >  1:2476073/1‑93 (MQ=255)
                                      ggTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCt                                                   >  1:4429298/1‑93 (MQ=255)
                                      ggTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCt                                                   >  2:1943607/1‑93 (MQ=255)
                                      ggTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCt                                                   >  2:117112/1‑93 (MQ=255)
                                      ggTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCt                                                   >  2:2167462/1‑93 (MQ=255)
                                      ggTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCt                                                   <  1:5217107/93‑1 (MQ=255)
                                       gTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCtt                                                  >  1:5842720/1‑93 (MQ=255)
                                       gTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCtt                                                  <  1:721233/93‑1 (MQ=255)
                                        tttCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTc                                                 >  1:100704/1‑93 (MQ=255)
                                         ttCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTcc                                                >  1:1484117/1‑93 (MQ=255)
                                         ttCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTcc                                                >  2:3444843/1‑93 (MQ=255)
                                          tCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCg                                               >  1:2839900/1‑93 (MQ=255)
                                            tGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCa                                             <  1:3619002/93‑1 (MQ=255)
                                            tGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCa                                             >  1:2000603/1‑93 (MQ=255)
                                            tGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCa                                             >  2:6190661/1‑93 (MQ=255)
                                            tGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCa                                             >  2:657853/1‑93 (MQ=255)
                                               aaCACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACgc                                          >  2:4655845/1‑93 (MQ=255)
                                               aaCACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACgc                                          <  1:987513/93‑1 (MQ=255)
                                                acacCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACgct                                         <  2:4802815/93‑1 (MQ=255)
                                                acacCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACgct                                         <  2:1337522/93‑1 (MQ=255)
                                                 cacCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACgctg                                        >  1:484487/1‑93 (MQ=255)
                                                  acCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACgctgc                                       <  1:1052690/93‑1 (MQ=255)
                                                   ccATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACgctgct                                      >  2:1634187/1‑93 (MQ=255)
                                                     aTGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTcc                                    >  1:321471/1‑93 (MQ=255)
                                                      tGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCg                                   <  2:2122762/93‑1 (MQ=255)
                                                       gCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCATAATTTCACCTTCCGCACGCTGCTCCgg                                  <  2:1592563/93‑1 (MQ=255)
                                                        cccACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCggg                                 <  1:2194037/93‑1 (MQ=255)
                                                          cACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTa                               >  2:1399244/1‑93 (MQ=255)
                                                          cACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTa                               >  2:2623319/1‑93 (MQ=255)
                                                          cACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTa                               >  1:3638447/1‑93 (MQ=255)
                                                            cTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACa                             >  2:2540755/1‑93 (MQ=255)
                                                            cTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACa                             >  2:4276564/1‑93 (MQ=255)
                                                            cTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACa                             <  2:3190990/93‑1 (MQ=255)
                                                                   acgcaGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCaaa                      >  2:4409485/1‑93 (MQ=255)
                                                                   acgcaGCAGTGCGATATCCTGAGGGGTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCaaa                      >  1:3105393/1‑93 (MQ=255)
                                                                    cgcaGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCAAAc                     <  1:3780739/93‑1 (MQ=255)
                                                                     gcagcaGTGCGATATTCTGAGGGTTTTTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTGCGCACGCTGCTCCGGGTACAGTTCAAACa                    >  2:2487589/1‑93 (MQ=255)
                                                                       agcagTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCAAACAtt                  <  1:1724496/93‑1 (MQ=255)
                                                                        gcagTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCAAACAttt                 <  1:4711053/93‑1 (MQ=255)
                                                                        gcagTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCAAACAttt                 <  2:5752859/93‑1 (MQ=255)
                                                                        gcagTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCAAACAttt                 >  1:5194610/1‑93 (MQ=255)
                                                                        gcagTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCAAACAttt                 >  2:2074597/1‑93 (MQ=255)
                                                                         cagTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCTGGGTACAGTTCAAACAtttt                <  2:5143014/93‑1 (MQ=255)
                                                                         cagTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCAAACAtttt                <  2:5170832/93‑1 (MQ=255)
                                                                          agTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCAAACATTttg               >  1:779878/1‑93 (MQ=255)
                                                                           gTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCAAACATTttgt              <  1:3460735/93‑1 (MQ=255)
                                                                           gTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCAAACATTttgt              <  1:1859169/93‑1 (MQ=255)
                                                                              ctatatCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCAAACATTTTGTTGa           <  1:3028120/91‑1 (MQ=255)
                                                                              cGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCAAACATTTTGTTGa           <  2:2615924/93‑1 (MQ=255)
                                                                              cGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCAAACATTTTGTTGa           >  1:262687/1‑93 (MQ=255)
                                                                               gATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCAAACATTTTGTTGaa          <  2:5098137/93‑1 (MQ=255)
                                                                                atatCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCAAACATTTTGTTGAAg         >  1:5318969/1‑93 (MQ=255)
                                                                                 tatCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCAAACATTTTGTTGAAgg        >  1:2923673/1‑93 (MQ=255)
                                                                                 tatCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCAAACATTTTGTTGAAgg        <  2:6101581/93‑1 (MQ=255)
                                                                                  atCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCAAACATTTTGTTGAAGGt       <  1:2504248/93‑1 (MQ=255)
                                                                                  atCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCAAACATTTTGTTGAAGGt       >  1:117771/1‑93 (MQ=255)
                                                                                  atCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCAAACATTTTGTTGAAGGt       <  2:3079270/93‑1 (MQ=255)
                                                                                  atCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCAAACATTTTGTTGAAGGt       <  2:2346334/93‑1 (MQ=255)
                                                                                   tCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCAAACATTTTGTTGAAGGTa      >  2:3630056/1‑93 (MQ=255)
                                                                                   tCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCAAACATTTTGTTGAAGGTa      <  1:1567644/93‑1 (MQ=255)
                                                                                   tCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCAAACATTTTGTTGAAGGTa      >  2:4910185/1‑93 (MQ=255)
                                                                                    ccTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCAAACATTTTGTTGAAGGTAc     <  1:5680727/93‑1 (MQ=255)
                                                                                      tGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCAAACATTTTGTTGAAGGTACgc   >  1:3327921/1‑93 (MQ=255)
                                                                                          |                                                                                         
AAGCGATGTTGTCGTAGATGGACATCGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGAGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCAAACATTTTGTTGAAGGTACGCA  >  NC_000913/3905658‑3905837

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