Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,941,514 T→A 5.4% E15V (GAA→GTA)  PP_3476 ← type II secretion system protein G

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,941,5140TA 5.4% 229.4 / 2.4 74E15V (GAA→GTA) PP_3476type II secretion system protein G
Reads supporting (aligned to +/- strand):  ref base T (36/34);  new base A (2/2);  total (38/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.76e-01

GGTATCACCGTAGTAGTGGTCCAGCGCCTCGCGCAGCACTGCCAGGCTCTGACGCAGGGTGGCCTCGCGCGAGCTTTCCAGGCTGTTGAAGTAGCGCGGCATGGCGATGGTCATCAGCGTGGCGATGATCGCCATCACCACCAGCAGTTCGATCAGGGTGAAGCCTTTGCTGCGTTTCATGGCTGTCACCATTGCCCGTAAGGGATATTGTTCAGGCCTTTGCCGCGGGCCCTGGAATACACGTCGAACACGTCTTCGCCCTCGCGCGGGTTGTCGGGGCTGCTGTCGTAGG  >  NC_002947/3941366‑3941657
                                                                                                                                                    |                                                                                                                                               
ggTATCACCGTAGTAGTGGTCCAGCGCCTCGCGCAGCACTGCCAGGCTCTGACGCAGGGTGGCCTCGCGCGAGCTTTCCAGGCTGTTGAAGTAGCGCGGCATGGCGATGGTCATCAGCGTGGCGATGATCGCCATCACCACCAGCAGtt                                                                                                                                                 <  5:174196/149‑1 (MQ=255)
      aCCGTAGTAGTGGTCCAGCGCCTCGCGCAGCACTGCCAGGCTCTGACGCAGGGTGGCCTCGCGCGAGCTTTCCAGGCTGTTGAAGTAGCGCGGCATGGCGATGGTCATCAGCGTGGCGATGATCGCCATCACCACCAGCAGTTCGAt                                                                                                                                             >  7:34399/1‑147 (MQ=255)
      aCCGTAGTAGTGGTCCAGCGCCTCGCGCAGCACTGCCAGGCTCTGACGCAGGGTGGCCTCGCGCGAGCTTTCCAGGCTGTTGAAGTAGCGCGGCATGGCGATGGTCATCAGCGTGGCGATGATCGCCATCACCACCAGCAGTTCGAt                                                                                                                                             >  4:141020/1‑147 (MQ=255)
       ccGTAGTAGTGGTCCAGCGCCTCGCGCAGCACTGCCAGGCTCTGACGCAGGGTGGCCTCGCGCGAGCTTTCCAGGCTGTTGAAGTAGCGCGGCATGGCGATGGTCATCAGCGTGGCGATGATCGCCATCACCACCAGCAGTTCGATCAg                                                                                                                                          <  3:204727/149‑1 (MQ=255)
        cGTAGTAGTGGTCCAGCGCCTCGCGCAGCACTGCCAGGCTCTGACGCAGGGTGGCCTCGCGCGAGCTTTCCAGGCTGTTGAAGTAGCGCGGCATGGCGATGGTCATCAGCGTGGCGATGATCGCCATCACCACCAGCAGTTCGATCAgg                                                                                                                                         <  6:1984/149‑1 (MQ=255)
            gtagtGGTCCAGCGCCTCGCGCAGCACTGCCAGGCTCTGACGCAGGGTGGCCTCGCGCGAGCTTTCCAGGCTGTTGAAGTAGCGCGGCATGGCGATGGTCATCAGCGTGGCGATGATCGCCATCACCACCAGCAGTTCGATCAg                                                                                                                                          <  6:79563/144‑1 (MQ=255)
            gtagtGGTCCAGCGCCTCGCGCAGCACTGCCAGGCTCTGACGCAGGGTGGCCTCGCGCGAGCTTTCCAGGCTGTTGAAGTAGCGCGGCATGGCGATGGTCATCAGCGTGGCGATGATCGCCATCACCACCAGCAGTTCGATCAg                                                                                                                                          >  5:79563/1‑144 (MQ=255)
             tagtGGTCCAGCGCCTCGCGCTGCACTGCCAGGCTCTGACGCAGGGTGGCCTCGCGCGAGCTTTCCAGGCTGTTGAAGTAGCGCGGCATGGCGATGGTCATCAGCGTGGCGATGATCGCCATCACCACCAGCAGTTCGATCAGTGTGaa                                                                                                                                    <  6:185867/149‑1 (MQ=255)
             tagtGGTCCAGCGCCTCGCGCAGCACTGCCAGGCTCTTACGCAGGGTGGCCTCGCGCGAGCTTTCCAGGCTGTTGAAGTAGCGCGGCATGGCGATGGTCATCAGCGTGGCGATGATCGCCATCACCACCAGCAGTTCGATCAGGGTGaa                                                                                                                                    <  1:62504/149‑1 (MQ=255)
                tGGTCCAGCGCCTCGCGCAGCACTGCCAGGCTCTGACGCAGGGTGGCCTCGCGCGAGCTTTCCAGGCTGTTGAAGTAGCGCGGCATGGCGATGGTCATCAGCGTGGCGATGATCGCCATCACCACCAGCAGTTCGATCAGGGTGAAGcc                                                                                                                                 >  8:103351/1‑149 (MQ=255)
                tGGTCCAGCGCCTCGCGCAGCACTGCCAGGCTCTGACGCAGGGTGGCCTCGCGCGAGCTTTCCAGGCTGTTGAAGTAGCGCGGCATGGCGATGGTCATCAGCGTGGCGATGATCGCCATCACCACCAGCAGTTCGATCAGGGTGAAGcc                                                                                                                                 >  5:4192/1‑149 (MQ=255)
                  gTCCAGCGCCTCGCGCAGCACTGCCAGGCTCTGACGCAGGGTGGCCTCGCGCGAGCTTTCCAGGCTGTTGAAGTAGCGCGGCATGGCGATGGTCATCAGCGTGGCGATGATCGCCATCACCACCAGCAGTTCGATCAGGGTGAAGCCtt                                                                                                                               >  7:286357/1‑149 (MQ=255)
                  gTCCAGCGCCTCGCGCAGCACTGCCAGGCTCTGACGCAGGGTGGCCTCGCGCGAGCTTTCCAGGCTGTTGAAGTAGCGCGGCATGGCGATGGTCATCAGCGTGGCGATGATCGCCATCACCACCAGCAGTTCGATCAGGGTGAAGCCtt                                                                                                                               <  8:286357/149‑1 (MQ=255)
                    ccAGCGCCTCGCGCAGCACTGCCAGGCTCTGACGCAGGGTGGCCTCGCGCGAGCTTTCCAGGCTGTTGAAGTAGCGCGGCATGGCGATGGTCATCAGCGTGGCGATGATCGCCATCACCACCAGCAGTTCGATCAGGGTGAAGCCTTtg                                                                                                                             >  5:287323/1‑149 (MQ=255)
                           cTCGCGCAGCACTGCCAGGCTCTGACGCAGGGTGGCCTCGCGCGAGCTTTCCAGGCTGTTGAAGTAGCGCGGCATGGCGATGGTCATCAGCGTGGCGATGATCGCCATCACCACCAGCAGTTCGATCAGGGTGAAGCCTTTGCTGCGtt                                                                                                                      <  6:287323/149‑1 (MQ=255)
                                         ccAGGCTCTGACGCAGGGTGGCCTCGCGCGAGCTTTCCAGGCTGTTGAAGTAGCGCGGCATGGCGATGGTCATCAGCGTGGCGATGATCGCCATCACCACCAGCAGTTCGATCAGGGTGAAGCCTTTGCTGCGTTTCATGGCTGTCAc                                                                                                         <  7:103351/148‑1 (MQ=255)
                                         ccAGGCTCTGACGCAGGGTGGCCTCGCGCGAGCTTTCCAGGCTGTTGAAGTAGCGCGGCATGGCGATGGTCATCAGCGTGGCGATGATCGCCATCACCACCAGCAGTTCGATCAGGGTGAAGCCTTTGCTGCGTTTCATGGCTGTCAc                                                                                                         >  5:255168/1‑148 (MQ=255)
                                         ccAGGCTCTGACGCAGGGTGGCCTCGCGCGAGCTTTCCAGGCTGTTGAAGTAGCGCGGCATGGCGATGGTCATCAGCGTGGCGATGATCGCCATCACCACCAGCAGTTCGATCAGGGTGAAGCCTTTGCTGCGTTTCATGGCTGTCAc                                                                                                         <  6:84318/148‑1 (MQ=255)
                                          cAGGCTCTGACGCAGGGTGGCCTCGCGCGAGCTTTCCAGGCTGTTGAAGTAGCGCGGCATGGCGATGGTCATCAGCGTGGCGATGATCGCCATCACCACCAGCAGTTCGATCAGGGTGAAGCCTTTGCTGCGTTTCATGGCTGTCACCa                                                                                                       >  1:57259/1‑149 (MQ=255)
                                            ggCTCTGACGCAGGGTGGCCTCGCGCGAGCTTTCCAGGCTGTTGAAGTAGCGCGGCATGGCGATGGTCATCAGCGTGGCGATGATCGCCATCACCACAAGCAGTTCGATCAGGGTGAAGCAGTTGATGCGTTTCATGGCTGTCACCAtt                                                                                                     >  5:60000/1‑149 (MQ=255)
                                              ctctGACGCAGGGTGGCCTCGCGCGAGCTTTCCAGGCTGTTGAAGTAGCGCGGCATGGCGATGGTCATCAGCGTGGCGATGATCGCCATCACCACCAGCAGTTCGATCAGGGTGAAGCCTTTGCTGCGTTTCATGGCTGTCACCATTGc                                                                                                   <  1:81406/149‑1 (MQ=255)
                                               tctGACGCAGGGTGGCCTCGCGCGAGCTTTCCAGGCTGTTGAAGTAGCGCGGCATGGCGATGGTCATCAGCGTGGCGATGATCGCAATCACCACCAGCAGTTCGATCAGGGTGAAGCCTTTGCTGCGTTTCATGGCTGTCACCATTGcc                                                                                                  >  1:340273/1‑149 (MQ=255)
                                                   aCGCAGGGTGGCCTCGCGCGAGCTTTCCAGGCTGTTGAAGTAGCGCGGCATGGCGATGGTCATCAGCGTGGCGATGATCGCCATCACCACCAGCAGTTCGATCAGGGTGAAGCCTTTGCTGCGTTTCATGGCTGTCACCATTGccc                                                                                                 <  3:141020/146‑1 (MQ=255)
                                                     gCAGGGTGGCCTCGCGCGAGCTTTCCAGGCTGTTGAAGTAGCGCGGCATGGCGATGGTCATCAGCGTGGCGATGATCGCCATCACCACCAGCAGTTCGATCAGGGTGAAGCCTTTGCTGCGTTTCATGGCTGTCACCATTGCCCGTAAg                                                                                            >  8:267779/1‑149 (MQ=255)
                                                       aGGGTGGCCTCGCGCGAGCTTTCCAGGCTGTTGAAGTAGCGCGGCATGGCGATGGTCATCAGCGTGGCGATGATCGCCATCACCACCAGCAGTTCGATCAGGGTGAAGCCTTTGCTGCGTTTCATGGCTGTCACCATTGCCCGTAAggg                                                                                          <  6:255168/149‑1 (MQ=255)
                                                       aGGGTGGCCTCGCGCGAGCTTTCCAGGCTGTTGAAGTAGCGCGGCATGGCGATGGTCATCAGCGTGGCGATGATCGCCATCACCACCAGCAGTTCGATCAGGGTGAAGCCTTTGCTGCGTTTCATGGCTGTCACCATTGCCCGTAAggg                                                                                          <  7:267779/149‑1 (MQ=255)
                                                        gggTGGCCTCGCGCGAGCTTTCCAGGCTGTTGAAGTAGCGCGGCATGGCGATGGTCATCAGCGTGGCGATGATCGCCATCACCACCAGCAGTTCGATCAGGGTGAAGCCTTTGCTGCGTTTCATGGCTGTCACCATTGCCCGTAAggg                                                                                          >  3:35158/1‑148 (MQ=255)
                                                        gggTGGCCTCGCGCGAGCTTTCCAGGCTGTTGAAGTAGCGCGGCATGGCGATGGTCATCAGCGTGGCGATGATCGCCATCACCACCAGCAGTTCGATCAGGGTGAAGCCTTTGCTGCGTTTCATGGCTGTCACCATTGCCCGTAAggg                                                                                          <  4:35158/148‑1 (MQ=255)
                                                         ggTGGCCTCGCGCGAGCTTTCCAGGCTGTTGAAGTAGCGCGGCATGGCGATGGTCATCAGCGTGGCGATGATCGCCATCACCACCAGCAGTTCGATCAGGGTGAAGCCTTTGCTGCGTTTCATGGCTGTCACCATTGCCCGTAAGGGat                                                                                        <  4:71869/149‑1 (MQ=255)
                                                         ggTGGCCTCGCGCGAGCTTTCCAGGCTGTTGAAGTAGCGCGGCATGGCGATGGTCATCAGCGTGGCGATGATCGCCATCACCACCAGCAGTTCGATCAGGGTGAAGCCTTTGCTGCGTTTCATGGCTGTCACCATTGCCCGTAAGGGa                                                                                         >  2:323658/1‑148 (MQ=255)
                                                             gCCTCGCGCGAGCTTTCCAGGCTGTTGAAGTAGCGCGGCATGGCGATGGTCATCAGCGTGGCGATGATCGCCATCACCACCAGCAGTTCGATCAGGGTGAAGCCTTTGCTGCGTTTCATGGCTGTCACCATTGCCCGTAAGGGATATTg                                                                                    >  7:160865/1‑149 (MQ=255)
                                                                  gcgcgAGCTTTCCAGGCTGTTGAAGTAGCGCGGCATGGCGATGGTCATCAGCGTGGCGATGATCGCCATCACCACCAGCAGTTCGATCAGGGTGAAGCCTTTGCTGCGTTTCATGGCTGTCACCATTGCCCGTAAGGGATATTGTTCa                                                                                <  8:34399/148‑1 (MQ=255)
                                                                  gcgcgAGCTTTCCAGGCTGTTGAAGTAGCGCGGCATGGCGATGGTCATCAGCGTGGCGATGATCGCCATCACCACCAGCAGTTCGATCAGGGGGAAGCCTTTGCTGCGTATCATGGCTGTCACCATTGCCCGTAAGGGATATTGTTCAg                                                                               <  4:120250/149‑1 (MQ=255)
                                                                  gcgcgAGCTTTCCAGGCTGTTGAAGTAGCGCGGCATGGCGATGGTCATCAGCGTGGCGATGATCGCCATCACCACCAGCAGTACGATCAGGGTGAAGCCTTTGCTGCGTTTCATGGCTGTCACCATTGCCCGTAAGGGATATTGTTCAg                                                                               >  7:85825/1‑149 (MQ=255)
                                                                  gcgcgAGCTTTCCAGGCTGTTGAAGTAGCGCGGCATGGCGATGGTCATCAGCGTGGCGATGATCGCCATCACCACCAGCAGTACGATCAGGGTGAAGCCTTTGCTGCGTTTCATGGCTGTCACCATTGCCCGTAAGGGATATTGTTCAg                                                                               >  7:81255/1‑149 (MQ=255)
                                                                  gcgcgAGCTTTCAAGGCGGGTGAAGTAGCGCGGCATGGCGCTGGTCATCAGCGTGGCGATGATCGCCATCACCACCAGCAGTTCGATCAGGGTGAAGCCTTTGCTGCGTTTCATGGCTGTCACCATTGCCCGTAAGGGATATTGTTCAg                                                                               <  8:160865/149‑1 (MQ=255)
                                                                      gAGCTTTCCAGGCTGTTGAAGTAGCGCGGCATGGCGATTGTCATCAGCGTGGCGATGATCGCCATCACCACCAGCAGTTCGATCAGGGTGAAGCCTTTGCTGCGTTTCATGGCTGTCACCATTGCCCGTAAGGGATATTGTTCAGGCCt                                                                           <  5:234208/149‑1 (MQ=255)
                                                                           ttCCAGGCTGTTGAAGTAGCGCGGCATGGCGATGGTCATCAGCGTGGCGATGATCGCCATCACCACCAGCAGTTCGATCAGGGTGAAGCCTTTGCTGCGTTTCATGGCTGTCACCATTGCCCGTAAGGGATATTGTTCAGGCCttt                                                                         <  5:245085/146‑1 (MQ=255)
                                                                           ttCCAGGCTGTTGAAGTAGCGCGGCATGGCGATGGTCATCAGCGTGGCGATGATCGCCATCACCACCAGCAGTTCGATCAGGGTGAAGCCTTTGCTGCGTTTCATGGCTGTCACCATTGCCCGTAAGGGATATTGTTCAGGCCttt                                                                         >  6:245085/1‑146 (MQ=255)
                                                                                 gCTGTTGAAGTTGCGCGGGATGGCGATGGTCATCAGCGTGGCTATGATCGCCATAACCTCCAGCAGTTCGATCAGGGTGAAGCCTTTGCTGCGTTTCATGGCTGTCACCATTGGCCGTAAGGTATATTGTTCAGGCCTTTGCCGCGGGc                                                                <  6:60000/149‑1 (MQ=255)
                                                                                 gCTGTTGAAGTAGCGCGGCATGGCGATGGTCATCAGCGTGGCGATGATCGCCATCACCACCAGCAGTTCGATCAGGGTGAAGCCTTTGCTGCGTTTCATGGCTGTCACCATTGCCCGTAAg                                                                                            >  6:194786/1‑121 (MQ=255)
                                                                                 gCTGTTGAAGTAGCGCGGCATGGCGATGGTCATCAGCGTGGCGATGATCGCCATCACCACCAGCAGTTCGATCAGGGTGAAGCCTTTGCTGCGTTTCATGGCTGTCACCATTGCCCGTAAg                                                                                            <  5:194786/121‑1 (MQ=255)
                                                                                 gCTGTTGAAGTAGCGCGGCATGGCGATGGTCATCAGCGTGGCGATGATCGCCATCACCACCAGCAGTTCGATCAGGGTGAAGCCTTTGCTGCGTTTCATGGCTGTCACCATTGCCCGTAAGGGATATTGTTCAGGCCTTTGCCGCGGGc                                                                >  7:218988/1‑149 (MQ=255)
                                                                                  cTGTTGAAGTAGCGCGGCATGGCGATGGTCATCAGCGTGGCGATGATCGCCATCACCACCAGCAGTTCGATCAGGGTGAAGCCTTTGCTGCGTTTCATGGCTGTCACCATTGCCCGTAAGGGATATTGTTCAGGCCTTTGCCGCGGGcc                                                               >  6:229585/1‑149 (MQ=255)
                                                                                  cTGTTGAAGTAGCGCGGCATGGCGATGGTCATCAGCGTGGCGATGATCGCCATCACCACCAGCAGTTCGATCAGGGTGAAGCCTTTGCTGCGTTTCATGGCTGTCACCATTGCCCGTAAGGGATATTGTTCAGGCCTTTGCCGCGGGcc                                                               <  2:57259/149‑1 (MQ=255)
                                                                                     ttGAAGTAGCGCGGCATGGCGATGGTCATCAGCGTGGCGATGATCGCCATCACCACCAGCAGTTCGATCAGGGTGAAGCCTTTGCTGCGTTTCATGGCTGTCTCCATTGCCCGTAAGGGATATTGTTCAGGCCTTTGCCGCGGGCCCt                                                             <  6:4192/148‑1 (MQ=255)
                                                                                       gAAGTAGCGCGGCATGGCGATGGTCATCAGCGTGGCGATGATCGCCATCACCACCAGCAGTTCGATCAGGGTGAAGCCTTTGCTGCGTTTCATGGCTGTCACCATTGCCCGTAAGGGATATTGTTCAGGCCTTTGCCGCGGGCCCTgg                                                           >  5:30173/1‑148 (MQ=255)
                                                                                       gAAGTAGCGCGGCATGGCGATGGTCATCAGCGTGGCGATGATCGCCATCACCACCAGCAGTTCGATCAGGGTGAAGCCTTTGCTGCGTTTCATGGCTGTCACCATTGCCCGTAAGGGATATTGTTCAGGCCTTTGCCGCGGGCCCTgg                                                           >  1:178197/1‑148 (MQ=255)
                                                                                        aaGTAGCGCGGCATGGCGATGGTCATCAGCGTGGCGATGATCGCCATCACCACCAGCAGTTCGATCAGGGTGAAGCCTTTGCTGCGTTTCATGGCTGTCACCATTGCCCGTAAGGGATATTGTTCAGGCCTTTGCCGCGGGCCCTGGaa                                                         >  3:116339/1‑149 (MQ=255)
                                                                                              cgcgTCATGGCGATGGTCATCAGCGTGGCGATGATCGCCATCACCACCAGCAGTTCGATCAGGGTGAAGCCTTTGCTGCGTTTCATGGCTGTCACCATTGCCCGTAAGGGATATTGTTCAGGCCTTTGCCGCGGGCCCTGGAATACACg                                                   >  8:265809/1‑149 (MQ=255)
                                                                                               tcgGCATGGCGATGGTCATCAGCGTGGCGATGATCGCCATCACCACCAGCAGTTCGATCAGGGTGAAGCCTTTGCTGCGTTTCATGGCTGTCACCATTGCCCGTAAGGGATATTGTTCAGGCCTTTGCCGAGGGCCCTGCAATACACGt                                                  >  8:181433/2‑149 (MQ=255)
                                                                                               gcgGCATGGCGATGGTCATCAGCGTGGCGATGATCGCCATCACCACCAGCAGTTCGATCAGGGTGAAGCCTTTGCTGCGTTTCATGGCTGTCACCATTGCCCGTAAGGGATATTGTTCAGGCCTTTGCCGCGGGCCCTGGAATACACGt                                                  >  4:54058/1‑149 (MQ=255)
                                                                                                     tGGCGATGGTCATCAGCGTGGCGATGATCGCCATCACCACCAGCAGTTCGATCAGGGTGAAGCCTTTGCTGCGTTTCATGGCTGTCACCATTGCCCGTGAGGGATATTGTTCAGGCCTTTGCCGCGGGCCCTGGAATACACGTCGAaca                                            >  8:24949/1‑149 (MQ=255)
                                                                                                       gCGATGGTCATCAGCGTGGCGATGATCGCCATCACCACCAGCAGTTCGATCAGGGTGAAGCCTTTGCTGCGTTTCATGGCTGTCACCATTGCCCGTAAGGGATATTGTTCAGGCCTTTGCCGCGGGCCCTGGAATACACGTCGAacac                                           <  7:144550/148‑1 (MQ=255)
                                                                                                       gCGATGGTCATCAGCGTGGCGATGATCGCCATCACCACCAGCAGTTCGATCAGGGTGAAGCCTTTGCTGCGTTTCATGGCTGTCACCATTGCCCGTAAGGGATATTGTTCAGGCCTTTGCCGCGGGCCCTGGAATACACGTCGAacac                                           <  5:229585/148‑1 (MQ=255)
                                                                                                        cGATGGTCATCAGCGTGGCGATGATCGCCATCACCACCAGCAGTTCGATCAGGGTGa                                                                                                                                     <  2:186266/57‑1 (MQ=255)
                                                                                                        cGATGGTCATCAGCGTGGCGATGATCGCCATCACCACCAGCAGTTCGATCAGGGTGa                                                                                                                                     >  1:186266/1‑57 (MQ=255)
                                                                                                              tcatcaGCGTGGCGATGATCGCCATCACCACCAGCAGTTCGATCAGGGTGAAGCCTTTGCTGCGTTTCATGGCTGTCACCATTGCCCGTAAGGGATATTGTTCAGGCCTTTGCCGCGGGCCCTGGAATACACGTCGAACACGTCTTCg                                    <  8:195595/148‑1 (MQ=255)
                                                                                                              tcatcaGCGTGGCGATGATCGCCATCACCACCAGCAGTTCGATCAGGGTGAAGCCTTTGCTGCGTTTCATGGCTGTCACCATTGCCCGTAAGGGATATTGTTCAGGCCTTTGCCGCGGGCCCTGGAATACACGTCGAACACGTCTTCg                                    >  7:195595/1‑148 (MQ=255)
                                                                                                                 tcaGCGTGGCGATGATCGCCATCACCACCAGCAGTTCGATCAGGGTGAAGCCTTTGCTGCGTTTCATGGCTGTCACCATTGCCCGTAAGGGATATTGTTCAGGCCTTTGCCGCGGGCCCTGGAATACACGTCGAACACGTCTTCGCCCt                                >  2:103023/1‑149 (MQ=255)
                                                                                                                 tcaGCGTGGCGATGATCGCCATCACCACCAGCAGTACGATCAGGGTGAAGCCTTTGCTGCGTTTCATGGCTGTCACCATTGCCCGTAAGGGATATTGTTCAGGCCTTTGCCGCGGGCCCTGGAATACACGTCGAACACGTCTTCGCCCt                                <  8:81255/149‑1 (MQ=255)
                                                                                                                 tcaGCGTGGCGATGATCGCCATCACCACCAGCAGTACGATCAGGGTGAAGCCTTTGCTGCGTTTCATGGCTGTCACCATTGCCCGTAAGGGATATTGTTCAGGCCTTTGCCGCGGGCCCTGGAATACACGTCGAACACGTCTTCGCCCt                                <  8:85825/149‑1 (MQ=255)
                                                                                                                  caGCGTGGCGATGATCGCCATCACCCCCAGCAGTTCGATCAGGGTGAAGCCTTTGCTGCGTTTCATGGCTGTCACCATTGCCCGTAAGGGATATTGTTCAGGCCTTTGCCGCGGGCCCTGGAATACACGTCGAACACGTCTTCGCCCTc                               >  5:330350/1‑149 (MQ=255)
                                                                                                                      gTGGCGATGATCGCCATCACCACCAGCAGTTCGATCAGGGTGAAGCCTTTGCTGCGTTTCATGGCTGTCACCATTGCCCGTGAGGGATATTGTTCAGGCCTTTGCCGCGGGCCCTGGAATACACGTCGAACACGTCTTCGCCCTcgcgc                           <  7:24949/149‑1 (MQ=255)
                                                                                                                      gTGGCGATGATCGCCATCACCACCAGCAGTTCGATCAGGGTGAAGCCTTTGCTGCGTTTCATGGCTGTCACCATTGCCCGTAAGGGATATTGTTCAGGCCTTTGCCGCGGGCCCTGGAATACACGTCGAACACGTCTTCGCCCTcgcgc                           >  6:330247/1‑149 (MQ=255)
                                                                                                                         gCGATGATCGCCATCACCACCAGCAGTTCGATCAGGGTGAAGCCTTTGCTGCGTTTCATGGCTGTCACCATTGCCCGTAAGGGATATTGTTCAGGCCTTTGCCGCGGGCCCTGGAATACACGTCGAACACGTCTTCGCCCTCGCGCggg                        <  2:340273/149‑1 (MQ=255)
                                                                                                                         gCGATGATCGCCATCACCACCAGCAGTTCGATCAGGGTGAAGCCTTTGCTGCGTTTCATGGCTGTCACCATTGCCCGTAAGGGATATTGTTCAGGCCTTTGCCGCGGGCCCTGGAATACACGTCGAACACGTCTTCGCCCTCGCGCgg                         <  7:265809/148‑1 (MQ=255)
                                                                                                                                 cGCCATCACCACCAGCAGTTCGATCAGGGTGAAGCCTTTGCTGCGTTTCATGGCTGTCACCATTGCCAGTAAGGGATATTGTTCATGCATTTGCCGCTGGCCCTTGAATACACGTCGAACACGTATTCGCCCTCGCGATGGTTGTCggg                >  1:330403/1‑149 (MQ=255)
                                                                                                                                   ccATCACCACCAGCAGTTCGATCAGGGTGAAGCCTTTGCTGCGTTTCATGGCTGTCACCATTGCCCGTAAGGGATATTGTTCAGGCCTTTGCCGCGGGCCCTGGAATACACGTCGAacac                                           <  4:209119/120‑1 (MQ=255)
                                                                                                                                   ccATCACCACCAGCAGTTCGATCAGGGTGAAGCCTTTGCTGCGTTTCATGGCTGTCACCATTGCCCGTAAGGGATATTGTTCAGGCCTTTGCCGCGGGCCCTGGAATACACGTCGAacac                                           >  3:209119/1‑120 (MQ=255)
                                                                                                                                          caccaGCAGTTCGATCAGGGTGAAGCCTTTGCTGCGTTTCATGGCTGTCACCATTGCCCGTAAGGGATATTGTTCAGGCCTTTGCCGCGGGCCCTGGAATACACGTCGAACACGTCTTCGCCCTCGCGCGGGTTGTCGGGGCTGCTGTc       >  7:230306/1‑149 (MQ=255)
                                                                                                                                           accaGCAGTTCGATCAGGGTGAAGCCTTTGCTGCGTTTCATGGCTGTCACCATTGCCCGTAAGGGATATTGTTCAGGCCTTTGCCGCGGGCCCTGGAATACACGTCGAACACGTCTTCGCCCTCGCGCGGGTTGTCGGGGCTGCTGTCg      >  5:90513/1‑149 (MQ=255)
                                                                                                                                           accaGCAGTTCGATCAGGGTGAAGCCTTTGCTGCGTTTCATGGCTGTCACCATTGCCCGTAAGGGATATTGTTCAGGCCTTTGCCGCGGGCCCTGGAATACACGTCGAACACGTCTTCGCCCTCGCGCGGGTTGTCGGGGCTGCTGTCg      >  5:196231/1‑149 (MQ=255)
                                                                                                                                             cagcagTTCGATCAGGGTGAAGCCTTTGCTGCGTTTCATGGCTGTCACCATTGCCCGTAAGGGATATTGTTCAGGCCTTTGCCGCGGGCCCTGGAATACACGTCGAACACGTCTTCGCCCTCGCGCGGGTTGTCGGGGCTGCTGTCGTa    <  1:323658/149‑1 (MQ=255)
                                                                                                                                               gcagTTCGATCAGGGTGAAGCCTTTGCTGCGTTTCATGGCTGTCACCATTGCCCGTAAGGGATATTGTTCAGGCCTTTGCCGCGGGCCCTGGAATACACGTCGAACACGTCTTCGCCCTCGCGCGGGTTGTCGGGGCTGCTGTCGTAgg  <  6:30173/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                               
GGTATCACCGTAGTAGTGGTCCAGCGCCTCGCGCAGCACTGCCAGGCTCTGACGCAGGGTGGCCTCGCGCGAGCTTTCCAGGCTGTTGAAGTAGCGCGGCATGGCGATGGTCATCAGCGTGGCGATGATCGCCATCACCACCAGCAGTTCGATCAGGGTGAAGCCTTTGCTGCGTTTCATGGCTGTCACCATTGCCCGTAAGGGATATTGTTCAGGCCTTTGCCGCGGGCCCTGGAATACACGTCGAACACGTCTTCGCCCTCGCGCGGGTTGTCGGGGCTGCTGTCGTAGG  >  NC_002947/3941366‑3941657

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