Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,218,235 T→A 8.7% S95T (TCG→ACG)  PP_4598 → membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,218,2350TA 8.7% 180.7 / 5.2 65S95T (TCG→ACG) PP_4598membrane protein
Reads supporting (aligned to +/- strand):  ref base T (27/31);  new base A (4/2);  total (31/34)
Fisher's exact test for biased strand distribution p-value = 4.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGGCGTATTTCTACCCGGTCAGCGCCACGCCGTTCACCTTGCTGGGCCTGTCGCTGTCGATCTTCATGAGCTTTCGCAACAACGCGTGCTATGACCGTTGGTGGGAGGGGCGCAAGGCCTGGGGCAAACTGATTATCGAAACCCGCTCGTTCGTGCGGGAAAGCGTGGTGATTACCGACGAAACACTGCGCACCGAACTGCTGCGCAACCTGTGTGGTTTCGCCCATGCACTCAATGCCCGGCTGCGTAACGAGGACGAACGGGCAGCTGCCAGGCCGTGGCTGGCCCGGGCC  >  NC_002947/5218089‑5218381
                                                                                                                                                   |                                                                                                                                                  
cGGCGTATTTCTACCCGGTCAGCGCCACGCCGTTCACCTTGCTGGGCCTGTCGCTGTCGATCTTCATGAGCTTTCGCAACAACGCGTGCTATGACCGTTGGTGGGAGGGGCGCAAGGCCTGGGGCAAACTGATTATCGAAACCCGCtcg                                                                                                                                                  >  2:87423/1‑149 (MQ=255)
 ggCGTATTTCTACCCGGTCAGCGCCACGCCGTTCACCTTGCTGGGCCTGTCGCTGTCGATCTTCATGAGCTTTCGCAACAACGCGTGCTATGACCGTTGGTGGGAGGGGCGCAAGGCCTGGGGCAAACTGATTATCGAAACCCGCtc                                                                                                                                                   <  8:225711/147‑1 (MQ=255)
 ggCGTATTTCTACCCGGTCAGCGCCACGCCGTTCACCTTGCTGGGCCTGTCGCTGTCGATCTTCATGAGCTTTCGCAACAACGCGTGCTATGACCGTTGGTGGGAGGGGCGCAAGGCCTGGGGCAAACTGATTATCGAAACCCGCtc                                                                                                                                                   >  7:225711/1‑147 (MQ=255)
   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)
                        ccACGCCGTTCACCTTGCTGGGCCTGTCGCTGTCGATCTTCATGAGCTTTCGCAACAACGCGTGCTATGACCGTTGGTGGGAGGGGCGCAAGGCCTGGGGCAAACTGATTATCGAAACCCGCtc                                                                                                                                                   >  1:366276/1‑124 (MQ=255)
                        ccACGCCGTTCACCTTGCTGGGCCTGTCGCTGTCGATCTTCATGAGCTTTCGCAACAACGCGTGCTATGACCGTTGGTGGGAGGGGCGCAAGGCCTGGGGCAAACTGATTATCGAAACCCGCtc                                                                                                                                                   <  2:366276/124‑1 (MQ=255)
                           cGCCGTTCACCTTGCTGGGCCTGTCGCTGTCGATCTTCATGAGCTTTCGCAACAACGCGTGCTATGACCGTTGGTGGGAGGGGCGCAAGGCCTGGGGCAAACTGATTATCGAAACCCGCTCGTTCGTGCGGGAAAGCGTGGTGATTAcc                                                                                                                       <  5:86230/149‑1 (MQ=255)
                            gCCGTTCACCTTGCTGGGCCTGTCGCTGTCGATCTTCATGAGCTTTCGCAACAACGCGTGCTATGACCGTTGGTGGGAGGGGCGCAAGGCCTGGGGCAAACTGATTATCGAAACCCGCTCGTTCGTGCGGGAAAGCGTGGTGATTACcg                                                                                                                      <  1:174210/149‑1 (MQ=255)
                               gTTCACCTTGCTGGGCCTGTCGCTGTCGATCTTCATGAGCTTTCGCAACAACGCGTGCTATGACCGTTGGTGGGAGGGGCGCAAGGCCTGGGGCAAACTGATTATCGAAACCCGCACGAACGAGCggg                                                                                                                                        >  8:322921/1‑128 (MQ=255)
                               gTTCACCTTGCTGGGCCTGTCGCTGTCGATCTTCATGAGCTTTCGCAACAACGCGTGCTATGACCGTTGGTGGGAGGGGCGCAAGGCCTGGGGCAAACTGATTATCGAAACCCGCACGAACGAGCggg                                                                                                                                        <  7:322921/128‑1 (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                                                                                                                  <  3:328535/108‑1 (MQ=255)
                                                                         tCGCAACAACGCGTGCTATGACCGTTGGTGGGAGGGGCGCAAGGCCTGGGGCAAACTGATTATCGAAACCCGCTCGTTCGTGCGGGAAAGCGTGGTGATTACcgacga                                                                                                                  >  4:328535/1‑108 (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                                                                                           >  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                        <  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                        <  1:330449/141‑1 (MQ=255)
                                                                                                                                   gATTATCGAAACCCGCTCGTTCGTGCGGGAAAGCGTGGTGATTACCGACGAAACACTGCGCACCGAACTGCTGCGCAACCTGTGTGGTTTCGCCCATGCACTCAATGCCCGGCTGCGTAACGAGGACGAACGGGCAGCTGc                        >  2:330449/1‑141 (MQ=255)
                                                                                                                                        tCGAAACCCGCTCGTTCGTGCGGGAAAGCGTGGTGATTACCGACGaaa                                                                                                                <  8:148988/48‑1 (MQ=255)
                                                                                                                                        tCGAAACCCGCTCGTTCGTGCGGGAAAGCGTGGTGATTACCGACGaaa                                                                                                                >  7:148988/1‑48 (MQ=255)
                                                                                                                                        tCGAAACCCGCTCGTTCGTGCGGGAAAGCGTGGTGATTACCGACGAAACACTGCGCACCGAACTGCTGCGCAACCTGTGTGGTTTCGCCCATGCACTCAATGCCCGGCTGCGTAACGAGGACGAACGGGCAGCTGCCAGGCCGtggct            <  3:186406/148‑1 (MQ=255)
                                                                                                                                        cgctttcccGCACGAACGAGCGGGAAAGCGTGGTGATTACCGACGAAACACTGCGCACCGAACTGCTGCGCAACCTGTGTGGTTTCGCCCATGCACTCAATGCCCGGCTGCGTAACGAGGACGAACGGGCAGCTGCCAGGCCGtggctg           >  5:179792/7‑149 (MQ=255)
                                                                                                                                                  cTCGTTCGTGCGGGAAAGCGTGGTGATTACCGACGAAACACTGCGCACCGAACTGCTGCGCAACCTGTGTGGTTTCGCCCATGCACTCAATGCCCGGCTGCGTAACGAGGACGAACGGGCAGCTGCCAGGCCGTGGCTGGCCCGGGcc  >  7:165070/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CGGCGTATTTCTACCCGGTCAGCGCCACGCCGTTCACCTTGCTGGGCCTGTCGCTGTCGATCTTCATGAGCTTTCGCAACAACGCGTGCTATGACCGTTGGTGGGAGGGGCGCAAGGCCTGGGGCAAACTGATTATCGAAACCCGCTCGTTCGTGCGGGAAAGCGTGGTGATTACCGACGAAACACTGCGCACCGAACTGCTGCGCAACCTGTGTGGTTTCGCCCATGCACTCAATGCCCGGCTGCGTAACGAGGACGAACGGGCAGCTGCCAGGCCGTGGCTGGCCCGGGCC  >  NC_002947/5218089‑5218381

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: