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,7480AG99.6% 893.6 / ‑4.8 277S83P (TCT→CCT) pstBATP‑binding component of high‑affinity phosphate‑specific transport system
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (140/136);  total (140/137)
Fisher's exact test for biased strand distribution p-value = 4.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.39e-01

Only 100 of 288 total aligned reads displayed.

AAGCGATGTTGTCGTAGATGGACATCGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGAGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCAAACATTTTGTTGAAGGTACGCA  >  NC_000913/3905658‑3905837
                                                                                          |                                                                                         
aaGCGATGTTGTCGTAGATGGACATCGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGt                                                                                         <  1:4206513/93‑1 (MQ=255)
 aGCGATGTTGTCGTAGATGGACATCGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGtt                                                                                        >  1:3869937/1‑93 (MQ=255)
  gCGATGTTGTCGTAGATGGACATCGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTg                                                                                       >  1:2927026/1‑93 (MQ=255)
   cGATGTTGTCGTAGATGGACATCGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTgg                                                                                      <  1:584698/93‑1 (MQ=255)
   cGATGTTGTCGTAGATGGACATCGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTgg                                                                                      <  1:1940517/93‑1 (MQ=255)
   cGATGTTGTCGTAGATGGACATCGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTgg                                                                                      <  2:4903050/93‑1 (MQ=255)
     aTGTTGTCGTAGATGGACATCGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTc                                                                                    <  2:2074144/93‑1 (MQ=255)
     aTGTTGTCGTAGATGGACATCGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTc                                                                                    >  2:5282995/1‑93 (MQ=255)
       gttgtCGTAGATGGACATCGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAg                                                                                  >  1:3728845/1‑93 (MQ=255)
       gttgtCGTAGATGGACATCGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAg                                                                                  >  2:3489440/1‑93 (MQ=255)
         tgtCGTAGATGGACATCGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGa                                                                                <  1:3980306/93‑1 (MQ=255)
          gtCGTAGATGGACATCGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGAt                                                                               <  2:365617/93‑1 (MQ=255)
           tCGTAGATGGACATCGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGAtg                                                                              >  1:2632663/1‑93 (MQ=255)
            cGTAGATGGACATCGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGAtgt                                                                             <  1:5966404/93‑1 (MQ=255)
             gTAGATTGACATCGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGAtgtt                                                                            <  2:5374460/93‑1 (MQ=255)
             gTAGATGGACATCGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGAtgtt                                                                            <  2:4510708/93‑1 (MQ=255)
             gTAGATGGACATCGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGAtgtt                                                                            >  1:3292951/1‑93 (MQ=255)
             gTAGATGGACATCGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGAtgtt                                                                            >  1:767799/1‑93 (MQ=255)
              tAGATGGACATCGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGAtgttg                                                                           >  1:1604513/1‑93 (MQ=255)
              tAGATGGACATCGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGAtgttg                                                                           <  2:3259038/93‑1 (MQ=255)
               aGATGGACATCGGAACCGGCGTCGGGTTCTGGAACACCATGCCCACTTTCGCACGCACCAGTGCGATATCCTGAGCGTTGGTCAGGAtgttgt                                                                          >  2:1037123/1‑93 (MQ=255)
               aGATGGACATCGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGAtgttgt                                                                          <  1:2804212/93‑1 (MQ=255)
                gATGGACATTGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTc                                                                         >  2:1496252/1‑93 (MQ=255)
                gATGGACATCGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTc                                                                         >  1:1827166/1‑93 (MQ=255)
                gATGGACATCGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTc                                                                         <  2:1834371/93‑1 (MQ=255)
                 aTGGACATCGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCg                                                                        <  2:3963322/93‑1 (MQ=255)
                   ggACATCGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGcc                                                                      <  2:5063016/93‑1 (MQ=255)
                   ggACATCGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGcc                                                                      >  2:406129/1‑93 (MQ=255)
                   ggACATCGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGcc                                                                      <  2:6666242/93‑1 (MQ=255)
                   ggACATCGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGGTGGTCAGGATGTTGTCGcc                                                                      >  1:6124098/1‑93 (MQ=255)
                     aCATCGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCAt                                                                    >  1:3146214/1‑93 (MQ=255)
                      cATCGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATc                                                                   >  2:1851718/1‑93 (MQ=255)
                      cATCGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCACTGCGATATCCTGAGGGCTGGCCAGGATGTTGTCGCCATc                                                                   >  2:3033730/1‑93 (MQ=255)
                       aTCGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCa                                                                  >  2:1261608/1‑93 (MQ=255)
                        taggACACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCaa                                                                 <  1:743000/91‑1 (MQ=255)
                        tCGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCaa                                                                 >  1:6162310/1‑93 (MQ=255)
                        tCGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCaa                                                                 >  1:2276288/1‑93 (MQ=255)
                         cGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAg                                                                <  1:6622530/93‑1 (MQ=255)
                           gAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCa                                                              <  1:1791841/93‑1 (MQ=255)
                            aaaCGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAg                                                             <  2:850435/93‑1 (MQ=255)
                             aaCGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGa                                                            >  2:6105224/1‑93 (MQ=255)
                              aCGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGaa                                                           >  1:1189448/1‑93 (MQ=255)
                              aCGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGaa                                                           >  1:1963368/1‑93 (MQ=255)
                                ggCGTCGGTTTCTGGAACACTATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAAtt                                                         <  2:6255746/93‑1 (MQ=255)
                                ggCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAAtt                                                         <  2:4452928/93‑1 (MQ=255)
                                ggCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAAtt                                                         <  1:2227447/93‑1 (MQ=255)
                                ggCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAAtt                                                         <  2:2909414/93‑1 (MQ=255)
                                 gCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAAttt                                                        <  1:5503266/93‑1 (MQ=255)
                                     cGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCAcc                                                    >  2:3077617/1‑93 (MQ=255)
                                     cGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCAcc                                                    >  1:4984132/1‑93 (MQ=255)
                                     cGGTTTCTGGAACACCATGCCCAATTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCATAATTTCAcc                                                    <  2:6723718/93‑1 (MQ=255)
                                      ggTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCt                                                   <  2:1201381/93‑1 (MQ=255)
                                      ggTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCt                                                   <  1:2505688/93‑1 (MQ=255)
                                       gTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCtt                                                  <  2:5499095/93‑1 (MQ=255)
                                       gTTTCTAGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCtt                                                  >  2:3003112/1‑93 (MQ=255)
                                        tttCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTc                                                 >  2:5609818/1‑93 (MQ=255)
                                        tttCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTc                                                 >  1:1019031/1‑93 (MQ=255)
                                         ttCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTcc                                                <  1:4150359/93‑1 (MQ=255)
                                              gAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACg                                           <  1:4673990/93‑1 (MQ=255)
                                               aaCACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACgc                                          <  1:5022533/93‑1 (MQ=255)
                                                acacCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACgct                                         >  2:2437737/1‑93 (MQ=255)
                                                  acCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACgctgc                                       <  1:551436/93‑1 (MQ=255)
                                                  acCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACgctgc                                       >  2:2311529/1‑93 (MQ=255)
                                                   ccATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACgctgct                                      >  1:6294009/1‑93 (MQ=255)
                                                    cATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTc                                     >  1:4365659/1‑93 (MQ=255)
                                                    cATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTc                                     <  1:5526907/93‑1 (MQ=255)
                                                     aTGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTcc                                    <  1:1736245/93‑1 (MQ=255)
                                                        cccACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCggg                                 <  2:23504/93‑1 (MQ=255)
                                                        cccACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCggg                                 <  2:5919451/93‑1 (MQ=255)
                                                        cccACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCggg                                 <  1:886165/93‑1 (MQ=255)
                                                          cACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTa                               <  1:601536/93‑1 (MQ=255)
                                                          cACCTTCGCACGCAGCATTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGACATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTa                               <  1:4745854/93‑1 (MQ=255)
                                                           ccTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTAc                              <  2:1781490/92‑1 (MQ=255)
                                                           aCTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTAc                              >  2:6565411/1‑93 (MQ=255)
                                                             tttCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAg                            <  2:4597193/93‑1 (MQ=255)
                                                                  cacgcaGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCaa                       >  1:1304196/1‑93 (MQ=255)
                                                                   acgcaGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCaaa                      <  2:2138561/93‑1 (MQ=255)
                                                                    cgcaGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCAAAc                     <  1:6713550/93‑1 (MQ=255)
                                                                     gcagcaGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCAAACa                    <  1:1218158/93‑1 (MQ=255)
                                                                      cagcagTGCGATATCCTGAGGGTTGGTCAGGATTTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCAAACAt                   <  2:1595458/93‑1 (MQ=255)
                                                                       agcagTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCAAACAtt                  >  1:5429168/1‑93 (MQ=255)
                                                                        gcagTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCAAACAttt                 >  1:357506/1‑93 (MQ=255)
                                                                          agTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCAAACATTttg               <  2:4376296/93‑1 (MQ=255)
                                                                            tGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCAAACATTttgtt             >  1:339708/1‑93 (MQ=255)
                                                                            tGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCAAACATTttgtt             <  1:1580014/93‑1 (MQ=255)
                                                                            tGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCAAACATTttgtt             >  2:916085/1‑93 (MQ=255)
                                                                             gCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCAAACATTttgttg            >  2:6319752/1‑93 (MQ=255)
                                                                               gATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCAAACATTTTGTTGaa          <  1:2041138/93‑1 (MQ=255)
                                                                                atatCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCAAACATTTTGTTGAAg         >  1:4596587/1‑93 (MQ=255)
                                                                                 tatCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCAAACATTTTGTTGAAgg        >  2:518968/1‑93 (MQ=255)
                                                                                 tatCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCAAACATTTTGTTGAAgg        >  1:3665311/1‑93 (MQ=255)
                                                                                 tatCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCAAACATTTTGTTGAAgg        <  2:1723341/93‑1 (MQ=255)
                                                                                 tatCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCAAACATTTTGTTGAAgg        >  2:2809447/1‑93 (MQ=255)
                                                                                  atCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCAAACATTTTGTTGAAGGt       <  2:6045350/93‑1 (MQ=255)
                                                                                  atCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCAAACATTTTGTTGAAGGt       <  1:5249611/93‑1 (MQ=255)
                                                                                    ccTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCAAACATTTTGTTGAAGGTAc     <  1:6504096/93‑1 (MQ=255)
                                                                                     cTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCAAACATTTTGTTGAAGGTACg    >  1:1437232/1‑93 (MQ=255)
                                                                                      tGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCAAACATTTTGTTGAAGGTACgc   <  1:4806288/93‑1 (MQ=255)
                                                                                      tGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCAAACATTTTGTTGAAGGTACgc   >  1:4285667/1‑93 (MQ=255)
                                                                                       gaggGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCAAACATTTTGTTGAAGGTACgca  <  1:5386776/93‑1 (MQ=255)
                                                                                          |                                                                                         
AAGCGATGTTGTCGTAGATGGACATCGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGAGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCAAACATTTTGTTGAAGGTACGCA  >  NC_000913/3905658‑3905837

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