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% 846.2 / ‑5.1 266Y365H (TAT→CAT) intCputative prophage Sf6‑like integrase
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (139/126);  total (140/126)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.62e-01

Only 100 of 269 total aligned reads displayed.

AACGAGCACGAATGGCCTGCTGACGCTATTGAAGTGCAACTGGCACATGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGTATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCGGACTGGCTTGATGAGAAGGTGGAGTGAGCGACCTTAACAACTATCGAATAG  >  NC_000913/2465567‑2465746
                                                                                          |                                                                                         
 aCGAGCACGAATGGCCTGCTGACGCTATTGAAGTGCAACTGGCACATGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCAtc                                                                                        >  1:6620681/1‑93 (MQ=255)
 aCGAGCACGAATGGCCTGCTGACGCTATTGAAGTGCAACTGGCACATGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCAtc                                                                                        >  1:3953633/1‑93 (MQ=255)
  cGAGCACGAATGGCCTGCTGACGCTATTGAAGTGCAACTGGCACATGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCAtct                                                                                       <  2:2667923/93‑1 (MQ=255)
  cGAGCACGAATGGCCTGCTGACGCTATTGAAGTGCAACTGGCACATGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCAtct                                                                                       <  1:345692/93‑1 (MQ=255)
   gAGCACGAATGGCCTGCTGACGCTATTGAAGTGCAACTGGCACATGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCAtctc                                                                                      >  2:6963380/1‑93 (MQ=255)
   gAGCACGAATGGCCTGCTGACGCTATTGAAGTGCAACTGGCACATGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCAtctc                                                                                      <  1:1445351/93‑1 (MQ=255)
   gAGCACGAATGGCCTGCTGACGCTATTGAAGTGCAACTGGCACATGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCAtctc                                                                                      <  2:1012871/93‑1 (MQ=255)
     gCACGAATGGCCTGCTGACGCTATTGAAGTGCAACTGGCACATGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGa                                                                                    >  1:700947/1‑93 (MQ=255)
      cACGAATGGCCTGCTGACGCTATTGAAGTGCAACTGGCACATGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGAt                                                                                   >  2:1283074/1‑93 (MQ=255)
       aCGAATGGCCTGCTGACGCTATTGAAGTGCAACTGGCACATGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATa                                                                                  >  1:6417892/1‑93 (MQ=255)
        cGAATGGCCTGCTGACGCTATTGAAGTGCAACTGGCACATGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATaa                                                                                 >  1:1642656/1‑93 (MQ=255)
        cGAATGGCCTGCTGACGCTATTGAAGTGCAACTGGCACATGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATaa                                                                                 <  1:2382478/93‑1 (MQ=255)
          aaTGGCCTGCTGACGCTATTGAAGTGCAACTGGCACATGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAAc                                                                               <  2:2918066/93‑1 (MQ=255)
          aaTGGCCTGCTGACGCTATTGAAGTGCAACTGGCACATGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAAc                                                                               >  1:2446663/1‑93 (MQ=255)
            tGGCCTGCTGACGCTATTGAAGTGCAACTGGCACATGCAAACGGCGGAGCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGc                                                                             >  1:1192909/1‑93 (MQ=255)
             ggCCTGCTGACGCTATTGAAGTGCAACTGGCACATGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCa                                                                            >  2:6751662/1‑93 (MQ=255)
               ccTGCTGACGCTATTGAAGTGCAACTGGCACATGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCaga                                                                          <  1:504697/93‑1 (MQ=255)
                ctgctgACGCTATTGAAGTGCAACTGGCACATGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCagag                                                                         >  2:3105537/1‑93 (MQ=255)
                ctgctgACGCTATTGAAGTGCAACTGGCACATGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCagag                                                                         >  2:4677967/1‑93 (MQ=255)
                 tgctgACGCTATTGAAGTGCAACTGGCACATGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAgaga                                                                        >  2:4738399/1‑93 (MQ=255)
                 tgctgACGCTATTGAAGTGCAACTGGCACATGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAgaga                                                                        >  1:283662/1‑93 (MQ=255)
                 tgctgACGCTATTGAAGTGCAACTGGCACATGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAgaga                                                                        >  1:7446572/1‑93 (MQ=255)
                 tgctgACGCTATTGAAGTGCAACTGGCACATGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAgaga                                                                        >  1:1805744/1‑93 (MQ=255)
                 tgctgACGCTATTGAAGTGCAACTGGCACATGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAgaga                                                                        >  1:3860281/1‑93 (MQ=255)
                   ctgACGCTATTGAAGTGCAACTGGCACATGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGaaa                                                                      <  2:1966421/93‑1 (MQ=255)
                   ctgACGCTATTGAAGTGCAACTGGCACATGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGaaa                                                                      >  2:4978216/1‑93 (MQ=255)
                   ctgACGCTATTGAAGTGCAACTGGCACATGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGaaa                                                                      <  1:1078960/93‑1 (MQ=255)
                    tgACGCTATTGAAGTGCAACTGGCACATGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGAAat                                                                     <  1:5344551/93‑1 (MQ=255)
                    tgACGCTATTGAAGTGCAACTGGCACATGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGAAat                                                                     <  2:4173901/93‑1 (MQ=255)
                      aCGCTATTGAAGTGCAACTGGCACATGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGAAatga                                                                   >  1:6456002/1‑93 (MQ=255)
                       cGCTATTGAAGTGCAACTGGCACATGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGAAatgat                                                                  <  1:2552906/93‑1 (MQ=255)
                        gCTATTGAAGTGCAACTGGCACATGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGAAatgatg                                                                 >  2:5135524/1‑93 (MQ=255)
                        gCTATTGAAGTGCAACTGGCACATGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGAAatgatg                                                                 >  2:3974339/1‑93 (MQ=255)
                         cTATTGAAGTGCAACTGGCACATGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGc                                                                >  1:4210966/1‑93 (MQ=255)
                         cTATTGAAGTGCAACTGGCACATGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGc                                                                >  2:3697417/1‑93 (MQ=255)
                          tATTGAAGTGCAACTGGCACATGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGNTAAACGCAGAGAANTGATGCa                                                               >  1:4699824/1‑93 (MQ=255)
                          tATTGAAGTGCAACTGGCACATGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGCa                                                               >  2:3902104/1‑93 (MQ=255)
                          tATTGAAGTGCAACTGGCACATGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGCa                                                               >  2:3337543/1‑93 (MQ=255)
                           aTTGAAGTGCAACTGGCACATGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGCaa                                                              >  1:3620046/1‑93 (MQ=255)
                            ttGAAGTGCAACTGGCACATGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGCAAt                                                             >  2:4116171/1‑93 (MQ=255)
                             tGAAGTGCAACTGGCACATGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGCAAtg                                                            >  1:7343710/1‑93 (MQ=255)
                             tGAAGTGCAACTGGCACATGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGCAAtg                                                            <  2:1762306/93‑1 (MQ=255)
                             tGAAGTGCAACTGGCACATGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGCAAtg                                                            <  1:936307/93‑1 (MQ=255)
                             tGAAGTGCAACTGGCACATGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGCAAtg                                                            <  2:2805931/93‑1 (MQ=255)
                             tGAAGTGCAACTGGCACATGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGCAAtg                                                            <  2:3787727/93‑1 (MQ=255)
                              gAAGTGCAACTGGCACATGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGCAAtgg                                                           <  1:7526069/93‑1 (MQ=255)
                                aGTGCAACTGGCACATGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGCAAtggtg                                                         >  2:2671362/1‑93 (MQ=255)
                                aGTGCAACTGGCACATGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGCAAtggtg                                                         >  2:520007/1‑93 (MQ=255)
                                     aaCTGGCACATGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCgg                                                    <  1:6245469/93‑1 (MQ=255)
                                     aaCTGGCACATGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCgg                                                    >  1:387850/1‑93 (MQ=255)
                                     aaCTGGCACATGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCgg                                                    >  1:3581247/1‑93 (MQ=255)
                                      aCTGGCACATGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCGGa                                                   <  2:1279467/93‑1 (MQ=255)
                                        tGGCACATGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCGGGCt                                                 >  2:923802/1‑93 (MQ=255)
                                        tGGCACATGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCGGACt                                                 <  2:5438584/93‑1 (MQ=255)
                                        tGGCACATGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCGGACt                                                 >  1:6168131/1‑93 (MQ=255)
                                         ggCACATGCAACCGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCGGACTg                                                <  2:7180038/93‑1 (MQ=255)
                                         ggCACATGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCGGACTg                                                <  1:7763010/93‑1 (MQ=255)
                                          gCACATGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCGGACTgg                                               <  2:2431658/93‑1 (MQ=255)
                                             caTGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGGGGACTGGCtt                                            >  2:2168619/1‑93 (MQ=255)
                                             caTGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCGGACTGGCtt                                            <  1:5688118/93‑1 (MQ=255)
                                             caTGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCGGACTGGCtt                                            <  1:2243559/93‑1 (MQ=255)
                                               tGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCGGACTGGCttta                                          >  2:6384512/1‑91 (MQ=255)
                                               tGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCGGACTGGCTtgg                                          >  2:3819483/1‑92 (MQ=255)
                                                 cAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCGGACTGGCTtgatg                                        >  2:570716/1‑93 (MQ=255)
                                                  aaaCGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCGGACTGGCTtggtga                                       >  1:1910039/1‑93 (MQ=255)
                                                  aaaCGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCGGACTGGCTtgatga                                       >  2:6168723/1‑93 (MQ=255)
                                                        cggATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCGGACTGGCTTGATGAGAAGGt                                 <  1:1716327/93‑1 (MQ=255)
                                                        cggATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCGGACTGGCTTGATGAGAAGGt                                 <  2:3511987/93‑1 (MQ=255)
                                                         ggATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCGGACTGGCTTTATGAGAAgggg                                >  2:5787505/1‑91 (MQ=255)
                                                         ggATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCGGACTGGCTTGATGAGAAGGTg                                <  2:2372085/93‑1 (MQ=255)
                                                          gATCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCGGACTGGCTTGATGAGAAGGTgg                               <  1:692084/93‑1 (MQ=255)
                                                            tCTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCGGACTGGCTTGATGAGAAGGTGGAg                             <  1:6820136/93‑1 (MQ=255)
                                                             cTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCGGACTGGCTTGATGAGAAGGTGGAGt                            <  2:221481/93‑1 (MQ=255)
                                                             cTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCGGACTGGCTTGATGAGAAGGTGGAGt                            <  1:7639930/93‑1 (MQ=255)
                                                             cTGTGCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCGGACTGGCTTGATGAGAAGGGGGAgg                            >  2:6387880/1‑92 (MQ=255)
                                                              tgtgCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCGGACTGGCTTGATGAGAAGGGGGAGTg                           >  1:5418580/1‑93 (MQ=255)
                                                              tgtgCGTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCGGACTGGCTTGATGAGAAGGGGGAGTg                           >  1:4858694/1‑93 (MQ=255)
                                                                   gTGGGATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCGGACTGGCTTGATGAGAAAGTGGAGGGAGCGa                      >  2:3064495/1‑93 (MQ=255)
                                                                     gggATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCGGACTGGCTTGATGAGAAGGTGGAGTGAGCGAcc                    >  1:3649876/1‑93 (MQ=255)
                                                                      ggATTTACAACCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCGGACTGGCTTGATGAGAAGGTGGAGTGAGCGACCt                   <  1:5637694/93‑1 (MQ=255)
                                                                         tttACAACCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCGGACTGGCTTGATGAGAAGGTGGAGTGAGCGACCTTaa                <  2:1657373/93‑1 (MQ=255)
                                                                            taaaCCATGCTCAGCATTTCAAAAAACGAAGAAAAATGATGCAATGGTGGGCGGACTGGCTTGATGAGAAGGTGGAGTGAGCGACCTTaacaa             <  1:6550438/91‑1 (MQ=255)
                                                                              aaCCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCGGACTGGCTTGATGAGAAGGTGGAGTGAGCGACCTTaaccaat           >  1:5428485/1‑89 (MQ=255)
                                                                              aaCCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCGGACTGGCTTGATGAGAAGGTGGAGTGAGCGACCTTAACAACt           >  2:7890494/1‑93 (MQ=255)
                                                                              aaCCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCGGACTGGCTTGATGAGAAGGGGGGGGGAGCGACCTTAACAACt           >  1:1823617/1‑93 (MQ=255)
                                                                               aCCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCGGACTGGCTTGATGAGAAGGTGGGGTGAGCGACCTTAACAACTa          >  1:7261356/1‑93 (MQ=255)
                                                                               aCCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCGGACTGGCTTGATGAGAAGGTGGAGTGAGCGACCTTAACAACTa          >  2:7712516/1‑93 (MQ=255)
                                                                               aCCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCGGACTGGCTTGATGAGAAGGTGGAGTGAGCGACCTTAACAACTa          >  1:3314144/1‑93 (MQ=255)
                                                                               aCCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCGGACTGGCTTGATGAGAAGGTGGAGTGAGCGACCTTAACAACTa          <  1:3195435/93‑1 (MQ=255)
                                                                               aCCATGCTCAGCATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCGGACTGGCTTGATGAGAAGGTGGAGTGAGCGACCTTAACAACTa          <  2:443899/93‑1 (MQ=255)
                                                                                ccATGCTCAGCATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCGGACTGGCTTGATGAGAAGGTGGAGTGAGCGACCTTAACAACTAt         >  2:7343471/1‑93 (MQ=255)
                                                                                   tGCTCAGCATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCGGACTGGCTTGATGAGAAGGTGGAGTGAGCGACCTTAACAACTATCGa      <  2:7519819/93‑1 (MQ=255)
                                                                                   tGCTCAGCATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCGGACTGGCTTGATGAGAAAGGGGGGTTAGCGACCTTAACAACTATCGa      >  2:7606698/1‑93 (MQ=255)
                                                                                    gCTCAGCATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCGGACTGGCTTTATTAGAAGGGGGAGGGAGCGACCTTAACAACTATCGaa     >  2:5335118/1‑93 (MQ=255)
                                                                                    gCTCAGCATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCGGACTGGCTTGATGAGAAGGTGGAGTGAGCGACCTTAACAACTATCGaa     <  1:5755590/93‑1 (MQ=255)
                                                                                    gCTCAGCATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCGGACTGGCTTGATGAGAAGGTGGAGTGAGCGACCTTAACAACTATCGaa     <  2:1443516/93‑1 (MQ=255)
                                                                                      tCAGCATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCGGACTGGCTTGATGAGAAGGTGGAGTGAGCGACCTTAACAACTATCGAATa   <  2:6833890/93‑1 (MQ=255)
                                                                                      tCAGCATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCGGACTGGCTTGATGAGAAAGTGGAGTGAGCGACCTTAACAACTCTCGAATa   >  1:3734124/1‑93 (MQ=255)
                                                                                       cAGCATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGGGGACTGGCTTGGTGAGAAGGGGGGGGGAGCGACCTTAACAACTAACGAAtcg  >  2:6826683/1‑91 (MQ=255)
                                                                                       cAGCATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCGGACTGGCTTGATGAGAAGGTGGAGTGAGCGACCTTAACAACTATCGAATAg  <  1:4635510/93‑1 (MQ=255)
                                                                                          |                                                                                         
AACGAGCACGAATGGCCTGCTGACGCTATTGAAGTGCAACTGGCACATGCAAACGGCGGATCTGTGCGTGGGATTTACAACCATGCTCAGTATCTCGATAAACGCAGAGAAATGATGCAATGGTGGGCGGACTGGCTTGATGAGAAGGTGGAGTGAGCGACCTTAACAACTATCGAATAG  >  NC_000913/2465567‑2465746

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: