Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,087,337 G→A 5.8% R170R (CGC→CGT)  pgaC ← biofilm PGA synthase PgaCD, catalytic subunit; poly‑beta‑1,6‑N‑acetyl‑D‑glucosamine synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,087,3370GA 5.8% 216.8 / 2.2 70R170R (CGC→CGT) pgaCbiofilm PGA synthase PgaCD, catalytic subunit; poly‑beta‑1,6‑N‑acetyl‑D‑glucosamine synthase
Reads supporting (aligned to +/- strand):  ref base G (37/28);  new base A (2/2);  total (40/30)
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

TTCGCTTGATCAAACCAATAATTGAGGAATACTCGCCAACCTGAATTTTACCCACCAGGGTAGAACGTGTTCGAATACGAGGATTACCGGTTACGGCACCCACACGCGGGTTGTACAACATCGGTTCCACAATATATGCCGCCGCATCGCGGTCTAATAACGCATCGCCATCAATGCACACCAGATATTCACTTTTCGCCGCGGCAGCTCCGGTTTTAAGCGCAATGGCTTTCCCCTGGTTTTGCGCCAGATGAATGACCCGCAAATGGGGAATTTGTGCAGCCATGCGATCCAGGA  >  NC_000913/1087189‑1087485
                                                                                                                                                    |                                                                                                                                                    
