Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,597,756 T→G 6.2% K685Q (AAG→CAG) ‡ PP_4071 ← membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,597,7560TG 6.2% 200.7 / 4.5 65K685Q (AAG→CAG) ‡PP_4071membrane protein
Reads supporting (aligned to +/- strand):  ref base T (28/33);  new base G (2/2);  total (30/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.83e-01

GCCTTTGGCAGTGAGCATTGGCGCCAGCAAAACACCATTGTCGCCCTGGCTCGA‑GCCGGAAGACAGCTGGTAAACGACATCGAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTCTGCTGCACGTGGCTGACACGCTGGCGATACTGCGGCAGGTCGGTATCCGCGTAGGTCAGCGCATACTTCAGGTGCCGCATGAGGTCTCGTTGCAGCT  >  NC_002947/4597610‑4597899
                                                                                                                                                   |                                                                                                                                               
gCCTTTGGCAGTGAGCATTGGCGCCAGCAAAACACCATTGTCGCCCTGGCTCGA‑GCCGGAAGACAGCTGGTAAACGACATCGAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCtt                                                                                                                                                 >  2:89410/1‑147 (MQ=255)
 ccTTTGGCAGTGAGCATTGGCGCCAGCAAAACACCATTGTCGCCCTGGCTCTA‑GCCGTAAGACATCTGGTAAACGACATCGAACGCTGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAg                                                                                                                                              <  4:21638/149‑1 (MQ=255)
   tttGGCAGTGAGCATTGGCGCCAGCAAAACACCATTGTCGCCCTGGCTCGA‑GCCGGAAGACAGCTGGTAAACGACATCGAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGc                                                                                                                                            >  3:302477/1‑149 (MQ=255)
      ggCAGTGAGCATTGGCGCCAGCAAAACACCATTGTCGCCCTGGCTCGATGCCGGATGACAGCTGGTAAACGACATCGAACGCAGGCCCACCCGGATGACCCAGGTCCAGACCGGAATGCAACTGCTCCTGGGCCTGCTGCTTCAGGccc                                                                                                                                          >  1:274218/1‑149 (MQ=255)
        cAGTGAGCATTGGCGCCAGCAAAACACCATTGTCGCCCTGGCTCGA‑GCCGGAAGACAGCTGGTAAACGACATCGAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCa                                                                                                                                       >  4:209840/1‑149 (MQ=255)
         aGTGAGCATTGGCGCCAGCAAAACACCATTGTCGCCCTGGCTCGA‑GCCGGAAGACAGCTGGTAAACGACATCGAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCAt                                                                                                                                      <  6:35141/149‑1 (MQ=255)
         aGTGAGCATTGGCGCCAGCAAAACACCATTGTCGCCCTGGCTCGA‑GCCGGAAGACAGCTGGTAAACGACATCGAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGAAGCAGCAGGcccagga                                                                                                                                      <  2:279680/149‑5 (MQ=255)
           tGAGCATTGGCGCCAGCAAAACACCATTGTCGCCCTGGCTCGA‑GCCGGAAGACAGCTGGTAAACGACATCGAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATa                                                                                                                                     >  1:368929/1‑148 (MQ=255)
             aGCATTGGCGCCAGCAAAACACCATTGTCGCCCTGGCTCGA‑GCCGGAAGACAGCTGGTAAACGACATCGAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATaca                                                                                                                                   <  2:117999/148‑1 (MQ=255)
                aTTGGCGCCAGCAAAACACCATTGTCGCCCTGGCTCGA‑GCCGGAAGACAGCTGGTAAACGACATCGAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACAcgcg                                                                                                                               >  1:315283/1‑149 (MQ=255)
                 ttGGCGCCAGCAAAACACCATTGTCGCCCTGGCTCGA‑GCCGGAAGACAGCTGGTAAACGACATCGAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGt                                                                                                                              >  3:180996/1‑149 (MQ=255)
                     ccccAGCAAAACACCCTTGGCCCCCTGGCTCGA‑GCCGGAAGACAGCTGGTAAACGACATCGAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTtgtg                                                                                                                          <  6:361247/147‑1 (MQ=255)
                      gcCAGCAAAACACCATTGTCGCCCTGGCTCGA‑GCCGGAAGACAGCTGGTAAACGACATCGAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTgg                                                                                                                         <  1:199410/149‑1 (MQ=255)
                      gcCAGCAAAACACCATTGTCGCCCTGGCTCGA‑GCCGGAAGACAGCTGGTAAACGACATCGAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCCGGGCCTGAAGCAGCAGGcccag                                                                                                                                        >  2:2362/1‑132 (MQ=255)
                      gcCAGCAAAACACCATTGTCGCCCTGGCTCGA‑GCCGGAAGACAGCTGGTAAACGACATCGAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCCGGGCCTGAAGCAGCAGGcccag                                                                                                                                        <  1:2362/134‑3 (MQ=255)
                            aaaaCACCATTGTCGCCCTGGCTCGA‑GCCGGAAGACAGCTGGTAAACGACATCGAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACgg                                                                                                                   <  7:84059/149‑1 (MQ=255)
                             aaaCACCATTGTCGCCCTGGCTCGA‑GCCGGAAGACAGCTGGTAAACGACATCGAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGc                                                                                                                  <  1:89410/149‑1 (MQ=255)
                               agacCATTGTCGCCCTGGCTCGA‑GCCGGAAGACAGCTGGTAAACGACATCGAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCAc                                                                                                                <  2:274218/147‑1 (MQ=255)
                               acacCATTGTCGCCCTGGCTCGA‑GCCGGAAGACAGCTGGTAAACGACATCGAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCAc                                                                                                                <  7:122587/149‑1 (MQ=255)
                                         gccacTGGCTCGA‑GCGGGAAGACAGGAGGTATACGAGATCGAACGCAGGCCGCTCCGTATGACGCAGGTCCAGTCCGGTATGCTCCTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAAt                                                                                                      <  8:342163/145‑1 (MQ=255)
                                           cccTGGCTCGA‑GCCGGAAGACAGCTGGTAAACGACATCGAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGt                                                                                                    <  8:353204/149‑1 (MQ=255)
                                            ccTGGCTCGA‑GCCGGAAGACAGCTGGTAAACGACATCGAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGGTCCTGGGACTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTc                                                                                                   <  5:291866/149‑1 (MQ=255)
                                             cTGGCTCGA‑GCCGGAAGACAGCTGGTAAACGACATCGAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTct                                                                                                  <  7:233248/149‑1 (MQ=255)
                                                    gA‑GCCGGAAGACAGCTGGTAAACGACATCGAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCg                                                                                                                                         >  1:96635/1‑103 (MQ=255)
                                                    gA‑GCCGGAAGACAGCTGGTAAACGACATCGAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCg                                                                                                                                         <  2:96635/103‑1 (MQ=255)
                                                     a‑GCCGGAAGACAGCTGGTAAACGACATCGAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTCTGCTGCACg                                                                                          >  5:154126/1‑149 (MQ=255)
                                                       gCCGGAAGACAGCTGGTAAACGACATCGAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTCTGCTGCACGt                                                                                         <  8:28268/149‑1 (MQ=255)
                                                           gAAGACAGCTGGTAAACGACATCGAACGCAGGCCCCACCTGATGACGCAGTTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTCTGCTGCACGTGGCt                                                                                     <  7:272834/149‑1 (MQ=255)
                                                                   cTGGTAAACGACATCGAACGCAGGCCCCACCTGATGACGCAGGTCCAGCCCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGGATACACGCGTTGTGGCACCGGCACCTTACGCAATGTCTGCTGCACGTGGATGACCCGCt                                                                             >  8:344747/1‑149 (MQ=255)
                                                                         aaCGACATCGAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTCTGCTGCACGTGGCTGACACGCTGGCGa                                                                        >  1:133795/1‑148 (MQ=255)
                                                                           cGACATCGAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTCTGCTGCACGTGGCTGACACGCTGGCGATAc                                                                     >  2:64910/1‑149 (MQ=255)
                                                                              cATCGAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTCTGCTGCACGTGGCTGACACGCTGGCGATACTGc                                                                  <  2:83825/149‑1 (MQ=255)
                                                                                tCGAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTCTGCTGCACGTGGCTGACACGCTGGCGATACTGCgg                                                                >  2:38645/1‑149 (MQ=255)
                                                                                tCGAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTCTGCTGCACGTGGCTGACACGCTGGCGATACTGCgg                                                                <  5:220103/149‑1 (MQ=255)
                                                                                 cGAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTCTGCTGCACGTGGCTGACACGCTGGCGATa                                                                      >  7:300360/1‑142 (MQ=255)
                                                                                 cGAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTCTGCTGCACGTGGCTGACACGCTGGCGATa                                                                      <  8:300360/142‑1 (MQ=255)
                                                                                  gAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTCTGCTGCACGTGGCTGACACGCTGGCGATACTGCGGc                                                               <  4:149064/148‑1 (MQ=255)
                                                                                     cGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGGGTTGTGGCAACGGCACCTTACGCAATGTCTGCTGCACGTGGCTGACACGCTGGCGATACTGCGGCAGGt                                                           <  3:209840/149‑1 (MQ=255)
                                                                                          gCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTCTGCTGCACGTGGCTGACACGCTGGCGATACTGCGGCAGGTCGGTa                                                      <  1:96236/149‑1 (MQ=255)
                                                                                            cccACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTCTGCTGCACGTGGCTGACACGCTGGCGATACTGCGGCAGGTCGGTATc                                                    >  2:344313/1‑149 (MQ=255)
                                                                                                    atgaCGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTCTGCTGCACGTGGCTGACACGCTGGCGATACTGCGGCAGGTCGGTATCCGCGTa                                              <  1:344313/147‑1 (MQ=255)
                                                                                                    atgaCGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTCTGCTGCACGTGGCTGACACGCTGGCGATACTGCGGCAGGTCGGTATCCGCGTAgg                                            >  1:42805/1‑149 (MQ=255)
                                                                                                        cGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTCTGCTGCACGTGGCTGACACGCTGGCGATACTGCGGCAGGTCGGTATCCGCGTAGGTCAg                                        >  6:230347/1‑149 (MQ=255)
                                                                                                           agttCCAGACCGGTATGCAACTGCTCCTGTGCCTGCTGCTTCAGGCCCGCATACAAGCGTTGTGGCAACGGCACCTTACGCAATGTCTGCTGCACGTGGCTGACACGCTGTCTATAATGCGGAAGGTCGTTATCAGCGTAGGTAAgcg                                      >  4:291251/4‑148 (MQ=255)
                                                                                                               ccAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTCTGCTGCACGTGGCTGACACGCTGGCGATACTGCGGCAGGTCGGTATCCGCGTAGGTCAg                                        <  6:205319/142‑1 (MQ=255)
                                                                                                               ccAGACCGGTATGCAACTGCTCATGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTCTGCTGCACGTGGCTGACACGCTGGCGATACTGCGGCAGGTCGGTATCCGCGTAGGTCAg                                        >  5:205319/1‑142 (MQ=255)
                                                                                                                cAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTCTGCTGCACGTGGCTGACACGCTGGCGATACTGCGGCAGGTCGGTATCCGCGTAGGTCAGCGCATAc                                 <  4:180996/148‑1 (MQ=255)
                                                                                                                  gACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCACGCATACACGCGTTGTGGAAACGGCACCTTACGCAATGTCTGCTGCACGTGGCTGACACGCTGTCGATACTGCGGCAGGTCGGTATCCGCGTAGGTCAGCGCATACTTc                              <  3:291251/149‑1 (MQ=255)
                                                                                                                   aCCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTCTGCTGCACGTGGCTGACACGCTGGCGATACTGCGGCAGGTCGGTATCCGCGTAGGTCAGCGCATACTTCa                             <  5:230347/149‑1 (MQ=255)
                                                                                                                     cGGTATGAAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTCTGCTGCACGTGGCTGACACGCTGGCGATACTGCGGCAGGTCGGTATCCGCGTAGGTCAGCGCATACTTCa                             <  7:344747/147‑1 (MQ=255)
                                                                                                                         aTGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTCTGCTGCACGTGGCTGACACGCTGGCGATACTGCGGCAGGTCGGTATCCGCGTAGGTCAGCGCATACTTCAGGTGcc                       >  2:136022/1‑149 (MQ=255)
                                                                                                                         aTGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTCTGCTGCACGTGGCTGACACGCTGGCGATACTGCGGCAGGTCGGTATCCGCGTAGGTCAGCGCATACTTCAGGTGcc                       >  2:141659/1‑149 (MQ=255)
                                                                                                                           gCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTCTGCTGCACGTGGCTGACACGCTGGCGATACTGCGGCAGGTCGGTATCCGCGTAGGTCAGCGCATACTTCAg                            >  6:248131/1‑142 (MQ=255)
                                                                                                                           gCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTCTGCTGCACGTGGCTGACACGCTGGCGATACTGCGGCAGGTCGGTATCCGCGTAGGTCAGCGCATACTTCAg                            <  5:248131/142‑1 (MQ=255)
                                                                                                                              aCTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTCTGCTGCACGTGGCTGACACGCTGGCGATACTGCGGCAGGTCGGTATCCGCGTAGGTCAGCGCATACTTCAGGTGCCGCAt                   >  8:133220/1‑148 (MQ=255)
                                                                                                                                 gtatgcgggCCTGAAGCAGCAGGCCCACATACACGCGTTGTGGCAACGGCACCTTACGCAATGTCTGCTGCACGTGGCTGACACGCTGGCGATACTGCGGCAGGTCGGTATCCGCGTAGGTCAGCGCATACTTCAGGTGCCGCATGAgg               >  5:117674/7‑149 (MQ=255)
                                                                                                                                  cTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTCTGCTGCACGTGGCTGACACGCTGGCGATACTGCGGCAGGTCGGTATCCGCGTAGGTCAGCGCATACTTCAGGTGCCGCATGAg                <  1:64910/147‑1 (MQ=255)
                                                                                                                                  cTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTCTGCTGCACGTGGCTGACACGCTGGCGATACTGCGGCAGGTCGGTATCCGCGTAGGTCAGCGCATACTTCAGGTGCCGCATGAGGt              >  5:77520/1‑149 (MQ=255)
                                                                                                                                   tCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTCTGCTGCACGTGGCTGACACGCTGGCGATACTGCGGCAGGTCGGTATCCGCGTAGGTCAGCGCATACTTCAGGTGCCGCATGAGGtc             >  4:127963/1‑149 (MQ=255)
                                                                                                                                      tGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTCTGCTGCACGTGGCTGACACGCTGGCGATACTGCGGCAGGTCGGTATCCGCGTAGGTCAGCGCATACTTCAGGTGCCGCATGAGGTCTCg          >  5:84735/1‑149 (MQ=255)
                                                                                                                                       gggCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTCTGCTGCACGTGGCTGACACGCTGGCGATACTGCGGCAGGTCGGTATCCGCGTAGGTCAGCGCATACTTCAGGTGCCGCATGAGGTCTCGt         >  5:378012/1‑149 (MQ=255)
                                                                                                                                        ggCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTCTGCTGCACGTGGCTGACACGCTGGCGATACTGCGGCAGGTCGGTATCCGCGTAGGTCAGCGCATACTTCAGGTGCCGCATGAGGTCTCg          <  1:38645/147‑1 (MQ=255)
                                                                                                                                          ccTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGGCTGCTGCACGTGGCTGACACGCTGGCGATACTGCGGCAGGTCGGTATCCGCGTAGGTCAGCGCATACTTCAGGTGCCGCATGAGGTCTCGTTGc      <  3:255670/149‑1 (MQ=255)
                                                                                                                                              ctgctTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTCTGCTGCACGTGGCTGACACGCTGGCGATACTGCGGCAGGTCGGTATCCGCGTAGGTCAGCGCATACTTCAGGTGCCGCATGAGGTCTCGTTGCAGCt  >  8:57782/1‑149 (MQ=255)
                                                                                                                                                gctTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTCTGCTGCACGTGGCTGACACGCTGGCGATACTGCGGCAGGTCGGTATCcgcg                                                <  7:113390/101‑1 (MQ=255)
                                                                                                                                                gctTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTCTGCTGCACGTGGCTGACACGCTGGCGATACTGCGGCAGGTCGGTATCcgcg                                                >  8:113390/1‑101 (MQ=255)
                                                                                                                                                   |                                                                                                                                               
GCCTTTGGCAGTGAGCATTGGCGCCAGCAAAACACCATTGTCGCCCTGGCTCGA‑GCCGGAAGACAGCTGGTAAACGACATCGAACGCAGGCCCCACCTGATGACGCAGGTCCAGACCGGTATGCAACTGCTCCTGGGCCTGCTGCTTCAGGCCCGCATACACGCGTTGTGGCAACGGCACCTTACGCAATGTCTGCTGCACGTGGCTGACACGCTGGCGATACTGCGGCAGGTCGGTATCCGCGTAGGTCAGCGCATACTTCAGGTGCCGCATGAGGTCTCGTTGCAGCT  >  NC_002947/4597610‑4597899

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