Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,757,897 A→T 5.1% Q390H (CAA→CAT) ‡ ompQ → outer membrane pyoverdine efflux protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,757,8970AT 5.1% 249.0 / 2.1 78Q390H (CAA→CAT) ‡ompQouter membrane pyoverdine efflux protein
Reads supporting (aligned to +/- strand):  ref base A (39/35);  new base T (2/2);  total (41/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.50e-01

GCCTGGTCGCCCCGGTGTTCAACCGAGGCCGTCTGCGCGCCGCGCAGCAACTGGCCGAGGCCGAACAGGAAGAATTGCTGGAGGCTTACCGCACCAGCATCCTGGCCGCTTTCGCCGATGTAGAAAAAGCACTGAACGCCACCCAAGGGGTGGACCGTCAACGCCAGTGGCAGGATCAGGAAGTGGAACAATCACGCATCGCCTTCCAGCTGGCCGAGCGCCGCTACCGCGCCGGCGCCGAGACCTTGCTGACCGTGCTCGACACCCAGCGCAGCCTGTACCAGGCCCAGGAC  >  NC_002947/4757752‑4758044
                                                                                                                                                 |                                                                                                                                                   
gCCTGGTCGCCCCGGTGTTCAACCGAGGCCGTCTGCGCGCCGCGCAGCAACTGGCCGAGGCCGAACAGGAAGAATTGGTGGAGGCTTACCGCACCAGCATCCTGGCCGCTTTCGCCGATGTAGAAAAAGCACTGAACGCCACCCAAggg                                                                                                                                                  >  1:132422/1‑149 (MQ=255)
gCCTGGTCGCCCCGGTGTTCAACCGAGGCCGTCTGCGCGCCGCGCAGCAACTGGCCGAGGCCGAACAGGAAGAATTGCTGGAGGCTTACCGCACCAGCATCCTGGCCGCTTTCGCCGATGTAGAAAAAGCACTGAACGCCACCCAAggg                                                                                                                                                  >  7:316854/1‑149 (MQ=255)
gCCTGGTCGCCCCGGTGTTCAACCGAGGCCGTCTGCGCGCCGCGCAGCAACTGGCCGAGGCCGAACAGGAAGAATTGCTGGAGGCTTACCGCACCAGCATCCTGGCCGCTTTCGCCGATGTAGAAAAAGCACTGAACGCCACCCAAgg                                                                                                                                                   <  4:279757/148‑1 (MQ=255)
 ccTGGTCGCCCCGGTGTTCAACCGAGGCCGTCTGCGCGCCGCGCAGCAACTGGCCGAGGCCGAACAGGAAGAATTGCTGGAGGCTTACCGCACCAGCATCCTGGCCGCTTTCGCCGATGTAGAAAAAGCACTGAACGCCACCCAAgggg                                                                                                                                                 >  1:62309/1‑149 (MQ=255)
      tCGCCCCGGTGTTCAACCGAGGCCGTCTGCGCGCCGCGCAGCAACTGGCCGAGGCCGAACAGGAAGAATTGCTGGAGGCTTACCGCACCAGCATCCTGGCCGCTTTCGCCGATGTAGAAAAAGCACTGAACGCCACCCAAGGTGTGGAc                                                                                                                                            <  2:135595/149‑1 (MQ=255)
        gCCCCGGTGTTCAACCGAGGCCGTCTGCGCGCCGCGCAGCAACTGGCCGAGGCCGAACAGGAAGAATTGCTGGAGGCTTACCGCACCAGCATCCTGGCCGCTTTCGCCGATGTAGAAAAAGCACTGAACGCCACCCCTTGGGTggcg                                                                                                                                            <  2:98855/147‑3 (MQ=255)
        gCCCCGGTGTTCAACCGAGGCCGTCTGCGCGCCGCGCAGCAACTGGCCGAGGCCGAACAGGAAGAATTGCTGGAGGCTTACCGCACCAGCATCCTGGCCGCTTTCGCCGATGTAGAAAAAGCACTGAACGCCACCCCTTGGGTggcg                                                                                                                                            >  1:98855/1‑145 (MQ=255)
        gCCCCGGTGTTCAACCGAGGCCGTCTGCGCGCCGCGCAGCAACTGGCCGAGGCCGAACAGGAAGAATTGCTGGAGGCTTACCGCACCAGCATCCTGGCCGCTTTCGCCGATGTAGAAAAAGCACTGAACGCCACCCAAGGGGTGGAcc                                                                                                                                           <  4:267482/148‑1 (MQ=255)
          cccGGTGTTCAACCGAGGCCGTCTGCGCGCCGCGCAGCAACTGGCCGAGGCCGAACAGGAAGAATTGCTGGAGGCTTACCGCACCAGCATCCTGGCCGCTTTCGCCGATGTAGAAAAAGCACTGAACGCCACCCAAGGGGTGGACCGt                                                                                                                                         >  6:316019/1‑148 (MQ=255)
             ggTGTTCAACCGAGGCCGTCTGCGCGCCGCGCAGCAACTGGCCGAGGCCGAACAGGAAGAATTGCTGGAGGCTTACCGCACCAGCATCCTGGCCGCTTTCGCCGATGTAGAAAAAGCAGTGAACGCCACCCAAGGGGTGGACCGTCAAc                                                                                                                                     <  1:33250/149‑1 (MQ=255)
             ggTGTTCAACCGAGGCCGTCTGCGCGCCGCGCAGCAACTGGCCGAGGCCGAACAGGAAGAATTGCTGGAGGCTTACCGCACCAGCATCCTGGCCGCTTTCGCCGATGTAGAAAAAGCACTGAACGCCACCCAAGGGGTGGACCGTCAAc                                                                                                                                     <  7:306179/149‑1 (MQ=255)
                 ttCAACCGAGGCCGTCTGCGCGCCGCGCAGCAACTGGCCGAGGCCGAACAGGAAGAATTGCTGGAGGCTTACCGCACCAGCATCCTGGCCGCTTTCGCCGATGTAGAAAAAGCACTGAACGCCACCCAAGGGGTGGACCGTCAACGc                                                                                                                                   >  5:303805/1‑147 (MQ=255)
                 ttCAACCGAGGCCGTCTGCGCGCCGCGCAGCAACTGGCCGAGGCCGAACAGGAAGAATTGCTGGAGGCTTACCGCACCAGCATCCTGGCCGCTTTCGCCGATGTAGAAAAAGCACTGAACGCCACCCAAGGGGTGGACCGTCAACGc                                                                                                                                   >  4:352854/1‑147 (MQ=255)
                   cAACCGAGGCCGTCTGCGCGCCGCGCAGCAACTGGCCGAGGCCGAACAGGAAGAATTGCTGGAGGCTTACCGCACCAGCATCCTGGCCGCTTTCGCCGATGTAGAAAAAGCACTGAACGCCACCCCTTGGGTggcgtt                                                                                                                                          >  1:103597/1‑134 (MQ=255)
                   cAACCGAGGCCGTCTGCGCGCCGCGCAGCAACTGGCCGAGGCCGAACAGGAAGAATTGCTGGAGGCTTACCGCACCAGCATCCTGGCCGCTTTCGCCGATGTAGAAAAAGCACTGAACGCCACCCCTTGGGTggcgtt                                                                                                                                          <  2:103597/138‑5 (MQ=255)
                      ccGAGGCCGTCTGCGCGCCGCGCAGCAACTGGCCGAGGCCGAACAGGAAGAATTGCTGGAGGCTTACCGCACCAGCATCCTGGCCGCTTTCGCCGATGTAGAAAAAGCACTGAACGCCACCCAAGGGGTGGACCGTCAACGCCAGTgg                                                                                                                             >  2:41433/1‑148 (MQ=255)
                          ggCCGTCTGCGCGCCGCGCAGCAACTGGCCGAGGCCGAACAGGAAGAATTGCTGGAGGCTTACCGCACCAGCATCCTGGCCGCTTTCGCCGATGTAGAAAAAGCACTGAACGCCACCCAAGGGGTGGACCGTCAACGCCAGTGGCAGGa                                                                                                                        >  5:378289/1‑149 (MQ=255)
                           gCCGTCTGCGCGCCGCGCAGCAACTGGCCGAGGCCGAACAGGAAGAATTGCTGGAGGCTTACCGCACCAGCATCCTGGCCGCTTTCGCCGATGTAGAAAAAGCACTGAACGCCACCCAAGGGGTGGACCGTCAACGCCAGTGGCAGGAt                                                                                                                       <  8:316854/149‑1 (MQ=255)
                               tCTGCGCGCCGCGCAGCAACTGGCCGAGGCCGAACAGGAAGAATTGCTGGAGGCTTACCGCACCAGCATCCTGGCCGCTTTCGCCGATGTAGAAAAAGCACTGAACGCCACCCAAGGGGTGGACCGTCAACGCCAGTGGCAGGATCAgg                                                                                                                   <  4:215036/149‑1 (MQ=255)
                                   cgcgccgcgcAGCAACTGGCCGAGGCCGAACAGGAAGAATTGCTGGAGGCTTACCGCACCAGCATCCTGGCCGCTTTCGCCGATGTAGAAAAAGCACTGAACGCCACCCAAGGGGTGGACCGTCAACGCCAGTGGCAGGATCAGGAAg                                                                                                                <  3:363314/148‑1 (MQ=255)
                                    gcgccgcgcAGCAACTGGCCGAGGCCGAACAGGAAGAATTGCTGGAGGCTTACCGCACCAGCATCCTGGCCGCTTTCGCCGATGTAGAAAAAGCACTGAACGCCACCCAAGGGGTGGACCGTCAACGCCAGTGGCAGGATCAGGAAGt                                                                                                               <  6:303805/148‑1 (MQ=255)
                                       ccgcgcAGCAACTGGCCGAGGCCGAACAGGAAGAATTGCTGGAGGCTTACCGCACCAGCATCCTGGCCGCTTTCGCCGATGTAGAAAAAGCACTGAACGCCACCCAAGGGGTGGACCGTCAACGCCAGTGGc                                                                                                                            >  5:318668/1‑132 (MQ=255)
                                       ccgcgcAGCAACTGGCCGAGGCCGAACAGGAAGAATTGCTGGAGGCTTACCGCACCAGCATCCTGGCCGCTTTCGCCGATGTAGAAAAAGCACTGAACGCCACCCAAGGGGTGGACCGTCAACGCCAGTGGc                                                                                                                            <  6:318668/132‑1 (MQ=255)
                                       ccgcgcAGCAACTGGCCGAGGCCGAACAGGAAGAATTGCTGGAGGCTTACCGCACCAGCATCCTGGCCGCTTTCGCCGATGTAGAAAAAGCACTGAACGCCACCCAAGGGGTGGACCGTCAACGCCAGTGGCAGGATCAGGAAGTGGaa                                                                                                           <  3:395894/149‑1 (MQ=255)
                                        cgcgcAGCAACTGGCCGAGGCCGAACAGGAAGAATTGCTGGAGGCTTACCGCACCAGCATCCTGGCCGCTTTCGCCGATGTAGAAAAAGCACTGAACGCCACCCAAGGGGTGGACCGTCAACGCCAGTGGCAGGATCAGGAAGTGGaa                                                                                                           >  4:155820/1‑148 (MQ=255)
                                        cgcgcAGCAACTGGCCGAGGCCGAACAGGAAGAATTGCTGGAGGCTTACCGCACCAGCATCCTGGCCGCTTTCGCCGATGTAGAAAAAGCACTGAACGCCACCCAAGGGGTGGACCGTCAACGCCAGTGGCAGGATCAGGAAGTGGAAc                                                                                                          <  4:175991/149‑1 (MQ=255)
                                            cagcaACTGGCCGAGGCCGAACAGGAAGAATTGCTGGAGGCTTACCGCACCAGCATCCTGGCCGCTTTCGCCGATGTAGAAAAAGCACTGAACGCCACCCAAGGGGTGGACCGTCAACGCCAGTGGCAGGATCAGGAAGTGGAACAAt                                                                                                       >  2:306665/1‑148 (MQ=255)
                                             agcaACTGGCCGAGGCCGAACAGGAAGAATTGCTGGAGGCTTACCGCACCAGCATCCTGGCCGCTTTCGCCGATGTAGAAAAAGCACTGAACGCCACCCAAGGGGTGGACCGTCAACGCCAGTGGCAGGATCAGGAAGTGGAACAATCa                                                                                                     <  4:367258/149‑1 (MQ=255)
                                              gcaACTGGCCGAGGCCGAACAGGAAGAATTGCTGGAGGCTTACCGCACCAGCATCCTGGCCGCTTTCGCCGATGTAGAAAAAGCACTGAACGCCACCCAAGGGGTGGACCGTCAACGCCAGTGGCAGGATCAGGAAGTGGAACAATCAc                                                                                                    <  6:268438/149‑1 (MQ=255)
                                               caACTGGCCGAGGCCGAACAGGAAGAATTGCTGGAGGCTTACCGCACCAGCATCCTGGCCGCTTTCGCCGATGTAGAAAAAGCACTGAACGCCACCCAAGGGGTGGACCGTCAACGCCAGTGGCAGGATCAGGAAGTGGAACAATCACg                                                                                                   >  1:276767/1‑149 (MQ=255)
                                                       cgaggccgaACAGGAAGAATTGCTGGAGGCTTACCGCACCAGCATCCTGGCCGCTTTCGCCGATGTAGAAAAAGCACTGAACGCCACCCAAGGGGTGGACCGTCAACGCCAGTGGCAGGATCAGGAAGTGGAACAATCACGCATCGCCt                                                                                           >  8:66459/1‑149 (MQ=255)
                                                          ggccgaACAGGAAGAATTGCTGGAGGCTTACCGCACCAGCATCCTGGCCGCTTTCGCCGATGTAGAAAAAGCACTGAACGCCACCCAAGGGGTGGACCGTCAACGCCAGTGGCAGGATCAGGAAGTGGAACAATCACGCATCGCCTTcc                                                                                        <  5:316019/149‑1 (MQ=255)
                                                          ggccgaACAGGAAGAATTGCTGGAGGCTTACCGCACCAGCATCCTGGCCGCTTTCGCCGATGTAGAAAAAGCACTGAACGCCACCCAAGGGGTGGACCGTCAACGCCAGTGGCAGGATCAGGAAGTGGAACAATCACGCATCGCCTTc                                                                                         >  6:166131/1‑148 (MQ=255)
                                                           gccgaACAGGAAGAATTGCTGGAGGCTTACCGCACCAGCATCCTGGCCGCTTTCGCCGATGTAGAAAAAGCTCTGAACGCCACCCAAGGGGTGGACCGTCAACGCCAGTGGCAGGATCAGGAAGt                                                                                                               >  7:58823/1‑125 (MQ=255)
                                                                    gaagaaTTGCTGGAGGCTTACCGCACCAGCATCCTGGCCGCTTTCGCCGATGTATAAAAAGGCCTGAACGCCACCCAAGGGGTGGGCCGTCAACGCCAGGGGCAGGCTCAGGAAGTGGGACAATCACGCATCGCCGTCCCGCTGGCCGa                                                                              >  1:292272/1‑149 (MQ=255)
                                                                    gaagaaTTGCTGGAGGCTTACCGCACCAGCATCCTGGCCGCTTTCGCCGATGTAGAAAAAGCACTGAACGCCACCCAAGGGGTGGACCGTCAACGCCAGTGGCAGGATCAGGAAGTGGAACAATCACGCATc                                                                                               <  5:42592/132‑1 (MQ=255)
                                                                    gaagaaTTGCTGGAGGCTTACCGCACCAGCATCCTGGCCGCTTTCGCCGATGTAGAAAAAGCACTGAACGCCACCCAAGGGGTGGACCGTCAACGCCAGTGGCAGGATCAGGAAGTGGAACAATCACGCATc                                                                                               >  6:42592/1‑132 (MQ=255)
                                                                    gaagaaTTGCTGGAGGCTTACCGCACCAGCATCCTGGCCGCTTTCGCCGATGTAGAAAAAGCACTGAACGCCACCCAAGGGGTGGACCGTCAACGCCAGTGGCAGGATCAGGAAGTGGAACAATCACGCATc                                                                                               <  1:262887/132‑1 (MQ=255)
                                                                    gaagaaTTGCTGGAGGCTTACCGCACCAGCATCCTGGCCGCTTTCGCCGATGTAGAAAAAGCACTGAACGCCACCCAAGGGGTGGACCGTCAACGCCAGTGGCAGGATCAGGAAGTGGAACAATCACGCATc                                                                                               >  2:262887/1‑132 (MQ=255)
                                                                    gaagaaTTGCTGGAGGCTTACCGCACCAGCATCCTGGCCGCTTTCGCCGATGTAGAAAAAGCACTGAACGCCACCCAAGGGGTGGACCGTCAACGCCAGTGGCAGGATCAGGAAGTGGAACAATCACGCATCGCCTTCCAGCTGGcc                                                                                >  7:321150/1‑147 (MQ=255)
                                                                        aaTTGCTGGAGGCTTACCGCACCAGCATCCTGGCCGCTTTCGCCGATGTAGAAAAAGCACTGAACGCCACCCAAGGGGTGGACCGTCAACGCCAGTGGCAGGATCAGGAAGTGGAACAATCACGCATCGCCTTCCAGCTGGCCGAgcgc                                                                          <  8:321150/149‑1 (MQ=255)
                                                                             cTGGAGGCTTACCGCACCAGCATCCTGGCCGCTTTCGCCGATGTAGAAAAAGCACTGAACGCCACCCAAGGGGTGGACCGTCAACGCCAGTGGCAGGATCAGGAAGTGGAACAATCACGCATCGCCTTCCAGCTGGCCGAGCGCCGCTa                                                                     >  6:153495/1‑149 (MQ=255)
                                                                               ggAGGCTTACCGCACCAGCATCCTGGCCGCTTTCGCCGATGTAGAAAAAGCACTGAACGCCACCCAAGGGGTGGACCGTCAACGCCAGTGGCAGGATCAGGAAGTGGAACAATCACGCATCGCCTTCCAGCTGGCCGAGCGCCGCTAcc                                                                   >  1:392042/1‑149 (MQ=255)
                                                                                            aCCAGCATCCTGGCCGCTTTCGCCGATGTAGAAAAAGCACTGAACGCCACCCAAGGGGTGGACCGTCAACGCCAGTGGCAGGATCAGGAAGTGGAACAATCACGCATCGCCTTCCAGCTGGCCGAGCGCCGCTACCGCGCCGGCGCCga                                                      <  5:153495/149‑1 (MQ=255)
                                                                                              cAGCATCCTGGCCGCTTTCGCCGATGTAGAAAAAGCACTGAACGCCACCCAAGGGGTGGACCGTCAACGCCAGTGGCAGGATCAGGAAGTGGAACAATCACGCATCGCCTTCCAGCTGGCCGAGCGCCGCTACCGCGCCGGCGCCgaga                                                    >  1:271082/1‑149 (MQ=255)
                                                                                              cAGCATCCTGGCCGCTTTCGCCGATGTAGAAAAAGCACTGAACGCCACCCAAGGGGTGGACCGTCAACGCCAGTGGCAGGATCAGGAAGTGGAACAATCACGCATCGCCTTCCAGCTGGCCGAGCGCCGCTACCGCGCCGGCGCCgaga                                                    >  5:183780/1‑149 (MQ=255)
                                                                                              cAGCATCCTGGCCGCTTTCGCCGATGTAGAAAAAGCACTGAACGCCACCCAAGGGGTGGACCGTCAACGCCAGTGGCAGGATCAGGAAGTGGAACAATCACGCATCGCCTTCCAGCTGGCCGAGCGCCGCTACCGCGCCGGCGCCgaga                                                    >  8:85805/1‑149 (MQ=255)
                                                                                                gCATCCTGGCCGCTTTCGCCGATGTAGAAAAAGCACTGAACGCCACCCAAGGGGTGGACCGTCAACGCCAGTGGCAGGATCAGGAAGTGGAACAATCTCGCATCGCCTTCCAGCTGGCCGAGCGCCGCTACCGCGCCGGCGCCGAGAcc                                                  <  1:306665/149‑1 (MQ=255)
                                                                                                    ccTGGCCGCTTTCGCCGATGTAGAAAAAGCACTGAACGCCACCCAAGGGGTGGACCGTCAACGCCAGTGGCAGGATCAGGAAGTGGAACAATCACGCATCGCCTTCCAGCTGGCCGAGCGCCGCTACCGCGCCGGCGCCGAGACCTTGc                                              >  8:4563/1‑149 (MQ=255)
                                                                                                         ccGCTTTCGCCGATGTAGAAAAAGCACTGAACGCCACCCAAGGGGTGGACCGTCAACGCCAGTGGCAGGATCAGGAAGTGGAACAATCACGCATCGCCTTCCAGCTGGCCGAGCGCCGCTACCGCGCCGGCGCCGAGACCTTGCTGAcc                                         <  7:85805/149‑1 (MQ=255)
                                                                                                         ccGCTTTCGCCGATGTAGAAAAAGCACTGAACGCCACCCAAGGGGTGGACCGTCAACGCCAGTGGCAGGATCAGGAAGTGGAACAATCACGCATCGCCTTCCAGCTGGCCGAGCGCCGCTACCGCGCCGGCGCCGAGACCTTGCTGAcc                                         <  7:66459/149‑1 (MQ=255)
                                                                                                         ccGCTTTCGCCGATGTAGAAAAAGCACTGAACGCCACCCAAGGGGTGGACCGTCAACGCCAGTGGCAGGATCAGGAAGTGGAACAATCACGCATCGCCTTCCAGCTGGCCGAGCGCCGCTACCGCGCCGGCGCCGAGACCTTGCTGAcc                                         <  1:41433/149‑1 (MQ=255)
                                                                                                           gCTTTCGCCGATGTAGAAAAAGCACTGAACGCCACCCAAGGGGTGGACCGTCAACGCCAGTGGCAGGATCAGGAAGTGGAACAATCACGCATCGCCTTCCAGCTGGCCGAGCGCCGCTACCGCGCCGGCGCCGAGACCTTGCTGACCGt                                       >  8:297629/1‑149 (MQ=255)
                                                                                                              ttCGCCGATGTAGAAAAAGCACTGAACGCCACCCAAGGGGTGGACCGTCAACGCCAGTGGCAGGATCAGGAAGTGGAACAATCAc                                                                                                    <  2:264846/85‑1 (MQ=255)
                                                                                                              ttCGCCGATGTAGAAAAAGCACTGAACGCCACCCAAGGGGTGGACCGTCAACGCCAGTGGCAGGATCAGGAAGTGGAACAATCACGCATCGCCTTCCAGCTGGCCGAGCGCCGCTACCGCGCCGGCGCCGAGACCTTGCTGACCGTGCt                                    >  8:199206/1‑149 (MQ=255)
                                                                                                                 gCCGATGTAGAAAAAGCACTGAACGCCACCCAAGGGGTGGACCGTCAACGCCAGTGGCAGGATCAGGAAGTGGAACAATCa                                                                                                     <  4:341733/81‑1 (MQ=255)
                                                                                                                 gCCGATGTAGAAAAAGCACTGAACGCCACCCAAGGGGTGGACCGTCAACGCCAGTGGCAGGATCAGGAAGTGGAACAATCa                                                                                                     >  3:341733/1‑81 (MQ=255)
                                                                                                                  ccGATGTAGAAAAAGCACTGAACGCCACCCAAGGGGTGGACCGTCAACGCCAGTGGCAGGATCAGGAAGTGGAACAATCACGCATCGCCTTCCAGCTGGCCGAGCGCCGCTACCGCGCCGGCGCCGAGACCTTGCTGACCGTGCTCGa                                 >  4:326389/1‑148 (MQ=255)
                                                                                                                   cGATGTAGAAAAAGCACTGAACGCCACCAAAGGGGTGGACCGTCAACGCAAGTGGCATGATCAGGAAGTGGAACAATCACGCATCGCCTTCCCGCTGGCCCAGCGCCCCAACCCCGCCGGCGCCGACACCATGCTGACACTGCTCGcca                               >  6:268734/1‑149 (MQ=255)
                                                                                                                     aTGTAGAAAAAGCACTGAACGCCACCCAAGGGGTGGACCGTCAACGCCAGTGGCAGGATCAGGAAGTGGAACAATCACGCATCGCCTTCCAGCTGGCCGAGCGCCGCTACCGCGCCGGCGCCGAGACCTTGCTGACCGTGCTCGACAcc                             >  2:364749/1‑149 (MQ=255)
                                                                                                                      tGTAGAAAAAGCACTGAACGCCACCCAAGGGGTGGACCGTCAACGCCAGTGGCAGGATCAGGAAGTGGAACAATCACGCATCGCCTTCCAGCTGGCCGAGCGCCGCTACCGCGCCGGCGCCGAGACCTTGCTGACCGTGCTCGACAccc                            <  2:276767/149‑1 (MQ=255)
                                                                                                                        tAGAAAAAGCACTGAACGCCACCCAAGGGGTGGACCGTCAACGCCAGTGGCAGGATCAGGAAGTGGAACAATCACGCATCGCCTTCCAGCTGGCCGAGCGCCGCTACCGCGCCGGCGCCGAGACCTTGCTGACCGTGCTCGACACCCAg                          >  2:2661/1‑149 (MQ=255)
                                                                                                                         aGAAAAAGCACTGAACGCCACCCAAGGGGTGGACCGTCAACGCCAGTGGCAGGATCAGGAAGTGGAACAATCACGCATCGCCTTCCAGCTGGCCGAGCGCCGCTACCGCGCCGGCGCCGAGACCTTGCTGACCGTGCTCGACACCCAgc                         >  1:219065/1‑149 (MQ=255)
                                                                                                                          gAAAAAGCACTGAACGCCACCCAAGGGGTGGACCGTCAACGCCAGTGGCAGGATCAGGAAGTGGAACAATCACGCATCGCCTTCCAGCTGGCCGAGCGCCGCTACCGCGCCGGCGCCGAGACCTTGCTGACCGTGCTCGACACCCAgcg                        >  1:348514/1‑149 (MQ=255)
                                                                                                                            aaaaGCACTGAACGCCACCCAAGGGGTGGACCGTCAACGCCAGTGGCAGGATCAGGAAGTGGAACAATCACGCATCGCCTTCCAGCTGGCCGAGCGCCGCTACCGCGCCGGCGCCGAGACCTTGCTGACCGTGCTCGACACCCAgcg                        >  1:98070/1‑147 (MQ=255)
                                                                                                                              aaGCACTGAACGCCACCCAAGGGGTGGACCGTCAACGCCAGTGGCAGGATCAGGAAGTGGAACAATCACGCATCGCCTTCCAGCTGGCCGAGCGCCGCTACCGCGCCGGCGCCGAGACCTTGCTGACCGTGCTCGACACCCAGCGCAGc                    <  7:4563/149‑1 (MQ=255)
                                                                                                                              aaGCACTGAACGCCACCCAAGGGGTGGACCGTCAACGCCAGTGGCAGGATCAGGAAGTGGAACAATCACGCATCGCCTTCCAGCTGGCCGAGCGCCGCTACCGCGCCGGCGCCGAGACCTTGCTGACCGTGCTCGACACCCAGCGCAGc                    <  3:352854/149‑1 (MQ=255)
                                                                                                                               aGCACTGAACGCCACCCAAGGGGTGGACCGTCAACGCCAGTGGCAGGATCAGGAAGTGGAAc                                                                                                          >  6:373267/1‑62 (MQ=255)
                                                                                                                               aGCACTGAACGCCACCCAAGGGGTGGACCGTCAACGCCAGTGGCAGGATCAGGAAGTGGAAc                                                                                                          <  5:373267/62‑1 (MQ=255)
                                                                                                                                    tGAACGCCACCCAAGGGGTGGACCGTCAACGCCAGTGGCAGGATCAGGAAGTGGAACAATCACGCATCGCCTTCCAGCTGGCCGAGCGCCGCTACCGCGCCGGCGCCGAGACCTTGCTGACCGTGCTCGACACCCAGCGCAGCCTGTAc              >  2:394668/1‑149 (MQ=255)
                                                                                                                                      aaCGCCACCCAAGGGGTGGACCGTCAACGCCAGTGGCAGGATCAGGAAGTGGAACAATCACGCATCGCCTTCCAGCTGGCCGAGCGCCGCTACCGCGCCGGCGCCGAGACCTTGCTGACCGTGCTCGACACCCAGCGCAGCCTGTAcc             <  2:292272/148‑1 (MQ=255)
                                                                                                                                          ccACCCAAGGGGTGGACCGTCAACGCCAGTGGCAGGATCAGGAAGTGGAACAATCACGCATCGCCTTCCAGCTGGCCGAGCGCCGCTACCGCGCCGGCGCCGAGACCTTGCTGACCGTGCTCGACACCCAGCGCAGCCTGTACCAGGcc        <  2:271082/149‑1 (MQ=255)
                                                                                                                                           cACCCAAGGGGTGGACCGTCAACGCCAGTGGCAGGATCAGGAAGTGGAACAATCACGCATCGCCTTCCAGCTGGCCGAGCGCCGCTACCGCGCCGGCGCCGAGACCTTGCTGACCGTGCTCGACACCCAGCGCAGCCTGTACCAGGccc       <  2:132422/149‑1 (MQ=255)
                                                                                                                                           cACCCAAGGGGTGGACCGTCAACGCCAGTGGCAGGATCAGGAAGTGGAACAATCACGCATCGCCTTCCAGCTGGCCGAGCGCCGCTACCGCGCCGGCGCCGAGACCTTGCTGACCGTGCTCGACACCCAGCGCAGCCTGTACCAGGcc        >  5:339399/1‑148 (MQ=255)
                                                                                                                                              ccAAGGGGTGGACCGTCAACGCCAGTGGCAGGATCAGGAAGTGGAACAATCACGCATCGCCTTCCAGCTGGCCGAGCGCCGCTACCGCGCCGGCGCCGAGACCTTGCTGACCGTGCTCGACACCCAGCGCAGCCTGTACCAGGCCCAgg    <  3:155820/149‑1 (MQ=255)
                                                                                                                                               cAAGGGGTGGACCGTCAACGCCAGTGGCAGGATCAGGAAGTGGAACAATCACGCATCGCCTTCCAGCTGGCCGAGCGCCGCTACCGCGCCGGCGCCGAGACCTTg                                               <  1:173974/105‑1 (MQ=255)
                                                                                                                                               cAAGGGGTGGACCGTCAACGCCAGTGGCAGGATCAGGAAGTGGAACAATCACGCATCGCCTTCCAGCTGGCCGAGCGCCGCTACCGCGCCGGCGCCGAGACCTTg                                               >  2:173974/1‑105 (MQ=255)
                                                                                                                                               cAAGGGGTGGACCGTCAACGCCAGTGGCAGGATCAGGAAGTGGAACAATCACGCATCGCCTTCCAGCTGGCCGAGCGCCGCTACCGCGCCGGCGCCGAGACCTTGCTGACCGTGCTCGACACCCAGCGCAGCCTGTACCAGGCCCAGGa   <  1:364749/149‑1 (MQ=255)
                                                                                                                                                aaGGGGTGGACCGTCAACGCCAGTGGCAGGATCAGGAAGTGGAACAATCACGCATCGCCTTCCAGCTGGCCGAGCGCCGCTACCGCGCCGGCGCCGAGACCTTGCTGACCGTGCTCGACACCCAGCGCAGCCTGTACCAGGCCCAGGAc  <  3:326389/149‑1 (MQ=255)
                                                                                                                                                aaGGGGTGGACCGTCAACGCCAGTGGCAGGATCAGGAAGTGGAACAATCACGCATCGCCTTCCAGCTGGCCGAGCGCCGCTACCGCGCCGGCGCCGAGACCTTGCTGACCGTGCTCGACACCCAGCGCAGCCTGTACCAGGCCCAGGAc  >  5:371719/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
GCCTGGTCGCCCCGGTGTTCAACCGAGGCCGTCTGCGCGCCGCGCAGCAACTGGCCGAGGCCGAACAGGAAGAATTGCTGGAGGCTTACCGCACCAGCATCCTGGCCGCTTTCGCCGATGTAGAAAAAGCACTGAACGCCACCCAAGGGGTGGACCGTCAACGCCAGTGGCAGGATCAGGAAGTGGAACAATCACGCATCGCCTTCCAGCTGGCCGAGCGCCGCTACCGCGCCGGCGCCGAGACCTTGCTGACCGTGCTCGACACCCAGCGCAGCCTGTACCAGGCCCAGGAC  >  NC_002947/4757752‑4758044

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