Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 966,595 A→T 8.7% L41* (TTG→TAG) ‡ PP_0828 ← membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947966,5950AT 8.7% 224.2 / 8.8 77L41* (TTG→TAG) ‡PP_0828membrane protein
Reads supporting (aligned to +/- strand):  ref base A (29/41);  new base T (4/3);  total (33/44)
Fisher's exact test for biased strand distribution p-value = 4.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCTGCGCGAACTGTTCACCAGCGCGCCCCACAGGGTGCGAGCCTCCCACCAGCGGCTGTAGGCGCTGCTGTTGCGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGT‑‑TCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTCACCACCAGGCGGTACTTGCGGGCGATGACGTCTTTCACGCGAGGGTACCTCGAAACGCGGGG  >  NC_002947/966447‑966740
                                                                                                                                                      |                                                                                                                                                 
gCTGCGCGAACTGTTCACCAGCGCGCCCCACAGGGTGCGAGCTTCCCACCAGCGGCTGTAGGCGCTGCTGTTGCGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGT‑‑TCa                                                                                                                                                   <  1:125699/149‑1 (MQ=255)
     gcgAACTGTTCACCAGCGCGCCCCACAGGGTGCGAGCCTCCCACCAGCGGCTGTAGGCGCTGCTGTTGCGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGT‑‑TCAAGAAc                                                                                                                                              <  5:177650/149‑1 (MQ=255)
      cgAACTGTTCACCAGCGCGCCCCACAGGGTGCGAGCCTCCCACCAGCGGCTGTAGGCGCTGCTGTTGCGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGT‑‑TCAAGAACa                                                                                                                                             <  1:55122/149‑1 (MQ=255)
           tGTTCACCAGCGCGCCCCACAGGGTGCGAGCCTCCCACCAGCGGCTGTAGGCGCTGCTGTTGCGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGT‑‑TCTTGAACAGc                                                                                                                                           >  8:20099/1‑146 (MQ=255)
             ttCACCAGCGCGCCCCACAGGGTGCGAGCCTCCCACCAGCGGCTGTAGGCGCTGCTGTTGCGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGT‑‑TCTTGAACAGCaa                                                                                                                                         <  7:20099/146‑2 (MQ=255)
                aCCAGCGCGCCCCACAGGGTGCGAGCCTCCCACCAGCGGCTGTAGGCGCTGCTGTTGCGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGT‑‑TCAAGAACAGCATGAAATc                                                                                                                                   >  1:1296/1‑149 (MQ=255)
                 ccAGCGCGCCCCACAGGGTGCGAGCCTCCCACCAGCGGCTGTAGGCGCTGCTGTTGCGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGT‑‑TCAAGAACAGCATGAAATc                                                                                                                                   >  5:332724/1‑148 (MQ=255)
                 ccAGCGCGCCCCACAGGGTGCGAGCCTCCCACCAGCGGCTGTAGGCGCTGCTGTTGCGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGT‑‑TCAAGAACAGCATGAAATc                                                                                                                                   >  6:261108/1‑148 (MQ=255)
                  cAGCGCGTCCCACAGGGTGCGAGCCTCCCACCAGCGGCTGTAGGCGCTGCTGTTGCGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGT‑‑TCAAGAACAGCATGAAATCGa                                                                                                                                 >  6:224394/1‑149 (MQ=255)
                   gcttcgcCCCACAGGGTGCGAGCCTCACACCAGCGGCTGTAGGCGCTGCTGTTGCGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGT‑‑TCAAGAACAGCATGAAATCGAt                                                                                                                                <  6:214194/145‑1 (MQ=255)
                   aGCGCGCCCCACAGGGTGCGAGCCTCCCACCAGCGGCTGTAGGCGCTGCTGTTGCGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGT‑‑TCAAGAACAGCATGAAATCGAt                                                                                                                                <  7:75738/149‑1 (MQ=255)
                   aGCGCGCCCCACAGGGTGCGAGCCTCCCACCAGCGGCTGTAGGCGCTGCTGTTGCGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGT‑‑TCAAGAACAGCATGAAATCGAt                                                                                                                                <  4:28709/149‑1 (MQ=255)
                    gcgcgcCCCACAGGGTGCGAGCCTCCCACCAGCGGCTGTAGGCGCTGCTGTTGCGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGT‑‑TCAAGAACAGCATGAAATCGATg                                                                                                                               <  3:214499/149‑1 (MQ=255)
                       cgcCCCACAGGGTGCGAGCCTCCCACCAGCGGCTGTAGGCGCTGCTGTTGCGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGT‑‑TCAAGAACAGCATGAAATCGATGGTg                                                                                                                            <  6:248186/149‑1 (MQ=255)
                        gcCCCACAGGGTGCGAGCCTCCCACCAGCGGCTGTAGGCGCTGCTGTTGCGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGT‑‑TCAAGAACAGCATGAAATCGATGGTGa                                                                                                                           >  1:80496/1‑149 (MQ=255)
                         ccccACAGGGTGCGAGCCTCCCACCAGCGGCTGTAGGCGCTGCTGTTGCGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGT‑‑TCAAGAACAGCATGAAATCGATGGTGAc                                                                                                                          >  5:365129/1‑149 (MQ=255)
                          cccACAGGGTGCGAGCCTCCCACCAGCGGCTGTAGGCGCTGCTGTTGCGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGT‑‑TCAAGAACAGCATGAAATCGATGGTGAcc                                                                                                                         <  6:1214/149‑1 (MQ=255)
                            cacaGGGTGCGAGCCTCCCACCAGCGGCTGTAGGCGCTGCTGTTGCGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGT‑‑TCAAGAACAGCATGAAATCGATGGTGACCAg                                                                                                                       >  2:237809/1‑149 (MQ=255)
                             acaGGGTGCGAGCCTCCCACCAGCGGCTGTAGGCGCTGCTGTTGCGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGT‑‑TCAAGAACAGCATGAAATCGATGGTGACCAGc                                                                                                                      <  6:365129/149‑1 (MQ=255)
                                 ggTGCGAGCCTCCCACCAGCGGCTGTAGGCGCTGCTGTTGCGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTAAACTTGCTGT‑‑TGAAGAACAGCATGAAATCGATGGTGACCAGCGCAt                                                                                                                  >  7:166968/1‑149 (MQ=255)
                                  gTGCGAGCCTCCCACCAGCGGCTGTAGGCGCTGCTGTTGCGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGT‑‑TCAAGAACAGCATGAAATCGATGGTGACCAGCACATc                                                                                                                 <  6:317149/149‑1 (MQ=255)
                                     cGAGCCTCCCACCAGCGGCTGTAGGCGCTGCTGTTGCGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGT‑‑TCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCa                                                                                                              <  5:224394/149‑1 (MQ=255)
                                       aGCCTCCCACCAGCGGCTGTAGGCGCTGCTGTTGCGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGT‑‑TCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGa                                                                                                            <  5:261108/149‑1 (MQ=255)
                                        gCCTCCCCCCAGCGGCTGTAGGCGCTGCTGTTTCGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGT‑‑TCACGAACAGCAGGAAATCGATGGTGACCAGCACATCCCAGag                                                                                                           >  8:332038/1‑148 (MQ=255)
                                              caccaGCGGCTGTAGGCGCTGCTGTTGCGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGT‑‑TCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATcagcag                                                                                                     >  2:6215/1‑149 (MQ=255)
                                               acctGCGGCTGTAGGCGCTGCTGTTGCGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGTGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGT‑‑TCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCagcagc                                                                                                    <  7:142404/149‑1 (MQ=255)
                                               accaGCGGCTGTAGGCGCTGCTGTTGCGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGT‑‑TCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCagcagc                                                                                                    <  1:6215/149‑1 (MQ=255)
                                                     ggCGGTAGGCGCTGCTGTTGCGAAAACTCACCCGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGT‑‑TCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAAc                                                                                              <  8:274108/149‑1 (MQ=255)
                                                     cgCTGTAGGCGCTTGTGTTGCGAAAACTCACCCGCACCACCAGCGCTGAGCCCAGCATGGGCAGCGGGATCAGTGGCAGGGTGAACTTGCTGT‑‑TCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAAc                                                                                              <  8:274104/148‑1 (MQ=255)
                                                      gCTGTAGGCGCTGCTGTTGCGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGT‑‑TCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACa                                                                                             <  2:353330/149‑1 (MQ=255)
                                                          tAGGCGCTGCTGTTGCGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGT‑‑TCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGa                                                                                         <  2:165728/149‑1 (MQ=255)
                                                               gctgctGTTGCGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGT‑‑TCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGc                                                                                    >  4:71364/1‑149 (MQ=255)
                                                                   ctgTTGCGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGT‑‑TCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGa                                                                                         <  3:60967/140‑1 (MQ=255)
                                                                   ctgTTGCGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGT‑‑TCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGa                                                                                         >  4:60967/1‑140 (MQ=255)
                                                                   ctgTTGCGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGT‑‑TCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGa                                                                                 <  3:39148/148‑1 (MQ=255)
                                                                          gAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGT‑‑TCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACa                                                                                             <  3:222031/129‑1 (MQ=255)
                                                                          gAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGT‑‑TCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACa                                                                                             >  4:222031/1‑129 (MQ=255)
                                                                           aaaaCTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGT‑‑TCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGa                                                                        >  8:142814/1‑149 (MQ=255)
                                                                                tCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGT‑‑TCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATca                                                                                                         >  5:13944/1‑111 (MQ=255)
                                                                                tCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGT‑‑TCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATca                                                                                                         <  6:13944/111‑1 (MQ=255)
                                                                                  aCCAGCACCACCAGCGCTGTTCCCAGCAGGGTCAGCGGGATCAGTGTCAGGGTGAACTTGCTGT‑‑TCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGATCAGCGACCAGCCGATATAGCCGATGGTc                                                                   <  3:54815/147‑1 (MQ=255)
                                                                                  aCCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGT‑‑TCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATgg                                                                     <  6:234167/145‑1 (MQ=255)
                                                                                  aCCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGT‑‑TCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATTGtag                                                                  >  5:234167/1‑146 (MQ=255)
                                                                                  aCCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGT‑‑TCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTc                                                                   <  6:332724/147‑1 (MQ=255)
                                                                                  aCCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGT‑‑TCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTc                                                                   <  2:1296/147‑1 (MQ=255)
                                                                                  aCCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGT‑‑TCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCt                                                                  <  1:237809/148‑1 (MQ=255)
                                                                                   ccAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGTAGGCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGCCCAGCCGATATAGCCGATGGtat                                                                  >  1:112772/1‑147 (MQ=255)
                                                                                   ccAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGT‑‑TCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCtt                                                                 >  8:301806/1‑148 (MQ=255)
                                                                                    cAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGT‑‑TCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTca                                                               >  2:8397/1‑149 (MQ=255)
                                                                                    cAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGT‑‑TCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTca                                                               <  1:275375/149‑1 (MQ=255)
                                                                                    cAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGT‑‑TCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTca                                                               >  4:23083/1‑149 (MQ=255)
                                                                                    cAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGT‑‑TCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTca                                                               <  7:142814/149‑1 (MQ=255)
                                                                                     aGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGT‑‑TCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTcac                                                              >  8:109036/1‑149 (MQ=255)
                                                                                       caccacCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGT‑‑TCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATcagca                                                                                                      >  3:324769/1‑107 (MQ=255)
                                                                                       caccacCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGT‑‑TCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATcagca                                                                                                      <  4:324769/107‑1 (MQ=255)
                                                                                                     aGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGT‑‑TCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTCACCACCAGGCGGTACtt                                               <  3:249626/148‑1 (MQ=255)
                                                                                                          agcagGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGT‑‑TCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTCACCACCAGGCAGTACTTGCGGGc                                         >  1:129889/1‑149 (MQ=255)
                                                                                                            cagGGTCAGCGGGATCAGTGGCAGGGTAAACTTGCTGT‑‑TCAAGAACAGCATGAAATCGATGGTGACCAGCGCATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTCACCACCAGGCGGTACTTGCGGGCGa                                       <  8:166968/149‑1 (MQ=255)
                                                                                                             agtgTCAGCGGTATCAGTGGCAGGGTGAACTTGCTGT‑‑TCTATAACAGCATTAAATCGATGGTGACCAGAACATCCCATATCAGCAGCCAGAACAGCGACCATACGATATAGCCGATGGTCTTCACCACCAGTCGGTACTTGCGGGCGAt                                      <  2:294554/146‑1 (MQ=255)
                                                                                                                gTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGT‑‑TCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTCACCACCAGGCGGTACTTGCGGGCGATGAc                                   >  2:312177/1‑149 (MQ=255)
                                                                                                                gTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGT‑‑TCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTCACCACCAGGCGGTACTTGCGGGCGATGAc                                   <  3:71364/149‑1 (MQ=255)
                                                                                                                             gTGGCAGGGTGAACTTGCTGT‑‑TCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTCACCACCAGGCGGTACTTGCGGGCGATGACGTCTTTCAcgc                        >  3:365661/1‑147 (MQ=255)
                                                                                                                             gTGGCAGGGTGAACTTGCTGT‑‑TCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTCACCACCAGGCGGTACTTGCGGGCGATGACGTCTTTCACGCGa                      <  1:8397/149‑1 (MQ=255)
                                                                                                                             gTGGCAGGGTGAACTTGCTGT‑‑TCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTCACCACCAGGCGGTACTTGCGGGCGATGACGTCTTTCACGCGa                      <  3:181377/149‑1 (MQ=255)
                                                                                                                              tGGCAGGGTGAACTTGCTGT‑‑TCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTCACCACCAGGCGGTACTTGCGGGCGATGACGTCTTTCACGCGAg                     <  2:80496/149‑1 (MQ=255)
                                                                                                                                gCAGGGTGAACTTGCTGT‑‑TCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTCACCACCAGGCGGTACTTGCGGGCGATGACGTCTTTCACGCGAggg                   <  2:129889/149‑1 (MQ=255)
                                                                                                                                  aGGGTGAACTTGCTGT‑‑TCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGa                                                                        <  8:203402/94‑1 (MQ=255)
                                                                                                                                     gTGAACTTGCTGT‑‑TCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTCACCACCAGGCGGTACTTGCGGGCGATGACGTCTTTCACGCGAGGGTACCt              <  2:112772/149‑1 (MQ=255)
                                                                                                                                       tttcatGCTGT‑‑TCTTGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTCACCACCAGGCGGTACTTGCGGGCGATGACGTCTTTCACGCGAGGGTACCTc             >  3:6686/6‑148 (MQ=255)
                                                                                                                                        aaCTTGCTGT‑‑TCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTCACCACCAGGCGGTACTTGCGGGCGATGACGTCTTTCACGCGAggg                   >  4:163290/1‑141 (MQ=255)
                                                                                                                                        aaCTTGCTGT‑‑TCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTCACCACCAGGCGGTACTTGCGGGCGATGACGTCTTTCACGCGAggg                   <  3:163290/141‑1 (MQ=255)
                                                                                                                                        aaCTTGCTGT‑‑TCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTCACCACCAGGCGGGACTTGCGGGCGATGACGTCTTTCACGCGAGGGTACCTCg            >  7:348965/1‑148 (MQ=255)
                                                                                                                                          catGCTGT‑‑TCTTCAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTCACCACCAGGCGGTACTTGCGGGCGATGACGTCTTTCACGCGAGGGTACCTCGaaa         >  7:176727/3‑149 (MQ=255)
                                                                                                                                           atGCTGT‑‑TCTTGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTCACCACCAGGCGGTACTTGCGGGCGATGACGTCttt                             >  5:183103/2‑128 (MQ=255)
                                                                                                                                           atGCTGT‑‑TCTTGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCACCCGATATAGCCGATGGTCTTCACCACCAGGCGGTACTTGCGGGCGATGACGTCttt                             <  6:183103/127‑1 (MQ=255)
                                                                                                                                                    tCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGAAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTCACCACCAGGCGGTACTTGCGGGCGATGACGTCTTTCACGCGAGGGTACCTCGAAACGCgggg  >  7:199690/1‑148 (MQ=255)
                                                                                                                                                     cAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTCACCACCAGGCg                                                     <  5:30132/96‑1 (MQ=255)
                                                                                                                                                     cAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTCACCACCAGGCg                                                     >  6:30132/1‑96 (MQ=255)
                                                                                                                                                      aaGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTCACCACCAGGCGGTACTTGCGGGCGATGACGTCTTTCACGCGAgg                    >  4:23/1‑128 (MQ=255)
                                                                                                                                                      aaGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTCACCACCAGGCGGTACTTGCGGGCGATGACGTCTTTCACGCGAgg                    <  3:23/128‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                 
GCTGCGCGAACTGTTCACCAGCGCGCCCCACAGGGTGCGAGCCTCCCACCAGCGGCTGTAGGCGCTGCTGTTGCGAAAACTCACCAGCACCACCAGCGCTGAGCCCAGCAGGGTCAGCGGGATCAGTGGCAGGGTGAACTTGCTGT‑‑TCAAGAACAGCATGAAATCGATGGTGACCAGCACATCCCAGATCAGCAGCCAGAACAGCGACCAGCCGATATAGCCGATGGTCTTCACCACCAGGCGGTACTTGCGGGCGATGACGTCTTTCACGCGAGGGTACCTCGAAACGCGGGG  >  NC_002947/966447‑966740

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: