Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,650,396 T→C intergenic (‑89/‑118) b2520 ← / → b2521 hypothetical protein/putative thiosulfate sulfur transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,650,3960TC87.2% 624.3 / 72.4 227intergenic (‑89/‑118)b2520/b2521hypothetical protein/putative thiosulfate sulfur transferase
Reads supporting (aligned to +/- strand):  ref base T (8/21);  new base C (107/91);  total (115/112)
Fisher's exact test for biased strand distribution p-value = 9.38e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.47e-01

Only 100 of 233 total aligned reads displayed.

ATTTTTCATCCCTGCCACAATGGCCCGTTAGCAACGTCGAATAATTATTATGCGTGAGAATTTACGTTATGTAATTCAGTATCACGCCCATACAATCCGAATTTTAGTGAGAATTGTGGCTCAGACAGCCAAATCAGCCTCTTTTGCTTGTCGCTTTTCACTAAACATCCGACAATTTAGCC  >  NC_000913/2650306‑2650487
                                                                                          |                                                                                           
atttttcatCCCTGCCACAATGGTCCGTTAGCAACGTCGAATAATTATTATGCGTGAGAATTTACGTTATGTAATTCAGTATCACGCCCacac                                                                                           >  2:5243047/1‑93 (MQ=255)
atttttcatCCCTGCCACAATGGCCCGTTAGCAACGTCGAATAATTATTATGCGTGAGAATTTACGTTATGTAATTCAGTATCACGCCCacac                                                                                           >  2:2027021/1‑93 (MQ=255)
        tCCCTGCCACAATGGCCCGTTAGCAACGTCGAATAATTATTATGCGTGAGAATTTACGTTATGTAATTCAGTATCACGCCCACACAATCCGaa                                                                                   >  2:4138946/1‑93 (MQ=255)
         cccTGCCACAATGGCCCGTTAGCAACGTCGAATAATTATTATGGGTGAAAATTTACGTTATGTAATTCAGTATCACGCCCACACAATCCGAAt                                                                                  <  1:3907738/93‑1 (MQ=255)
           cTGCCACAATGGCCCGTTAGCAACGTCGAATAATTATTATGCGTGAGAATTTACGTTATGTAATTCAGTATCACGCCCACACAATCCGAAttt                                                                                >  2:146899/1‑93 (MQ=255)
            tGCCACAATGGCCCGTTAGCAACGTCGAATAATTATTATGCGTGAGAATTTACGTTATGTAATTCAGTATCACGCCCACACAATCCGAAtttt                                                                               >  2:900307/1‑93 (MQ=255)
             gCCACAATGGCCCGTTAGCAACGTCGAATAATTATTATGCGTGAGAATTTACGTTATGTAATTCAGTATCACGCCCACACAATCCGAATTTTa                                                                              >  1:5182829/1‑93 (MQ=255)
               cacaATGGCCCGTTAGCAACGTCGAATAATTATTATGCGTGAGAATTTACGTTATGTAATTCAGTATCACACCCACACAATCAGAATTTTAGt                                                                            >  2:3433303/1‑93 (MQ=255)
                acaATGGCCCGTTAGCAACGTCGAATAATTATTATGCGTGAGAATTTACGTTATGTAATTCAGTATCACGCCCACACAATCCGAATTTTAGTg                                                                           <  1:1485470/93‑1 (MQ=255)
                acaATGGCCCGTTAGCAACGTCGAATAATTATTATGCGTGAGAATTTACGTTATGTAATTCAGTATCACGCCCACACAATCCGAATTTTAGTg                                                                           >  2:506032/1‑93 (MQ=255)
                 caATGGCCCGTTAGCAACGTCGAATAATTATTATGCGTGAGAATTTACGTTATGTAATTCAGTATCACGCCCATACAATCCGAATTTTAGTga                                                                          <  1:872787/93‑1 (MQ=255)
                  aaTGGCCCGTTAGCAACGTCGAATAATTATTATGCGTGAGAATTTACGTTATGTAATTCAGTATCATGCCCACACAATCCGAATTTTAGTgag                                                                         >  1:1171712/1‑93 (MQ=255)
                  aaTGGCCCGTTAGCAACGTCGAATAATTATTATGCGTGAGAATTTACGTTATGTAATTCAGTATCACGCCCATACAATCCGAATTTTAGTgag                                                                         >  1:1011293/1‑93 (MQ=255)
                  aaTGGCCCGTTAGCAACGTCGAATAATTATTATGCGTGAGAATTTACGTTATGTAATTCAGTATCACGCCCACACAATCCGAATTTTAGTgag                                                                         >  1:6043292/1‑93 (MQ=255)
                  aaTGGCCCGTTAGCAACGTCGAATAATTATTATGCGTGAGAATTTACGTTATGTAATTCAGTATCACGCCCACACAATCCGAATTTTAGTgag                                                                         <  2:3496729/93‑1 (MQ=255)
                  aaTGGCCCGTTAGCAACGTCGAATAATTATTATGCGTGAGAATTTACGTTATGTAATTCAGTATCACGCCCACACAATCCGAATTTTAGTgag                                                                         >  1:3543202/1‑93 (MQ=255)
                  aaTGGCCCGTTAGCAACGTCGAATAATTATTATGCGTGAGAATTTACGTTATGTAATTCAGTATCACGCCCACACAATCCGAATTTTAGTgag                                                                         >  2:5913243/1‑93 (MQ=255)
                  aaTGGCCCGTTAGCAACGTCGAATAATTATTATGCGTGAGAATTTACGTTATGTAATTCAGTATCACGCCCACACAATCCGAATTTTAGTgag                                                                         >  2:5990679/1‑93 (MQ=255)
                   aTGGCCCGTTAGCAACGTCGAATAATTATTATGCGTGAGAATTTACGTTATGTAATTCAGTATCACGCCCACACAATCCGAATTTTAGTgaga                                                                        >  1:301927/1‑93 (MQ=255)
                   aTGGCCCGTTAGCAACGTCGAATAATTATTATGCGTGAGAATTTACGTTATGTAATTCAGTATCACGCCCACACAATCCGAATTTTAGTgaga                                                                        >  1:3360394/1‑93 (MQ=255)
                    tGGCCCGTTAGCAACGTCGAATAATTATTATGCGTGAGAATTTACGTTATGTAATTCAGTATCACGCCCACACAATCCGAATTTTAGTGAGaa                                                                       >  1:5268214/1‑93 (MQ=255)
                     ggCCCGTTAGCAACGTCGAATAATTATTATGCGTGAGAATTTACGTTATGTAATTCAGTATCACGCCCACACAATCCGAATTTTAGTGAGAAt                                                                      <  1:185886/93‑1 (MQ=255)
                     ggCCCGTTAGCAACGTCGAATAATTATTATGCGTGAGAATTTACGTTATGTAATTCAGTATCACGCCCACACAATCCGAATTTTAGTGAGAAt                                                                      >  1:1958885/1‑93 (MQ=255)
                     ggCCCGTTAGCAACGTCGAATAATTATTATGCGTGAGAATTTACGTTATGTAATTAAGTATCACGCCCACACAATCCGAATTTTAGTGAGAAt                                                                      <  2:4007052/93‑1 (MQ=255)
                       cccGTTAGCAACGTCGAATAATTATTATGCGTGAGAATTTACGTTATGTAATTCAGTATCACGCCCACACAATCCGAATTTTAGTGAGAATtg                                                                    >  1:3673268/1‑93 (MQ=255)
                        ccGTTAGCAACGTCGAATAATTATTATGCGTGAGAATTTACGTTATGTAATTCAGTATCACGCCCACACAATCCGAATTTTAGTGAGAATtgt                                                                   <  2:1665973/93‑1 (MQ=255)
                          gTTAGCAACGTCGAATAATTATTATGCGTGAGAATTTACGTTATGTAATTCAGTATCACGCCCATACAATCCGAATTTTAGTGAGAATTGTgg                                                                 <  1:5256865/93‑1 (MQ=255)
                          gTTAGCAACGTCGAATAATTATTATGCGTGAGAATTTACGTTATGTAATTCAGTATCACGCCCACACAATCCGAATTTTAGTGAGAATTGTgg                                                                 <  2:2617718/93‑1 (MQ=255)
                          gTTAGCAACGTCGAATAATTATTATGCGTGAGAATTTACGTTATGTAATTCAGTATCACGCCCACACAATCCGAATTTTAGTGAGAATTGTgg                                                                 <  1:4726843/93‑1 (MQ=255)
                             aGCAACGTCGAATAATTATTATGCGTGAGAATTTACGTTATGTAATTCAGTATCACGCCCACACAATCCGAATTTTAGTGAGAATTGTGGCTc                                                              <  1:4388950/93‑1 (MQ=255)
                              gCAACGTCGAATAATTATTATGCGTGAGAATTTACGTTATGTAATTCAGTATCACGCCCATACAATCCGAATTTTAGTGAGAATTGTGGCTCa                                                             <  1:2785228/93‑1 (MQ=255)
                              gCAACGTCGAATAATTATTATGCGTGAGAATTTACGTTATGTAATTCAGTATCACGCCCACACAATCCGAATTTTAGTGAGAATTGTGGCTCa                                                             <  2:6225338/93‑1 (MQ=255)
                              gCAACGTCGAATAATTATTATGCGTGAGAATTTACGTTATGTAATTCAGTATCACGCCCACACAATCCGAATTTTAGTGAGAATTGTGGCTCa                                                             >  2:1269516/1‑93 (MQ=255)
                                aaCGTTGAATAATTATTATGCGTGAGAATTTACGTTATGTAATTCAGTATCACGCCCATACACTCCGAATTTTAGTGAGAATTGTGGCTCAGa                                                           >  1:1369815/1‑93 (MQ=255)
                                aaCGTCGAATAATTATTATGCGTGAGAATTTACGTTATGTAATTCAGTATCACGCCCACACAATCCGAATTTTAGTGAGAATTGTGGCTCAGa                                                           <  1:2382957/93‑1 (MQ=255)
                                 aCGTCGAATAATTATTATGCGTGAGAATTTACGTTATGTAATTCAGTATCACGCCCACACAATCCGAATTTTAGTGAGAATTGTGGCTCAGAc                                                          <  2:6545202/93‑1 (MQ=255)
                                    tCGAATAATTATTATGCGTGAGAATTTACGTTATGTAATTCAGTATCACGCCCACACAATCCGAATTTTAGTGAGAATTGTGGCTCAGACAGc                                                       >  2:2203187/1‑93 (MQ=255)
                                    tCGAATAATTATTATGCGTGAGAATTTACGTTATGTAATTCAGTATCACGCCCACACAATCCGAATTTTAGTGAGAATTGTGGCTCAGACAGc                                                       <  1:958136/93‑1 (MQ=255)
                                      gAATAATTATTATGCGTGAGAATTTACGTTATGTAATTCAGTATCACGCCCACACAATCCGAATTTTAGTGAGAATTGTGGCTCAGACAGCCa                                                     >  2:5509031/1‑93 (MQ=255)
                                      gAATAATTATTATGCGTGAGAATTTACGTTATGTAATTCAGTATCACGCCCACACAATCCGAATTTTAGTGAGAATTGTGGCTCAGACAGCCa                                                     >  1:1886672/1‑93 (MQ=255)
                                       aataatTATTATGCGTGAGAATTTACGTTATGTAATTCAGTATCACGCCCACACAATCCGAATTTTAGTGAGAATTGTGGCTCAGACAGCCaa                                                    >  2:4840783/1‑93 (MQ=255)
                                        ataatTATTATGCGTGAGAATTTACGTTATGTAATTCAGTATCACGCCCACACAATCCGAATTTTAGTGAGAATTGTGGCTCAGACAGCCaaa                                                   >  2:5929002/1‑93 (MQ=255)
                                        ataatTATTATGCGTGAGAATTTACGTTATGTAATTCAGTATCACGCCCACACAATCCGAATTTTAGTGAGAATTGTGGCTCAGACAGCCaaa                                                   >  1:1109405/1‑93 (MQ=255)
                                         taatTATTATGCGTGAGAATTTACGTTATGTAATTCAGTATCACGCCCACACAATCCGAATTTTAGTGAGAATTGTGGCTCAGACAGCCAAAt                                                  >  2:2275496/1‑93 (MQ=255)
                                         taatTATTATGCGTGAGAATTTACGTTATGTAATTCAGTATCACGCCCACACAATCCGAATTTTAGTGAGAATTGTGGCTCAGACAGCCAAAt                                                  <  1:211384/93‑1 (MQ=255)
                                          aatTATTATGCGTGAGAATTTACGTTATGTAATTCAGTATCACGCCCATACAATCCGAATTTTAGTGAGAATTGTGGCTCAGACAGCCAAATc                                                 <  1:3799064/93‑1 (MQ=255)
                                          aatTATTATGCGTGAGAATTTACGTTATGTAATTCAGTATCACGCCCACACAATCCGAATTTTAGTGAGAATTGTGGCTCAGACAGCCAAATc                                                 <  2:4675426/93‑1 (MQ=255)
                                          aatTATTATGCGTGAGAATTTACGTTATGTAATTCAGTATCACGCCCACACAATCCGAATTTTAGTGAGAATTGTGGCTCAGACAGCCAAATc                                                 <  2:4436216/93‑1 (MQ=255)
                                            ttattaTGCGTGAGAATTTACGTTATGTAATTCAGTATCACGCCCACACAATCCGAATTTTAGTGAGAATGTGGCTCAGACAGCCAAATCAGc                                              >  2:3750312/1‑93 (MQ=255)
                                              attatGCGTGAGAATTTACGTTATGTAATTCAGTATCACGCCCACACAATCCGAATTTTAGTGAGAATTGTGGCTCAGACAGCCAAATCAGcc                                             <  1:3245418/93‑1 (MQ=255)
                                               ttatGCGTGAGAATTTACGTTATGTAATTCAGTATCACGCCCACACAATCCGAATTTTAGTGAGAATTGTGGCTCAGACAGCCAAATCAGCct                                            >  1:2529529/1‑93 (MQ=255)
                                               ttatGCGTGAGAATTTACGTTATGTAATTCAGTATCACGCCCACACAATCCGAATTTTAGTGAGAATTGTGGCTCAGACAGCCAAATCAGCct                                            >  1:4352173/1‑93 (MQ=255)
                                               ttatGCGTGAGAATTTACGTTATGTAATTCAGTATCACGCCCACACAATCCGAATTTTAGTGAGAATTGTGGCTCAGACAGCCAAATCAGCct                                            >  2:502113/1‑93 (MQ=255)
                                                tatGCGTGAGAATTTACGTTATGTAATTCAGTATCACGCCCACACAATCCGAATTTTAGTGAGAATTGTGGCTCAGACAGCCAAATCAGCctc                                           >  2:5526713/1‑93 (MQ=255)
                                                   gCGTGAGAATTTACGTTATGTAATTCAGTATCACGCCCATACAATCCGAATTTTAGTGAGAATTGTGGCTCAGACAGCCAAATCAGCCTCttt                                        <  2:2791744/93‑1 (MQ=255)
                                                   gCGTGAGAATTTACGTTATGTAATTCAGTATCACGCCCACACAATCCGAATTTTAGTGAGAATTGTGGCTCAGACAGCCAAATCAGCCTCttt                                        >  1:4033002/1‑93 (MQ=255)
                                                    cGTGAGAATTTACGTTATGTAATTCAGTATCACGCCCACACAATCCGAATTTTAGTGAGAATTGTGGCTCAGACAGCCAAATCAGCCTCtttt                                       <  2:5110234/93‑1 (MQ=255)
                                                    cGTGAGAATTTACGTTATGTAATTCAGTATCACGCCCACACAATCCGAATTTTAGTGAGAATTGTGGCTCAGACAGCCAAATCAGCCTCtttt                                       <  1:915269/93‑1 (MQ=255)
                                                      tGAGAATTTACGTTATGTAATTCAGTATCACGCCCACACAATCCGAATTTTAGTGAGAATTGTGGCTCAGACAGCCAAATCAGCCTCTTTTGc                                     <  2:357353/93‑1 (MQ=255)
                                                       gagaATTTACGTTATGTAATTCAGTATCACGCCCATACAATCCGAATTTTAGTGAGAATTGTGGCTCAGACAGCCAAATCAGCCTCTTTTGCt                                    <  1:4770671/93‑1 (MQ=255)
                                                       gagaATTTACGTTATGTAATTCAGTATCACGCCCACACAATCCGAATTTTAGTGAGAATTGTGGCTCAGACAGCCAAATCAGCCTCTTTTGCt                                    >  2:6059622/1‑93 (MQ=255)
                                                       gagaATTTACGTTATGTAATTCAGTATCACGCCCACACAATCCGAATTTTAGTGAGAATTGTGGCTCAGACAGCCAAATCAGCCTCTTTTGCt                                    >  2:3299516/1‑93 (MQ=255)
                                                       gagaATTTACGTTATGTAATTCAGTATCACGCCCACACAATCCGAATTTTAGTGAGAATTGTGGCTCAGACAGCCAAATCAGCCTCTTTTGCt                                    <  1:160395/93‑1 (MQ=255)
                                                       gagaATTTACGTTATGTAATTCAGTATCACGCCCACACAATCCGAATTTTAGTGAGAATTGTGGCTCAGACAGCCAAATCAGCCTCTTTTGCt                                    <  1:6220946/93‑1 (MQ=255)
                                                        agaATTTACGTTATGTAATTCAGTATCACGCCCACACAATCCGAATTTTAGTGAGAATTGTGGCTCAGACAGCCAAATCAGCCTCTTTTGCtt                                   <  1:2044727/93‑1 (MQ=255)
                                                           aTTTACGTTATGTAATTCAGTATCACGCCCACACAATCCGAATTTTAGTGAGAATTGTGGCTCAGACAGCCAAATCAGCCTCTTTTGCTTGTc                                >  2:5223456/1‑93 (MQ=255)
                                                            tttACGTTATGTAATTCAGTATCACGCCCACACAATCCGAATTTTAGTGAGAATTGTGGCTCAGACAGCCAAATCAGCCTCTTTTGCTTGTCg                               >  2:3156099/1‑93 (MQ=255)
                                                              tACGTTATGTAATTCAGTATCACGCCCACACAATCCGAATTTTAGTGAGAATTGTGGCTCAGACAGCCAAATCAGCCTCTTTTGCTTGTCGCt                             >  1:5774022/1‑93 (MQ=255)
                                                               aCGTTATGTAATTCAGTATCACGCCCACACAATCCGAATTTTAGTGAGAATTGTGGCTCAGACAGCCAAATCAGCCTCTTTTGCTTGTCGCtt                            <  2:1887604/93‑1 (MQ=255)
                                                                cGTTATGTAATTCAGTATCACGCCCACACAATCCGAATTTTAGTGAGAATTGTGGCTCAGACAGCCAAATCAGCCTCTTTTGCTTGTCGCttt                           <  1:3078476/93‑1 (MQ=255)
                                                                 gTTATGTAATTCAGTATCACGCCCACACAATCCGAATTTTAGTGAGAATTGTGGCTCAGACAGCCAAATCAGCCTCTTTTGCTTGTCGCtttt                          >  2:4560749/1‑93 (MQ=255)
                                                                 gTTATGTAATTCAGTATCACGCCCACACAATCCGAATTTTAGTGAGAATTGTGGCTCAGACAGCCAAATCAGCCTCTTTTGCTTGTCGCtttt                          >  2:1538962/1‑93 (MQ=255)
                                                                 gTTATGTAATTCAGTATCACGCCCACACAATCCGAATTTTAGTGAGAATTGTGGCTCAGACAGCCAAATCAGCCTCTTTTGCTTGTCGCtttt                          >  1:3128666/1‑93 (MQ=255)
                                                                  ttATGTAATTCAGTATCACGCCCACACCATCCGAATTTTAGTGAGCATTGTGGCTCAGACAGCCAAATCAGCCTCTTTTGCTTGTCGCTTTTc                         >  2:2908845/1‑93 (MQ=255)
                                                                  ttATGTAATTCAGTATCACGCCCACACAATCCGAATTTTAGTGAGAATTGTGGCTCAGACAGCCAAATCAGCCTCTTTTGCTTGTCGCTTTTc                         <  1:3632270/93‑1 (MQ=255)
                                                                   tATGTAATTCAGTATCACGCCCACACAATCCGAATTTTAGTGAGAATTGTGGCTCAGACAGCCAAATCAGCCTCTTTTGCTTGTCGCTTTTCa                        >  2:2240655/1‑93 (MQ=255)
                                                                    aTGTAATTCAGTATCACGCCCATACAATCCGAATTTTAGTGAGAATTGTGGCTCAGACAGCCAAATCAGCCTCTTTTGCTTGTCGCTTTTCAc                       <  2:4237661/93‑1 (MQ=255)
                                                                    aTGTAATTCAGTATCACGCCCACACAATCCGAATTTTAGTGAGAATTGTGGCTCAGACAGCCAAATCAGCCTCTTTTGCTTGTCGCTTTTCAc                       <  2:2099903/93‑1 (MQ=255)
                                                                     tGTAATTCAGTATCACGCCCACACAATCCGAATTTTAGTGAGAATTGTGGCTCAGACAGCCAAATCAGCCTCTTTTGCTTGTCGCTTTTCACt                      <  1:4680641/93‑1 (MQ=255)
                                                                     tGTAATTCAGTATCACGCCCACACAATCCGAATTTTAGTGAGAATTGTGGCTCAGACAGCCAAATCAGCCTCTTTTGCTTGTCGCTTTTCACt                      >  2:5513876/1‑93 (MQ=255)
                                                                        aaTTCAGTATCACGCCCACACAATCCGAATTTTAGTGAGAATTGTGGCTCAGACAGCCAAATCAGCCTCTTTTGCTTGTCGCTTTTCACTaaa                   >  2:1609901/1‑93 (MQ=255)
                                                                         aTTCAGTATCACGCCCACACAATCCGAATTTTAGTGAGAATTGTGGCTCAGACAGCCAAATCAGCCTCTTTTGCTTGTCGCTTTTCACTAAAc                  >  2:1409895/1‑93 (MQ=255)
                                                                          ttCAGTATCACGCCCACACAATCCGAATTTTAGTGAGAATTGTGGCTCAGACAGCCAAATCAGCCTCTTTTGCTTGTCGCTTTTCACTAAACa                 >  2:5723231/1‑93 (MQ=255)
                                                                          ttCAGTATCACGCCCACACAATCCGAATTTTAGTGAGAATTGTGGCTCAGACAGCCAAATCAGCCTCTTTTGCTTGTCGCTTTTCACTAAACa                 <  1:4104213/93‑1 (MQ=255)
                                                                           tCAGTATCACGCCCATACAATCCGAATTTTAGTGAGAATTGTGGCTCAGACAGCCAAATCAGCCTCTTTTGCTTGTCGCTTTTCACTAAACAt                >  2:4635789/1‑93 (MQ=255)
                                                                            cAGTATCACGCCCACACAATCCGAATTTTAGTGAGAATTGTGGCTCAGACAGCCAAATCAGCCTCTTTTGCTTGTCGCTTTTCACTAAACATc               <  1:1783521/93‑1 (MQ=255)
                                                                            cAGTATCACGCCCACACAATCCGAATTTTAGTGAGAATTGTGGCTCAGACAGCCAAATCAGCCTCTTTTGCTTGTCGCTTTTCACTAAACATc               >  1:6320080/1‑93 (MQ=255)
                                                                             aGTATCACGCCCACACAATCCGAATTTTAGTGAGAATTGTGGCTCAGACAGCCAAATCGGCCTCTTTTGCTTGTCGCTTTTCACTAAACATcc              <  2:3119294/93‑1 (MQ=255)
                                                                             aGTATCACGCCCACACAATCCGAATTTTAGTGAGAATTGTGGCTCAGACAGCCAAATCAGCCTCTTTTGCTTGTCGCTTTTCACTAAACATcc              <  1:4289362/93‑1 (MQ=255)
                                                                             aGTATCACGCCCACACAATCCGAATTTTAGTGAGAATTGTGGCTCAGACAGCCAAATCAGCCTCTTTTGCTTGTCGCTTTTCACTAAACATcc              <  1:5411179/93‑1 (MQ=255)
                                                                                aTCACGCCCACACAATCCGAATTTTAGTGAGAATTGTGGCTCAGACAGCCAAATCAGCCTCTTTTGCTTGTCGCTTTTCACTAAACATCCGAc           >  1:2665034/1‑93 (MQ=255)
                                                                                aTCACGCCCACACAATCCGAATTTTAGTGAGAATTGTGGCTCAGACAGCCAAATCAGCCTCTTTTGCTTGTCGCTTTTCACTAAACATCCGAc           >  2:239002/1‑93 (MQ=255)
                                                                                aTCACGCCCACACAATCCGAATTTTAGTGAGAATTGTGGCTCAGACAGCCAAATCAGCCTCTTTTGCTTGTCGCTTTTCACTAAACATCCGAc           <  2:2647757/93‑1 (MQ=255)
                                                                                aTCACGCACACACAATCCGAATTTTAGTGAGAATTGTGGCTCAGACAGCCAAATCAGCCTCTTTTGCTTGTCGCTTTTCACTAAACATCCGAc           >  1:3693417/1‑93 (MQ=255)
                                                                                    cGCCCACACAATCCGAATTTTAGTGAGAATTGTGGCTCAGACAGCCAAATCAGCCTCTTTTGCTTGTCGCTTTTCACTAAACATCCGACAAtt       >  1:5685670/1‑93 (MQ=255)
                                                                                      cccACACAATCCGAATTTTAGTGAGAATTGTGGCTCAGACAGCCAAATCAGCCTCTTTTGCTTGTCGCTTTTCACTAAACATCCGACAATTTa     >  2:3932032/1‑93 (MQ=255)
                                                                                       ccACACAATCCGAATTTTAGTGAGAATTGTGGCTCAGACAGCCAAATCAGCCTCTTTTGCTTGTCGCTTTTCACTAAACATCCTACAATTTaa    >  2:5869521/1‑92 (MQ=255)
                                                                                       ccACACAATCCGAATTTTAGTGAGAATTGTGGCTCAGACAGCCAAATCAGCCTCTTTTGCTTGTCGCTTTTCACTAAACATCCGACAATTTAg    >  2:3956640/1‑93 (MQ=255)
                                                                                       ccACACAATCCGAATTTTAGTGAGAATTGTGGCTCAGACAGCCAAATCAGCCTCTTTTGCTTGTCGCTTTTCACTAAACATCCGACAATTTAg    >  1:4997261/1‑93 (MQ=255)
                                                                                       ccACACAATCCGAATTTTAGTGAGAATTGTGGCTCAGACAGCCAAATCAGCCTCTTTTGCTTGTCGCTTTTCACTAAACATCCGACAATTTAg    <  1:4758573/93‑1 (MQ=255)
                                                                                          |                                                                                           
ATTTTTCATCCCTGCCACAATGGCCCGTTAGCAACGTCGAATAATTATTATGCGTGAGAATTTACGTTATGTAATTCAGTATCACGCCCATACAATCCGAATTTTAGTGAGAATTGTGGCTCAGACAGCCAAATCAGCCTCTTTTGCTTGTCGCTTTTCACTAAACATCCGACAATTTAGCC  >  NC_000913/2650306‑2650487

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: