Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 448,768 C→T 15.5% W625* (TGG→TAG)  cyoB ← cytochrome o ubiquinol oxidase subunit I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913448,7680CT15.5% 372.7 / 46.0 148W625* (TGG→TAG) cyoBcytochrome o ubiquinol oxidase subunit I
Reads supporting (aligned to +/- strand):  ref base C (58/66);  new base T (12/11);  total (71/77)
Fisher's exact test for biased strand distribution p-value = 6.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 155 total aligned reads displayed.

CCTGCCTTAGTAATCTCATCGAAATGCTGGTTTTCCAGTTTTTCGATTTCTGCCACCGGCACGTAGTAATCCACGTCCTCGTCGAAGCTTTTCACGATCCAGGTGATGATCATGCCTGCGAAGCCAACAATCGCCAGCCACCAGATATGCCAGATCATGGCGAAACCGAAGATGGTGGAGAAAGCTGCAATGACGATAC  >  NC_000913/448669‑448867
                                                                                                   |                                                                                                   
ccTGCCTTAGTAATCTCATCGAAATGCTGGTTTTCCAGTTTTTCGATTTCTGCCACCGGCACGTAGTAATCCACGTCCTCGTCGAAGCTTTTCACGATcc                                                                                                     >  1:2729330/1‑100 (MQ=255)
 cTGCCTTAGTAATCTCATCGAAATGCTGGTTTTCCAGTTTTTCGATTTCTGCCACCGGCACGTAGTAATCCACGTCCTCGTCGAAGCTTTTAACGATCCa                                                                                                    >  1:905717/1‑100 (MQ=255)
 cTGCCTTAGTAATCTCATCGAAATGCTGGTTTGCCAGTTTTTCGATTTCTCCCCCCGGCACGTAGTAATCCACGTCCTCATCAAAGCTTTTCACGACCCa                                                                                                    >  2:901712/1‑100 (MQ=255)
   gCCTTAGTAATCTCATCGAAATGCTGGTTTTCCAGTTTTTCGATTTCTGCCACCGGCACGTAGTAATCCACGTCCTCGTCGAAGCTTTTCACGATCCAgg                                                                                                  <  1:1217139/100‑1 (MQ=255)
   gCCTTAGTAATCTCATCGAAATGCTGGTTTTCCAGTTTTTCGATTTCTGCCACCGGCACGTAGTAATCCACGTCCTCGTCGAAGCTTTTCACGATCCAgg                                                                                                  <  2:3317257/100‑1 (MQ=255)
   gCCTTAGTAATCTCATCGAAATGCTGGTTTTCCAGTTTTTCGATTTCTGCCACCGGCACGTAGTAATCCACGTCCTCGTCGAAGCTTTTCACGATCCAgg                                                                                                  <  2:2408656/100‑1 (MQ=255)
    ccTTAGTAATCTCATCGAAATGCTGGTTTTCCAGTTTTTCGATTTCTGCCACCGGCACGTAGTAATCCACGTCCTCGTCGAAGCTTTTCACGATCTAGGt                                                                                                 >  1:2925045/1‑100 (MQ=255)
      ttAGTAATCTCATCGAAATGCTGGTTTTCCAGTTTTTCGATTTCTGCCACCGGCACGTAGTAATCCACGTCCTCGTCGAAGCTTTTCACGATCCAGGtga                                                                                               >  2:3004149/1‑100 (MQ=255)
        aGTAATCTCATCGAAATGCTGGTTTTCCAGTTTTTCGATTTCTGCCACCGGCACGTAGTAATCCACGTCCCCGTCGAAGCTTTTCACGATCCAGGtgatg                                                                                             >  2:1473079/1‑100 (MQ=255)
         gTAATCTCATCGAAATGCTGGTTTTCCAGTTTTTCGATTTCTGCCACCGGCACGTAGTAATCCACGTCCTCGTCGTAGCTTTTCACGATCTAGGagatga                                                                                            >  1:1600587/1‑100 (MQ=255)
         gTAATCTCATCGAAATGCTGGTTTTCCAGTTTTTCGATTTCTGCCACCGGCACGTAGTAATCCACGTCCTCGTCGAAGCTTTTCACGATCCAGGtgatga                                                                                            <  2:1307443/100‑1 (MQ=255)
          tAATCTCATCGAAATGCTGGTTTTCCAGTTTTTCGATTTCTGCCACCGGCACGTAGTAATCCACGTCCTCGTCGAAGATTTTCACGTTCCAGGTgatgat                                                                                           <  1:2634545/100‑1 (MQ=255)
            aTCTCATCGAAATGTTGGTTTTCCAGTTTTTTTTTTTTTTTTTCCGGCAGGTAGTAATCCACGTCCTCGTCGTTGCTTTTCACGATCCAGGtgat                                                                                              >  1:2422210/1‑95 (MQ=255)
             tctcATCGAAATGCTGGTTTTCCAGTTTTTTTTTTTTTTTCACCGGCACGTAGTAATCCACGTCCTCGTCGAAGCTTTTCACGATCTAGGTGATGATCAt                                                                                        >  1:1616438/1‑100 (MQ=255)
             tctcATCGAAATGCTGGTTTTCCAGTTTTTCGATTTCTGCCACCGGCACGTAGTAATCCACGTCCTCGTCGAAGCTTTTCACGATCCAGGTGATGATCAt                                                                                        >  2:622338/1‑100 (MQ=255)
              ctcATCGAAATGCTGGTTTTCCAGTTTTTCGATTTCTGCCACCGGCACGTAGTAATCCACGTCCTCGTCGAAGCTTTTCACGATCCAGGTGATGATCATg                                                                                       <  1:1187149/100‑1 (MQ=255)
               tcATCGAAACGCTGTTTTTCCAGTTTTTCGATTTCTGCCACCGGCACGTAGTAATACCCGTCCTCCTCGAAGCTTTTCACGATCCACATGATGATCCTGc                                                                                      >  1:125738/1‑100 (MQ=255)
                cATCGAAATGCTGGTTTTCCAGTTTTTCGATTTCTGCCACCGGCACGTAGTAATCCACGTCCTCGTCGAAGCTTTTCACGATCTAGGTGATGATCATGcc                                                                                     >  2:1117055/1‑100 (MQ=255)
                cATCGAAATGCTGGTTTTCCAGTTTTTCGATTTCTGCCACCGGCACGTAGTAATCCACGTCCTCGTCGAAGCTTTTCACGATCCAGGTGATGATCATGcc                                                                                     <  1:780860/100‑1 (MQ=255)
                    gAAATGCTGGTTTTCCAGTTTTTCGATTTCTGCCACCGGCACGTAGTAATCCACGTCCTCGTCGAAGCTTTTCACGATCCAGGTGATGATCATGCCTGCg                                                                                 >  1:1916005/1‑100 (MQ=255)
                    gAAATGCTGGTTTTCCAGTTTTTCGATTTCTGCCACCGGCACGTAGTAATCCACGTCCTCGTCGAAGCTTTTCACGATCCAGGTGATGATCATGCCTGCg                                                                                 >  2:1543401/1‑100 (MQ=255)
                     aaaTGCTGGTTTTCCAGTTTTTTTTTTTTTTTTACCGGCACGTAGTAATCCCCGTCCTCGTCGTTGCTTTTCACGATCCAGGTGATGATCATGCCTGCGa                                                                                >  1:3763167/1‑100 (MQ=255)
                       aTGCTGGTTTTCCAGTTTTTCGATTTCTGCCACCGGCACGTAGTAATCCACGTCCTCGTCGAAGCTTTTCACGATCCAGGTGATGATCATGCCTGCGAAg                                                                              <  2:3609835/100‑1 (MQ=255)
                       aTGCTGGTTTTCCAGTTTTTCGATTTCTGCCACCGGCACGTAGTAATCCACGTCCTCGTCGAAGCTTTTCACGATCCAGGTGATGATCATGCCTGCGAAg                                                                              <  2:242869/100‑1 (MQ=255)
                         gCTGGTTTTCCAGTTTTTCGATTTCTGCCACCGGCACGTAGTAATCCACGTCCTCGTCGAAGCTTTTCACGATCTAGGTGATGATCATGCCTGCGAAGcc                                                                            >  2:2000503/1‑100 (MQ=255)
                         gCTGGTTTTCCAGTTTTTCGATTTCTGCCACCGGCACGTAGTAATCCACGTCCTCGTCGAAGCTTTTCACGATCCAGGTGATGATCATGCCTGCGAAGcc                                                                            >  1:1184285/1‑100 (MQ=255)
                          cTGGTTTTCCAGTTTTTCGATTTCTGCCACCGGCACGTAGTAATCCACGTCCTCGTCGAAGCTTTTCACGATCCAGGTGATGATCATGCCTGCGAAGccn                                                                           <  1:3294024/100‑2 (MQ=255)
                            ggTTTTCCAGTTTTTCGATTTCTGCCACCGGCACGTAGTAATCTACGTCCTCGTCGAAGCTTTTCACGATCCAGGTGATGATCATGCCTGCGAAGCcaac                                                                         <  2:3237893/100‑1 (MQ=255)
                             gTTTTCCAGTTTTTCGATTTCTGCCACCGGCACGTAGTAATCCACGTCCTCGTCGAAGCTTTTCACGATCCAGGTGATGATCATGCCTGCGAAGCcaaca                                                                        >  2:152909/1‑100 (MQ=255)
                                  ccAGTTTTTCGATTTCTGCCACCGGCACGTAGTAATCCACGTCCTCGTCGAAGCTTTTCACGATCTAGGTGATGATCATGCCTGCGAAGCCAACAATCgc                                                                   >  2:2683642/1‑100 (MQ=255)
                                   cAGTTTTTCGATTTCTGCCACCGGCACGTAGTAATCCACGTCCTCGTCGAAGCTTTTCACGATCCAGGTGATGATCATGCCTGCGAAGCCAACAATCgcc                                                                  >  2:2000727/1‑100 (MQ=255)
                                     gTTTTTCGATTTCTGCCACCGGCACGTAGTAATCCACGTCCTCGTCGAAGCTTTTCACGATCTAGGTGATGATCATGCCTGCGAAGCCAACAATCgccag                                                                >  2:1750692/1‑100 (MQ=255)
                                     gTTTTTCGATTTCTGCCACCGGCACGTAGTAATCCACGTCCTCGTCGAAGCTTTTCACGATCCAGGTGATGATCATGCCTGCGAAGCCAACAATCgccag                                                                >  2:3457996/1‑100 (MQ=255)
                                     gTTTTTCGATTTCTGCCACCGGCACGTAGTAATCCACGTCCTCGTCGAAGCTTTTCACGATCCAGGTGATGATCATGCCTGCGAAGCCAACAATCgccag                                                                <  2:2550169/100‑1 (MQ=255)
                                      tttttttttttCTTCCACCGGCACGTAGTAATCCACGTCCTCGTCGAAGCTTTTCACGATCCAGGTGATGATCATGCCTGCGAAGCCAACAATCgccagc                                                               >  1:3249748/9‑100 (MQ=255)
                                      tttttCGATTTCTGCCACCGGCACGTAGTAATCCACGTCCTCGTCGAAGCTTTTCACGATCCAGGTGATGATCATGCCTGCGAAGCCAACAATCgccagc                                                               <  1:629362/100‑1 (MQ=255)
                                         ntCGATTTCTGCCACCGGCACGTAGTAATCCACGTCCTCGTCGAAGCTTTTCACGATCCAGGTGATGATCATGCCTGCGAAGCCAACAATCGCCAGccac                                                            >  2:3744363/2‑100 (MQ=255)
                                           cGATTTCTGCCACCGGCACGTAGTAATCCACGTCCTCGTCGAAGCTTTTCACGATCCAGGTGATGATCATGCCTGCGAAGCCAACAATCGCCAGccacca                                                          >  1:79638/1‑100 (MQ=255)
                                            gATTTCTGCCACCTGCACGTAGTAATCCACGTCCTCGTCGAAGCTTTTCACGATCCAGGTGATGATCATGCCTGCGAAGCCAACAATCGCCAGCCACCAg                                                         <  1:1009373/100‑1 (MQ=255)
                                            gATTTCTGCCACCGGCACGTAGTAATCCACGTCCTCGTCGAAGCTTTTCACGATCCAGGTGATGATCATGCCTGCGAAGCCAACAATCGCCAGCCACCAg                                                         >  2:2867096/1‑100 (MQ=255)
                                                  tGCCACCGGCACGTAGTAATCCACGTCCTCGTCGAAGCTTTTCACGATCCAGGTGATGATCATGCCTGCGAAGCCAACAATCGCCAGCCACCAGATATGc                                                   <  2:231819/100‑1 (MQ=255)
                                                   gCCACCGGCACGTAGTAATCCACGTCCTCGTCGAAGCTTTTCACGATCCAGGTGATGATCATGCCTGCGAAGCCAACAATCGCCAGCCACCAGATATGcc                                                  >  1:3732830/1‑100 (MQ=255)
                                                     cACCGGCACGTAGTAATCCACGTCCTCGTCGAAGCTTTTCACGATCCAGGTGATGATCATGCCTGCGAAGCCAACAATCGCCAGCCACCAGATATGCCAg                                                <  1:256510/100‑1 (MQ=255)
                                                      aCCGGCACGTAGTAATCCACGTCCTCGTCGAAGCTTTTCACGATCCAGGTGATGATCATGCCTGCGAAGCCAACAATCGCCAGCCACCAGATATGCCAGa                                               <  2:2212697/100‑1 (MQ=255)
                                                      aCCGGCACGTAGTAATCCACGTCCTCGTCGAAGCTTTTCACGATCCAGGTGATGATCATGCCTGCGAAGCCAACAATCGCCAGCCACCAGATATGCCAGa                                               <  2:3336173/100‑1 (MQ=255)
                                                         ggCACGTAGTAATCCACGTCCTCGTCGAAGCTTTTCACGATCCAGGTGATGATCATGCCTGCGAAGCCAACAATCGCCAGCCACCAGATATGCCAGATCa                                            >  2:551327/1‑100 (MQ=255)
                                                         ggCACGTAGTAATCCACGTCCTCGTCGAAGCTTTTCACGATCCAGGTGATGATCATGCCTGCGAAGCCAACAATCGCCAGCCACCAGATATGCCAGATCa                                            <  1:2427274/100‑1 (MQ=255)
                                                           cACGTAGTAATCCACGTCCTCGTCGAAGCTTTTCACGATCTAGGTGATGATCATGCCTGCGAAGCCAACAATCGCCAGCCACCAGATATGCCCGATCATg                                          >  1:650721/1‑100 (MQ=255)
                                                           cACGTAGTAATCCACGTCCTCGTCGAAGCTTTTCACGATCTAGGTGATGATCATGCCTGCGAAGCCAACAATCGCCAGCCACCAGATATGCCAGATCATg                                          <  1:1912866/100‑1 (MQ=255)
                                                           cACGTAGTAATCCACGTCCTCGTCGAAGCTTTTCACGATCCAGGTGATGCTCATGCCTGCGAAGCCAACAATCGCCAGCCACCAGATATGCCAGATCa                                            <  1:1665125/98‑1 (MQ=255)
                                                           cACGTACTAATCCACGTCCTCGTCGAAGCTGTACACGATCCAGATGATGAACATGCCTGCCAAGCCAACAATCGCCAGCCACCAGATATGCCAGATCATg                                          <  1:1651717/100‑1 (MQ=255)
                                                            aCGTAGTAATCCACGTCCTCGTCGAAGCTTTTCACGATCTAGGTGATGATCATGCCTGCGAAGCCAACAATCGCCAGCCACCAGATATGCCAGATCATgg                                         <  1:875823/100‑1 (MQ=255)
                                                             cGTAGTAATCCACGTCCTCGTCGAAGCTTTTCACGATCCCGGTGATGATCATGCCTGCGCTGTCAACAATCGCCAGCAACCAGATCTGCCCGATCATGGc                                        >  1:3763086/1‑100 (MQ=255)
                                                              gtagtaATCCACGTCCTCGTCGAAGCTTTTCACGATCTAGGTGATGATCATGCCTGCGAAGCCAACAATCGCCAGCCACCAGATATGCCAGATCATGGCg                                       <  1:2268274/100‑1 (MQ=255)
                                                              gtagtaATCCACGTCCTCGTCGAAGCTTTTCACGATCCAGGTGATGATCATGCCTGCGAAGTCAACAATCGCCAGCCACCAGATATGCCAGATCATGGCg                                       <  1:2350336/100‑1 (MQ=255)
                                                              gtagtaATCCACGTCCTCGTCGAAGCTTTTCACGATCCAGGTGATGATCATGCCTGCGAAGCCAACAATCGCCAGCCACCAGATATGCCAGATCATGGCg                                       <  2:2986064/100‑1 (MQ=255)
                                                               tagtaATCCACGTCCTCGTCGAAGCTTTTCACGATCCAGGTGATGATCATGCCTGCGAAGCCAACAATCGCCAGCCACCAGATATGCCAGATCATGGCGa                                      >  2:2729433/1‑100 (MQ=255)
                                                                 gtaATCCACGTCCTCGTCGAAGCTTTTCACGATCCAGGTGATGATCATGCCTGCGAAGTCAACAATCGCCAGCCACCAGATATGCCAGATCATGGCGaaa                                    <  1:2079172/100‑1 (MQ=255)
                                                                 gtaATCCACGTCCTCGTCGAAGCTTTTCACGATCCAGGTGATGATCATGCCTGCGAAGCCAACAATCGCCAGCCACCAGATATGCCAGATCATGGCGaaa                                    >  2:2033238/1‑100 (MQ=255)
                                                                  taATCCACGTCCTCGTCGAAGCTTTTCACGATCCAGGTGATGATCATGCCTGCGAAGCCAACAATCGCCAGCCACCAGATATGCCAGATCATGGCGAAAc                                   >  2:1002002/1‑100 (MQ=255)
                                                                  taATCCACGTCCTCGTCGAAGCTTTTCACGATCCAGGTGATGATCATGCCTGCGAAGCCAACAATCGCCAGCCACCAGATATGCCAGATCATGGCGAAAc                                   <  1:254183/100‑1 (MQ=255)
                                                                    aTCCACGTCCTCGTCGAAGCTTTTCACGATCCAGGTGATGATCATGCCTGCGAAGCCAACAATCGCCAGCCACCAGATATGCCAGATCATGGCGAAACCg                                 >  2:3782142/1‑100 (MQ=255)
                                                                      ccACGTCCTCGTCGAAGCTTTTCACGATCCAGGTGATGATCATGCCTGCGAAGCCAACAATCGCCAGCCACCAGATATGCCAGATCATGGCGAAACCGaa                               >  2:1566594/1‑100 (MQ=255)
                                                                       cACGTCCTCGTCGAAGCTTTTCACGATCCAGGTGATGATCATGCCTGCGAAGACAAGAATCGCCAGAAACAAGATATGACAGCTCATGGCGAACCCGa                                >  2:1649527/1‑98 (MQ=255)
                                                                         cGTCCTCGTCGAAGCTTTTCACGATCCAGGTGATGATCATGCCTGCGAAGCCAACAATCGCCAGCCACCAGATATGCCAGATCATGGCGAAACCGAAGAt                            >  1:492935/1‑100 (MQ=255)
                                                                           tCCTCGTCGAAGCTTTTCACGATCTAGGGGATGATCATGCCTGCGAAGCCAACAATCGCCAGCCACCAGATATGCCAGATCATGGCGAAACCGAAGAtgg                          <  2:45381/100‑1 (MQ=255)
                                                                           tCCTCGTCGAAGCTTTTCACGATCCAGGTGATGATCATGCCTGCGAAGCCAACAATCGCCAGCCACCAGATATGCCAGATCATGGCGAAACCGAAGAtgg                          >  1:3682673/1‑100 (MQ=255)
                                                                            ccTCGTCGAAGCTTTTCACGATCCAGGTGATGATCATGCCTGCGAAGCCAACAATCGCCAGCCACCAGATATGCCAGATCATGGCGAAACCGAAGAtggt                         >  2:749745/1‑100 (MQ=255)
                                                                             cTCGTCGAAGCTTTTCACGATCCAGGTGATGATCATGCCTGCGAAGCCAACAATCGCCAGCCACCAGATATGCCAGATCATGGCGAAACCGAAGAtggtg                        >  1:150066/1‑100 (MQ=255)
                                                                             cTCGTCGAAGCTTTTCACGATCCAGGTGATGATCATGCCTGCGAAGCCAACAATCGCCAGCCACCAGATATGCCAGATCATGGCGAAACCGAAGAtggtg                        >  2:308868/1‑100 (MQ=255)
                                                                              tcgtcgAAGCTTTTCACGATCCAGGTGATGATCATGCCTGCGAAGCCAACAATCGCCAGCCACCAGATATGCCAGATCATGGCGAAACCGAAGAtggtgg                       >  2:1863698/1‑100 (MQ=255)
                                                                               cgtcgAAGCTTTTCACGATCTAGGTGATGATCATGCCTGCGAAGCCAACAATCGCCAGCCACCAGATATGCCAGATCATGGCGAAACCGAAGAtggtgg                       <  1:1518539/99‑1 (MQ=255)
                                                                               cgtcgAAGCTTTTCACGATCTAGGTGATGATCATGCCTGCGAAGCCAACAATCGCCAGCCACCAGATATGCCAGATCATGGCGAAACCGAAGATGGTGga                      >  2:1108114/1‑100 (MQ=255)
                                                                               cgtcgAAGCTTTTCACGATCCAGGTGATGATCATGCCTGCGAAGCCAACAATCGCCAGCCACCAGATATGCCAGATCATGGCGAAACCGAAGATGGTGga                      <  1:2727670/100‑1 (MQ=255)
                                                                                gtcgAAGCTTTTCACGATCCAGGTGATGATCATGCCTGCGAAGCCAACAATCGCCAGCCACCAGATATGCCAGATCATGGCGAAACCGAAGATGGTGgag                     >  1:2786130/1‑100 (MQ=255)
                                                                                gtcgAAGCTTTACCCGAACCAGCAGATGATAATGCCTGAGACGCCCAACATCGCCAGCCAACAGATATGACAGATCATGGCGAAACCGAAggtggttg                       <  1:2612495/98‑3 (MQ=255)
                                                                                  cgAAGCTTTTCACGATCCAGGTGATGATCATGCCTGCGAAGCCAACAATCGCCAGCCACCAGATATGCCAGATCATGGCGAAACCGAAGATGGTGGAGaa                   >  2:1634462/1‑100 (MQ=255)
                                                                                   gAAGCTTTTCACGATCTAGGTGATGATCATGCCTGCGAAGCCAACAATCGCCAGCCACCAGATATGCCAGATCATGGCGAAACCGAAGATGGTGGAGaaa                  <  2:2359926/100‑1 (MQ=255)
                                                                                   gAAGCTTTTCACGATCCAGGTGATGATCATGCCTGCGAAGCCAACAATCGCCAGCCACCAGATATGCCAGATCATGGCGAAACCGAAGATGGTGGAGaaa                  >  2:3651838/1‑100 (MQ=255)
                                                                                   gAAGCTTTTCACGATCCAGGTGATGATCATGCCTGCGAAGCCAACAATCGCCAGCCACCAGATATGCCAGATCATGGCGAAACCGAAGATGGTGGAGaaa                  >  2:718653/1‑100 (MQ=255)
                                                                                    aaGCTTTTCACGATCCAGGTGATGATCATGCCTGCGAAGCCAACAATCGCCAGCCACCAGATATGCCAGATCATGGCGAAACCGAAGATGGTGGAGAAAg                 >  1:2957532/1‑100 (MQ=255)
                                                                                     aGCTTTTCACGATCCAGGTGATGATCATGCCTGCGAAGCCAACAATCGCCAGCCACCAGATATGCCAGATCATGGCGAAACCGAAGATGGTGGAGAAAGc                >  2:3031159/1‑100 (MQ=255)
                                                                                       cTTTTCACGATCCAGGTGATGATCATGCCTGCGAAGCTAACAATCGCCAGCCACCAGTTATGCCAGATCATGGCGAAACCGAAGATGGTGGAGAAAGCt               <  1:2077633/99‑1 (MQ=255)
                                                                                       cTTTTCACGATCCAGGTGATGATCATGCCTGCGAAGCCAACAATCGCCAGCCACCAGATATGCCAGATCATGGCGAAACCGAAGATGGTGGAGAAAGCTg              >  2:1090104/1‑100 (MQ=255)
                                                                                         tttCACGATCCAGGTGATGATCATGCCTGCGAAGCCAACAATCGCCAGCCACCAGATATGCCAGATCATGGCGAAACCGAAGATGGTGGAGAAAGCTGCa            <  2:3721734/100‑1 (MQ=255)
                                                                                         tttCACGATCCAGGTGATGATCATGCCTGCGAAGCCAACAATCGCCAGCCACCAGATATGCCAGATCATGGCGAAACCGAAGATGGTGGAGAAAGCTGCa            <  1:2174928/100‑1 (MQ=255)
                                                                                         tttCACGATCCAGGTGATGATCATGCCTGCGAAGCCAACAATCGCCAGCCACCAGATATGCCAGATCATGGCGAAACCGAAGATGGTGGAGAAAGCTGCa            >  2:1187341/1‑100 (MQ=255)
                                                                                          ttCACGATCCAGGTGATGATCATGCCTGCGAAGCCAACAATCGCCAGCCACCAGATATGCCAGATCATGGCGAAACCGAAGATGGTGGAGAAAGCTGCaa           >  2:3503172/1‑100 (MQ=255)
                                                                                            cACGATCCAGGTGATGATCATGCCTGCGAAGCCAACAATCGCCAGCCACCAGATATGCCAGATCATGGCGAAACCGAAGATGGTGGAGAAAGCTGCAATg         <  2:492941/100‑1 (MQ=255)
                                                                                              cGATCTAGGTGATGATCATGCCTGCGAAGCCAACAATCGCCAGCCACCAGATATGCCAGATCATGGCGAAACCGAAGATGGTGGAGAAAGCTGCAATGAc       <  2:2604042/100‑1 (MQ=255)
                                                                                              cGATCCAGGTGATGATCATGCCTGCGAAGCCAACAATCGCCAGCCACCAGATATGCCAGATCATGGCGAAACCGAAGATGGTGGAGAAAGCTGCAATGAc       <  1:1358110/100‑1 (MQ=255)
                                                                                              cGATCCAGGTGATGATCATGCCTGCGAAGCCAACAATCGCCAGCCACCAGATATGCCAGATCATGGCGAAACCGAAGATGGTGGAGAAAGCTGCAATGAc       <  2:150069/100‑1 (MQ=255)
                                                                                              cGATCCAGGTGATGATCATGCCTGCGAAGCCAACAATCGCCAGCCACCAGATATGCCAGATCATGGCGAAACCGAAGATGGTGGAGAAAGCTGCAATGAc       <  2:1208981/100‑1 (MQ=255)
                                                                                               gATCCAGGTGATGATCATGCCTGCGAAGCCAACAATCGCCAGCCACCAGATATGCCAGATCATGGCGAAACCGAAGATGGTGGAGAAAGCTGCAATGACg      <  1:1051604/100‑1 (MQ=255)
                                                                                                aTCCAGGTGCTGATCATGCCTGCTAAGCCAACAATCGCCAGCCACCAGATATGCCAGATCATGGCGAAACCAAAGATGATGGAGAAAGCTGCAATGACGa     >  1:1644787/1‑100 (MQ=255)
                                                                                                  ccAGGTGATGATCATGCCTGCGAAGCCAACAATCGCCAGCCACCAGATATGCCAGATCATGGCGAAACCGAAGATGGTGGAGAAAGCTGCAATGACGATa   >  2:2915344/1‑100 (MQ=255)
                                                                                                  ccAGGTGATGATCATGCCTGCGAAGCCAACAATCGCCAGCCACCAGATATGCCAGATCATGGCGAAACCGAAGATGGTGGAGAAAGCTGCAATGACGATa   <  2:2057338/100‑1 (MQ=255)
                                                                                                  ccAGGTGATGATCATGCCTGCGAAGCCAACAATCGCCAGCCACCAGATATGCCAGATCATGGCGAAACCGAAGATGGTGGAGAAAGCTGCAATGACGATa   <  2:1023535/100‑1 (MQ=255)
                                                                                                   cAGGTGATGATCATGCCTGCGAAGCCAACAATCGCCAGCCACCAGATATGCCAGATCATGGCGAAACCGAAGATGGTGGAGAAAGCTGCAATGACGATAc  >  1:1448226/1‑100 (MQ=255)
                                                                                                   cAGGTGATGATCATGCCTGCGAAGCCAACAATCGCCAGCCACCAGATATGCCAGATCATGGCGAAACCGAAGATGGTGGAGAAAGCTGCAATGACGATAc  >  2:1396642/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
CCTGCCTTAGTAATCTCATCGAAATGCTGGTTTTCCAGTTTTTCGATTTCTGCCACCGGCACGTAGTAATCCACGTCCTCGTCGAAGCTTTTCACGATCCAGGTGATGATCATGCCTGCGAAGCCAACAATCGCCAGCCACCAGATATGCCAGATCATGGCGAAACCGAAGATGGTGGAGAAAGCTGCAATGACGATAC  >  NC_000913/448669‑448867

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