Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,218,242 T→A 7.7% V97E (GTG→GAG)  PP_4598 → membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,218,2420TA 7.7% 188.3 / 4.9 65V97E (GTG→GAG) 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.60e-01

TTTCTACCCGGTCAGCGCCACGCCGTTCACCTTGCTGGGCCTGTCGCTGTCGATCTTCATGAGCTTTCGCAACAACGCGTGCTATGACCGTTGGTGGGAGGGGCGCAAGGCCTGGGGCAAACTGATTATCGAAACCCGCTCGTTCGTGCGGGAAAGCGTGGTGATTACCGACGAAACACTGCGCACCGAACTGCTGCGCAACCTGTGTGGTTTCGCCCATGCACTCAATGCCCGGCTGCGTAACGAGGACGAACGGGCAGCTGCCAGGCCGTGGCTGGCCCGGGCCGATGCGATC  >  NC_002947/5218096‑5218390
                                                                                                                                                   |                                                                                                                                                    
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                                                                                                                       >  7:219036/1‑129 (MQ=255)
                                                 tcgATCTTCATGAGCTTTCGCAACAACGCGTGCTATGACCGTTGGTGGGAGGGGCGCAAGGCCTGGGGCAAACTGATTATCGAAACCCGCTCGTTCGTGCGGGAAAGCGTGGTGATTACCGACGAAaca                                                                                                                       <  8:219036/129‑1 (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                                                                                                                                                 <  4:29700/93‑1 (MQ=255)
                                                           tGAGCTTTCGCAACAACGCGTGCTATGACCGTTGGTGGGAGGGGCGCAAGGCCTGGGGCAAACTGATTATCGAAACCCGCACGAACGAGCggg                                                                                                                                                 >  3:29700/1‑93 (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                                                                                  >  4:287064/1‑127 (MQ=255)
                                                                                        cGTTGGTGGGAGGGGCGCAAGGCCTGGGGCAAACTGATTATCGAAACCCGCTCGTTCGTGCGGGAAAGCGTGGTGATTACCGACGAAACACTGCGCACCGAACTGCTGCGCAACCTGTGTGGTTTCg                                                                                  <  3:287064/127‑1 (MQ=255)
                                                                                         gTTGGTGGGAGGGGCGCAAGGCCTGGGGCAAACTGATTATCGAAACCCGCTCGTTCGTGCGGGAAAGCGTGGTGATTACCGACGAAACACTGCGCACCGAACTGCTGCACAACCtgtgt                                                                                         <  3:193189/119‑1 (MQ=255)
                                                                                         gTTGGTGGGAGGGGCGCAAGGCCTGGGGCAAACTGATTATCGAAACCCGCTCGTTCGTGCGGGAAAGCGTGGTGATTACCGACGAAACACTGCGCACCGAACTGCTGCACAACCtgtgt                                                                                         >  4:193189/1‑119 (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                                                                                                    >  4:185267/1‑89 (MQ=255)
                                                                                                            ggCCTGGGGCAAACTGATTATCGAAACCCGCTCGTTCGTGCGGGAAAGCGTGGTGATTACCGACGAAACACTGCGCACCGAACtgctgc                                                                                                    <  3:185267/89‑1 (MQ=255)
                                                                                                              ccTGGGGCAAACTGATTATCGAAACCCGCTCGTTCGTGCGGGAAAGCGTGGTGATTACCGACGAAACACTGCGCACCGAACTGCTGCGCAACCTGTGTGGTTTCGCCCATGCACTCAATGCCCGGCTGCGTAACGAGGACGAACggg                                        >  6:19793/1‑147 (MQ=255)
                                                                                                               cTGGGGCAAACTGATTATCGAAACCCGCTCGTTCGTGCGGGAAAGCGTGGTGATTACCGACGAAACACTGCGCACCGAACTGCTGCGCAACCTGTGTGGTTTCGCCCATGCACTCAATGCCCGGCTGCGTAACGAGGACGAACGGGCAg                                     <  2:230724/149‑1 (MQ=255)
                                                                                                                 ggggCAAACTGATTATCGAAACCCGCTCGTTCGTGCGGGAAAGCGTGGTGATTACCGACGAAACACTGCGCACCGAACTGCTGCGCAACCTGTGTGGTTTCGCCCATGCACTCAATGCCCGGCTGCGTAACGAGGACGAACGGGCAGCt                                   <  7:267084/149‑1 (MQ=255)
                                                                                                                         aCTGATTATCGAAACCCGCTCGTTCGTGCGGGAAAGCGTGGTGATTACCGACGAAACACTGCGCACCGAACTGCTGCGCAACCTGTGTGGTTTCGCCCATGCACTCAATGCCCGGCTGCGTAACGAGGACGAACGGGCAGCTGc                                 >  8:97557/1‑144 (MQ=255)
                                                                                                                         aCTGATTATCGAAACCCGCTCGTTCGTGCGGGAAAGCGTGGTGATTACCGACGAAACACTGCGCACCGAACTGCTGCGCAACCTGTGTGGTTTCGCCCATGCACTCAATGCCCGGCTGCGTAACGAGGACGAACGGGCAGCTGc                                 <  7:97557/144‑1 (MQ=255)
                                                                                                                           tGATTATCGAAACCCGCTCGTTCGTGCGGGAAAGCGTGGTGATTACCGACGAAACACTGCGCACCGAACTGCTGCGCAACCTGTGTGGTTTCGCCCATGCACTCAATGccc                                                                >  2:213905/1‑111 (MQ=255)
                                                                                                                           tGATTATCGAAACCCGCTCGTTCGTGCGGGAAAGCGTGGTGATTACCGACGAAACACTGCGCACCGAACTGCTGCGCAACCTCTGTGGTTTCGCCCATGCACTCAATGccc                                                                <  1:213905/111‑1 (MQ=255)
                                                                                                                            gATTATCGAAACCCGCTCGTTCGTGCGGGAAAGCGTGGTGATTACCGACGAAACACTGCGCACCGAACTGCTGCGCAACCTGTGTGGTTTCGCCCATGCACTCAATGCCCGGCTGCGTAACGAGGACGAACGGGCAGCTGc                                 <  1:330449/141‑1 (MQ=255)
                                                                                                                            gATTATCGAAACCCGCTCGTTCGTGCGGGAAAGCGTGGTGATTACCGACGAAACACTGCGCACCGAACTGCTGCGCAACCTGTGTGGTTTCGCCCATGCACTCAATGCCCGGCTGCGTAACGAGGACGAACGGGCAGCTGc                                 >  2:330449/1‑141 (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)
                                                                                                                                                  gtGCGGGAAAGCGTGGTGATTACCGACGAAACACTGCGCACCGAACTGCTGCGCAACCTGTGTGGTTTCGCCCAt                                                                             >  6:7498/1‑75 (MQ=255)
                                                                                                                                                  gtGCGGGAAAGCGTGGTGATTACCGACGAAACACTGCGCACCGAACTGCTGCGCAACCTGTGTGGTTTCGCCCAt                                                                             <  5:7498/75‑1 (MQ=255)
                                                                                                                                                  gtGCGGGAAAGCGTGGTGATTACCGACGAAACACTGCGCACCGAACTGCTGCGCAACCTGTGTGGTTTCGCCCATGCACTCAATGCCCGGCTGCGTAACGAGGACGAACGGGCAGCTGCCAGGCCGTGGCTGGCCCGGGCCGATGCGa    <  5:219196/148‑1 (MQ=255)
                                                                                                                                                   tGCGGGAAAGCGTGGTGATTACCGACGAAACACTGCGCACCGAACTGCTGCGCAACCTGTGTGGTTTCGCCCATGCACTCAATGCCCGGCTGCGTAACGAGGACGAACGGGCAGCTGCCAGGCCGTGGCTGGCCCGGGCCGATGCGATc  <  8:165070/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
TTTCTACCCGGTCAGCGCCACGCCGTTCACCTTGCTGGGCCTGTCGCTGTCGATCTTCATGAGCTTTCGCAACAACGCGTGCTATGACCGTTGGTGGGAGGGGCGCAAGGCCTGGGGCAAACTGATTATCGAAACCCGCTCGTTCGTGCGGGAAAGCGTGGTGATTACCGACGAAACACTGCGCACCGAACTGCTGCGCAACCTGTGTGGTTTCGCCCATGCACTCAATGCCCGGCTGCGTAACGAGGACGAACGGGCAGCTGCCAGGCCGTGGCTGGCCCGGGCCGATGCGATC  >  NC_002947/5218096‑5218390

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: