Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,218,239 T→A 7.6% F96Y (TTC→TAC) ‡ PP_4598 → membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,218,2390TA 7.6% 190.6 / 4.4 65F96Y (TTC→TAC) ‡PP_4598membrane protein
Reads supporting (aligned to +/- strand):  ref base T (28/32);  new base A (3/2);  total (31/34)
Fisher's exact test for biased strand distribution p-value = 6.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01

CGTATTTCTACCCGGTCAGCGCCACGCCGTTCACCTTGCTGGGCCTGTCGCTGTCGATCTTCATGAGCTTTCGCAACAACGCGTGCTATGACCGTTGGTGGGAGGGGCGCAAGGCCTGGGGCAAACTGATTATCGAAACCCGCTCGTTCGTGCGGGAAAGCGTGGTGATTACCGACGAAACACTGCGCACCGAACTGCTGCGCAACCTGTGTGGTTTCGCCCATGCACTCAATGCCCGGCTGCGTAACGAGGACGAACGGGCAGCTGCCAGGCCGTGGCTGGCCCGGGCCGAT  >  NC_002947/5218092‑5218384
                                                                                                                                                    |                                                                                                                                                 
cGTATTTCTACCCGGTCAGCGCCACGCCGTTCACCTTGCTGGGCCTGTCGCTGTCGATCTTCATGAGCTTTCGCAACAACGCGTGCTATGACCGTTGGTGGGAGGGGCGCAAGGCCTGGGGCAAACTGATTATCGAAACCCGCtcgttc                                                                                                                                                  <  2:220614/149‑1 (MQ=255)
    tttCTACCCGGTCAGCGCCACGCCGTTCACCTTGCTGGGCCTGTCGCTGTCGATCTTAATGAGCTTTCGCAACAACGCGTGCTATGACCGTTGGTGGGAGGGGCGCAAGGCCTGGGGCAAACTGATTATCGAAACCCGCTCGTTCGTGc                                                                                                                                              >  1:104616/1‑149 (MQ=255)
      gcTACCCGCTCAGCGCCACGACGTTCACCTTGCTGGGCCTGTCGCTGTCGATCTTCATGAGCTTTCGCAACAACGCGTGCTATGACCGTTGGTGGGAGGGGCGCAAGGCCTGGGGCAAACTGATTATCGAAACCCGCTCGTTCGTGCg                                                                                                                                             <  4:202948/147‑1 (MQ=255)
       cTACCCGGTCAGCGCCACGCCGTTCACCTTGCTGGGCCTGTCGCTGTCGATCTTCATGAGCTTTCGCAACAACGCGTGCTATGACCGTTGGTGGGAGGGGCGCAAGGCCTGGGGCAAACTGATTATCGAAACCCGCTCGTTCGTGCgg                                                                                                                                            <  2:262810/148‑1 (MQ=255)
          cccGGTCAGCGCCACGCCGTTCACCTTGCTGGGCCTGTCGCTGTCGATCTTCATGAGCTTTCGCAACAACGCGTGCTATGACCGTTGGTGGGAGGGGCGCAAGGCCTGGGGCAAACTGATTATCGAAACCCGCTCGTTCGTGCGGGa                                                                                                                                          >  8:267084/1‑147 (MQ=255)
                 agggcCACGCCGTTCACCTTGCTGGGCGTGTCGCGGTCGTGCTTCATGAGCTTTCGCAACAACGCGTCCTATGAGCGTTGGTGGGGGGGGCGCAAGGCCTGTGGCAAACTGATTATCGAAACCCGCTCGTTCGTGCGGGAAAGCgtggt                                                                                                                                 <  5:127084/146‑1 (MQ=255)
                   cgcCACGCCGTTCACCTTGCTGGGCCTGTCGCTGTCGATCTTCATGAGCTTTCGCAACAACGCGTGCTATGACCGTTGGTGGGAGGGGCGCAAGGCCTGGGGCAAACTGATTATCGAAACCCGCTCGTTCGTGCGGGAAAGCGTGGTGa                                                                                                                               <  2:171572/149‑1 (MQ=255)
                        cGCCGTTCACCTTGCTGGGCCTGTCGCTGTCGATCTTCATGAGCTTTCGCAACAACGCGTGCTATGACCGTTGGTGGGAGGGGCGCAAGGCCTGGGGCAAACTGATTATCGAAACCCGCTCGTTCGTGCGGGAAAGCGTGGTGATTAcc                                                                                                                          <  5:86230/149‑1 (MQ=255)
                         gCCGTTCACCTTGCTGGGCCTGTCGCTGTCGATCTTCATGAGCTTTCGCAACAACGCGTGCTATGACCGTTGGTGGGAGGGGCGCAAGGCCTGGGGCAAACTGATTATCGAAACCCGCTCGTTCGTGCGGGAAAGCGTGGTGATTACcg                                                                                                                         <  1:174210/149‑1 (MQ=255)
                            gTTCACCTTGCTGGGCCTGTCGCTGTCGATCTTCATGAGCTTTCGCAACAACGCGTGCTATGACCGTTGGTGGGAGGGGCGCAAGGCCTGGGGCAAACTGATTATCGAAACCCGCACGAACGAGCggg                                                                                                                                           <  7:322921/128‑1 (MQ=255)
                            gTTCACCTTGCTGGGCCTGTCGCTGTCGATCTTCATGAGCTTTCGCAACAACGCGTGCTATGACCGTTGGTGGGAGGGGCGCAAGGCCTGGGGCAAACTGATTATCGAAACCCGCACGAACGAGCggg                                                                                                                                           >  8:322921/1‑128 (MQ=255)
                                   ttGCTGGGCCTGTCGCTGTCGATCTTCATGAGCTTTCGCAACAACGCGTGCTATGACCGTTGGTGGGAGGGGCGCAAGGCCTGGGGCAAACTGATTATCGAAACCCGCTCGTTCGTGCGGGaaa                                                                                                                                        <  1:88322/124‑1 (MQ=255)
                                   ttGCTGGGCCTGTCGCTGTCGATCTTCATGAGCTTTCGCAACAACGCGTGCTATGACCGTTGGTGGGAGGGGCGCAAGGCCTGGGGCAAACTGATTATCGAAACCCGCTCGTTCGTGCGGGaaa                                                                                                                                        >  2:88322/1‑124 (MQ=255)
                                       tGGGCCTGTCGCTGTCGATCTTCATGAGCTTTCGCAACAACGCGTGCTATGACCGTTGGTGGGAGGGGCGCAAGGCCTGGGGCAAACTGATTATCGAAACCCGCTCGTTCGTGCGGGAAAGCGTGGTGATTACCGACGAAACACTgcg                                                                                                            >  3:198421/1‑148 (MQ=255)
                                                     tcgATCTTCATGAGCTTTCGCAACAACTCGTGCTATGACCGTTGGTGGGAGGGGCGCAAGGCCTGGGGCAAACTGATTATCGAAACCCGCTCGTTCGTGCGGGAAAGCGTGGTGATTACCGACGAAACACTGCGCACCGAACtgctgc                                                                                              <  2:89993/148‑1 (MQ=255)
                                                     tcgATCTTCATGAGCTTTCGCAACAACGCGTGCTATGACCGTTGGTGGGAGGGGCGCAAGGCCTGGGGCAAACTGATTATCGAAACCCGCTCGTTCGTGCGGGAAAGCGTGGTGATTACCGACGAAaca                                                                                                                 <  8:219036/129‑1 (MQ=255)
                                                     tcgATCTTCATGAGCTTTCGCAACAACGCGTGCTATGACCGTTGGTGGGAGGGGCGCAAGGCCTGGGGCAAACTGATTATCGAAACCCGCTCGTTCGTGCGGGAAAGCGTGGTGATTACCGACGAAaca                                                                                                                 >  7:219036/1‑129 (MQ=255)
                                                     tcgATCTTCATGAGCTTTCGCAACAACGCGTGCTATGACCGTTGGTGGGAGGGGCGCAAGGCCTGGGGCAAACTGATTATCGAAACCCGCTCGTTCGTGCGGGAAAGCGTGGTGATTACCGACGAAACACTGCGCACCGAACTGCTgcg                                                                                             >  1:231619/1‑149 (MQ=255)
                                                       gATCTTCATGAGCTTTCGCAACAACGCGTGCTATGACCGTTGGTGGGAGGGGCGCAAGGCCTGAGGGCAAACTTAATAAACAAAACCGCACGTTCGTGCGGGAAAGCGTGGTGATTACCTACGAAAAACTGCGCCCCCAActgcatcc                                                                                             >  1:303114/1‑144 (MQ=255)
                                                              aTGAGCTTTCGCAACAACGCGTGCTATGACCGTTGGTGGGAGGGGCGCAAGGCCTGGGGCAAACTGATTATCGAAACCCGCGCGTTCGTGCGGGAAAGCGTGGTGATTACCGACGAAACACTGCGCACCGAACTGCTGCGCAACCtgt                                                                                     <  8:25248/148‑1 (MQ=255)
                                                               tGAGCTTTCGCAACAACGCGTGCTATGACCGTTGGTGGGAGGGGCGCAAGGCCTGGGGCAAACTGATTATCGAAACCCGCTCGTTCGTGCGGGAAAGCGTGGTGATTACCGACGAAACACTGCGCACCGAACTGCTGCGCAACCtgtgt                                                                                   <  4:362999/149‑1 (MQ=255)
                                                               tGAGCTTTCGCAACAACGCGTGCTATGACCGTTGGTGGGAGGGGCGCAAGGCCTGGGGCAAACTGATTATCGAAACCCGCTCGTTCGTGCGGGAAAGCGTGGTGATTACCGACGAAACACTGCGCACCGAACTGCTGCGCAACCtgtgt                                                                                   <  4:198421/149‑1 (MQ=255)
                                                               tGAGCTTTCGCAACAACGCGTGCTATGACCGTTGGTGGGAGGGGCGCAAGGCCTGGGGCAAACTGATTATCGAAACCCGCACGAACGAGCggg                                                                                                                                           >  3:29700/1‑93 (MQ=255)
                                                               tGAGCTTTCGCAACAACGCGTGCTATGACCGTTGGTGGGAGGGGCGCAAGGCCTGGGGCAAACTGATTATCGAAACCCGCACGAACGAGCggg                                                                                                                                           <  4:29700/93‑1 (MQ=255)
                                                                gAGCTTTCGCAACAACGCGTGCTATGACCGTTGGTGGGAGGGGCGCAAGGCCTGGGGCAAACTGATTATCGAAACCCGCTCGTTCGTGCGGGAAAGCGTGGTGATTACCGACGAAACACTGCGCACCGAACTGCTGCGCAACCtgtgtg                                                                                  >  5:24671/1‑149 (MQ=255)
                                                                 aGCTTTCGCAACAACGCGTGCTATGACCGTTGGTGGGAGGGGCGCAAGGCCTGGGGCAAACTGATTATCGAAACCCGCTCGTTCGTGCGGGAAAGCGTGGTGATTACCGACGAAACACTGCGCACCGAACTGCTGCGCAACCtgtgtg                                                                                  <  6:24671/148‑1 (MQ=255)
                                                                 aGCTTTCGCAACAACGCGTGCTATGACCGTTGGTGGGAGGGGCGCAAGGCCTGGGGCAAACTGATTATCGAAACCCGCTCGTTCGTGCGGGAAAGCGTGGTGATTACCGACGAAACACTGCGCACCGAACTGCTGCGCAACCtgtgtg                                                                                  <  1:87423/148‑1 (MQ=255)
                                                                      tCGCAACAACGCGTGCTATGACCGTTGGTGGGAGGGGCGCAAGGCCTGGGGCAAACTGATTATCGAAACCCGCTCGTTCGTGCGGGAAAGCGTGGTGATTACcgacga                                                                                                                     >  4:328535/1‑108 (MQ=255)
                                                                      tCGCAACAACGCGTGCTATGACCGTTGGTGGGAGGGGCGCAAGGCCTGGGGCAAACTGATTATCGAAACCCGCTCGTTCGTGCGGGAAAGCGTGGTGATTACcgacga                                                                                                                     <  3:328535/108‑1 (MQ=255)
                                                                           acaacGCGTGCTATGACCGTTGGTGGGAGGGGCGCAAGGCCTGGGGCAAACTGATTATCGAAACCCGCTCGTTCGTGCGGGAAAGCgt                                                                                                                                    >  4:10077/1‑88 (MQ=255)
                                                                           acaacGCGTGCTATGACCGTTGGTGGGAGGGGCGCAAGGCCTGGGGCAAACTGATTATCGAAACCCGCTCGTTCGTGCGGGAAAGCgt                                                                                                                                    <  3:10077/88‑1 (MQ=255)
                                                                            caacGCGTGCTATGACCGTTGGTGGGAGGGGCGCAAGGCCTGGGGCAAACTGATTATCGAAACCCGCTCGTTCGTGCGGGAAAGCGTGGTGATTACCGACGAAACACTGCGCACCGAACTGCTGCGCAACCTGTGTGGTTTCGCCCATg                                                                      >  4:186406/1‑149 (MQ=255)
                                                                                     cTATGACCGTTGGTGGGAGGGGCGCAAGGCCTGGGGCAAACAGATTATCGAAACCCGCTCGTTCGTGCGGGAAAGCGTGGTGATTACCGACGAAACACTGCGCACCGAACTGCTGCGCAAc                                                                                         >  5:87722/1‑121 (MQ=255)
                                                                                     cTATGACCGTTGGTGGGAGGGGCGCAAGGCCTGGGGCAAACAGATTATCGAAACCCGCTCGTTCGTGCGGGAAAGCGTGGTGATTACCGACGAAACACTGCGCACCGAACTGCTGCGCAAc                                                                                         <  6:87722/121‑1 (MQ=255)
                                                                                      tATGACCGTTGGTGGGAGGGGCGCAAGGCCTGGGGCAAACTGATTATCGAAACCCGCTCGTTCGTGCGGGAAAGCGTGTGATTACCGACGAAACACTACGCACCGAACTGCTGCGCAACCTGTGTGGt                                                                                <  3:79101/128‑1 (MQ=255)
                                                                                      tATGACCGTTGGTGGGAGGGGCGCAAGGCCTGGGGCAAACTGATTATCGAAACCCGCTCGTTCGTGCGGGAAAGCGTGTGATTACCGACGAAACACTACGCACCGAACTGCTGCGCAACCTGTGTGGt                                                                                >  4:79101/1‑128 (MQ=255)
                                                                                         gACCGTTGGTGGGAGGGGCGCAAGGCCTGGGGCAAACTGATTATCGAAACCCGCTCGTTCGTGCGGGAAAGCGTGGTGATTACCGACGAAACACTGCGCACCGAACTGCTGCGCAACCAGTGTGGTTTCGCCCATGCACTCAATGCCCg                                                         >  6:208558/1‑149 (MQ=255)
                                                                                            cGTTGGTGGGAGGGGCGCAAGGCCTGGGGCAAACTGATTATCGAAACCCGCTCGTTCGTGCGGGAAAGCGTGGTGATTACCGACGAAACACTGCGCACCGAACTGCTGCGCAACCTGTGTGGTTTCg                                                                            <  3:287064/127‑1 (MQ=255)
                                                                                            cGTTGGTGGGAGGGGCGCAAGGCCTGGGGCAAACTGATTATCGAAACCCGCTCGTTCGTGCGGGAAAGCGTGGTGATTACCGACGAAACACTGCGCACCGAACTGCTGCGCAACCTGTGTGGTTTCg                                                                            >  4:287064/1‑127 (MQ=255)
                                                                                             gTTGGTGGGAGGGGCGCAAGGCCTGGGGCAAACTGATTATCGAAACCCGCTCGTTCGTGCGGGAAAGCGTGGTGATTACCGACGAAACACTGCGCACCGAACTGCTGCACAACCtgtgt                                                                                   >  4:193189/1‑119 (MQ=255)
                                                                                             gTTGGTGGGAGGGGCGCAAGGCCTGGGGCAAACTGATTATCGAAACCCGCTCGTTCGTGCGGGAAAGCGTGGTGATTACCGACGAAACACTGCGCACCGAACTGCTGCACAACCtgtgt                                                                                   <  3:193189/119‑1 (MQ=255)
                                                                                                  tggGAGGGGCGCAAGGCCTGGGGCAAACTGATTATCGAAACCCGCTCGTTCGTGCGGGAAAGCGTGGTGATTACCGACGAAACACTGCGCACCGAACTGCTGCGCAACCTGTGTGGTTTCGCCCATGCACTCAATGCCCGGCTGCGTaa                                                <  2:303114/149‑1 (MQ=255)
                                                                                                    ggAGGGGCGCAAGGCCTGGGGCAAACTGATTATCGAAACCCGCTCGTTCGTGCGGGAAAGCGTGGTGATTACCGACGAAACACTGCGCACCGAACTGCTGCGCAACCTGTGTGGTTTCGCCCATGCACTCAATGCCCGGCTGCGTAAc                                               <  2:104616/148‑1 (MQ=255)
                                                                                                      aGGGGCGCAAGGCCTGGGGCAAACTGATTATCGAAACCCGCTCGTTCGTGCGGGAAAGCGTGGTGATTACCGACGAAACACTGCGCACCGAACTGCTGCGCAACCTGTGTGGTTTCGCCCATGCACTCAATGCCCGGCTGCGTAACGAg                                            >  6:219196/1‑149 (MQ=255)
                                                                                                       ggggCGCAAGGCCTGGGGCAAACTGATTATCGAAACCCGCTCGTTCGTGCGGGAAAGCGTGGTGATTACCGACGAAACACTGCGCACCGAACTGCTGCGCAACCTGTGTGGTTTCGCCCATGCACTCAATGCCCGGCTGCGTAACGAgg                                           <  2:231619/149‑1 (MQ=255)
                                                                                                       ggggCGCAAGGCCTGGGGCAAACTGATTATCGAAACCCGCTCGTTCGTGCGGGAAAGCGTGGTGATTACCGACGAAACACTGCGCACCGAACTGCTGCGCAACCTGTGTGGTTTCGCCCATGCACTCAATGCCCGGCTGCGTAACGAg                                            >  1:230724/1‑148 (MQ=255)
                                                                                                        gggCGCAAGGCCTGGGGCAAACTGATTATCGAAACCCGCTCGTTCGTGCGGGAAAGCGTGGTGATTACCGACGAAACACTGCGCACCGAACTGCTGCGCAACCTGTGTGGTTTCGCCCATGCACTCAATGCCCGGCTGCGTAACGAGGa                                          >  7:347551/1‑149 (MQ=255)
                                                                                                         ggCGCAAGGCCTGGGGAAAACTTATTATCGAAACCCGCTCGTTCGTGCGGGAAAGTGTGGTGATTACCGACGAAACACTGCGCACCTAACTGCTGCGCAACCTGTGTGGTTTCg                                                                            <  7:56870/114‑1 (MQ=255)
                                                                                                         ggCGCAAGGCCTAAAACAAACTGATTATCGAAACCCACTCATTCATACAGGAAAGCGTGGTAATTACAGACAAAACAATGCACACCGAACTGCTGCGAAAACTGTGTGGTTTCg                                                                            >  8:56870/1‑114 (MQ=255)
                                                                                                          gcgcAAGGCCTGGGGCAAACTGATTATCGAAACCCGCTCGTTCGTGCGGGAAAGCGTGGTGATTACCGACGAAACACTGCGCACCGAACTGCTGCGCAACCTGTGTGGTTTCGCCCATGCACTCAATGCCCGGCTGCGTAACGAGGAc                                         >  4:226792/1‑148 (MQ=255)
                                                                                                                ggCCTGGGGCAAACTGATTATCGAAACCCGCTCGTTCGTGCGGGAAAGCGTGGTGATTACCGACGAAACACTGCGCACCGAACtgctgc                                                                                              <  3:185267/89‑1 (MQ=255)
                                                                                                                ggCCTGGGGCAAACTGATTATCGAAACCCGCTCGTTCGTGCGGGAAAGCGTGGTGATTACCGACGAAACACTGCGCACCGAACtgctgc                                                                                              >  4:185267/1‑89 (MQ=255)
                                                                                                                  ccTGGGGCAAACTGATTATCGAAACCCGCTCGTTCGTGCGGGAAAGCGTGGTGATTACCGACGAAACACTGCGCACCGAACTGCTGCGCAACCTGTGTGGTTTCGCCCATGCACTCAATGCCCGGCTGCGTAACGAGGACGAACggg                                  >  6:19793/1‑147 (MQ=255)
                                                                                                                   cTGGGGCAAACTGATTATCGAAACCCGCTCGTTCGTGCGGGAAAGCGTGGTGATTACCGACGAAACACTGCGCACCGAACTGCTGCGCAACCTGTGTGGTTTCGCCCATGCACTCAATGCCCGGCTGCGTAACGAGGACGAACGGGCAg                               <  2:230724/149‑1 (MQ=255)
                                                                                                                     ggggCAAACTGATTATCGAAACCCGCTCGTTCGTGCGGGAAAGCGTGGTGATTACCGACGAAACACTGCGCACCGAACTGCTGCGCAACCTGTGTGGTTTCGCCCATGCACTCAATGCCCGGCTGCGTAACGAGGACGAACGGGCAGCt                             <  7:267084/149‑1 (MQ=255)
                                                                                                                             aCTGATTATCGAAACCCGCTCGTTCGTGCGGGAAAGCGTGGTGATTACCGACGAAACACTGCGCACCGAACTGCTGCGCAACCTGTGTGGTTTCGCCCATGCACTCAATGCCCGGCTGCGTAACGAGGACGAACGGGCAGCTGc                           <  7:97557/144‑1 (MQ=255)
                                                                                                                             aCTGATTATCGAAACCCGCTCGTTCGTGCGGGAAAGCGTGGTGATTACCGACGAAACACTGCGCACCGAACTGCTGCGCAACCTGTGTGGTTTCGCCCATGCACTCAATGCCCGGCTGCGTAACGAGGACGAACGGGCAGCTGc                           >  8:97557/1‑144 (MQ=255)
                                                                                                                               tGATTATCGAAACCCGCTCGTTCGTGCGGGAAAGCGTGGTGATTACCGACGAAACACTGCGCACCGAACTGCTGCGCAACCTGTGTGGTTTCGCCCATGCACTCAATGccc                                                          >  2:213905/1‑111 (MQ=255)
                                                                                                                               tGATTATCGAAACCCGCTCGTTCGTGCGGGAAAGCGTGGTGATTACCGACGAAACACTGCGCACCGAACTGCTGCGCAACCTCTGTGGTTTCGCCCATGCACTCAATGccc                                                          <  1:213905/111‑1 (MQ=255)
                                                                                                                                gATTATCGAAACCCGCTCGTTCGTGCGGGAAAGCGTGGTGATTACCGACGAAACACTGCGCACCGAACTGCTGCGCAACCTGTGTGGTTTCGCCCATGCACTCAATGCCCGGCTGCGTAACGAGGACGAACGGGCAGCTGc                           >  2:330449/1‑141 (MQ=255)
                                                                                                                                gATTATCGAAACCCGCTCGTTCGTGCGGGAAAGCGTGGTGATTACCGACGAAACACTGCGCACCGAACTGCTGCGCAACCTGTGTGGTTTCGCCCATGCACTCAATGCCCGGCTGCGTAACGAGGACGAACGGGCAGCTGc                           <  1:330449/141‑1 (MQ=255)
                                                                                                                                     tCGAAACCCGCTCGTTCGTGCGGGAAAGCGTGGTGATTACCGACGaaa                                                                                                                   >  7:148988/1‑48 (MQ=255)
                                                                                                                                     tCGAAACCCGCTCGTTCGTGCGGGAAAGCGTGGTGATTACCGACGaaa                                                                                                                   <  8:148988/48‑1 (MQ=255)
                                                                                                                                     tCGAAACCCGCTCGTTCGTGCGGGAAAGCGTGGTGATTACCGACGAAACACTGCGCACCGAACTGCTGCGCAACCTGTGTGGTTTCGCCCATGCACTCAATGCCCGGCTGCGTAACGAGGACGAACGGGCAGCTGCCAGGCCGtggct               <  3:186406/148‑1 (MQ=255)
                                                                                                                                     cgctttcccGCACGAACGAGCGGGAAAGCGTGGTGATTACCGACGAAACACTGCGCACCGAACTGCTGCGCAACCTGTGTGGTTTCGCCCATGCACTCAATGCCCGGCTGCGTAACGAGGACGAACGGGCAGCTGCCAGGCCGtggctg              >  5:179792/7‑149 (MQ=255)
                                                                                                                                               cTCGTTCGTGCGGGAAAGCGTGGTGATTACCGACGAAACACTGCGCACCGAACTGCTGCGCAACCTGTGTGGTTTCGCCCATGCACTCAATGCCCGGCTGCGTAACGAGGACGAACGGGCAGCTGCCAGGCCGTGGCTGGCCCGGGcc     >  7:165070/1‑148 (MQ=255)
                                                                                                                                                 cgttcgtGCGGGAAAGCGTGGTGATTACCGACGAAACACTGCGCACCGAACTGCTGCGCAACCTGTGTGGTTTCGCCCATGCACTCAATGCCCGGCTGCGTAACGAGGACGAACGGGCAGCTGCCAGGCCGTGGCTGGCCCGGGCCGAt  <  3:226792/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
CGTATTTCTACCCGGTCAGCGCCACGCCGTTCACCTTGCTGGGCCTGTCGCTGTCGATCTTCATGAGCTTTCGCAACAACGCGTGCTATGACCGTTGGTGGGAGGGGCGCAAGGCCTGGGGCAAACTGATTATCGAAACCCGCTCGTTCGTGCGGGAAAGCGTGGTGATTACCGACGAAACACTGCGCACCGAACTGCTGCGCAACCTGTGTGGTTTCGCCCATGCACTCAATGCCCGGCTGCGTAACGAGGACGAACGGGCAGCTGCCAGGCCGTGGCTGGCCCGGGCCGAT  >  NC_002947/5218092‑5218384

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