Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 219,763 T→C 5.6% N78S (AAT→AGT)  rcsF ← putative outer membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913219,7630TC 5.6% 224.1 / 2.1 72N78S (AAT→AGT) rcsFputative outer membrane protein
Reads supporting (aligned to +/- strand):  ref base T (36/31);  new base C (2/2);  total (38/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AACCGATACATACAGCCTGACGATAGCAGCCTGGCGTACCGCTGGTGACTTCGCAGCTATGCAGTAATACAGCATTGGCTTTCATTTTAGAGGCGTTGATTTGCATCCGCTTACGTGCGGTTGGAATGCTCGGCGGAGAGTCCTGATTAGAGGCCTGGCAAGAGTCGCCACTGACTTCACCGAGATCGCGGAACGGTTTGCCGACTAATTCTTCTGCATTGGTATAAATTCGGACCGGCGTGGCGCGCGGCGCTTTCGGTTTTGCAGGCTCCGCTTTCGGCTGGGGTGCAGT  >  NC_000913/219617‑219908
                                                                                                                                                  |                                                                                                                                                 
aaCCGATACATACAGCCTGACGATAGCAGCCTGGCGTACCGCTGGTGACTTCGCAGCTATGCAGTAATACAGCATTGGGTTTGATTTTAGAGGCGTTGATTTGCCTCCGCTAAGGTGGGTTTGGATTGCTCGGGTGAGAGTCGTGAt                                                                                                                                                   >  4:342390/1‑147 (MQ=255)
 aCCGATACATACAGCCTGACGATAGCAGCCTGGCGTACCGCTGGTGACTTCGCAGCTATGCAGTAATACAGCATTGGCTTTCATTTTAGAGGCGTTGATTTGCATCCGCTTACGTGCGGTTGGAATGCTCGGCGGAGAGTCCTGAt                                                                                                                                                   >  7:406574/1‑146 (MQ=255)
  ccGATACATACAGCCTGACGATAGCAGCCTGGCGTACCGCTGGTGACTTCGCAGCTATGCAGTAATACAGCATTGGCTTTCATTTTAGAGGCGTTGATTTGCATCCGCTTACGTGCGGTTGGAATGCTCGGCGGAGAGTCCTGATTag                                                                                                                                                <  4:180409/148‑1 (MQ=255)
   cGATACATACAGCCTGACGATAGCAGCCTGGCGTACCGCTGGTGACTTCGCAGCTATGCAGTAATACAGCATTGGCTTTCATTTTAGAGGCGTTGATTTGCATCCGCTTACGTGCGGTTGGAATGCTCGGCGGAGAGTCCTGATTaga                                                                                                                                               <  5:277312/148‑1 (MQ=255)
    gATACATACAGCCTGAAGAGAGCAGCCTGGCGTACCGCTGGTGACTTCGCAGCTATGCAGTAATACAGCATTGGCTTTCATTTTAGAGGCGTAGATTTGCATCCGCATACGTGCGGATGGAATGCTCGGCGGAGAGTCCTGATTagag                                                                                                                                              >  4:5206/1‑148 (MQ=255)
      tacatacaGCCTGCCGATAGCAGCCTGGCGTACCGCTGGTGACTTCGCAGCTATGCAGTAATACAGCATTGGCTTTCATTTTAGAGGCGTTGATTTGCATCCGCTTACGTGCGGTTGGAATGCTCGGCGGAGAGTCCTGATTAGAGGc                                                                                                                                            <  6:264567/148‑1 (MQ=255)
                   aagctAGCAGCCTGCCCGCCCGCCGGGGACTTCGCAGATAGGCAGAAATACAGCATGGGCTTTCATTATAGAGGCGGTGATTAGCCTCCGCTAACCGGCGGGTGGACTGCTCGGCGGAGAGTCCTGATTAGAGGCCGGGCAAGAGt                                                                                                                                 <  4:258689/142‑1 (MQ=255)
                   aCGATAGCAGCCTGGCGTACCGCTGGTGACTTCGCAGCTATGCAGTAATACAGCATTGGCTTTCATTTTAGAGGCGTTGATTTGCATCCGCTTACGTGCGGTTGGAATGCTCGGCGGAGAGTCCTGATTAGAGGCCTGGCAAGAGTCg                                                                                                                               <  8:406574/148‑1 (MQ=255)
                   aCGATAGCAGCCTGGCGTACCGCTGGTGACTTCGCAGCTATGCAGTAATACAGCATTGGCTTTCATTTTAGAGGCGTTGATTTGCATCCGCTTACGTGCGGTTGGAATGCTCGGCGGAGAGTCCTGATTAGAGGCCTGGCAAGAGTCg                                                                                                                               >  1:199543/1‑148 (MQ=255)
                    cGATAGCAGCCTGGCGTACCGCTGGTGACTTCGCAGCTATGCAGTAATACAGCATTGGCTTTCATTTTAGAGGCGTTGATTTGCATCCGCTTACGTGCGGTTGGAATGCTCGGCGGAGAGTCCTGATTAGAGCCCTGGCAAGAGTCGc                                                                                                                              >  2:22282/1‑148 (MQ=255)
                     gATAGCAGCCTGGCGTACCGCTGGTGACTTCGCAGCTATGCTGTAATACAGCATTGGCTTTCATTTTAGAGGCGTTGATTTGCATCCGCTTACGTGCGGTTGGAATGCTCGGCGGAGAGTCCTGATTAGAGGCCTGGCAAGAGTCGcc                                                                                                                             <  2:252081/148‑1 (MQ=255)
                         gcagcCTGGCGTACCGCTGGTGACTTCGCAGCTATGCAGTAATACAGCATTGGCTTTCATTTTAGAGGCGTTGATTTGCATCCGCTTACGTGCGGTTGGAATGCTCGGCGGAGAGTCCTGATTAGAGGCCTGGCAAGAGTCGCCACTg                                                                                                                         <  5:244707/148‑1 (MQ=255)
                            gcCTGGCGTACCGCTGGTGACTTCGCAGCTATGCAGTAATACAGCATTGGCTTTCATTTTAGAGGCGTTGATTTGCATCCGCTTACGTGCGGTTGGAATGCTCGGCGGAGAGTCCTGATTAGAGGCCTGGCAAGAGTCGCCACTGACt                                                                                                                      >  5:112869/1‑148 (MQ=255)
                              cTGGCGTACCGCTGGTGACTTCGCAGCTATGCAGTAATACAGCATTGGCTTTCATTTTAGAGGCGTTGATTTGCATCCGCTTACGTGCGGTTGGAATGCTCGGCGGAGAGTCCTGATTAGAGGCCTGGCAAGAGTCGCCACTGACtt                                                                                                                     >  7:44715/1‑147 (MQ=255)
                                        gCTGGTGACTTCGCAGCTATGCAGTAATACAGCATTGGCTTTCATTTTAGAGGCGTTGATTTGCATCCGCTTACGTGCGGTTGGAATGCTCGGCGGAGAGTCCTGATTAGAGGCCTTGCAAGAGTCGCCACTGACTTCACCGAGATcg                                                                                                          >  3:338840/1‑148 (MQ=255)
                                          tGGTGACTTCGCAGCTATGCAGTAATACAGCATTGGCTTTCATTTTAGAGGCGTTGATTTGCATCCGCTTACGTGCGGTTGGAATGCTCGGCGGAGAGTCCTGACTAGAGGCCTGGCAAGAGTCGCCACTGACTTCACCGAGATcgcg                                                                                                        >  2:111856/1‑148 (MQ=255)
                                          tGGTGACTTCGCAGCTATGAAGTAATACAGCATTGGCTTTCATTTTAGAGTCGTTGATTTGCATACGCTTACGTGCGGTTGGAATGCTCGGCGGAGAGTCCTGACTAGAGGCCTGGCAAGAGTCGCCACTGACTTCACCGAGATcgcg                                                                                                        >  2:105911/1‑148 (MQ=255)
                                          tGGTCACTTCGCAGCAATGCAGTAATACAGCATTGGCTTTCATTTTAGAGGCGTTGATTTGCATCCGCATACGTGCGGTTGGAATGATCGGCTGAGAGTCCTGATTATAGGCATGGCAATAGTCTCCACTGACTTCACCGAGATcgcg                                                                                                        >  4:360908/1‑148 (MQ=255)
                                           ggTGACTTCGCAGCTATGCAGTAATACAGCATTGGCTTTCATTTTAGAGGCGTTGATTTGCATCCGCTTACGTGCGGTTGGAATGCTCGGCGGAGAGTCCTGATTAGAGGCCTGGCAAGAGTCGCCACTGACTTCACCGAGATcgcg                                                                                                        >  5:130983/1‑147 (MQ=255)
                                               aCTTCGCAGCTATGCAGTAATACAGCATTGGCTTTCATTTTAGAGGCGTTGATTTGCATCCGCTTACGTGCGGTTGGAATGCTCGGCGGAGAGTCCTGATTAGAGGCCTGGCAAGAGTCGCCACTGACTTCACCGAGATCGCGGAAc                                                                                                    <  6:112869/147‑1 (MQ=255)
                                                  tCGCAGCTATGCAGTAATACAGCATTGGCTTTCATTTTAGAGGCGTTGATTTGCATCCGCTTACGTGCGGTTGGAATGCTCGGCGGAGAGTCCTGATTAGAGGCCTGGCAAGAGTCGCCACTGACTTCACCGAGATCGCGGAACGGtt                                                                                                >  3:157741/1‑148 (MQ=255)
                                                  tCGCAGCTATGCAGTAATACAGCATTGGCTTTCATTTTAGAGGCGTTGATTTGCATCCGCTTACGTGCGGTTGGAATGCTCGGCGGAGAGTCCTGATTAGAGGCCTGGCAAGAGTCGCCACTGACTTCACCGAGATCGCGGAACGGtt                                                                                                >  8:321847/1‑148 (MQ=255)
                                                   cGCAGCTATGCAGTAATACAGCATTGGCTTTCATTTTAGAGGCGTTGATTTGCATCCGCTTACGTGCGGTTGGAATGCTCGGCGGAGAGTCCTGATTAGAGGCCTGGCAAGAGTCGCCACTGACTTCACCGAGATCGCGGAACGGt                                                                                                 <  4:338840/146‑1 (MQ=255)
                                                      aGCTATGCAGTAATTCAGCATTGGCTTTCATTTTAGTGGCGTTGATTTGCATCGGCTTACGTGCGGTTGGAATGCTCGGCGGAGAGTCCTGATTAGTGGCCTGGCAAGAGTCGCCAGTGACTTCACCGAGATCGCGGAACGGttt                                                                                               <  1:9566/145‑1 (MQ=255)
                                                      aGCTATGCAGTAATACAGCATTGGCTTTCATTTTAGAGGCGTTGATTTGCATCCGCTTACGTGCGGTTGGAATGCTCGGCGGAGAGTCCTGATTAGAGGCCTGGCAAGAGTCGCCACTGACTTCACCGAGATCGCGGAACGGTTTGcc                                                                                            <  2:199543/148‑1 (MQ=255)
                                                         tATGCAGTAATACAGCATTGGCTTTCATTTTAGAGGCGTTGATTTGCATCCGCTTACGTGCGGTTGGAATGCTCGGCGGAGAGTCCTGATTAGAGGCCTGGCAAGAGTCGCCACTGACTTCACCGAGATCGCGGAACGGTTTGCCGAc                                                                                         >  8:5799/1‑148 (MQ=255)
                                                         tATGCAGTAATACAGCATTGGCTTTCATTTTAGAGGCGTTGATTTGCATCCGCTTACGTGCGGTTGGAATGCTCGGCGGAGAGTCCTGATTAGAGGCCTGGCAAGAGTCGCCACTGACTTCACCGAGATCGCGGAACGGTTTGCCGAc                                                                                         >  4:313296/1‑148 (MQ=255)
                                                            gCAGTAATACAGCATTGGCTTTCATTTTAGAGGCGTTGATTTGCATCCGCTTACGTGCGGTTGGAATGCTCGGCGGAGAGTCCTGATTAGAGGCCTGGCAAGAGTCGCCACTGACTTCACCGAGATCGCGGAACGGTTTGCCGACTaa                                                                                      >  3:218711/1‑148 (MQ=255)
                                                                  aTACAGCATTGGCTTTCATTTTAGAGGCGTTGATTTGCATCCGCTTACGTGCGGTTGGAATGCTCGGCGGAGAGTCCTGATTAGAGGCCTGGCAAGAGTCGCCACTGACTTCACCGAGATCGCGGAACGGTTTGCCGACTAAttcttc                                                                                <  7:114936/148‑1 (MQ=255)
                                                                   tACAGCATAGGCTTTCATTTTAGAGGCGTTGATTTGCATCCGCTTACGTGCGGTTGGAATGCTCGGCGGAGAGTCCTGATTAGAGGCCTGGCAAGAGTCGCCACTGACTTCCCCGAGATAGCGGAACGGATTGCCGACTAATacttct                                                                               >  2:156563/1‑148 (MQ=255)
                                                                    aCAGCATTGGCTTTCATTTTAGAGGCGTTGATTTGCATCCGCTTACGTGCGGTTGGAATGCTCGGCGGAGAGTCCTGATTAGAGGCCTGGCAAGAGTCGCCACTGACTTCACCGAGATCGCGGAACGGTTTGCCGACTAATTCTTCTg                                                                              >  1:280116/1‑148 (MQ=255)
                                                                       gCATTGGCTTTCATTTTAGAGGCGTTGATTTGCATCCGCTTACGTGCGGTTGGAATGCTCGGCGGAGAGTCCTGATTAGAGGCCTGGCAAGAGTCGCCACTGACTTCACCGAGATCGCGGAACGGTTTGCCGACTAAttct                                                                                  <  7:121626/141‑1 (MQ=255)
                                                                       gCATTGGCTTTCATTTTAGAGGCGTTGATTTGCATCCGCTTACGTGCGGTTGGAATGCTCGGCGGAGAGTCCTGATTAGAGGCCTGGCAAGAGTCGCCACTGACTTCACCGAGATCGCGGAACGGTTTGCCGACTAAttct                                                                                  >  8:121626/1‑141 (MQ=255)
                                                                        cATTGGCTTTCATTTTAGAGGCGTTGATTTGCATCCGCTTACGTGCGGTTGGAATGCTCGGCGGAGAGTCCTGATTAGAGGCCTGGCAAGAGTCGCCACTGACTTCACCGAGATCGCGGAACGGTTTGCCGACTAATTCTTCTGCAtt                                                                          <  3:67691/148‑1 (MQ=255)
                                                                        cATTGGCTTTCATTTTAGAGGCGTTGATTTGCATCCGCTTACGTGCGGTTGGAATGCTCGGCGGAGAGTCCTGATTAGAGGCCTGGCAAGAGTCGCCACTGACTTCACCGAGATCGCGGAACGGTTTGCCGACTAATTCTTCTGCAtt                                                                          <  4:218711/148‑1 (MQ=255)
                                                                          ttGGCTTTCATTTTAGAGGCGTTGATTTGCATCCGCTTACGTGCGGTTGGAATTCTCGGCGGAGAGTCCTGATTAGAGGCCTGGCAAGAGTCGACACTGACTTCAACGATATAGCGTAACGGTTGGCCGACTAATTCTTCTGCATTgg                                                                        >  2:325758/1‑148 (MQ=255)
                                                                           tGGCTTTCATTTTAGAGGCGTTGATTTGCATCCGCTTACGTGCGGTTGGAATGCTCGGCGGAGAGTCCTGATTAGAGGCCTGGCAAGAGTCGCCACTGACTTCACCGAGATCGCGGAACGGTTTGCCGACTAATTCTTCTGCATTgg                                                                        >  7:169443/1‑147 (MQ=255)
                                                                                ttCATTTTAGAGGCGTTGATTTGCATCCGCTTACGTGCGGTTGGAATGCTCGGCGGAGAGTCCTGAtt                                                                                                                                                  <  4:357861/68‑1 (MQ=255)
                                                                                ttCATTTTAGAGGCGTTGATTTGCATCCGCTTACGTGCGGTTGGAATGCTCGGCGGAGAGTCCTGAtt                                                                                                                                                  >  3:357861/1‑68 (MQ=255)
                                                                                ttCATTTTAGAGGCGTTGATTTGCATCCGCTTACGTGCGGTTGGAATGCTCGGCGGAGAGTCCTGATTAGAGGCCTGGCAAGAGTCGCCACTGACTTCACCGAGATCGCGGAACGGTTTGCCGACTAATTCTTCTGCATTGGTATaaa                                                                  >  5:159089/1‑148 (MQ=255)
                                                                                 tCATTTTAGAGGCGTTGATTTGCATCCGCTTACGTGCGGTTGGAATGCTCGGCGGAGAGTCCTGATTAGAGGCCTGGCAAGAGTCGCCACTGACTTCACCGAGATCGCGGAACGGTTTGCCGACTAATTCTTCTGCATTGGTATaaa                                                                  <  7:321847/147‑1 (MQ=255)
                                                                                 tCATTTTAGAGGCGTTGATTTGCATCCGCTTACGTGCGGTTGGAATGCTCGGCGGAGAGTCCTGATTAGAGGCCTGGCAAGAGTCGCCACTGACTTCACCGAGATCGCGGAACGGTTTGCCGACTAATTCTTCTGCATTGGTATAAAt                                                                 >  4:375132/1‑148 (MQ=255)
                                                                                   aTTTTAGAGGCGTTGATTTGCATCCGCTTACGTGCGGTTGGAATGCTCGGCGGAGAGTCCTGATTAGAGGCCTGGCAAGAGTCGCCACTGACTTCACCGAGATCGCGGAACGGTTTGCCGACTAATTCTTCTGCATTGGTATAAAtt                                                                <  1:22282/147‑1 (MQ=255)
                                                                                    ttttAGAGGCGTTGATTTGCATCCGCTTACGTGCGGTTGGAATGCTCGGCGGAGAGTCCTGATTAGAGGCCTGGCAAGAGTCGCCACTGACTTCACCGAGATCGCGGAACGGTTTGCCGACTAATTCTTCTGCATTGGTATAAATTCg                                                              >  5:119155/1‑148 (MQ=255)
                                                                                         gagGCGTTGATTTGCATCCGCTTACGTGCTGTTGGAATGCTCGGCGTAGAGTCCTGATTAGATGACTGGCAAGAGTCTCCACTGACTTCAACGAGATCGCGGAACTGTTTGCCGACTAATTCTTCTGCATTGGTATAAATTCGGAcc                                                          <  1:325758/147‑1 (MQ=255)
                                                                                         gagGCGTTGATTTGCATCCGCTTACGTGCGGTTGGAATGCTCGGCGGAGAGTCCTGATTAGAGGCCTGGCAAGAGTCGCCACTGACTTCACCGAGATCGCGGAACGGTTTGCCGACTAATTCTTCTGCATTGGTATAAATTCGGAcc                                                          <  8:169443/147‑1 (MQ=255)
                                                                                         gagGCGTTGATTTGCATCCGCTTACGTGCGGTTGGAATGCTCGGCGGAGAGTCCTGATTAGAGGCCTGGCAAGAGTCGCCACTGACTTCACCGAGATCGCGGAACGGTTTGCCGACTAATTCTTCTGCATTGGTATAAATTCGGACCg                                                         >  3:141090/1‑148 (MQ=255)
                                                                                             cGTTGATTTGCATCCGCTTAGGTGCGGTTGGAATGCTCGGCGGAGAGTCCTGATTAGAGGCCTGGCAAGAGTCGCCACTGACTTCACCGAGATCGCGGAACGGTTTGCCGACTAATTCTTCTGCATTGGTATAAATTCGGACCGGCGt                                                     >  3:254351/1‑148 (MQ=255)
                                                                                              gTTGATTTGCATCCGCTTACGTGCGGTTGGAATGCTCGGCGGAGAGTCCTGATTAGAGGCCTGGCAAGAGTCGCCACTGACTTCACCGAGATCGCGGAACGGTTTGCCGACTAATTCTTCTGCAt                                                                           >  4:395240/1‑125 (MQ=255)
                                                                                              gTTGATTTGCATCCGCTTACGTGCGGTTGGAATGCTCGGCGGAGAGTCCTGATTAGAGGCCTGGCAAGAGTCGCCACTGACTTCACCGAGATCGCGGAACGGTTTGCCGACTAATTCTTCTGCAt                                                                           <  3:395240/125‑1 (MQ=255)
                                                                                                tGATTTGCATCCGCTTACGTGCGGTTGGAATGCTCGGCGGAGAGTCCTGATTAGAGGCCTGGCAAGAGTCGCCACTGACTTCACCGAGATCGCGGAACGGTTTGCCGACTAATTCTTCTGCATTGGTATAAATTCGGACCGGCGTGgc                                                  >  1:72410/1‑148 (MQ=255)
                                                                                                          ccGCTTACGTGCGGTTGGAATGCTCGGCGGAGAGTCCTGATTAGAGGCCTGGCAAGAGTCGCCACTGACTTCACCGAGATCGCGGAACGGTTTGCCGACTAATTCTTCTGCATTGGTATAAATTCGGACCGGCGTGGCGCGCGGCGCt                                        >  3:342143/1‑148 (MQ=255)
                                                                                                            gCTTACGTGCGGTTGGAATGCTCGGCGGAGAGTCCTGATTAGAGGCCTGGCAAGAGTCGCCACTGACTTCACCGAGATCGCGGAACGGTTTGCCGACTAATTCTTCTGCATTGGTATAAATTCGGACCGGCGTGGCGCGCGGCGCttt                                      >  4:123169/1‑148 (MQ=255)
                                                                                                             cTTACGTGCGGTTGGAATGCTCGGCGGAGAGTCCTGATTAGAGGCCTGGCAAGAGTCGCCACTGACTTCACCGAGATCGCGGAACGGTTTGCCGACTAATTCTTCTGCATTGGTATAAATTCGGACCGGCGTGGCGCGCGGCGCTTTc                                     >  4:55243/1‑148 (MQ=255)
                                                                                                             cTTACGTCCGGTTGGACTGCTCTGCGTATACTACTCAGGAGACGACTGGCAATAGTATCAACTGTCTTCACCGAGATCGCGGAACGGTTTGCCGACTAATTCTTCTGCATTGGTATAAATTCGGACCGGCGTGGCGCGCGGCGCTTTc                                     <  4:264223/148‑1 (MQ=255)
                                                                                                                aCGTGCGGTTGGAATGCTCGGCGGAGAGTCCTGATTAGAGGCCTGGCAAGAGTCGCCACTGACTTCACCGAGATCGCGGAACGGTTTGCCGACTAATTCTTCTGCATTGGTATAAATTCGGACCGGCGTGGCGCGCGGCGCTTTCGGt                                  >  1:208500/1‑148 (MQ=255)
                                                                                                                     cGGTTGGAATGCTCGGCGGAGAGTCCTGATTAGAGGCCTGGCAAGAGTCGCCACTGACTTCACCGAGATCGCGGAACGGTTTGCCGACTAATTCTTCTGCATTGGTATAAATTCGGACCGGCGTGGCGCGCGGCGCTTTCGGTTTTg                              <  6:119155/147‑1 (MQ=255)
                                                                                                                     cGGTTGGAATGCTCGGCGGAGAGTCCTGATTAGAGGCCTGGCAAGAGTCGCCACTGACTTCACCGAGATCGCGGAACGGTTTGCCGACTAATTCTTCTGCATTGGTATAAATTCGGACCGGCGTGGCGCGCGGCGCTTTCGGTTTTg                              <  2:72410/147‑1 (MQ=255)
                                                                                                                     cGGTTGGAATGCTCGGCGGAGAGTCCTGATTAGAGGCCTGGCAAGAGTCGCCACTGACTTCACCGAGATCGCGGAACGGTTTGCCGACTAATTCTTCTGCATTGGTATAAATTCGGACCGGCGTGGCGCGCGGCGCTTTCGGTTTTGc                             <  3:375132/148‑1 (MQ=255)
                                                                                                                     cGGTTGGAATGCTCGGCGGAGAGTCCTGATTAGAGGCCTGGCAAGAGTCGCCACTGACTTCACCGAGATCGCGGAACGGTTTGCAGACTAATTCTTCTGCATTGGTATAAATTCGGAACGGCGTGTCGCGCGGCGCTTTCGGTTTTg                              >  6:237222/1‑147 (MQ=255)
                                                                                                                       gTTGGAATGCTCGGCGGAGAGTCCTGATTAGAGGCCTGGCAAGAGTCGCCACTGACTTCACCGAGATCGCGGAACGGTTTGCCGAGTAATTCTTCTGCATTGGTATAAATTCGGACCGGCGTGGCGCGCGGCCCTTTCGGTTTTGCAg                           >  6:360803/1‑148 (MQ=255)
                                                                                                                             aTGCTCGGCGTAGATTCCTGATTAGAGGCCTGGCAAGAGTCTCCACTGACTTCACCGAGATCGCGGAACGGTTTGCCGACTAATTCTTCTGCATTGGTATAAATTCGGACCGGCGTGGCGCGCGGCGCTTTCGGTGTTGCAGGCTCCg                     <  5:371140/148‑1 (MQ=255)
                                                                                                                             aTGCTCGGCGGAGAGTCCTGATTAGAGGCCTGGCAAGAGTCGCCACTGACTTCACCGAGATCGCGGAACGGTTTGCCGACTAATTCTTCTGCATTGGTATAAATTCGGACCGGCGTGGCGCGCGGCGCTTTCGGTTTTGCAGGCTCCg                     >  4:233892/1‑148 (MQ=255)
                                                                                                                               gCTCGGCGGAGAGTCCTGATTAGAGGCCTGGCAAGAGTCGCCACTGAATTCACCGAGATCGCGGAACGGTTTGCCGACTAATTCTTCTGCATTGGTATAAATTCGGACCGGCGTGGCGCGCGGCGCTTTCGGTTTTGCAGGCTCCGCt                   >  3:57992/1‑148 (MQ=255)
                                                                                                                                cgcggcggAGAGTCCTGATTAGAGGCCTGGCAAGAGTCGCCACTGACTTCACCGAGATCGCGGAACGGTTTGCCGACTAATTCTTCTGCATTGGTATAAATTCGGACCGGCGTGGCGCGCGGCGCTTTCGGTTTTGCAGGCTCCGCtt                  <  7:342632/146‑1 (MQ=255)
                                                                                                                                  cggcggAGAGTCCTGATTAGAGGCCTGGCAAGAGTCGCCACTGACTTCACCGAGATCGCGGAACGGTTTGCCGACTAATTCTTCTGCATTGGTATAAATTCGGACCGGCGTGGCGCGCGGCGCTTTCGGTTTTGCAGGCTCCGCTTTc                >  7:265101/1‑148 (MQ=255)
                                                                                                                                  cggcggAGAGTCCTGATTAGAGGCCTGGCAAGAGTCGCCACTGACTTCACCGAGATCGCGGAACGGTTTGCCGACTAATTCTTCTGCATTGGTATAAATTCGGACCGGCGTGGCGCGCGGCGCTTTCGGTTTTGCAGGCTCCGCTTTc                >  3:86521/1‑148 (MQ=255)
                                                                                                                                     cggAGAGTCCTGATTAGAGGCCTGGCAAGAGTCGCCACTGACTTCACCGAGATCGCGGAACGGTTTGCCGACTAATTCTTCTGCATTGGTATAAATTCGGACCGGCGTGGCGCGCGGCGCTTTCGGTTTTGCAGGCTCCGCTTTCGGc             <  4:149394/148‑1 (MQ=255)
                                                                                                                                     cggAGAGTCCTGATTAGAGGCCTGGCAAGAGTCGCCACTGACTTCACCGAGATCGCGGAACGGTTTGCCGACTAATTCTTCTGCATTGGTATAAATTCGGACCGGCGTGGCGCGCGGCGCTTTCGGTTTTGCAGGCTCCGCTTTCGGc             <  6:130983/148‑1 (MQ=255)
                                                                                                                                       gagaGTCCTGATTAGAGGCCTGGCAAGAGTCGCCACTGACTTCACCGAGATCGCGGAACGGTTTGCCGACTAATTCTTCTGCATTGGTATAAATTCGGACCGGCGTGGCGCGCGGCGCTTCGGTTTTGCAGGCTCCGCTTTCGGCTgg          <  2:280116/148‑1 (MQ=255)
                                                                                                                                        agagTCCTGATTAGAGGCCTGGCAAGAGTAGCCACTGACTTCACCGAGATCGCGGAACGGTTTGCCGACTAATTCTTCTGCATTGGTATAAATTCGGCCCGGCGTGGCGCGCGGCGCTTTCGGTTTTGCAGGCTCCGCTTTCGGCTgg          <  7:407905/148‑1 (MQ=255)
                                                                                                                                             ccTGACTAGAGGCCTGGCAAGAGTCGCCACTGACTTCACCGAGATCGCGGAACGGTTTGCCGACTAATTCTTCTGCATTGGTATAAATTCGGACCGGCGTGGCGCGCGGCGCTTTCGGTTTTGCAGGCTCCGCTTTCGGCTGGGGTGc     <  1:105911/148‑1 (MQ=255)
                                                                                                                                              cTGATTAGAGGCCTGGCAAGAGTCGCCACTGACTTCACCGAGATCGCGGAACGGTTTGCCGACTAATTCTTCTGCATTGGTATAAATTCGGACCGGCGTGGCGCGCGGCGCTTTCGGTTTTGCAGGCTCCGCTTTCGGCTGGGGTGc     <  3:313296/147‑1 (MQ=255)
                                                                                                                                              cTGATTAGAGGCCTGGCAAGAGTCGCCACTGACTTCACCGAGATCGCGGAACGGTTTGCCGACTAATTCTTCTGCATTGGTATAAATTCGGACCGGCGTGGCGCGCGGCGCTTTCGGTTTTGCAGGCGCCGCTTTCGGCTGGGGTGc     <  1:156563/147‑1 (MQ=255)
                                                                                                                                              cTGACTAGAGGCCTGGCAAGAGTCGCCACTGACTTCACCGAGATCGCGGAACGGTTTGCCGACTAATTCTTCTGCATTGGTATAAATTCGGACCGGCGTGGCGCGCGGCGCTTTCGGTTTTGCAGGCTCCGCTTTCGGCTGGGGTGc     <  1:111856/147‑1 (MQ=255)
                                                                                                                                                gATTAGAGGCCTGGCAAGAGTCGCCACTGACTTCACCGAGATCGCGGAACGGTTTGCCGACTAATTCTTCTGCATTGGTATAAATTCGGACCGGCGTGGCGCGCGGCGCTTTCGGTTTTGCAGGCTCCGCTTTCGGCTGGGGTGCAGt  >  8:201646/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
AACCGATACATACAGCCTGACGATAGCAGCCTGGCGTACCGCTGGTGACTTCGCAGCTATGCAGTAATACAGCATTGGCTTTCATTTTAGAGGCGTTGATTTGCATCCGCTTACGTGCGGTTGGAATGCTCGGCGGAGAGTCCTGATTAGAGGCCTGGCAAGAGTCGCCACTGACTTCACCGAGATCGCGGAACGGTTTGCCGACTAATTCTTCTGCATTGGTATAAATTCGGACCGGCGTGGCGCGCGGCGCTTTCGGTTTTGCAGGCTCCGCTTTCGGCTGGGGTGCAGT  >  NC_000913/219617‑219908

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