Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,561,348 A→G C49C (TGT→TGC)  b4563 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,561,3480AG93.8% 1074.8 / 31.0 392C49C (TGT→TGC) b4563hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (14/10);  new base G (165/200);  total (181/211)
Fisher's exact test for biased strand distribution p-value = 2.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.36e-01

Only 100 of 410 total aligned reads displayed.

ATGTGCAATGCAGCGCCAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCC‑GGGGGGACAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTCATGCCCGCTGCCTGCAACGCTTT  >  NC_000913/4561258‑4561440
                                                                                           |                                                                                            
aTGTGCAATGCAGCGCCAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCC‑GGGGGGGCa                                                                                            <  1:4601948/93‑1 (MQ=255)
 tgtgCAATGCAGCGCCAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCC‑GGGGGGGCAg                                                                                           <  1:3722273/93‑1 (MQ=255)
   tgCAATGCAGCGCCAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCC‑GGGGGGGCAGGt                                                                                         <  1:5842133/93‑1 (MQ=255)
   tgCAATGCAGCGCCAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCC‑GGGGGGGCAGGt                                                                                         <  1:4340039/93‑1 (MQ=255)
   tgCAATGCAGCGCCAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCC‑GGGGGGGCAGGt                                                                                         <  2:2802433/93‑1 (MQ=255)
       aTGCAGCCCCAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGGGGCAGGTCAg                                                                                      <  2:8831923/93‑1 (MQ=255)
        tGCAGCGCCAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCC‑GGGGGGGCAGGTCAGcc                                                                                    >  1:4464864/1‑93 (MQ=255)
        tGCAGCGCCAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCC‑GGGGGGGCAGGTCAGcc                                                                                    >  1:8574795/1‑93 (MQ=255)
          cAGCGCCAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGTATCACCCATTGTATTTCCTCGCC‑GGGGGGGCAGGTCAGCCaa                                                                                  <  2:809289/93‑1 (MQ=255)
           aGCGCCAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGGGGCAGGTCAGCCaa                                                                                  >  2:3796817/1‑93 (MQ=255)
           aGCGCCAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCC‑GGGGGGGCAGGTCAGCCaaa                                                                                 >  2:9812434/1‑93 (MQ=255)
           aGCGCCAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCC‑GGGGGGGCAGGTCAGCCaaa                                                                                 <  1:2618481/93‑1 (MQ=255)
            gcgcCAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCC‑GGGGGGGCAGGTCAGCCAAAt                                                                                <  1:8076545/93‑1 (MQ=255)
            gcgcCAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCC‑GGGGGGACAGGTCAGCCAAAt                                                                                <  2:4311055/93‑1 (MQ=255)
               ccAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCC‑GGGGGGGCAGGTCAGCCAAATGCa                                                                             <  2:7372326/93‑1 (MQ=255)
               ccAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCC‑GGGGGGGCAGGTCAGCCAAATGCa                                                                             <  1:20997/93‑1 (MQ=255)
                cAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGGGGCAGGTCAGCCAAATGCa                                                                             >  2:1993181/1‑93 (MQ=255)
                cAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCC‑GGGGGGGCAGGTCAGCCAAATGCAt                                                                            <  2:2490041/93‑1 (MQ=255)
                cAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCC‑GGGGGGGCAGGTCAGCCAAATGCAt                                                                            >  2:9564238/1‑93 (MQ=255)
                cAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCC‑GGGGGGGCAGGTCAGCCAAATGCAt                                                                            >  1:5157795/1‑93 (MQ=255)
                cAGACACCATCATGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCAATTGTATTTCCTCGCC‑GGGGGGGCAGGTCAGCCAAATGCAt                                                                            >  2:4644867/1‑93 (MQ=255)
                 cgACACCATCCTGCTGGCAATAAATTGACTCAGTCAGGCCAGGGATCACCCATTGTATTTCCTCGCC‑GGGGGGACAGGTCAGCCAAATGCATa                                                                           <  1:3393218/92‑1 (MQ=255)
                 aGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCC‑GGGGGGGCAGGTCAGCCAAATGCATa                                                                           <  2:7739490/93‑1 (MQ=255)
                  gACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCC‑GGGGGGGCAGGTCAGCCAAATGCATaa                                                                          >  1:8514427/1‑93 (MQ=255)
                  gACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCC‑GGGGGGGCAGGTCAGCCAAATGCATaa                                                                          <  1:5110328/93‑1 (MQ=255)
                   acacCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCTCC‑GGGGGGGCAGGTCAGCCAAATGCATAAc                                                                         <  2:4115471/93‑1 (MQ=255)
                   acacCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCC‑GGGGGGGCAGGTCAGCCAAATGCATAAc                                                                         >  2:4008287/1‑93 (MQ=255)
                    cacCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCC‑GGGGGGGCAGGTCAGCCAAATGCATAAcc                                                                        <  1:965032/93‑1 (MQ=255)
                     acCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCCGGGGGGGGCAGGTCAGCCAAATGCATAAcc                                                                        <  2:4498547/93‑1 (MQ=255)
                      ccATCCTGCTGGTAATAAATTATCTCAGTCGGCCCAGGTATCCCCCATTGTTTTTCCTCGCA‑GGGGAGTCAGGTCAGCCAAATGCATAACCCg                                                                      <  2:7454421/93‑1 (MQ=255)
                      ccATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCC‑GGGGGGGCAGGTCAGCCAAATGCATAACCCg                                                                      >  2:5350117/1‑93 (MQ=255)
                      ccATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCC‑GGGGGG‑CAGGTCAGCCTAATGCATAACCCGc                                                                     >  2:4451135/1‑93 (MQ=255)
                       cATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCC‑GGGGGGGCAGGTCAGCCAAATGCATAACCCGc                                                                     <  1:4900210/93‑1 (MQ=255)
                       cATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCC‑GGGGGGGCAGGTCAGCCAAATGCATAACCCGc                                                                     >  1:70670/1‑93 (MQ=255)
                        aTCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCC‑GGGGGGGCAGGTCAGCCAAATGCATAACCCGCa                                                                    <  2:6618114/93‑1 (MQ=255)
                         tCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCC‑GGGGGGGCAGGTCAGCCAAATGCATAACCCGCAg                                                                   <  2:1321051/93‑1 (MQ=255)
                          ccTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCC‑GGGGGGGCAGGTCAGCCAAATGCATAACCCGCAgc                                                                  <  1:6173750/93‑1 (MQ=255)
                           ctgctgGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCC‑GGGGGGGCAGGTCAGCCAAATGCATAACCCGCAgcc                                                                 >  1:3503386/1‑93 (MQ=255)
                             gctgGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCC‑GGGGGGGCAGGTCAGCCAAATGCATAACCCGCAgccgc                                                               >  1:3964903/1‑93 (MQ=255)
                             gctgGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCC‑GGGGGGGCAGGTCAGCCAAATGCATAACCCGCAgccgc                                                               >  1:5998218/1‑93 (MQ=255)
                              gtgGCAATAAATTGACTCAGTCAGCCCAGGGATCCCCCATTGTATTTCCTCGCC‑GGGGGGGCAGGTCAGCCAAATGCATAACCCGCAgccgcc                                                              <  1:2282066/92‑1 (MQ=255)
                              ctgGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCC‑GGGGGGGCAGGTCAGCCAAATGCATAACCCGCAgccgcc                                                              <  2:1811976/93‑1 (MQ=255)
                                ggCAATAAATTGACTCACTCACCCCAGCGATCACCAATTGTATTGCCTCGCC‑GGGGGGACAGGTCAGCCAAATGCATAACCCGCAGCcggcgc                                                            <  1:9177759/93‑1 (MQ=255)
                                 gCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCC‑GGGGGG‑CAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGt                                                          <  2:2949223/93‑1 (MQ=255)
                                     tAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCC‑GGGGGGGCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTaaa                                                       <  2:5699865/93‑1 (MQ=255)
                                     tAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCC‑GGGGGGGCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTaaa                                                       >  2:7248382/1‑93 (MQ=255)
                                      aaaTTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCC‑GGGGGGGCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGAGTAAAt                                                      >  2:10014099/1‑93 (MQ=255)
                                       aaTTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCC‑GGGGGGGCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATc                                                     <  1:1374730/93‑1 (MQ=255)
                                       aaTTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCC‑GGGGGGGCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATc                                                     <  1:2419528/93‑1 (MQ=255)
                                        aTTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCC‑GGGGGGGCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCa                                                    <  2:5837638/93‑1 (MQ=255)
                                         ttGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCC‑GGGGGGGCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCAc                                                   <  2:6741758/93‑1 (MQ=255)
                                         gtGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCC‑GGGGGGGCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCAc                                                   <  1:7405149/92‑1 (MQ=255)
                                           gACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCC‑GGGGGGGCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACgg                                                 <  2:151325/93‑1 (MQ=255)
                                            gcTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCC‑GGGGGGGCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACggg                                                <  2:6109342/92‑1 (MQ=255)
                                             cTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCC‑GGGGGGGCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACgggg                                               >  2:4781060/1‑93 (MQ=255)
                                               cagtcagCCCAGGGATCACCCATTGTATTTCCTCGCC‑GGGGGGGCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGaa                                             <  2:8258850/93‑1 (MQ=255)
                                                agtcagCCCAGGGATCACCCATTGTATTTCCTCGCC‑GGGGGGGCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAat                                            >  2:2201288/1‑93 (MQ=255)
                                                agtcagCCCAGGGATCACCCATTGTATTTCCTCGCC‑GGGGGGGCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAat                                            >  2:4914929/1‑93 (MQ=255)
                                                agtcagCCCAGGGATCACCCATTGTATTTCCTCGCC‑GGGGGGGCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAat                                            >  1:1445286/1‑93 (MQ=255)
                                                  tcagCCCAGGGATCACCCATTGTATTTCCTCGCC‑GGGGGGGCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAatat                                          <  1:1959361/93‑1 (MQ=255)
                                                  tcagCCCAGGGATCACCCATTGTATTTCCTCGCC‑GGGGGGGCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAatat                                          >  2:8633844/1‑93 (MQ=255)
                                                   cagCCCAGGGATCACCCATTGTATTTCCTCGCC‑GGGGGGGCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATc                                         >  2:9346214/1‑93 (MQ=255)
                                                     gCCCAGGGATCACCCATTGTATTTCCTCGCC‑GGGGGGGCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGc                                       >  1:3881045/1‑93 (MQ=255)
                                                      cccAGGGATCACCCATTGTATTTCCTCGCC‑GGGGGGGCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCt                                      <  2:9448050/93‑1 (MQ=255)
                                                       ccAGGGATCACCCATTGTATTTCCTCGCC‑GGGGGGGCAGGTCAGCCAAATGCCTAACCCGCAGCCGCCGCGTAAATCACGGGCAATATCGCTg                                     >  1:2890545/1‑93 (MQ=255)
                                                        cAGGGATCACCCATTGTATTTCCTCGCC‑GGGGGGGCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTgc                                    >  1:1028136/1‑92 (MQ=255)
                                                        cAGGGATCACCCATTGTATTTCCTCGCC‑GGGGGGGCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGa                                    >  2:5629490/1‑93 (MQ=255)
                                                         aGGGATCACCCATTGTATTTCCTCGCC‑GGGGGGGCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGAc                                   <  2:6485474/93‑1 (MQ=255)
                                                         aGGGATCACCCATTGTATTTCCTCGCC‑GGGGGGGCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGAc                                   >  1:408722/1‑93 (MQ=255)
                                                          gggATCACCCATTGTATTTCCTCGCC‑GGGGGGACAGGTGAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGAcg                                  >  2:6249533/1‑93 (MQ=255)
                                                           ggATCACCCATTGTATTTCCTCGCC‑GGGGGGGCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGAcgc                                 <  1:9525907/93‑1 (MQ=255)
                                                              tCACCCATTGTATTTCCTCGCC‑GGGGGGGCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGaa                              >  2:6820517/1‑93 (MQ=255)
                                                               cACCCATTGTATTTCCTCGCC‑GGGGGGGCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGaaa                             <  2:3289189/93‑1 (MQ=255)
                                                                tcccATTGTATTTCCTCGCC‑GGGGGGGCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGaaaa                            <  2:393654/92‑1 (MQ=255)
                                                                 cccATTGTATTTCCTCGCC‑GGGGGGGCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAg                           >  1:1726544/1‑93 (MQ=255)
                                                                 cccATTGTATTTCCTCGCC‑GGGGGGGCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAg                           <  2:3506611/93‑1 (MQ=255)
                                                                   cATTGTATTTCCTCGCC‑GGGGGGGCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGTGAATATCGCTGACGCGAAAAGTc                         >  2:5546599/1‑93 (MQ=255)
                                                                    aTTGTATTTCCTCGCC‑GGGGGGGCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTCa                        <  2:9911657/93‑1 (MQ=255)
                                                                    aTTGTATTTCCTCGCC‑GGGGGGGCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTCa                        <  1:4998757/93‑1 (MQ=255)
                                                                    aTTGTATTTCCTCGCC‑GGGGGGGCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTCa                        <  1:2839631/93‑1 (MQ=255)
                                                                    aTTGTATTTCCTCGCC‑GGGGGGGCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTCa                        >  1:7611066/1‑93 (MQ=255)
                                                                     ttGTATTTCCTCGCC‑GGGGGGGCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTCAt                       >  1:728606/1‑93 (MQ=255)
                                                                     ttGTATTTCCTCGCC‑GGGGGGGCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTCAt                       >  1:7762800/1‑93 (MQ=255)
                                                                        tATTTCCTCGCC‑GGGGGGGCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGAGGAATATCGCTGACGCGAAAAGTCATGcc                    >  1:1513851/1‑93 (MQ=255)
                                                                          tttCCTCGCC‑GGGGGGGCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTCATGCCCg                  <  2:2993233/93‑1 (MQ=255)
                                                                               tCGCC‑GGGGGGGCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTCATGCCCGctgcc             >  1:551758/1‑93 (MQ=255)
                                                                               tCGCC‑GGGGGGGCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTCATGCCCGctgcc             >  1:9898046/1‑93 (MQ=255)
                                                                                cGCC‑GGGGGGGCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTCATGCCCGctgcct            <  2:834328/93‑1 (MQ=255)
                                                                                cGCC‑GGGGGGGCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTCATGCCCGctgcct            >  1:4225560/1‑93 (MQ=255)
                                                                                 gCC‑GGGGGGGCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTCATGCCCGctgcctg           >  1:6521199/1‑93 (MQ=255)
                                                                                 gCC‑GGGGGGGCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTCATGCCCGctgcctg           >  1:8234701/1‑93 (MQ=255)
                                                                                  cc‑GGGGGGGCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTCATGCCCGctgcctgc          >  2:1249245/1‑93 (MQ=255)
                                                                                   cgggggggGCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTCATGCCCGctgcctgc          >  1:6861732/3‑93 (MQ=255)
                                                                                     ggggggGCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTCATGCCCGCTGCCTGCac        >  1:7144929/1‑92 (MQ=255)
                                                                                     ggggggGCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTCATGCCCGCTGCCTGCaa        <  1:9389383/93‑1 (MQ=255)
                                                                                      ggggggCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTCATGCCCGCTGCCTGCAAc       >  1:9004589/1‑93 (MQ=255)
                                                                                       ggggGCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTCATGCCCGCTGCCTGCAACg      <  1:5606345/93‑1 (MQ=255)
                                                                                       ggggGCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTCATGCCCGCTGCCTGCAACg      >  1:8849400/1‑93 (MQ=255)
                                                                                        ggggCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTCATGCCCGCTGCCTGCAACGc     <  1:9229741/93‑1 (MQ=255)
                                                                                           aCAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTCATGCCCGCTGCCTGCAACGCttt  >  1:3070443/1‑93 (MQ=255)
                                                                                           |                                                                                            
ATGTGCAATGCAGCGCCAGACACCATCCTGCTGGCAATAAATTGACTCAGTCAGCCCAGGGATCACCCATTGTATTTCCTCGCC‑GGGGGGACAGGTCAGCCAAATGCATAACCCGCAGCCGCCGCGTAAATCACGGGGAATATCGCTGACGCGAAAAGTCATGCCCGCTGCCTGCAACGCTTT  >  NC_000913/4561258‑4561440

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 19 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