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,6570TC98.5% 636.1 / ‑2.1 197Y365H (TAT→CAT) intCputative prophage Sf6‑like integrase
Reads supporting (aligned to +/- strand):  ref base T (2/1);  new base C (98/96);  total (100/97)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 200 total aligned reads displayed.

AACGAGCACGAATGGCCTGCTGACGCTATTGAAGTGCAACTGGCACATGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGTATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCGGACTGGCTTGATGAGAAGGTGGAGTGAGCGACCTTAACAACTATCGAATAG  >  NC_000913/2465567‑2465746
                                                                                          |                                                                                         
 aCGAGCACGAATGGCCTGCTGACGCTATTGAAGTGCAACTGGCACATGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCAtc                                                                                        <  2:2808901/93‑1 (MQ=255)
 aCGAGCACGAATGGCCTGCTGACGCTATTGAAGTGCAACTGGCACATGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCAtc                                                                                        >  1:3706242/1‑93 (MQ=255)
   gAGCACGAATGGCCTGCTGACGCTATTGAAGTGCAACTGGCACATGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCAtctc                                                                                      >  2:4342274/1‑93 (MQ=255)
     gCACGAATGGCCTGCTGACGCTATTGAAGTGCAACTGGCACATGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGa                                                                                    <  2:1047021/93‑1 (MQ=255)
      cACGAATGGCCTGCTGACGCTATTGAAGTGCAACTGGCACATGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGAt                                                                                   <  2:64919/93‑1 (MQ=255)
      cACGAATGGCCTGCTGACGCTATTGAAGTGCAACTGGCACATGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGAt                                                                                   <  2:3199185/93‑1 (MQ=255)
       aCGAATGGCCTGCTGACGCTATTGAAGTGCAACTGGCACATGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGTATCTCGATa                                                                                  >  1:3841563/1‑93 (MQ=255)
       aCGAATGGCCTGCTGACGCTATTGAAGTGCAACTGGCACACGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATa                                                                                  <  1:1865671/93‑1 (MQ=255)
        cGAATGGCCTGCTGACGCTATTGAAGTGCAACTGGCACATGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATaa                                                                                 <  2:415234/93‑1 (MQ=255)
          aaTGGCCTGCTGACGCTATTGAAGTGCAACTGGCACATGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAAc                                                                               >  1:72555/1‑93 (MQ=255)
          aaTGGCCTGCTGACGCTATTGAAGTGCAACTGGCACATGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAAc                                                                               >  2:1384242/1‑93 (MQ=255)
          aaTGGCCTGCTGACGCTATTGAAGTGCAACTGGCACATGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAAc                                                                               >  2:2045389/1‑93 (MQ=255)
            tGGCCTGCTGACGCTATTGAAGTGCAACTGTCACATGCAAATGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGc                                                                             >  2:3730323/1‑93 (MQ=255)
             ggCCTGCTGACGCTATTGAAGTGCAACTGGCACATGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCa                                                                            <  2:1677504/93‑1 (MQ=255)
              gCCTGCTGACGCTATTGAAGTGCAACTGGCACATGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCag                                                                           <  2:1086674/93‑1 (MQ=255)
              gCCTGCTGACGCTATTGAAGTGCAACTGGCACATGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCag                                                                           >  2:3550112/1‑93 (MQ=255)
               ccTGCTGACGCTATTGAAGTGCAACTGGCACATGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCaga                                                                          <  1:4107017/93‑1 (MQ=255)
                ctgctgACGCTATTGAAGTGCAACTGGCACATGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCagag                                                                         >  2:1308064/1‑93 (MQ=255)
                 tgctgACGCTATTGAAGTGCAACTGGCACATGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAgaga                                                                        <  1:738681/93‑1 (MQ=255)
                 tgctgACGCTATTGAAGTGCAACTGGCACATGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAgaga                                                                        >  1:271781/1‑93 (MQ=255)
                 tgctgACGCTATTGAAGTGCAACTGGCACATGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAgaga                                                                        >  1:27895/1‑93 (MQ=255)
                  gctgACGCTATTGAAGTGCAACTGGCACATGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGaa                                                                       <  1:3530050/93‑1 (MQ=255)
                   ctgACGCTATTGAAGTGCAACTGGCACATGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGaaa                                                                      <  1:865591/93‑1 (MQ=255)
                    tgACGCTATTGAAGTGCAACTGGCACATGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGAAat                                                                     <  2:855299/93‑1 (MQ=255)
                    tgACGCTATTGAAGTGCAACTGGCACATGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGAAat                                                                     >  1:1766499/1‑93 (MQ=255)
                    tgACGCTATTGAAGTGCAACTGGCACATGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGAAat                                                                     >  1:2363727/1‑93 (MQ=255)
                     gACGCTATTGAAGTGCAACTGGCACATGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGAAatg                                                                    >  2:1229627/1‑93 (MQ=255)
                      aCGCTATTGAAGTGCAACTGGCACATGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGAAatga                                                                   >  2:352053/1‑93 (MQ=255)
                       cGCTATTGAAGTGCAACTGGCACAGGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGAAatgat                                                                  <  2:175352/93‑1 (MQ=255)
                         cTATTGAAGTGCAACTGGCACATGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGc                                                                <  1:4324693/93‑1 (MQ=255)
                         cTATTGAAGTGCAACTGGCACATGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGc                                                                >  2:3581576/1‑93 (MQ=255)
                          tATTGAAGTGCAACTGGCACATGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGCa                                                               >  1:1837343/1‑93 (MQ=255)
                           aTTGAAGTGCAACTGGCACATGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGCaa                                                              <  2:4525191/93‑1 (MQ=255)
                            ttGAAGTGCAACTGGCACATGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGCAAt                                                             >  2:686960/1‑93 (MQ=255)
                            ttGAAGTGCAACTGGCACATGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGCAAt                                                             <  2:2454248/93‑1 (MQ=255)
                            ttGAAGTGCAACTGGCACATGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGCAAt                                                             >  2:2082158/1‑93 (MQ=255)
                             tGAAGTGCAACTGGCACATGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGCAAtg                                                            <  1:819826/93‑1 (MQ=255)
                               aaGTGCAACTGGCACATGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGCAAtggt                                                          >  2:3765454/1‑93 (MQ=255)
                               aaGTGCAACTGGCACATGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCACCATCTCGATAAACGCAGAGAAATGATGCAAtggg                                                          >  1:1240466/1‑92 (MQ=255)
                                 gTGCAACTGGCACATGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGCAAtggtgg                                                        >  2:1594304/1‑93 (MQ=255)
                                   gCAACTGGCACATGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGCAATGGTgggg                                                      >  1:1938652/1‑92 (MQ=255)
                                     aaCTGGCACATGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCgg                                                    <  1:2428379/93‑1 (MQ=255)
                                      aCTGGCACATGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCGGa                                                   >  1:4431270/1‑93 (MQ=255)
                                       cTGGCACATGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCGGAc                                                  >  1:2097142/1‑93 (MQ=255)
                                       cTGGCACATGCAAACGGCGGATCTGTGCGAAGGATTTACAACCATGCTCAGTATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCGGAc                                                  >  2:3029844/1‑93 (MQ=255)
                                       atgacacaTGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGATCAGCATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCGGAc                                                  <  1:3726869/89‑1 (MQ=255)
                                           cacaTGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCGGAGTggt                                              >  2:1016308/1‑92 (MQ=255)
                                           cacaTGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCGGACTGGc                                              <  1:4095990/93‑1 (MQ=255)
                                            acaTGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCGGACTGGCt                                             <  1:3004243/93‑1 (MQ=255)
                                            acaTGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCGGACTGGCt                                             >  1:2757570/1‑93 (MQ=255)
                                             caTGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCGGACTGGCtt                                            >  1:2548726/1‑93 (MQ=255)
                                             caTGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCGGACTGGCtt                                            <  1:2752081/93‑1 (MQ=255)
                                              aTGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCGGACTGGCTtg                                           <  2:2724604/93‑1 (MQ=255)
                                              aTGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCGGACTGGCTtg                                           <  1:1030182/93‑1 (MQ=255)
                                              aTGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCGGACTGGCTtg                                           >  2:144771/1‑93 (MQ=255)
                                               tGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCGGACTGGCTtga                                          <  2:3666728/93‑1 (MQ=255)
                                                gCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCGGACTGGCTtgat                                         <  2:1831800/93‑1 (MQ=255)
                                                 cAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCGGACTGGCTtgatg                                        <  1:1672909/93‑1 (MQ=255)
                                                  aaaCGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCGGACTGGCTtgatga                                       <  1:2028185/93‑1 (MQ=255)
                                                  aaaCGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCGGACTGGCCtggtgg                                       >  1:2955685/1‑92 (MQ=255)
                                                   aaCGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCGGACTGGCTTTATgag                                      >  2:1471705/1‑93 (MQ=255)
                                                       gcggATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCGGACTGGCTTGATGAGAAgg                                  <  2:1298325/93‑1 (MQ=255)
                                                           aTCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCGGACTGGCTTGATGAGAAGGTggg                              >  2:4416932/1‑92 (MQ=255)
                                                           aTCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCGGACTGGCTTGATGAGAAGGTGGa                              <  2:528860/93‑1 (MQ=255)
                                                           aTCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCGGACTGGCTTGATGAGAAGGTGGa                              <  2:3289300/93‑1 (MQ=255)
                                                           aTCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCGGACTGGCTTGATGAGAAGGTGGa                              <  2:3174393/93‑1 (MQ=255)
                                                            tCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCGGACTGGCTTGATGAGAAGgggggg                             >  2:1145554/1‑91 (MQ=255)
                                                            tCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCGGACTGGCTTGATGAGAAGGGGGAg                             >  1:1593688/1‑93 (MQ=255)
                                                             cTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCGGACTGGCTTTATGAGGAGGGGGAgg                            >  2:2280664/1‑92 (MQ=255)
                                                             cTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCGGACTGGCTTGATGAGAAGGTGGAGt                            <  1:1874960/93‑1 (MQ=255)
                                                              tgtgCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCGGACTGGCTTGATGAGAAGGGGGAggg                           >  1:3176417/1‑91 (MQ=255)
                                                               gtgCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCGGACTGGCTTGATGAGAAGGTGGAGTGa                          <  1:3370408/93‑1 (MQ=255)
                                                                tgCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCGGACTGGCTTGATGAGAAGGTGGAGTGAg                         <  1:779147/93‑1 (MQ=255)
                                                                ggCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCGGACTGGCTTGATGAGAAGGTGGAGTGAg                         <  1:2005878/92‑1 (MQ=255)
                                                                 gCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCGGACTGGCTTGATGAGAAGGTGGAGTGAGc                        <  2:757796/93‑1 (MQ=255)
                                                                 gCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCGGACTGGCTTGATGAGAAGGTGGAGTGAGc                        <  2:3359203/93‑1 (MQ=255)
                                                                  cGTGGGATTTACAACCATGCTCAGCAACTCGATAAACGCAGAGAAATGATGCAATGGTGGGCGGACTGGCTTGATGAGAAGGTGGAGTGAGCg                       <  2:4216516/93‑1 (MQ=255)
                                                                   gTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCGGACTGGCTTGATGAGAAGGGGGGGGGAGCgc                      >  1:1420388/1‑92 (MQ=255)
                                                                    tGGTATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCGGACTGGCTTGATGAGAAGGTGGAGTGAGCGAc                     <  1:30384/93‑1 (MQ=255)
                                                                    tGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCGGACTGGCTTGATGAGAAGGTGGAGTGAGCGAc                     <  2:4476462/93‑1 (MQ=255)
                                                                    tGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCGGACTGGCTTGATGAGAAGGGGGAGGGAGAGAc                     >  2:2169100/1‑93 (MQ=255)
                                                                        aTTTACAACCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCGGACTGGCTTGATGAGAAGGTGGAGTGAGCGACCTTa                 <  1:3997040/93‑1 (MQ=255)
                                                                         tttACAACCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCGGACTGGCTTGATGAGAAGGGGGGGGGAGCGACCTTaa                >  1:1794772/1‑93 (MQ=255)
                                                                         tttACAACCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCGGACTGGCTTGATGAGAAGGGGGAGTGAGCGACCTTaa                >  2:4023867/1‑93 (MQ=255)
                                                                          ttACAACCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCGGACTGGCTTGATGAGAAGGTGGAGTGAGCGACCTTaac               <  1:2232740/93‑1 (MQ=255)
                                                                            aCAACCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCGGACTGGCTTGATGAGAAGGTGGAGTGAGCGACCTTaacaa             <  1:4604446/93‑1 (MQ=255)
                                                                            aCAACCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCGGACTGGCTTGATGAGAAGGGGGAGTGGGCGACCTTaacaa             >  1:55876/1‑93 (MQ=255)
                                                                              aaCCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCGGACTGGCTTGATGAGAAGGTGGAGTGAGCGACCTTAACAACt           <  2:3898362/93‑1 (MQ=255)
                                                                               aCCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCGGACTGGCTTGATGAGAAGGTGGAGTGAGCGACCTTAACAACTa          >  1:287080/1‑93 (MQ=255)
                                                                               aCCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCGGACTGGCTTGATGAGAAGGTGGAGTGAGCGACCTTAACAACTa          >  1:3601249/1‑93 (MQ=255)
                                                                                 cATGCTCAGCATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCGGACTGGCTTGATGAGAAGGTGGAGTGAGCGACCTTAACAACTATc        >  1:1678229/1‑93 (MQ=255)
                                                                                  aTGCTCAGCATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCGGACTGGCTTGATGAGGAGGGGGAGTGAGCGACCTTAACAACTATCg       >  1:500786/1‑93 (MQ=255)
                                                                                  aTGCTCAGCATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCGGACTGGCTTGATGAGAAGGTGGAGGGAGCGACCTTAAAAACTATCg       >  2:1952637/1‑93 (MQ=255)
                                                                                   tGCTCAGCATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCGGACTGGCTTGATGAGAAGGTGGAGTGAGCGACCTTAACAACTATCGa      >  1:1881216/1‑93 (MQ=255)
                                                                                     cTCAGTCTCTCGAGAAACGCAGAGAAATTGGGCAATGGGGGGCGGACTGGCTTGATGAGAAGGTTGAGTGAGCGACCTTAACAACTATCGAAt    <  1:1171606/93‑1 (MQ=255)
                                                                                      tCAGCATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCGGACTGGCTTGATGAGAAGGTGGAGTGAGCGACCTTAACAACTATCGAATa   >  2:4246237/1‑93 (MQ=255)
                                                                                      tCAGCATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCGGACTGGCTTGATGAGAAGGGGGAGGGAGCGACCTTAACAACTATCGAATa   >  1:1305611/1‑93 (MQ=255)
                                                                                       cAGCATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCGGACTGGCTTGGTGAGAAGGTGGAGGGAGAGACCTTAACAACTATTGGATAg  >  2:223774/1‑93 (MQ=255)
                                                                                       cAGCATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCGGACTGGCTTGATGAGAAGGTGGAGTGAGCGACCTTAACAACTATCGAATAg  <  1:3772353/93‑1 (MQ=255)
                                                                                       cAGCATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCGGACTGGCTTGATGAGAAGGTGGAGTGAGCGACCTTAACAACTATCGAATAg  <  2:2477028/93‑1 (MQ=255)
                                                                                          |                                                                                         
AACGAGCACGAATGGCCTGCTGACGCTATTGAAGTGCAACTGGCACATGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGTATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCGGACTGGCTTGATGAGAAGGTGGAGTGAGCGACCTTAACAACTATCGAATAG  >  NC_000913/2465567‑2465746

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