Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,465,657 T→C Y365H (TAT→CAT)  intC → putative prophage Sf6‑like integrase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,465,6570TC99.6% 743.2 / NA 227Y365H (TAT→CAT) intCputative prophage Sf6‑like integrase
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (116/110);  total (117/110)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 228 total aligned reads displayed.

AACGAGCACGAATGGCCTGCTGACGCTATTGAAGTGCAACTGGCACATGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGTATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCGGACTGGCTTGATGAGAAGGTGGAGTGAGCGACCTTAACAACTATCGAATAG  >  NC_000913/2465567‑2465746
                                                                                          |                                                                                         
aaCGAGCACGAATGGCCTGCTGACGCTATTGAAGTGCAACTGGCACATGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCAt                                                                                         >  1:1653750/1‑93 (MQ=255)
    aGCACGAATGGCCTGCTGACGCTATTGAAGTGCAACTGGCACATGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCg                                                                                     >  1:1133445/1‑93 (MQ=255)
    aGCACGAATGGCCTGCTGACGCTATTGAAGTGCAACTGGCACATGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCg                                                                                     >  2:1812242/1‑93 (MQ=255)
      cACGAATGGCCTGCTGACGCTATTGAAGTGCAACTGGCACATGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGAt                                                                                   >  1:4805916/1‑93 (MQ=255)
        cGAATGGCCTGCTGACGCTATTGAAGTGCAACTGGCACATGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATaa                                                                                 >  2:1579507/1‑93 (MQ=255)
        cGAATGGCCTGCTGACGCTATTGAAGTGCAACTGGCACATGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATaa                                                                                 >  1:3155953/1‑93 (MQ=255)
         gAATGGCCTGCTGACGCTATTGAAGTGCAACTGGCACATGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATaaa                                                                                >  1:2691037/1‑93 (MQ=255)
          aaTGGCCTGCTGACGCTATTGAAGTGCAACTGGCACATGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAAc                                                                               >  1:1483117/1‑93 (MQ=255)
          aaTGGCCTGCTGACGCTATTGAAGTGCAACTGGCACATGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAAc                                                                               <  2:4887604/93‑1 (MQ=255)
           aTGGCCTGCTGACGCTATTGAAGTGCAACTGGCACATGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACg                                                                              >  2:2177508/1‑93 (MQ=255)
             ggCCTGCTGACGCTATTGAAGTGCAACTGGCACATGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCa                                                                            >  2:5622598/1‑93 (MQ=255)
             ggCCTGCTGACGCTATTGAAGTGCAACTGGCACATGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCa                                                                            <  2:879768/93‑1 (MQ=255)
             ggCCTGCTGACGCTATTGAAGTGCAACTGGCACATGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCa                                                                            <  1:2406832/93‑1 (MQ=255)
               ccTGCTGACGCTATTGAAGTGCAACTGGCACATGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCaga                                                                          <  1:917690/93‑1 (MQ=255)
               ccTGCTGACGCTATTGAAGTGCAACTGGCACATGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCaga                                                                          <  1:2112567/93‑1 (MQ=255)
                 tgctgACGCTATTGAAGTGCAACTGGCACATGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAgaga                                                                        >  1:5986101/1‑93 (MQ=255)
                 tgctgACGCTATTGAAGTGCAACTGGCACATGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAgaga                                                                        >  2:5811154/1‑93 (MQ=255)
                 tgctgACGCTATTGAAGTGCAACTGGCACATGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAgaga                                                                        >  2:4742540/1‑93 (MQ=255)
                  gctgACGCTATTGAAGTGCAACTGGCACATGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGaa                                                                       >  2:3293383/1‑93 (MQ=255)
                   ctgACGCTATTGAAGTGCAACTGGCACATGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGaaa                                                                      <  2:172234/93‑1 (MQ=255)
                    tgACGCTATTGAAGTGCAACTGGCACATGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGAAat                                                                     >  2:3487108/1‑93 (MQ=255)
                    tgACGCTATTGAAGTGCAACTGGCACATGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGAAat                                                                     <  2:5413241/93‑1 (MQ=255)
                     gACGCTATTGAAGTGCAACTGGCACATGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGAAatg                                                                    <  2:364382/93‑1 (MQ=255)
                     gACGCTATTGAAGTGCAACTGGCACATGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGAAatg                                                                    <  1:2015625/93‑1 (MQ=255)
                     gACGCTATTGAAGTGCAACTGGCACATGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGAAatg                                                                    >  2:1470296/1‑93 (MQ=255)
                     gACGCTATTGAAGTGCAACTGGCACATGCAAACGGCGGATCTGTGCGTGGGATTTAAAACCATGCTCAGCATCTCGATAAACGCAGAGAAatg                                                                    >  1:3694951/1‑93 (MQ=255)
                      aCGCTATTGAAGTGCAACTGGCACATGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGAAatga                                                                   <  1:2295638/93‑1 (MQ=255)
                       cgttATTTAAGTGCAACTGTCTCATGCAAACGGCGGATCTGTGCGTGGGATTTACACCCATGCTCAGCATCTCGATAAACGCAGAGAAatgat                                                                  <  2:4417729/90‑1 (MQ=255)
                       cGCTATTGAAGTGCAACTGGCACATGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGAAatgat                                                                  >  2:5029280/1‑93 (MQ=255)
                        gCTATTGAAGTGCAACTGGCACATGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGAAatgatg                                                                 <  1:3856030/93‑1 (MQ=255)
                        gCTATTGAAGTGCAACTGGCACATGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGAAatgatg                                                                 <  1:4099231/93‑1 (MQ=255)
                         cTATTGAAGTGCAACTGGCACATGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGc                                                                >  1:1143424/1‑93 (MQ=255)
                         cTATTGAAGTGCAACTGGCACATGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGc                                                                <  2:3548207/93‑1 (MQ=255)
                          tATTGAAGTGCAACTGGCACATGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGCa                                                               >  2:659433/1‑93 (MQ=255)
                          tATTGAAGTGCAACTGGCACATGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGCa                                                               >  1:4997827/1‑93 (MQ=255)
                          tATTGAAGTGCAACTGGCACATGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGAGAAACGCAGAGAAATGATGCa                                                               >  1:4203613/1‑93 (MQ=255)
                           aTTGAAGTGCAACTGGCACATGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGCaa                                                              >  2:4230932/1‑93 (MQ=255)
                            ttGAAGTGCAACTGGCACATGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTGGATAAACGCAGAGAAATGATGCAAt                                                             <  2:1154085/93‑1 (MQ=255)
                              gAAGTGCAACTGGCACATGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGCAAtgg                                                           <  1:5332803/93‑1 (MQ=255)
                                aGTGCAACTGGCACATGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGCAAtggtg                                                         <  2:854627/93‑1 (MQ=255)
                                aGTGCAACTGGCACATGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGCAAtggtg                                                         <  1:3010009/93‑1 (MQ=255)
                                aGTGCAACTGGCACATGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGCAAtggtg                                                         <  1:2177707/93‑1 (MQ=255)
                                  tGCAACTGGCACATGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGCAATGGTggg                                                       <  1:2740574/93‑1 (MQ=255)
                                   gCAACTGGCACATGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGc                                                      >  1:3382660/1‑93 (MQ=255)
                                    cAACTGGCACATGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCg                                                     >  1:5504805/1‑93 (MQ=255)
                                      aCTGGCACATGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCGGa                                                   <  1:2540457/93‑1 (MQ=255)
                                        tGGCACATGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCGGACt                                                 <  1:259148/93‑1 (MQ=255)
                                         ggCACATGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATACTCAGCATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCGGACTg                                                <  1:3611843/93‑1 (MQ=255)
                                           cacaTGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCGGACTGGc                                              >  1:1630224/1‑93 (MQ=255)
                                            acaTGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCGGAGTgggt                                             >  1:4031684/1‑91 (MQ=255)
                                              aTGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCGGACTGGGttt                                           >  2:4396564/1‑92 (MQ=255)
                                               tGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCGGATTGGCTtga                                          <  2:6130332/93‑1 (MQ=255)
                                               tGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCGGACTGGCTtgt                                          >  2:6189066/1‑92 (MQ=255)
                                               tGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCGGACTGGCTtga                                          >  2:3794429/1‑93 (MQ=255)
                                                gCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCGGACTGGCTtgat                                         <  2:3452033/93‑1 (MQ=255)
                                                 cAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCGGACTGGCTtgatg                                        <  2:1065241/93‑1 (MQ=255)
                                                 cAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCGGACTGGCTtgatg                                        >  2:455221/1‑93 (MQ=255)
                                                 cAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCGGACTGGCTtgatg                                        <  1:4675500/93‑1 (MQ=255)
                                                   aaCGGCGGATGTGTGCGTGGGATTTACAACCATGCTCAGCATCTTGATAAACGCAGAGAAATGATGCAATGGTGGTCAGACTGGCTtgatggg                                      <  2:6260347/93‑3 (MQ=255)
                                                     cggcggATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCGGACTGGCTTGATGAGaa                                    <  1:846878/93‑1 (MQ=255)
                                                     cggcggATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCGGACTGGCTTGATGAGaa                                    <  2:1941405/93‑1 (MQ=255)
                                                      ggcggATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCGGCATGGCTTTGTGAGAAg                                   >  2:2411511/1‑93 (MQ=255)
                                                      ggcggATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCGGACTGGCTTGATGAGAAg                                   <  1:482539/93‑1 (MQ=255)
                                                       gcggATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCGGACTGGCTTGATGAGAAgg                                  >  2:763438/1‑93 (MQ=255)
                                                       gcggATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCGGACTGGCTTGATGAGAAgg                                  <  1:2800123/93‑1 (MQ=255)
                                                          gATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCGGACTGGCTTGATGAGAAGgggg                               >  2:397588/1‑93 (MQ=255)
                                                          gATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCGGACTGGCTTGATGAGAAGGTgg                               <  1:6167160/93‑1 (MQ=255)
                                                            tCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCGGACTGGCTTGATGAGAAGgggggg                             >  2:4089642/1‑91 (MQ=255)
                                                              tgtgCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCGGACTGGCTTGATGAGAAGGTGGAGTg                           <  2:3683727/93‑1 (MQ=255)
                                                              tgtgCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCGGACTGGCTTGATGAGAAGGTGGAGTg                           <  2:5517317/93‑1 (MQ=255)
                                                                tgCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCGGACTGGCTTGATGAGAAGGGggagggag                         >  1:893561/1‑93 (MQ=255)
                                                                 gCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCGGACTGGCTTGATGAGAAGGTGGAGTGAGc                        <  1:1417631/93‑1 (MQ=255)
                                                                 gCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCGGACTGGCTTGATGAGAAGGTGGAGTGAGc                        <  1:2834793/93‑1 (MQ=255)
                                                                  cGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCGGACTGGCTTGATGAGAAGGGggagggaggg                       >  1:1543569/1‑91 (MQ=255)
                                                                    tGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCGGACTGGCTTGATGAGAAGGTGGGGGGAGCGAc                     >  1:5130686/1‑93 (MQ=255)
                                                                    tGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCGGACTGGCTTGATGAGAAGGTGGAGTGAGCGAc                     <  2:5638293/93‑1 (MQ=255)
                                                                    tGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCGGACTGGCTTGATGAGAAGGGGGAGTGAGCGAc                     >  1:725795/1‑93 (MQ=255)
                                                                     gggATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCGGACTGGCTTGATGAGAAGGTGGGGTGAGCGAcc                    >  1:1800203/1‑93 (MQ=255)
                                                                     gggATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCGGACTGGCTTGATGAGAAGGTGGAGTGGGCGAcc                    >  2:4281704/1‑93 (MQ=255)
                                                                     gggATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCGGACTGGCTTGATGAGAAGGTGGAGTGAGCGAcc                    <  2:3195551/93‑1 (MQ=255)
                                                                     gggATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCGGACTGGCTTGATGAGAAGGGGGAGGGAGCGAcc                    >  1:5575919/1‑93 (MQ=255)
                                                                      ggATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCGGACTGGCTTGATGAGAAAGTGGAGTGAGCGACCt                   >  1:3931311/1‑93 (MQ=255)
                                                                         tttACAACCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCGGACTGGCTTGGTCAGTAGGGGGAGGGAGCGTCCttta                >  2:1264160/1‑91 (MQ=255)
                                                                         tttACAACCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCGGACTGGCTTGATGAGGAGGTGGAGGGAGAGACCTTaa                >  1:5648426/1‑93 (MQ=255)
                                                                         tttACAACCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCGGACTGGCTTGATGAGAAGGTGGAGTGAGCGACCTTaa                <  1:5173108/93‑1 (MQ=255)
                                                                          ttACAACCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCGGACTGGCTTGATGAGAAGGTGGAGTGAGCGACCTTaac               <  1:4250662/93‑1 (MQ=255)
                                                                            aCAACCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCGGACTGGCTTGATGAGAAAGTGGAGTGAGCGACCTTaacaa             >  2:2686815/1‑93 (MQ=255)
                                                                              aaCCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCGGACTGGCTTGATGAGAAGGTGGAGTGAGCGACCTTAACAACt           <  2:4105580/93‑1 (MQ=255)
                                                                              aaCCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCGGACTGGCTTGATGAGAAGGTGGAGTGAGCGAACTTAACAACt           >  2:5256080/1‑93 (MQ=255)
                                                                              aaCCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCGGACTGGCTTGATGAGAAGGTGGAGGGAGCGACCTTAACAACt           >  1:440194/1‑93 (MQ=255)
                                                                                  aTGCTCAGCATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCGGACTGGCTTGATGAGAAGGTGGAGTGAGCGACCTTAACAACTATCg       <  1:3291839/93‑1 (MQ=255)
                                                                                    gCTCAGCATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCGGACTGTCTTGATGAGAAGGTGGAGTGAGCGACCTTAACAACTATCGaa     <  1:4437164/93‑1 (MQ=255)
                                                                                    gCTCAGCATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCGGACTGGCTTGATGAGAAGGTGGAGTGAGCGACCTTAACAACTATCGaa     >  1:4982760/1‑93 (MQ=255)
                                                                                    gCTCAGCATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCGAACTGGCTTGATGAGAAGGTGGAGTGAGCGACCTTAACAACTATCGaa     <  2:4778607/93‑1 (MQ=255)
                                                                                     cTCAGCATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCGGACTGGCTTGATGAGAAGGTGGAGTGAGCGACCTTAACAACTATCGAAt    <  2:5639667/93‑1 (MQ=255)
                                                                                     cTCAGCATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCGGACTGGCTTGATGAGAAGGTGGAGTGAGCGACCTTAACAACTATCGAAt    <  2:5879691/93‑1 (MQ=255)
                                                                                      tCAGCATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCGGACTGGCTTGATGAGAAGGTGGAGTGAGCGACCTTAACAACTATCGAATa   <  2:5732520/93‑1 (MQ=255)
                                                                                      tCAGCATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCGGACTGGCTTGATGAGAAGGTGGAGTGAGCGACCTTAACAACTATCGAATa   <  1:4890910/93‑1 (MQ=255)
                                                                                       cAGCATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCGGACTGGCTTGATGAGAAGGTGGAGTGAGCGACCTTAACAACTATCGAATaa  >  1:4293263/1‑92 (MQ=255)
                                                                                       cAGCATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCGGACTGGCTTGATGAGAAGGTGGAGTGAGCGACCTTAACAACTATCGAATAg  <  1:3665748/93‑1 (MQ=255)
                                                                                          |                                                                                         
AACGAGCACGAATGGCCTGCTGACGCTATTGAAGTGCAACTGGCACATGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGTATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCGGACTGGCTTGATGAGAAGGTGGAGTGAGCGACCTTAACAACTATCGAATAG  >  NC_000913/2465567‑2465746

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