Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,380,890 G→T 6.1% A188S (GCC→TCC)  tnpT‑I → cointegrate resolution protein T

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,380,8900GT 6.1% 209.3 / 3.7 66A188S (GCC→TCC) tnpT‑Icointegrate resolution protein T
Reads supporting (aligned to +/- strand):  ref base G (31/31);  new base T (2/2);  total (33/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.95e-01

GCCAACAGACCGAGGTTGAGAATGCCCGCTTGGCGCAGGCTAACGGCGATCAGGAGGT‑GCGCCTGCAGGATCGCGACAGCCAAATTCGCTCGCTGGAGGAAAAGCATCAACATGCCCGTGATGCTTTGGAGCACTATCGTCAGGCCAGCAAGGAGCAGCGCGAGCAGGAACAGCGTCGACACGAGTCGCAGGTGCAACAACTGCAACTGGAACTGCGGCAACTGCAGCAGACCTTGATAATCAAGCAGGACGAACTGACGCAGCTGAACCGCGACAATGCACGCTTGCT  >  NC_002947/3380747‑3381035
                                                                                                                                                |                                                                                                                                                 
gCCAACAGACCGAGGTTGAGAATGCCCGCTTGGCGCAGGCTAACGGCGATCAGGAGGT‑GCGCCTGCAGGATCGCGACAGCCAAATTCGCTCGCTGGAGGAAAAGCATCAACATGCCCGTGATGCTTTGGAGCACTATCGTCAGGCCAg                                                                                                                                               >  1:154552/1‑148 (MQ=255)
  cAACAGACCGAGGTTGAGAATGCCCGCTTGGCGCAGGCTAACGGCGATCAGGAGGTAGCGCCTGCAGGATCGCGACAGCCAAATTCGCTCGCTGGAGGAAAAGCATCAACATGCCCGTGATGCTTTGGAGCACTATCGTCAGGCCAGc                                                                                                                                              >  1:280588/1‑148 (MQ=255)
          cGAGGTTGAGAATGCCCGCTTGGCGCAGGCTAACGGCGATCAGGAGGT‑GCGCCTGCAGGATCGCGACAGCCAAATTCGCTCGCTGGAGGAAAAGCATCAACATGCCCGTGATGCTTTGGAGCACTATCGTCAGGCCAGCAAGGagcagc                                                                                                                                    >  5:103063/1‑149 (MQ=255)
          cGAGGTTGAGAATGCCCGCTTGGCGCAGGCTAACGGCGATCAGGAGGT‑GCGCCTGCAGGATCGCGACAGCCAAATTCGCTCGCTGGAGGAAAAGCATCAACATGCCCGTGATGCTTTGGAGCACTATCGTCAGGCCAGCAAGGagcagc                                                                                                                                    >  5:140630/1‑149 (MQ=255)
            aGGTTGAGAATGCCCGCTTGGCGCAGGCTAACGGCGATCAGGAGGT‑GCGCCTGCAGGATCGCGACAGCCAAATTCGCTCGCTGGAGGAAAAGCATCAACATGCCCGTGATGCTTTGGAGCACTATCGTCAGGCCAGCAAGGAGCAgcgc                                                                                                                                  <  4:52743/149‑1 (MQ=255)
            aGGTTGAGAATGCCCGCTTGGCGCAGGCTAACGGCGATCAGGAGGT‑GCGCCTGCAGGATCGCGACAGCCAAATTCGCTCGCTGGAGGAAAAGCATCAACATGCCCGTGATGCTTTGGAGCACTATCGTCAGGCCAGCAAGGAGCAgcgc                                                                                                                                  <  6:140630/149‑1 (MQ=255)
                        cccGCTTGGCGCAGGCTAACGGCGATCAGGAGGT‑GCGCCTGCAGGATCGCGACAGCCAAATTCGCTCGCTGGAGGAAAAGCATCAACATGCCCGTGATGCTTTGGAGCACTATCGTCAGGCCAGCAAGGAGCAGCGCGAGCAGGAACAg                                                                                                                      >  5:195654/1‑149 (MQ=255)
                         ccGCTTGGCGCAGGCTAACGGCGATCAGGAGGT‑GCGCCTGCAGGATCGCGACAGCCAAATTCGCTCGCTGGAGGAAAAGCATCAACATGCCCGTGATGCTTTGGAGCACTATCGTCAGGCCAGCAAGGAGCAGCGCGAGCAg                                                                                                                            <  1:292143/142‑1 (MQ=255)
                         ccGCTTGGCGCAGGCTAACGGCGATCAGGAGGT‑GCGCCTGCAGGATCGCGACAGCCAAATTCGCTCGCTGGAGGAAAAGCATCAACATGCCCGTGATGCTTTGGAGCACTATCGTCAGGCCAGCAAGGAGCAGCGCGAGCAg                                                                                                                            >  2:292143/1‑142 (MQ=255)
                           gCTTGGCGCAGGCTAACGGCGATCAGGAGGT‑GCGCCTGCAGGATCGCGACAGCCAAATTCGCTCGCTGGAGGAAAAGCATCAACATGCCCGTGATGCTTTGGAGCACTATCGTCAGGCCAGCAAGGAGCAGCGCGAGCAGGAACAGCGt                                                                                                                   <  2:154552/149‑1 (MQ=255)
                            cTTGGCGCAGGCTAACGGCGATCAGGAGGT‑GCGCCTGCAGGATCGCGACAGCCAAATTCGCTCGCTGGAGGAAAAGCATCAACATGCCCGTGATGCTTTGGAGCACTATCGTCAGGCCAGCAAGGAGCAGCGCGAGCAGGAACAGCGTc                                                                                                                  <  2:280588/149‑1 (MQ=255)
                               ggCGCAGGCTAACGGCGATCAGGAGGT‑GCGCCTGCAGGATCGCGACAGCCAAATTCGCTCGCTGGAGGAAAAGCATCAACATGCCCGTGATGCTTTGGAGCACTATCGTCAGGCCAGCAAGGAGCAGCGCGAGCAGGAACAGCGTCGa                                                                                                                <  7:249605/148‑1 (MQ=255)
                                        tAACGGCGATCAGGAGGT‑GCGCCTGCAGGATCGCGACAGCCAAATTCGCTCGCTGGAGGAAAAGCATCAACATGCCCGTGATGCTTTGGAGCACTATCGTCAGGCCAGCAAGGAGCAGCGCGAGCAg                                                                                                                            >  8:207869/1‑127 (MQ=255)
                                        tAACGGCGATCAGGAGGT‑GCGCCTGCAGGATCGCGACAGCCAAATTCGCTCGCTGGAGGAAAAGCATCAACATGCCCGTGATGCTTTGGAGCACTATCGTCAGGCCAGCAAGGAGCAGCGCGAGCAg                                                                                                                            <  7:207869/127‑1 (MQ=255)
                                          aCGGCGATCAGGAGGT‑GCGCCTGCAGGATCGCGACAGCCAAATTCGCTCGCTGGAGGAAAAGCATCAACATGCCCGTGATGCTTTGGAGCACTATCGTCAGGCCAGCAAGGAGCAGCGCGAGCAGGAACAGCGTCGACACGAGTCGCa                                                                                                     <  2:128892/148‑1 (MQ=255)
                                            ggCGATCAGGAGGT‑GCGCCTGCAGGATCGCGACAGCCAAATTCGCTCGCTGGAGGAAAAGCATCAACATGCCCGTGATGCTTTGGAGCACTATCGTCAGGCCAGAAAGGAGCAGCGCGAGCACGACCAGCGTCTCCAAGAGTCGCAgct                                                                                                  >  5:358836/1‑147 (MQ=255)
                                              cGATCAGGAGGT‑GCGCCTGCAGGATCGCGACAGCCAAATTCGCTCGCTGGAGGAAAAGCATCAACATGCCCGTGATGCTTTGGAGCACTATCGTCAGGCCAGCAAGGAGCAGCGCGAGCAGGAACAGCGTCGACACGAGTCGCAGGTg                                                                                                 >  5:169928/1‑148 (MQ=255)
                                              cGATCAGGAGGT‑GCGCCTGCAGGATCGCGACAGCCAAATTCGCTCGCTGGAGGAAAAGCATCAACATGCCCGTGATGCTTTGGAGCACTATCGTCAGGCCAGCAAGGAGCAGCGCGAGCAGGAACAGCGTCGACACGAGTCGCAGGTGc                                                                                                <  5:117728/149‑1 (MQ=255)
                                                 tCAGGAGGT‑GCGCCTGCAGGATCGCGACAGCCAAATTCGCTCGCTGGAGGAAAAGCATCAACATGCCCGTGATGCTTTGGAGCACTATCGTCAGGCCAGCAAGGAGCAGCGCGAGCAGGAACAGCGTCGACACGAGTCGCAGGTGcaa                                                                                              >  5:271155/1‑148 (MQ=255)
                                                   aggaggT‑GCGCCTGCAGGATCGCGACAGCCAAATTCGCTCGCTGGAGGAAAAGCATCAACATGCCCGTGATGCTTTGGAGCACTATCGTCAGGCCAGCAAGGAGCAGCGCGAGCAGGAACAGCGTCGACACGAGTCGCAGGTGcaacaa                                                                                           <  6:195654/149‑1 (MQ=255)
                                                     gaggT‑GCGCCTGCAGGATCGCGACAGCCAAATTCGCTCGCTGGAGGAAAAGCATCAACATGCCCGTGATGCTTTGGAGCACTATCGTCAGGCCAGCAAGGAGCAGCGCGAGCAGGAACAGCGTCGACACGAGTCGCAGGTGCAAcaact                                                                                         <  6:103063/149‑1 (MQ=255)
                                                      aggT‑GCGCCTGCAGGATCGCGACAGCCAAATTCGCTCGCTGGAGGAAAAGCATCAACATGCCCGTGATGCTTTGGAGCACTATCGTCAGGCCAGCAAGGAGCAGCGCGAGCAGGAACAGCGTCGACACGAGTCGCAGGTGCAAcaactg                                                                                        >  5:51140/1‑149 (MQ=255)
                                                               cTGCAGGATCGCGACAGCCAAATTCGCTCGCTGGAGGAAAAGCATCAACATGCCCGTGATGCTTTGGAGCACTATCGTCAGGCCAGCAAGGAGCAGCGCGAGCAGGAACAGCGTCGACACGAGTCGCAGGTGCAACAACTGCAACTgg                                                                                 >  7:329649/1‑148 (MQ=255)
                                                                           gACAGCCAAATTCGCTCGCTGGAGGAAAAGCATCAACATGCCCGTGATGCTTTGGAGCACTATCGTCAGGCCAGCAAGGAGCAGCGCGAGCAGGAACAGCGTCGACACGAGTCGCAGGTGCAACAACTGCAACTGGAACTGCGGCAACt                                                                    <  6:169928/149‑1 (MQ=255)
                                                                              aGCCAAATTCGCTCGCTGGAGGAAAAGCATCAACATGCCCGTGATGCTTTGGAGCACTATCGTCAGGCCAGCAAGGAGCAGCGCGAGCAGGAACAGCGTCGACACGAGTCGCAGGTGCAACAACTGCAACTGGAACTGCGGCAACTg                                                                   >  2:88061/1‑147 (MQ=255)
                                                                                            gctGGAGGAAAAGCATCAACATGCCCGTGATGCTTTGGAGCACTATCGTCAGGCCAGCAAGGAGCAGCGCGAGCAGGAACAGCGTCGACACGAGTCGCAGGTGCAACAACTGCAACTGGAACTGCGGCAACTGCAGCAGACCTTGATaa                                                   >  6:262591/1‑149 (MQ=255)
                                                                                            gctGGAGGAAAAGCATCAACATGCCCGTGATGCTTTGGAGCACTATCGTCAGGCCAGCAAGGAGCAGCGCGAGCAGGAACAGCGTCGACACGAGTCGCAGGTGCAACAACTGCAACTGGAACTGCGGCAACTGCAGCAGACCTTGATa                                                    <  8:329649/148‑1 (MQ=255)
                                                                                                gaggaAAAGCATCAACATGCCCGTGATGCTTTGGAGCACTATCGTCAGGCCAGCAAGGAGCAGCGCGAGCAGGAACAGCGTCGACACGAGTCGCAGGTGCAACAACTGCAACTGGAACTGCTGCAGTTGccgcagttccag                                                       <  1:174801/141‑7 (MQ=255)
                                                                                                gaggaAAAGCATCAACATGCCCGTGATGCTTTGGAGCACTATCGTCAGGCCAGCAAGGAGCAGCGCGAGCAGGAACAGCGTCGACACGAGTCGCAGGTGCAACAACTGCAACTGGAACTGCTGCAGTTGccgcagttccag                                                       >  2:174801/1‑135 (MQ=38)
                                                                                                  ggaAAAGCATCAACATGCCCGTGATGCTTTGGAGCCCTATCGTCAGGCCAGCAAGGAGCAGCGCGAGCAGGAACAGCGACGACAAGACTCGCAGGTGCAACAACTGCAACTGGAACTGCGGCAACAGCAGCAGACCTTGATAATCAAg                                              >  5:135217/1‑148 (MQ=255)
                                                                                                  ggaAAAGCATCAACATGCCCGTGATGCTTTGGAGCACTATCGTCAGGCCAGCAAGGAGCAGCGCGAGCAGGAACAGCGTCGACACGAGTCGCAGGTGCAACAACTGCAACTGGAACTGCGGCAACTGCAGCAGACCTTGATAATCAAg                                              >  6:180326/1‑148 (MQ=255)
                                                                                                  ggaAAAGCATCAACATGCCCGTGATGCTTTGGAGCACTATCGTCAGGCCAGCAAGCAGCAGCGCGAGCAGGAACAGCGTCGACACGAGTCGCAGGTGCAACAACTGCAACTGGAACTGCGGCAACTGCAGCAGACCTTGATAATCAAg                                              >  1:155998/1‑148 (MQ=255)
                                                                                                    aaaaGCATCAAGATCCCCGTGATGCTTTGGAGCACTATCGTCAGGCCAGCAATGAGCAGCGCGAGCAGGTACAGCGTCGACACGAGTCGCAGGTGCAACATCTGCAACTTGAACTGCGGCAACTGCAGGAGACCTTGATATTCAAGCAg                                           <  6:135217/149‑1 (MQ=255)
                                                                                                    aaaaGCATCAACATGCCCGTGATGCTTTGGAGCACTATCGTCAGGCCAGCAAGGAGCAGCGCGAGCAGGAACAGCGTCGACACGAGTCGCAGGTGCAACAACTGCAACTGGAACTGCGGCAACTGCAGCAGACCTTGATAATCAAGCa                                            <  5:262591/148‑1 (MQ=255)
                                                                                                     aaaGCATCAACATGCCCGTGATGCTTTGGAGCACTATCGTCAGTCCAGCAAGGAGCAGCGCGAGCAGGAACAGCGTCGACACGAGTCGCAGGTGCAACAACTGCAACTGGAACTGCGGCAACTGCAGCAGACCTTGATAATCAAGCAg                                           >  1:213403/1‑148 (MQ=255)
                                                                                                     aaaGCATCAACATGCCCGTGATGCTTTGGAGCACTATCGTCAGTCCAGCAAGGAGCAGCGCGAGCAGGAACAGCGTCGACACGAGTCGCAGGTGCAACAACTGCAACTGGAACTGCGGCAACTGCAGCAGACCTTGATAATCAAGCAg                                           >  1:161622/1‑148 (MQ=255)
                                                                                                     aaaGCATCAACATGCCCGTGATGCTTTGGAGCACTATCGTCAGGCCAGCAAGGAGCAGCGCGAGCAGGAACAGCGTCGACACGAGTCGCAGGTGCAACAACTGCAACTGGAACTGCGGCAACTgcag                                                                >  1:338513/1‑127 (MQ=255)
                                                                                                     aaaGCATCAACATGCCCGTGATGCTTTGGAGCACTATCGTCAGGCCAGCAAGGAGCAGCGCGAGCAGGAACAGCGTCGACACGAGTCGCAGGTGCAACAACTGCAACTGGAACTGCGGCAACTgcag                                                                <  2:338513/127‑1 (MQ=255)
                                                                                                     aaaGCATCAACATGCCCGTGATGCTTTGGAGCACTATCGTCAGGCCAGCAAGGAGCAGCGCGAGCAGGAACAGCGTCGACACGAGTCGCAGGTGCAACAACTGCAACTGGAACTGCGGCAACTGCAGCAGACCTTGATAATCAAGCAgg                                          <  6:98136/149‑1 (MQ=255)
                                                                                                        gCATCAATATGCCCGTGATGCTTTGGAGCACTATCGTCAGGCCAGCAAGGAGCAGCGCGAGCAGGAACAGCGTCGACACGAGTCGCAGGTGCAACAACTGCAACTGGAACTGCGTCAACTGCAGCAGACCTTGATAATCAAGCa                                            <  7:263109/144‑1 (MQ=255)
                                                                                                        gCATCAACATGCCCGTGATGCTTTGGAGCACTATCCTCAGGCCAGCAAGGAGCAGCGCGAGCAGGAACAGCGTCGACACGAGTCGCAGGTGCAACAACAGCAACTGGAACTGCGGCAACTGCAGCAGACGATGATAATCAAGCAgaga                                        >  8:263109/1‑145 (MQ=39)
                                                                                                           tCAACATGCCCGTGATGCTTTGGAGCACTATCGTCAGGCCAGCAAGGAGCAGCGCGAGCAGGAACAGCGTCGACACGAGTCGCAGGTGCAACAACTGCAACTGGAACTGCGGCAACTGCAGCAGACCTTGATAATCAAGCAGGACGAAc                                    >  2:228988/1‑149 (MQ=255)
                                                                                                           tCAACATGCCCGTGATGCTTTGGAGCACTATCGTCAGGCCAGCAAGGAGCAGCGCGAGCAGGAACAGCGTCGACACGAGTCGCAGGTGCAACAACTGCAACTGGAACTGCGGCAACTGCAGCAGACCTTGATAATCAAGCAGGACGAAc                                    <  6:51140/149‑1 (MQ=255)
                                                                                                               cATGCCCGTGATGCTTTGGAGCACTATCGTCAGGCAAGAAAGGAGCAGCGAGAGCAGGAACAGCGTCGACACGAGTCGCAGGTGCAACAACTGCAACTGGAACTGCGGCAACTGCATCATACCTTGATAATCAAGAAGGACGAACTGAc                                >  3:125812/1‑149 (MQ=255)
                                                                                                                aTGCCCGTGATGCTTTGGAGCACTATCGTCAGTCCAGCAAGGAGCAGCGCGAGCAGGAACAGCGTCGACACGAGTCGCAGGTGCAACAACTGCAACTGGAACTGCGGCAACTGCAGCAGACCTTGATAATCAAGCAGGACGAACTGACg                               <  2:213403/149‑1 (MQ=255)
                                                                                                                aTGCCCGTGATGCTTTGGAGCACTATCGTCAGTCCAGCAAGGAGCAGCGCGAGCAGGAACAGCGTCGACACGAGTCGCAGGTGCAACAACTGCAACTGGAACTGCGGCAACTGCAGCAGACCTTGATAATCAAGCAGGACGAACTGACg                               <  2:161622/149‑1 (MQ=255)
                                                                                                                      gTGATGCTTTGGAGCACTATCGTCAGGCCAGCAAGGAGCAGCGCGAGCAGGAACAGCGTCGACACGAGTCGCAGGTGCAACAACTGCAACTGGAACTGCGGCAACTGCAGCAGACCTTGATAATCAAGCAGGACGAACTGACGCAGCt                          >  8:208249/1‑148 (MQ=255)
                                                                                                                         aTGCTTTGGAGCACTATCGCCAGGCCAGCAAGGAGCAGCGCGAGCAGGAACAGCGTCGACACGAGTCGCAGGTGCAACAACTGCAACTGGAACTGCGGCAACTGCAGCAGACCTTGa                                                      <  7:317106/117‑1 (MQ=255)
                                                                                                                         aTGCTTTGGAGCACTATCGCCAGGCCAGCAAGGAGCAGCGCGAGCAGGAACAGCGTCGACACGAGTCGCAGGTGCAACAACTGCAACTGGAACTGCGGCAACTGCAGCAGACCTTGa                                                      >  8:317106/1‑117 (MQ=255)
                                                                                                                           gCTTTGGAGCACTATCGTCAGGCCAGCAAGGAGCAGCGCGAGCAGGAACAGCGTCGACACGAGTCGCAGGTGCAACAACTGCAACTGGAACTGCGGCAACTGCAGCAGACCTTGATAATCAAGCAGGACGAACTGACGCAGCTGAACcg                    >  5:313724/1‑149 (MQ=255)
                                                                                                                            cTTTGGAGCACTATCGTCAGGCCAGCAAGGAGCAGCGCGAGCAGGAACAGCGTCGACACGAGTCGCAGGTGCAACAACTGCAACTGGAACTGCGGCAACTGCAGCAGACCTTGATAATCAAGCAGGACGAACTGACGCAGCTGAACcgc                   <  6:313724/149‑1 (MQ=255)
                                                                                                                               tGGAGGACTATCGTCAGGCAAGCAATGAGCAGCGAGTGCATGAACAGCGTCGACACGAGTCGCATGTGCAACAACTGCAACTTGAACTGCGGCAACTGCAGCAGAGCTTGATAATCAAGCAGGACGATCTGACGCAGCTGAACcgcgt                 <  6:358836/148‑2 (MQ=255)
                                                                                                                               tGGAGCACTATCGTCAGGCCAGCAAGGAGCAGCGCGATCAGGAACAGCGTCGACACGAGTCGCAGGTGCAACAACTTCAACTGGAACTGCGGCAACTGCAGCAGACCTTGATAATCAAGCATGACGAACTGACGCAGCTGAACcgc                   <  4:125812/146‑1 (MQ=255)
                                                                                                                               tGGAGCACTATCGTCAGGCCAGCAAGGAGCAGCGCGAGCAGGAACAGCGTCGACACGAGTCGCAGGTGCAACAACTGCAACTGGAACTGCGGCAACTGCAGCAGACCTTGATAATCAAGCAGGACGAACTGACGCAGCTGAACCGCGa                 <  1:228988/148‑1 (MQ=255)
                                                                                                                               tGGAGCACTATCGTCAGGCCAGCAAGGAGCAGCGCGAGCAGGAACAGCGTCGACACGAGTCGCAGGTGCAACAACTGCAACTGGAACTGCGGCAACTGCAGCAGACCTTGATAATCAAGCAGGACGAACTGACGCAGCTGAACCGCGa                 <  7:208249/148‑1 (MQ=255)
                                                                                                                                 gAGCACTATCGTCAGGCCAGCAAGGAGCAGCGCGAGCAGGAACAGCGTCGACACGAGTCGCAGGTGCAACAATTGCAACTGGAACTGCGGCAACTGCAGCAGACCTTGATAATCAAGCAGGACGAACTGACGCAGCTGAACCGCGACaa              >  4:31385/1‑149 (MQ=255)
                                                                                                                                      cTATCGTCAGGCCAGCAAGGAGCAGCGCGAGCAGGAACAGCGTCGACACGAGTCGCAGGTGCAACAACTGCAACTGGAACTGCGGCAACTGCAGCAGACCTTGATAATCAAGCAGGACGAACTGACGCAGCTGAACCGCGACAATGCAc         <  2:155998/149‑1 (MQ=255)
                                                                                                                                      cTATCGTCAGGACAGCAAGGAGCAGCTCGAGCAGGAAGAGCGTCGACACGAGTCGCAGGTGCAACAACTGCAACTGGAACTGCGGCAACTGCAGCAGAACTTGATaa                                                   <  2:112096/107‑1 (MQ=255)
                                                                                                                                       tATCGTCAGGCCAGCAAGGAGCAGCGCGAGCAGGAACAGCGTCGACACGAGTCGCAGGTGCAACAACTGCAACTGGAACTGCGGCAACTGCAGCAGACCTTGATAATCAAGCAGGACGAACTGACGCAGCTGAACCGCGACAATGCACg        >  3:167325/1‑149 (MQ=255)
                                                                                                                                         tCGTCAGGCCAGCAAGGAGCAGCGCGAGCAGGAACAGCGTCGACACGAGTCGCAGGTGCAACAATTGCAACTGGAACTGCGGCAACTGCAGCAGACCTTGATAATCAAGCAGGACGAACTGACGCAGCTGAACCGCGACAATGCACGc       <  3:31385/148‑1 (MQ=255)
                                                                                                                                         tCGTCAGGCCAGCAAGGAGCAGCGCGAGCAGGAACAGCGTCGACACGAGTCGCAGGTGCAACAACTGCAACTGGAACTGCGGCAACTGCAGCAGACCTTGATAATCAAGCAGGACGAACTGACGCAGCTGAACCGCGACAATGCACGc       <  4:167325/148‑1 (MQ=255)
                                                                                                                                         tCGTCAGGCCAGCAAGGAGCAGCGCGAGCAGGAACAGCGTCGACACGAGTCGCAGGTGCAACAACTGCAACTGGAACTGCGGCAACTGCAGCAGACCTTGATAATCAAGCAGGACGAACTGACGCAGCTGAACCGCGACAATGCACGCt      >  8:111864/1‑149 (MQ=255)
                                                                                                                                          cGTCAGGCCAGCAAGGAGCAGCGCGAGCAGGAACAGCGTCGACACGAGTCGCAGGTGCAACAACTGCAACTGGAACTGCGGCAACTGCAGCAGACCTTGATAATCAAGCAGGACGAACTGACGCAGCTGAACCGCGACAATGCACGCtt     >  4:156358/1‑149 (MQ=255)
                                                                                                                                           gTCAGGCCAGCAAGGAGCAGCGCGAGCAGGAACAGCGTCGACACGAGTCGCAGGTGCAACAACTGCAACTGGAACTGCGGCAACTGCAGCAGACCTTGATAATCAAGCAGGACGAACTGACGCAGCTGAACCGCGACAATGCACGCTTg    <  2:74414/149‑1 (MQ=255)
                                                                                                                                            tCAGGCCAGCAAGGAGCAGCGCGAGCAGGAACAGCGTCGACACGAGTCGCAGGTGCAACAACTGCAACTGGAACTGCGGCAACTGCAGCAGACCTTGATAATCAAGCAGGACGAACTGACGCAGCTGAACCGCGACAATGCACGCTTGc   >  5:289126/1‑149 (MQ=255)
                                                                                                                                              aGGCCAGCAAGGAGCAGCGCGAGCAGGAACAGCGTCGACACGAGTCGCAGGTGCAACAACTGCAACTGGAACTGCGGCAACTGCAGCAGACCTTGATAATCAAGCAGGACGAACTGACGCAGCTGAACCGCGACAATGCACGCTTGCt  >  5:321171/1‑148 (MQ=255)
                                                                                                                                                |                                                                                                                                                 
GCCAACAGACCGAGGTTGAGAATGCCCGCTTGGCGCAGGCTAACGGCGATCAGGAGGT‑GCGCCTGCAGGATCGCGACAGCCAAATTCGCTCGCTGGAGGAAAAGCATCAACATGCCCGTGATGCTTTGGAGCACTATCGTCAGGCCAGCAAGGAGCAGCGCGAGCAGGAACAGCGTCGACACGAGTCGCAGGTGCAACAACTGCAACTGGAACTGCGGCAACTGCAGCAGACCTTGATAATCAAGCAGGACGAACTGACGCAGCTGAACCGCGACAATGCACGCTTGCT  >  NC_002947/3380747‑3381035

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: ‑