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,7480AG98.0% 637.7 / ‑3.0 203S83P (TCT→CCT) pstBATP‑binding component of high‑affinity phosphate‑specific transport system
Reads supporting (aligned to +/- strand):  ref base A (3/1);  new base G (110/89);  total (113/90)
Fisher's exact test for biased strand distribution p-value = 6.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 215 total aligned reads displayed.

AAGCGATGTTGTCGTAGATGGACATCGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGAGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCAAACATTTTGTTGAAGGTACGCAGCA  >  NC_000913/3905658‑3905840
                                                                                          |                                                                                            
aaGCGATGTTGTCGTAGATGGACATCGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGt                                                                                            >  2:2105436/1‑93 (MQ=255)
aaGCGATGTTGTCGTAGATGGACATCGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGt                                                                                            <  2:1188285/93‑1 (MQ=255)
aaGCGATGTTGTCGTAGATGGACATCGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGt                                                                                            >  1:2092486/1‑93 (MQ=255)
 aGCGATGTTGTCGTAGATGGACATCGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGtt                                                                                           <  2:4259231/93‑1 (MQ=255)
 aGCGATGTTGTCGTAGATGGACATCGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGtt                                                                                           >  2:3347223/1‑93 (MQ=255)
 aGCGATGTTGTCGTAGATGGACATCGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGtt                                                                                           >  2:2383865/1‑93 (MQ=255)
 aGCGATGTTGTCGTAGATGGACATCGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGtt                                                                                           <  1:4463276/93‑1 (MQ=255)
 aGCGATGTTGTCGTAGATGGACATCGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGtt                                                                                           <  1:3788303/93‑1 (MQ=255)
 aGCGATGTTGTCGTAGATGGACATCGGAAACGGCGTCGGTTTCTGGAAAACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGtt                                                                                           >  2:3975703/1‑93 (MQ=255)
  gCGATGTTGTCGTAGATGGACATCGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTg                                                                                          >  1:1647395/1‑93 (MQ=255)
  gCGATGTTGTCGTAGATGGACATCGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTg                                                                                          >  1:331272/1‑93 (MQ=255)
   cGATGTTGTCGTAGATGGACATCGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTgg                                                                                         <  2:3070101/93‑1 (MQ=255)
   cGATGTTGTCGTAGATGGACATCGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTgg                                                                                         >  1:1377690/1‑93 (MQ=255)
    gATGTTGTCGTAGATGGACATCGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGt                                                                                        >  2:889676/1‑93 (MQ=255)
    gATGTTGTCGTAGATGGACATCGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGt                                                                                        >  1:4119205/1‑93 (MQ=255)
     aTGTTGTCGTAGATGGACATCGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTc                                                                                       <  2:4176895/93‑1 (MQ=255)
      tgttgtCGTAGATGGACATCGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCa                                                                                      >  2:1604778/1‑93 (MQ=255)
       gttgtCGTAGATGGACATCGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAg                                                                                     >  2:4326290/1‑93 (MQ=255)
       gttgtCGTAGATGGACATCGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAg                                                                                     <  1:3833798/93‑1 (MQ=255)
          gtCGTAGATGGACATCGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGAt                                                                                  <  1:3719715/93‑1 (MQ=255)
           tCGTAGATGGACATCGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGAtg                                                                                 <  1:3620830/93‑1 (MQ=255)
            cGTAGATGGACATCGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGAtgt                                                                                <  2:3564340/93‑1 (MQ=255)
             gTAGATGGACATCGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGAtgtt                                                                               >  2:858712/1‑93 (MQ=255)
             gTAGATGGACATCGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGAtgtt                                                                               <  1:4478693/93‑1 (MQ=255)
              tAGATGGACATCGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGAtgttg                                                                              >  1:3909377/1‑93 (MQ=255)
               aGATGGACATCGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGAtgttgt                                                                             >  1:1979581/1‑93 (MQ=255)
                  tGGACATCGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGc                                                                          >  1:106515/1‑93 (MQ=255)
                     aCATCGGAAACGGCGTCGGTTTCTGGAACTCCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCAt                                                                       <  1:2204143/93‑1 (MQ=255)
                        tCGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCaa                                                                    >  1:1304917/1‑93 (MQ=255)
                        tCGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCaa                                                                    >  2:1950707/1‑93 (MQ=255)
                          ggAAACGGCGTCGGTTTCTGTAACACCATGCCCACGTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGc                                                                  <  1:533623/93‑1 (MQ=255)
                          ggAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAga                                                                  >  1:2534958/1‑92 (MQ=255)
                          ggAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGc                                                                  <  2:5643/93‑1 (MQ=255)
                           gAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCCTCAAGCa                                                                 >  1:2151670/1‑93 (MQ=255)
                           gAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCa                                                                 <  2:2917809/93‑1 (MQ=255)
                           gAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCa                                                                 >  2:2319415/1‑93 (MQ=255)
                           gAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCa                                                                 >  2:2513712/1‑93 (MQ=255)
                           gAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCa                                                                 >  2:2822701/1‑93 (MQ=255)
                            aaaCGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAg                                                                >  1:3233014/1‑93 (MQ=255)
                            aaaCGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAg                                                                <  1:1106805/93‑1 (MQ=255)
                             aaCGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGa                                                               >  2:1638244/1‑93 (MQ=255)
                              aCGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGaa                                                              >  2:4520860/1‑93 (MQ=255)
                                 gCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAAttt                                                           >  1:459310/1‑93 (MQ=255)
                                  cGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTc                                                          <  1:595361/93‑1 (MQ=255)
                                  cGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTc                                                          <  2:167106/93‑1 (MQ=255)
                                   gTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCa                                                         <  1:2695787/93‑1 (MQ=255)
                                   gTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCa                                                         >  1:781863/1‑93 (MQ=255)
                                      ggTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCt                                                      >  2:3873728/1‑93 (MQ=255)
                                      ggTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCt                                                      >  1:2666277/1‑93 (MQ=255)
                                         ttCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTcc                                                   <  2:2250323/93‑1 (MQ=255)
                                         ttCTGGAACACCATGCCCACTGTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACGTTcc                                                   <  2:4319604/93‑1 (MQ=255)
                                           cTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGc                                                 >  2:708015/1‑93 (MQ=255)
                                            tgtaaCACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCa                                                >  2:4134043/4‑93 (MQ=255)
                                             ggAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCAc                                               >  2:4491530/1‑93 (MQ=255)
                                             ggAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCAc                                               >  2:3672755/1‑93 (MQ=255)
                                             ggAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCAc                                               >  2:2662805/1‑93 (MQ=255)
                                               aaCACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACgc                                             <  1:2877027/93‑1 (MQ=255)
                                               aaCACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACgc                                             <  2:258989/93‑1 (MQ=255)
                                                acacCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACgct                                            <  2:3507304/93‑1 (MQ=255)
                                                acacCATGCCAACTTACGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACgct                                            <  1:1081206/93‑1 (MQ=255)
                                                  acCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACgctgc                                          <  1:914009/93‑1 (MQ=255)
                                                     aTGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTcc                                       >  2:4576968/1‑93 (MQ=255)
                                                     aTGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTcc                                       >  1:1519372/1‑93 (MQ=255)
                                                     aTGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTcc                                       >  1:1705002/1‑93 (MQ=255)
                                                     aTGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTcc                                       >  2:4424824/1‑93 (MQ=255)
                                                      tGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCg                                      >  2:2337922/1‑93 (MQ=255)
                                                       gCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCgg                                     >  2:49946/1‑93 (MQ=255)
                                                        cccACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCggg                                    <  2:3797926/93‑1 (MQ=255)
                                                        cccACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCggg                                    >  2:3288802/1‑93 (MQ=255)
                                                         ccACTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGt                                   >  1:470429/1‑93 (MQ=255)
                                                           aCTTTCGCACGCAGCAGTGCGATATCCTGAGAGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTAc                                 >  2:1585734/1‑93 (MQ=255)
                                                            cTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACa                                <  2:1138986/93‑1 (MQ=255)
                                                            cTTTCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACa                                >  2:1145482/1‑93 (MQ=255)
                                                             tttCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAg                               <  2:1590771/93‑1 (MQ=255)
                                                              ttCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGt                              <  1:688089/93‑1 (MQ=255)
                                                              ttCGCACGCAGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGt                              <  1:4528553/93‑1 (MQ=255)
                                                              ttCGCACGCAGAAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGt                              <  2:1330603/93‑1 (MQ=255)
                                                                  cacgcaGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCaa                          >  1:176686/1‑93 (MQ=255)
                                                                   acgcaGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCaaa                         >  2:275100/1‑93 (MQ=255)
                                                                   acgcaGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCaaa                         >  1:4477643/1‑93 (MQ=255)
                                                                    cgcaGCAGTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCAAAc                        <  1:3232011/93‑1 (MQ=255)
                                                                       agcagTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCAAACAtt                     >  2:4048424/1‑93 (MQ=255)
                                                                         cagTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCAAACAtttt                   <  2:1570000/93‑1 (MQ=255)
                                                                          agTGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCAAACATTttg                  >  2:4693152/1‑93 (MQ=255)
                                                                            tGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCAAACATTttgtt                <  1:1936231/93‑1 (MQ=255)
                                                                            tGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCAAACATTttgtt                >  2:672595/1‑93 (MQ=255)
                                                                            tGCGATATCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCAACTTCCGCACGCTGCTCCGGGTACAGTTCAAACATTttgtt                >  2:1124729/1‑93 (MQ=255)
                                                                               gATATCCTGAGGGTTGGTCAGGATGTTGTCGCCGTCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCAAACATTTTGTTGaa             <  1:747606/93‑1 (MQ=255)
                                                                               gATATCCTGAGGGTGGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCAAACATTTTGTTGaa             <  2:2219967/93‑1 (MQ=255)
                                                                                atatCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCAAACATTTTGTTGAAg            >  1:3365278/1‑93 (MQ=255)
                                                                                 tatCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCAAACATTTTGTTGAAgg           >  2:149396/1‑93 (MQ=255)
                                                                                  atCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCAAACATTTTGTTGAAGGt          >  1:1155475/1‑93 (MQ=255)
                                                                                  atCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCAAACATTTTGTTGAAGGt          <  1:4403199/93‑1 (MQ=255)
                                                                                   tCCTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCAAACATTTTGTTGAAGGTa         <  1:4032307/93‑1 (MQ=255)
                                                                                     cTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCAAACATTTTGTTGAAGGTACg       >  1:260667/1‑93 (MQ=255)
                                                                                     cTGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCAAACATTTTGTTGAAGGTACg       >  1:2758693/1‑93 (MQ=255)
                                                                                      tGAGGGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCAAACATTTTGTTGAAGGTACgc      <  2:599874/93‑1 (MQ=255)
                                                                                       gaggGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCAAACATTTTGTTGAAGGTACgca     >  1:2443668/1‑93 (MQ=255)
                                                                                          agTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCAAACATTTTGTTGAAGGTACgcagca  <  1:3128053/93‑1 (MQ=255)
                                                                                          agTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCAAACATTTTGTTGAAGGTACgcagca  >  2:1578363/1‑93 (MQ=255)
                                                                                          |                                                                                            
AAGCGATGTTGTCGTAGATGGACATCGGAAACGGCGTCGGTTTCTGGAACACCATGCCCACTTTCGCACGCAGCAGTGCGATATCCTGAGAGTTGGTCAGGATGTTGTCGCCATCAAGCAGAATTTCACCTTCCGCACGCTGCTCCGGGTACAGTTCAAACATTTTGTTGAAGGTACGCAGCA  >  NC_000913/3905658‑3905840

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: ‑