Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,337,024 G→A 100% S3L (TCG→TTG)  adiC ← arginine:agmatine antiporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,337,0240GA99.1% 935.6 / ‑3.2 264S3L (TCG→TTG) adiCarginine:agmatine antiporter
Reads supporting (aligned to +/- strand):  ref base G (1/2);  new base A (137/124);  total (138/126)
Fisher's exact test for biased strand distribution p-value = 6.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 267 total aligned reads displayed.

CAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCGG  >  NC_000913/4336927‑4337119
                                                                                                 |                                                                                               
cAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCaaa                                                                                               <  1:3349508/100‑1 (MQ=255)
 aGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCaaa                                                                                               <  1:2921065/99‑1 (MQ=255)
   gTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGAc                                                                                            <  1:1409201/100‑1 (MQ=255)
       gCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGTCCACTTTGTGAGCATCAGCATCCAAAGACATaa                                                                                        >  2:2032785/1‑100 (MQ=255)
        cAGGTAACAGAAAAACACCTGTCCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAAt                                                                                       <  2:1378343/100‑1 (MQ=255)
         aGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAAtt                                                                                      <  1:1144435/100‑1 (MQ=255)
         aGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAAtt                                                                                      <  2:2771895/100‑1 (MQ=255)
         aGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAAtt                                                                                      <  1:3405520/100‑1 (MQ=255)
          ggTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTa                                                                                     >  2:1919377/1‑100 (MQ=255)
          ggTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTa                                                                                     >  2:1231921/1‑100 (MQ=255)
          ggTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTa                                                                                     <  1:1673863/100‑1 (MQ=255)
          ggTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTa                                                                                     >  2:146226/1‑100 (MQ=255)
             aaCAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAAc                                                                                  <  2:908953/100‑1 (MQ=255)
             aaCAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAAc                                                                                  <  2:1744561/100‑1 (MQ=255)
              aCAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACt                                                                                 >  1:2409005/1‑100 (MQ=255)
               cAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTc                                                                                >  2:1645035/1‑100 (MQ=255)
               cAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTc                                                                                >  1:1626237/1‑100 (MQ=255)
                aGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTcc                                                                               >  2:2706307/1‑100 (MQ=255)
                 gAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCt                                                                              >  2:1204876/1‑100 (MQ=255)
                 gAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCt                                                                              >  2:350071/1‑100 (MQ=255)
                  aaaaaCACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTg                                                                             <  1:3211362/100‑1 (MQ=255)
                   aaaaCACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGc                                                                            >  2:3320243/1‑100 (MQ=255)
                    aaaCACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCg                                                                           <  2:1942364/100‑1 (MQ=255)
                    aaaCACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCg                                                                           <  2:553392/100‑1 (MQ=255)
                        acCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAgg                                                                       <  1:614197/100‑1 (MQ=255)
                        acCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAgg                                                                       <  1:2622186/100‑1 (MQ=255)
                           tGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGa                                                                    >  1:22387/1‑100 (MQ=255)
                           tGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGa                                                                    >  2:1763641/1‑100 (MQ=255)
                            gACCCCATAATATTCCCCGAGACCATCAGGGTGACGGGGATTAAGCCCACGTTGTGGGCATCAGCAGCCAAAGACTTAATTGAACTCCGGCGAAGGCGAg                                                                   <  1:565536/100‑1 (MQ=255)
                              ccccATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCt                                                                 >  1:326142/1‑100 (MQ=255)
                               cccATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCtt                                                                >  1:3029363/1‑100 (MQ=255)
                                  ataataTTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCTTCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGt                                                             <  1:1919361/100‑1 (MQ=255)
                                  ataataTTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGt                                                             <  2:3084915/100‑1 (MQ=255)
                                   taataTTCCCCGACACCATCAGGGTGACGGGGATTAAGTCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTg                                                            <  1:1882975/100‑1 (MQ=255)
                                   taataTTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTg                                                            <  2:2299883/100‑1 (MQ=255)
                                   taataTTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTg                                                            <  1:2460761/100‑1 (MQ=255)
                                   taataTTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTg                                                            <  1:180565/100‑1 (MQ=255)
                                    aataTTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGa                                                           <  2:649901/100‑1 (MQ=255)
                                     atatTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGAc                                                          >  1:719141/1‑100 (MQ=255)
                                     atatTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGAc                                                          <  2:2769115/100‑1 (MQ=255)
                                     atatTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGAc                                                          <  1:100531/100‑1 (MQ=255)
                                      tatTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACa                                                         <  1:1391502/100‑1 (MQ=255)
                                        ttCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAAt                                                       >  1:2799358/1‑100 (MQ=255)
                                        ttCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAAt                                                       <  2:1061653/100‑1 (MQ=255)
                                        ttCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAAt                                                       >  1:1503691/1‑100 (MQ=255)
                                              gACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTa                                                 <  2:2599282/100‑1 (MQ=255)
                                              gACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTa                                                 <  2:326145/100‑1 (MQ=255)
                                              gACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTa                                                 <  2:2936207/100‑1 (MQ=255)
                                               acacCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTAc                                                <  1:84488/100‑1 (MQ=255)
                                               acacCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTAc                                                <  1:3310492/100‑1 (MQ=255)
                                                cacCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACg                                               <  1:538352/100‑1 (MQ=255)
                                                cacCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACg                                               <  1:3416218/100‑1 (MQ=255)
                                                cacCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACg                                               <  1:2525232/100‑1 (MQ=255)
                                                  ccATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGtc                                             >  1:2289497/1‑100 (MQ=255)
                                                   cATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGtca                                            >  1:810193/1‑100 (MQ=255)
                                                   cATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGtca                                            >  2:3267033/1‑100 (MQ=255)
                                                    aTCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGtcat                                           <  2:752586/100‑1 (MQ=255)
                                                    aTCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGtcat                                           >  1:209520/1‑100 (MQ=255)
                                                     tCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGtcata                                          >  2:2132004/1‑100 (MQ=255)
                                                     tCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGtcata                                          <  1:3294784/100‑1 (MQ=255)
                                                        gggTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGtcatatca                                       >  2:600678/1‑100 (MQ=255)
                                                          gTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGtcatatcata                                     <  1:1231910/100‑1 (MQ=255)
                                                          gTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGtcatatcata                                     <  1:427412/100‑1 (MQ=255)
                                                          gTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGtcatatcata                                     >  1:879664/1‑100 (MQ=255)
                                                             aCGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGc                                  >  2:3358611/1‑100 (MQ=255)
                                                              cGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGcc                                 <  2:1948096/100‑1 (MQ=255)
                                                              cGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGcc                                 >  2:41602/1‑100 (MQ=255)
                                                              cGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGcc                                 >  1:2452550/1‑100 (MQ=255)
                                                              cGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGcc                                 >  2:2399366/1‑100 (MQ=255)
                                                              cGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGcc                                 >  2:1842317/1‑100 (MQ=255)
                                                              cGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGcc                                 <  2:1445197/100‑1 (MQ=255)
                                                               ngggATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCt                                >  1:2569433/2‑100 (MQ=255)
                                                               ggggATTAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCt                                >  2:1524423/1‑100 (MQ=255)
                                                                    ttAAGTCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCAtt                           <  1:2032768/100‑1 (MQ=255)
                                                                     tAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTc                          <  2:831767/100‑1 (MQ=255)
                                                                     tAAGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTAACAATAAATTACGTCATATCATACGCCTGCAtt                           <  1:892340/99‑1 (MQ=255)
                                                                      aggCCCACTTTGTAAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCg                         >  2:2343498/3‑100 (MQ=255)
                                                                      aaGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCg                         >  2:2865724/1‑100 (MQ=255)
                                                                       aGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGa                        <  1:1842301/100‑1 (MQ=255)
                                                                       aGCCCACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGa                        >  2:1975198/1‑100 (MQ=255)
                                                                          ccACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAg                     <  1:2172755/100‑1 (MQ=255)
                                                                          ccACTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAg                     <  1:449525/100‑1 (MQ=255)
                                                                             cTTTGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAGATTACGTCATATCATACGCCTGCATTCGATAGTTg                  <  1:3223932/100‑1 (MQ=255)
                                                                               ttGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGaa                >  1:2197369/1‑100 (MQ=255)
                                                                               ttGTGAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGaa                <  1:1216629/100‑1 (MQ=255)
                                                                                 gtgAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAatt              >  2:2576000/1‑100 (MQ=255)
                                                                                 gtgAGCATCAGCATCCAAAGACATAATCAAACTCCTGCGAAGGCGAGCTGCGAGACAATAAATTACGACATATCATACGCCTGCATTCGATAGTTGaact              >  1:2309947/1‑98 (MQ=255)
                                                                                   gAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAattat            >  1:1439350/1‑100 (MQ=255)
                                                                                   gAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAattat            <  1:3459210/100‑1 (MQ=255)
                                                                                   gAGCATCAGCATCCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAattat            >  1:2985382/1‑100 (MQ=255)
                                                                                     gcatcagcatcCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAAttatta          >  2:787323/1‑100 (MQ=255)
                                                                                     gcatcagcatcCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAAttatta          >  1:700402/1‑100 (MQ=255)
                                                                                      catcagcatcCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCCTTCGATAGTTGAAttctta          >  1:2741554/1‑99 (MQ=255)
                                                                                      catcagcatcCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATtattat         <  2:2210604/100‑1 (MQ=255)
                                                                                         cagcatcCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGa      >  1:2372621/1‑100 (MQ=255)
                                                                                          agcatcCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGAc     >  2:1882376/1‑100 (MQ=255)
                                                                                           gcatcCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGAcc    >  2:315614/1‑100 (MQ=255)
                                                                                           gcatcCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGAcc    <  1:1075253/100‑1 (MQ=255)
                                                                                            catcCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGAcc    >  1:2782394/1‑99 (MQ=255)
                                                                                             atcCAAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCgg  <  1:3073040/100‑1 (MQ=255)
                                                                                                 |                                                                                               
CAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCGG  >  NC_000913/4336927‑4337119

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