ttCGCTTGATCAAACCAATAATTGAGGAATACTCGCCAACCTGAATTTTACCCACCAGGGTAGAACGTGTTCGAATACGAGGATTACCGGTTACGGCACCCACACGCGGGTTGTACAACATCGGTTCCACAATATATGCCGCCGCATcg                                                                                                                                                      <  5:47292/149‑1 (MQ=255)
  cGCTTGATCAAACCAATAATTGAGGAATACTCGCCAACCTGAATTTTACCCACCAGGGTAGAACGTGTTCGAATACGAGGATTACCGGTTACGGCACCCACACGCGGGTTGTACAACATCGGTTCCACAATATATGCCGCCGCATcgc                                                                                                                                                     <  5:153604/148‑1 (MQ=255)
    cTTGATCAAACCAATAATTGAGGAATACTCGCCAACCTGAATTTTACCCACCAGGGTAGAACGGGTTCGAATACGAGGATTACCGGTTACGGCACCCACACGCTGGTTGTACAACATCGGTTCCACAATATATGCCGCCGCATCGCGGt                                                                                                                                                  <  5:273117/149‑1 (MQ=255)
    cTTCATCAACCCAATAATCGAGGAATACTCGCCAACCTGAATTTTACCCACCAGGGTAGAACGTGTTCGAATACGAGGATTACCGGTTACGGCACCCACACGCGGGTTGTACAACATCGGTTCCACAATATATGCCGCCGCATCGCGGt                                                                                                                                                  <  2:231107/149‑1 (MQ=255)
     ttGATCAAACCAATAATTGAGGAATACTCGCCAACCTGAATTTTACCCACCAGGGTAGAACGTGTTCGAATACGAGGATTACCGGTTACGGCACCCACACGCGGGTTGTACAACATCGGATCCACAATATATGCCGCCGCATCGCGGt                                                                                                                                                  >  5:271958/1‑148 (MQ=255)
       gATCAAACCAATAATTGAGGAATACTCGCCAACCTGAATTTTACCCACCAGGGTAGAACGTGTTCGAATACGAGGATTACCGGTTACGGCACCCACACGCGGGTTGTACAACATCGGATCCACAATATATGCCGCCGCATCGCGGt                                                                                                                                                  <  6:271958/146‑1 (MQ=255)
          cAAACCAATAATTGAGGAATACTCGCCAACCTGAATTTTACAAACCGGTGTAGAAAGTGTTCGAATACGAGGATTACCGGTTACTGAACCCAAACGCGTGTTGTACAACATCGGTTAAAAAAAATATGACGCCGCAAAGAGGTAtaat                                                                                                                                             >  2:34039/1‑148 (MQ=255)
           aacccAAGAATTGAGGAATACTCGCCAACCTGAATTTTACCCACCAGGGTAGAACGTGTTCGAATACGAGGATTACCGGTTACGGCACCCACACGCGGGTTGTCCAACATCGGTTCCACAATATATGCCGCCGCATCACGGTCtaataa                                                                                                                                           >  1:212216/4‑149 (MQ=255)
           aaaCCAATAATTGAGGAATACTCGCCAACCTGAATTTTACCCACCAGGGTAGAACGTGTTCGAATACGAGGATTACCGGTTACGGCACCCACACGCGGGTTGTACAACATCGGTTCCACAATATATGCCGCCGCATCACGGTCtaataa                                                                                                                                           >  1:212213/1‑149 (MQ=255)
                 ataatTGAGGAATACTCGCCAACCTGAATTTTACCCACCAGGGTAGAACGTGTTCGAATACGAGGCTTACCGGTTACGGCACACACACGCGGGTTGTACAACAGCGGTTCCACAATATATGCCGCCGCATCCCGGTCTAATACCGCagc                                                                                                                                     >  6:216983/1‑147 (MQ=255)
                  taatTGAGGAATACTCGCCAACCTGAATTTTACCCACCAGGGTAGAACGTGTTCGAATACGAGGATTACCGGTTACGGCACCCACACGCGGGTTGTACAACATCGGTTCCACAATATATGCCGCCGCATCGCGGTCTAATAACGCATCg                                                                                                                                    >  4:2144/1‑149 (MQ=255)
                     ttGAGGAATACTCGCCAACCTGAATTTTACCCACCAGGGTAGAACGTGTTCGAATACGAGGATTACCGGTTACGGCACCCACACGCGGGTTGTACAACATCGGTTCCACAATATATGCCGCCGCATCGCGGTCTAATAACGCATCGCCa                                                                                                                                 >  5:126657/1‑149 (MQ=255)
                             tccTCGCCAACCTGAATTTTACCCACCAGGGTAGAACGTGTTCGAATACGAGGATTACCGGTTACGGCACCCACACGCGGGTTGTACAACATCGGTTCCACAATATATGCCGCCGCATCGCGGGCTAATAACGCATCGCCATCAATGca                                                                                                                         <  8:15147/147‑1 (MQ=255)
                                     aaCCTGAATTTTACCCACCAGGGTAGAACGTGTTCGAATACGAGGATTACCGGTTACGGCACCCACACGCGGGTTGTACAACATAGGTTCCACAATATATGCCGCCGCATCGCGGTCTAATAACGCATCGCCATCAATGCACACCAGat                                                                                                                 <  1:248928/149‑1 (MQ=255)
                                      aCCTGAATTTTACCCACCAGGGTAGAACGTGTTCGAATACGAGGATTACCGGTTACGGCACCCACACGCGGGTTGTACAACATCGGTTCCACAATATATGCCGCCGCATCGCGGTCTAATAACGCATCGCCATCAATGCACACCAGata                                                                                                                >  7:193822/1‑149 (MQ=255)
                                       ccTGAATTTTACCCACCAGGGTAGAACGTGTTCGAATACGAGGATTACCGGTTACGGCACCCACACGCGGGTTGTACAACATCGGTTCCACAATATATGCCGCCGCATCGCGGTCTAATAACGCATCGCCATCAATGCACACCAGatat                                                                                                               >  2:15091/1‑149 (MQ=255)
                                         tGAATTTTACCCACCAGTGTAGAACGTGTTCGAATACGAGGATTACCGGTTACGGCACCCACACGCGGGTTGTACAACATCGGTTCCACAATATATGCCGCCGCATCGCGGTCTAATAACGCATCGCCATCa                                                                                                                              <  5:70613/132‑1 (MQ=255)
                                         tGAATTTTACCCACCAGGGTAGAACGTGTTCGAATACGAGGATTACCGGTTACGGCACCCACACGCGGGTTGTACAACATCGGTTCCACAATATATGCCGCCGCATCGCGGTCTAATAACGCATCGCCATCa                                                                                                                              >  6:70613/1‑132 (MQ=255)
                                           aaTTTTACCCACCAGTGTAGAACGTGTTCGAATACGAGGATTACCGGTTACGGCACCCACACGCGGGTTGTACAACATCGGTTCCACAATATATGCCGCCGCATCGCGGTCTAATAACGCACCGCCATCAATGCACACCAGATATTCAc                                                                                                           <  4:14216/149‑1 (MQ=255)
                                            aTTTTACCCTCCAGGGTAGAACGTGTTCGAATACGAGGATTACCGGTTACGGCACCCACACGCGGGTTGTACAACATCGGTTCCACAATATATGCCGCCGCATCGCGGTCTAATAACGCATCGCCATCAATGCACACCAGATATTCAc                                                                                                           <  5:161602/148‑1 (MQ=255)
                                             ttttACCCACCAGGGTAGAACGTGTTCGAATACGAGGATTACCGGTTACGGCACCCACACGCGGGTTGTACAACATCGGTTCCACAATATATGCCGCCGCATCGCGGTCTAATAACGCATCGCCATCAATGCACACCAGATATTCACtt                                                                                                         <  1:244370/149‑1 (MQ=255)
                                                     accaGGGTAGAACGTGTTCGAATACGAGGATTACCGGTTACGGCACCCACACGCGGGTTGTACAACATCGGTTCCACAATATATGCCGCCGCATCGCGGTCTAATAACGCATCGCCATCAATGCACACCAGATATTCACTTTTcgccgc                                                                                                 >  7:314887/1‑149 (MQ=255)
                                                     accaGGGTAGAACGTGTTCGAATACGAGGATTACAGGTTACGGCACCCACACGCGGGTTGTACAACATCGGTTCCACAATATATGCCGCCGCATCGCGGTCTAAGACCGCCTCGCCAACAATGCACACCAGATATTCACATTTcgccgc                                                                                                 >  5:28950/1‑149 (MQ=255)
                                                      ccaGGGTAGAACGTGTTCGAATACGAGGATTACCGGTTACGGCACCCACACGCGGGTTGTACAACATCGGTTCCACAATATATGCCGCCGCATCGCGGTCTAATAACGCATCGCCATCAATGCACACCAGATATTCACTTTTCGCcgcg                                                                                                >  8:131790/1‑149 (MQ=255)
                                                         gggTAGAACGTGTTCGAATACGAGGATTACCGGTTACGGCACCCACACGCGGGTTGTACAACATCGGTTCCACAATATATGCCGCCGCATCGCGGt                                                                                                                                                  >  7:80711/1‑96 (MQ=255)
                                                         gggTAGAACGTGTTCGAATACGAGGATTACCGGTTACGGCACCCACACGCGGGTTGTACAACATCGGTTCCACAATATATGCCGCCGCATCGCGGt                                                                                                                                                  <  8:80711/96‑1 (MQ=255)
                                                              gatcGTGTTCGAATACGAGGATTTCCTGTTACGGCACCCACACGCGGGTTGTACAACATCGGTGCCACAATATATGCCGCCGCATAGCTGTCTAATATCGCCTCGCCATCAATGCACACCAGATATTCACTTATCGCCGCGGCAGCTcc                                                                                        <  8:83174/146‑1 (MQ=255)
                                                                 cGTGTTCGAATACGATGATTACCGGTTACGGCACCCACACGCGGGTTGTACTACATCGGTTCCACAATATATGCCGCCGCATCGCGGTCTAATAACGCATCGCCATCAATGCACACCAGATATTCACTTTTCGCCGCGGCAGCTCCGGt                                                                                     <  5:150973/149‑1 (MQ=255)
                                                                    gtTCGAATACGAGGATTACCGGTTACGGCACCCACACGCGGGTTGTACAACATCGGTTCCACAATATATGCCGCCGCATCGCGGTCTAATAACGCATCGCCATCAATGCACGCCAGATATTCACTTTTCGCCGCGGCAGCTCCGGttt                                                                                   >  1:260683/1‑148 (MQ=255)
                                                                          aTACGAGGATTACCGGTTACGGCACCCACACGCGGGTTGTACAACATCGGTTCCACAATATATGCCGCCGCATCGCGGTCTAATAACGCATCGCCATCAATGCACACCAGATATTCACTTTTCGCCGCGGCAGCTCCGGTTTTAAgcg                                                                             >  3:310980/1‑148 (MQ=255)
                                                                                ggATTACCGGTTACGGCACCCACACGCGGGTTGTACAACATCGGTTCCACAATATATGCCGCCGCATCGCGGTCTAATAACGCATCGCCATCAATGCACACCAGATATTCACTTTTCGCCGCGGCAGCTCCGGTTTTAAGCGCAATGGc                                                                      >  3:325209/1‑149 (MQ=255)
                                                                                   ttACCGGTTACGGCACACACACGCGGGTTGTAAAACATCGGTTCCACAATATATGCCGCCGCATCGCGGTCTAATAACGCATCGCCATCAATGCACACCAGATATTCACTTTTCGCCGCGGCAGCTCCGGTTTTAAGCGCACTGGCttt                                                                   >  8:188152/1‑149 (MQ=255)
                                                                                          ttACGGCACCCACACGCGGGTTGTACAACATCGGTTCCACAATATATGCCGCCGCATCGCCGTCTAATAACGCATCGCCATCAATGCACACCAGATATTCACTTTTCGCCGCGGCAGCTCCGGTTTTAAGCGCAATGGCTTTCCCCTgg                                                            >  5:177455/1‑149 (MQ=255)
                                                                                            aCGGCACCCACACGCGGGTTGTACAACATCGGTTCCACAATATATGCCGCCGCATCGCGGTCTAATAACGCATCGCCATCAATGCACACCAGATATTCACTTTTCGCCGCGGCAGCTCCGGTGTTAAGCGCAATGGCTTTCCCCTGGtt                                                          >  8:352646/1‑149 (MQ=255)
                                                                                              ggCACCCACACGCGGGTTTTACATCATCGGTTCCACAATATATGCCGCCGCATCGCGGTCTAATAACGCATCGCCATCTATGCACACCAGATATTCACTTTTCGACGCGGCAGCTCCGGTTTTAAGCGCAATGGCTTTCCCCTgg                                                            <  5:13946/145‑1 (MQ=255)
                                                                                              ggCACCCACACGCGGGTTGTACATCATCGGTTCCACAATATATGCCGCCGCATCGCGGTCTAATAACGCATCGCCATCAATGCACACCAGATATTAACTTTTCTCCGCGGCAGCTCCGGTTTTAAGCGCAATGGATTTCCCCTgg                                                            >  6:13946/1‑145 (MQ=255)
                                                                                                      cacGCGGGTTGTACAACATCGGTTCCACAATATATGCCGCCGCATCGCGGTCTAATAACGCATCGCCATCAATGCACACCAGATATTCACTTTTCGCCGCGGCAGCTCCGGTTTTAAGCGCAATGGCTTTCCCCTGGTTTTGCGCCAGa                                                <  2:308226/149‑1 (MQ=255)
                                                                                                           gggTTGTACAACATCGGTTCCACAATATATGCCGCCGCATCGCGGTCTAATAACGCATCGCCATCAATGCACACCAGATATTCACTTTTCGCCGCGGCAGCTCCGGTTTTAAGCGCAATGGCTTTCCCCTGGTTTTGCGCCAGatgaat                                           <  8:217743/149‑1 (MQ=255)
                                                                                                             gTTGTACAACATCGGTTCCACAATATATGCCGCCGCATCGCGGTCTAATAACGCATCGCCATCAATGCACACCAGATATTCACTTTTCGCCGCGGCAGCTCCGGTTTTAAGCGCAATGGCTTTCCCCTGGTTTTGCGCCAGatgaatga                                         >  6:152254/1‑149 (MQ=255)
                                                                                                              ttGTACAACATCGGTTCTACAATATATGCCGCCGCATCGCGGTCTAATAACGCATCGCCATCAATGCACACCAGATATTCACTTTTCGCCGCGGCAGCTCCGGTTTTAAGCGCAATGGCTTTCCCCTGGTTTTGCGCCAGATGAATGAc                                        >  6:140441/1‑149 (MQ=255)
                                                                                                              ttGTACAACATCGGTTCCACAATATATGCCGCCGCATCACGGTCTAATAACGCATCGCCATCAATGCACCCCAGATCTTCACTTTTCGCCGCGGCAGCTCCGGGGTTACGCGCAATGGCTTTCCCCTGGTTTTGCGCCAGATGAATGAc                                        <  2:212216/149‑1 (MQ=255)
                                                                                                              ttGTACAACATCGGTTCCACAATATATGCCGCCGCATCACGGTCTAATAACGCATCGCCATCAATGCACACCAGATATTCACTTTTCGCCGCGGCAGCTCCGGTTTTAAGCGCAATGGCTTTCCCCTGGTTTTGCGCCAGATGAATGAc                                        <  2:212213/149‑1 (MQ=255)
                                                                                                                 tACAACATCGGTTCCACAATATATGCCGGCGCATCGCGGGCTAATAACGCATCTCCCTCAATGCACACCAGATATTCACTTTTCGCCGCGGCAGCTCCGGTTTTAAGCGCAATGGCTTTCCCCTGGTTTTGCGCCAGATGAATGAccc                                      <  6:28950/148‑1 (MQ=255)
                                                                                                                 tACAACATCGGTTCCACAATATATGCCGCCGCATCGCGGTCTAATAACGCATCGCCATCAATGCACACCAGATATTCACTTTTCGCCGCGGCAGCTCCGGTTTTAAGCGCAATGGCTTTCCCCTGGTTTTGCGCCAGATGAATGACCCg                                     >  8:324776/1‑149 (MQ=255)
                                                                                                                 tACAACATCGGGTCCACAATATATGCCGCCGCATCGCGGTCTAATAACGCATCGCCATCAATGCACACCAGATATTCACTTTTCGCCGCGGCAGCTCCGGTTTTAAGCGCAATGGCTTTCCCCTGGTTTTGCGCCAGATGAATGAccc                                      <  4:93007/148‑1 (MQ=255)
                                                                                                                 tACAACATCCGTACCACAATATATGCCGCCGCATCGCGGTCTAATAACGCATCGCCATCAATGCACACCAGATATTCACTTTTCGCCGCGGCAGCTCCGGTTTTAATCGCAATGGCTTTCCCCTGGTTTTGCGCCAGATGAATGAccc                                      <  3:29080/148‑1 (MQ=255)
                                                                                                                  acaacaTCGGTTCCACAATATATGCCGCCGCATCGCGGTCTAATAACGCATCGCCATCAATGCACACCAGATATTCACTTTTCGCCGCGGCAGCTGCGGTTTTAAGCGCAATGGCTTTCCCCTGGTTTTGCGACAGATGAATGACCCg                                     >  6:197640/1‑148 (MQ=255)
                                                                                                                           gTTCCACAATATATGCCGCCGCATCGCGGTCTAATAACGCATCGCCATCAATGCACACCAGATATTCACTTTTCGCCGCGGCAGCTCCGGTTTTAAGCGCAATGGCTTTCCCCTGGTTTTGCGCCAGATGAATGACCCGCAAATgggg                            >  6:279686/1‑148 (MQ=255)
                                                                                                                            ttCCACAATATATGCCGCCGCATCGCGGTCTAATAACGCATCGCCATCAATGCACACCAGATATTCACTTTTCGCCGCGGCAGCTCCGGTTTTAAGCGCAATGGCTTTCCCCTGGTTTTGCGCCAGatga                                             <  4:222198/130‑1 (MQ=255)
                                                                                                                            ttCCACAATATATGCCGCCGCATCGCGGTCTAATAACGCATCGCCATCAATGCACACCAGATATTCACTTTTCGCCGCGGCAGCTCCGGTTTTAAGCGCAATGGCTTTCCCCTGGTTTTGCGCCAGatga                                             >  3:222198/1‑130 (MQ=255)
                                                                                                                            ttCCACAATATATGCCGCCGCATCGCGGTCTAATAACGCATCGCCATCAATGCACACCAGATATTCACTTTTCGCCGCGGCAGCTCCGGTTTTAAGCGCAATGGCTTTCCCCTGGTTTTGCGCCAGATGAATGACCCGCAAATGGGGaa                          >  2:200975/1‑149 (MQ=255)
                                                                                                                            ttCCACAATATATGCCGCCGCATCGCGGTCTAATAACGCATCGCCATCAATGCACACCAGATATTCACTTTTCGCCGCGGCAGCTCCGGTTTTAAGCGCAATGGCTTTCCCCTGGTTTTGCGCCAGATGAATGACCCGCAAATGGGGaa                          <  8:193822/149‑1 (MQ=255)
                                                                                                                              ccACAATATATGCCGCCGCATCGCGGTCTAATAACGCATCGCCATCAATGCACACCAGATATTCACTTTTCGCCGCGGCAGCTCCGGTTTTAAGCGCAATGGCTTTCCCCTGGTTTTGCGCCAGATGAATGACCCGCAAATGGGGAAtt                        >  8:7813/1‑149 (MQ=255)
                                                                                                                               cacaATATATGCCGCCGCATCGCGGTCTAATAACGCATCGCCATCAATGCACACCAGATATTCACTTTTCGCCGCGGCAGCTCCGGTTTTAAGCGCAATGGCTTTCCCCTGGTTTTGCGCCAGATGAATGAc                                        >  6:223178/1‑132 (MQ=255)
                                                                                                                               cacaATATATGCCGCCGCATCGCGGTCTAATAACGCATCGCCATCAATGCACACCAGATATTCACTTTTCGCCGCGGCAGCTCCGGTTTTAAGCGCAATGGCTTTCCCCTGGTTTTGCGCCAGATGAATGAc                                        <  5:223178/132‑1 (MQ=255)
                                                                                                                               cacaATATATGCCGCCGCATCGCGGTCTAATAACGCATCGCCATCAATGCACACCAGATATTCACTTTTCGCCGCGGCAGCTCCGGTTTTAAGCGCAATGGCTTTCCCCTGGTTTTGCGCCAGATGAATGACCCGCAAATGGGGAAttt                       >  5:315924/1‑149 (MQ=255)
                                                                                                                               cacaATATATGCCGCCGCATCGCGGTCTAATAACGCATCGCCATCAATGCACACCAGATATTCACTTTTCGCCGCGGCAGCTCCGGTTTTAAGCGCAATGGCTTTCCCCTGGTTTTGCGCCAGATGAATGACCCGCAAATGGGGAAttt                       <  4:310980/149‑1 (MQ=255)
                                                                                                                               cacaATATATGCCGCCGCATCGCGGTCTAATAACGCATCGCCATCAATGCACACCAGATATTCACTTTTCGCCGCGGCAGCTCCGGTTTTAAGCGCAATGGCTTTCCCCTGGTTTTGCGCCAGATGAATGACCCGCAAATGGGGAAttt                       >  6:137815/1‑149 (MQ=255)
                                                                                                                                 caATATATGCCGCCGCATCGCGGTCTAATAACGCATCGCCATCAATGCACACCAGATATTCACTTTTCGCCGCGGCAGCTCCGGTTTTAAGCGCAATGGCTTTCCCCTGGTTTTGCGCCAGATGAATGACCCGCAAATGGGGAATTtgt                     >  3:26922/1‑149 (MQ=255)
                                                                                                                                  aaactctGCCGCCGCATCGCGGTCTAATAACCCATCGCCATCAATGCACACCAGATATTCACTTTTCGCCTCGGCAGCTCCGGTTTTAAGCGCAATGGCTTTCCCCTGGTTTTGCGCCATATGAATGACCCGCAAATg                               <  8:329115/132‑1 (MQ=255)
                                                                                                                                  aaTATATGCCGCCGCATCGCGGTCTAATAACGCATCGCCATCAATCCACACCAGATAATCCCTTATCGCCGCGGCAGCTCCGGTTTTAAGCGCAAAGGCTTTCCCCTGGTTTTGCGCCAGATGAATGACCCGCAAATg                               >  7:329115/1‑138 (MQ=255)
                                                                                                                                   atatatGCCGCCGCATCGCGGTCTAATAACGCATCGCCATCAATGCACACCAGATATTCACTTTTCGCCGCGGCAGCTCCGGTTTTAAGCGCAATGGCTTTCCCCTGGTTTTGCGCCAGATGAATGACCCGCAAATGGGGAATTTGTGc                   <  7:352646/149‑1 (MQ=255)
                                                                                                                                         gccgccGCATCGCGGTCTAATAACGCATCGCCATCAATGCACACCAGATATTCACTTTTCGCCGCGGCAGCTCCGGTTTTAAGCGCAATGGCTTTCCCCTGGTTTTGCGCCAGATGAATGACCCGCAAATGGGGAATTTGTGCAGCCAt             >  5:48406/1‑149 (MQ=255)
                                                                                                                                          ccgccgCATCGCGGTCTAATAACGCATCGCCATCAATGCACACCAGATCTTCACTTTTCGCCGCGGCAGCTCCGGTTTTAAGCGCAAGGGCTTTCCCCTGGGTTTGCGCCAGATGAATGACCCGCAAATGGGGAATTTGTGCAGCCATg            <  5:197640/149‑1 (MQ=255)
                                                                                                                                          ccgccgCATCGCGGTCTAATAACGCATCGCCATCAATGCACACCAGATATTCACTTTTCGCCGCGGCAGCTCCGGTTTTAAGCGCAATGGCTTTCCCCTGGTTTTGCGCCAGATGAATGACCCGCAAATGGGGAATTTGTGCAGCCATg            >  7:158657/1‑149 (MQ=255)
                                                                                                                                           cgccgcATCGCGGTCTAATAACGCATCGCCATCAATGCACACCAGATATTCACTTTTCGCCGCGGCAGCTCCGGTTTTAAGCGCAATGGCTTTCCCCTGGTTTTGCGCCAGATGAATGACCCGCAAATGGGGAATTTGTGCAGCCATGc           >  8:27853/1‑149 (MQ=255)
                                                                                                                                            gccgcATCGCGGTCTAATAACGCATCGCCATCAATGCACACCAGATATTCACTTTTCGCCGCGGCAGCTCCGGTTTTAAGCGCAATGGCTTTCCCCTGGTTTTGCGCCAGATGAATGACCCGCAAATgg                              <  2:224716/129‑1 (MQ=255)
                                                                                                                                            gccgcATCGCGGTCTAATAACGCATCGCCATCAATGCACACCAGATATTCACTTTTCGCCGCGGCAGCTCCGGTTTTAAGCGCAATGGCTTTCCCCTGGTTTTGCGCCAGATGAATGACCCGCAAATgg                              >  1:224716/1‑129 (MQ=255)
                                                                                                                                             ccgcATCGCGGTCTAATAACGCATCGCCATCAATGCACACCAGATATTCCCTTTTCGCCGCGGCAGCTCCGGTTTTAAGCGCAATGGCTTTCCCCTGGTTTTGCGCCAGATGAATGACCCGCAAATGGGGAATTTGTGCAGCCATGCGa         >  5:199034/1‑149 (MQ=255)
                                                                                                                                             ccgcATCGCGGTCTAATAACGCATCGCCATCAATGCACACCAGATATTCACTTTTCGCCGCGGCAGCTCCGGTTTTAAGCGCAATGGCTTTCCCCTGGTTTTGCGCCAGATGAATGACCCGCAAATGGGGAATTTGTGCAGCCATGCGa         >  7:216544/1‑149 (MQ=255)
                                                                                                                                                cATCGCGGTCTAATAACGCATCGCCATCAATGCACACCAGATATTCACTTTTCGCCGCGGCAGCTCCGGTTTTAAGCGCAATGGCTTTCCCCTGGTTTTGCGCCAGATGAATGACCCGCAAATGGGGAATTTGTGCAGCCATGCGATcc      <  4:325209/149‑1 (MQ=255)
                                                                                                                                                 aTCGCGGTCTAATAACGCATCGCCATCAATGCACACCAGATATTCACTTTTCGCCGCGGCAGCTCCGGTTTTAAGCGCAATGGCTTTCCCCTGGTTTTGCGCCAGATGAATGACCCGCAAATGGGGAATTTGTGCAGCCATGCGATCCa     <  1:15091/149‑1 (MQ=255)
                                                                                                                                                    gcgGTCTAATAACGCATCGCCATCAATGCACACCAGATATTCACTTTTCGCCGCGGCAGCTCCGGTTTTAAGCGCAATGGCTTTCCCCTGGTTTTGCGCCAGATGAATGACCCGCAAATGGGGAATTTGTGCAGCCATGCGATCCAGGa  <  5:152254/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TTCGCTTGATCAAACCAATAATTGAGGAATACTCGCCAACCTGAATTTTACCCACCAGGGTAGAACGTGTTCGAATACGAGGATTACCGGTTACGGCACCCACACGCGGGTTGTACAACATCGGTTCCACAATATATGCCGCCGCATCGCGGTCTAATAACGCATCGCCATCAATGCACACCAGATATTCACTTTTCGCCGCGGCAGCTCCGGTTTTAAGCGCAATGGCTTTCCCCTGGTTTTGCGCCAGATGAATGACCCGCAAATGGGGAATTTGTGCAGCCATGCGATCCAGGA  >  NC_000913/1087189‑1087485

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: ‑