Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,331,667 G→A 5.2% G123S (GGC→AGC)  PP_3803 → cation ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,331,6670GA 5.2% 298.3 / 3.2 95G123S (GGC→AGC) PP_3803cation ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base G (41/49);  new base A (3/2);  total (44/51)
Fisher's exact test for biased strand distribution p-value = 6.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.25e-01

TGCTGGGTGAGCCCTATACCGGCCCTTACGGCAGCCTGTTCGGCTACTGCCTGCTGTTCGGCATCCTGCTCAACTTCCTGCGCAACCGTACCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAGATCAACGTGCACATCCTCGAGAACGTGCTGTTCGGTTCGGTGCTGACCGTTAGTGCCCAGGACTTGCTGGTACTGGGCATCGTCGCGGTGCTGGTGCTTGCCTTGGCGCTGCCGCT  >  NC_002947/4331521‑4331815
                                                                                                                                                  |                                                                                                                                                    
tgctgGGTGAGCCCTATACCGGCCCTTACGGCAGCCTGTTCGGCTACTGCCTGCTGTTCGGCATCCTGCTCAACTTCCTGCGCAACCGTACCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCCGTGTCGCTGGCGCTTGgc                                                                                                                                                    >  3:348979/1‑149 (MQ=255)
      gTGAGCCCCATACCGGCCCTTACGGCAGCCTGTTCGGCTACTGCCTGCTGTTCGGCATCCTGCTCAACTTCCTGCGCAACCGTACCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGc                                                                                                                                              <  5:351969/149‑1 (MQ=255)
           cccTATACCGGCCCTTACGGCAGCCTGTTCGGCTACTGCCTGCTGTTCGGCATCCTGCTCAACTTCCTGCGCAACCGTACCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCctgct                                                                                                                                         >  6:76499/1‑149 (MQ=255)
           cccTATACCGGCCCTTACGGCAGCCTGTTCGGCTACTGCCTGCTGTTCGGCATCCTGCTCAACTTCCTGCGCAACCGTACCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCctgct                                                                                                                                         >  5:81993/1‑149 (MQ=255)
                    ggCCCTTACGGCAGCCTGTTCGGCTACTGCCTGCTGTTCGGCATCCTGCTCAACTTCCTGCGCAACCGTACCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCt                                                                                                                                <  6:81993/149‑1 (MQ=255)
                     gCCCTTACGGCAGCCTGTTCGGCTACTGCCTGCTGTTCGGCATCCTGCTCAACTTCCTGCGCAACCGTACCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTg                                                                                                                               <  5:76499/149‑1 (MQ=255)
                     gCCCTTACGGCAGCCTGTTCGGCTACTGCCTGCTGTTCGGCATCCTGCTCAACTTCCTGCGCAACCGTACCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTg                                                                                                                               <  1:99252/149‑1 (MQ=255)
                        cTTACGGCAGCCTGTTCGGCTACTGCCTGCTGTTCGGCATCCTGCTCAACTTCCTGCGCAACCGTACCGGGCTGTCGCCGGATACACTGATTTGCGTGTTACTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTggca                                                                                                                            >  6:316941/1‑149 (MQ=255)
                         ttACGGCAGCCTGTTCGGCTACTGCCTGCTGTTCGGCATCCTGCTCAACTTCCTGCGCAACCGTACCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTggcag                                                                                                                           >  8:289988/1‑149 (MQ=255)
                           aCGGCAGCCTGTTCGGCTACTGCCTGCTGTTCGGCATCCTGCTCAACTTCCTGCGCAACCGTACCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTggcag                                                                                                                           <  1:108181/147‑1 (MQ=255)
                             ggCAGCCTGTTCGGCTACTGCCTGCTGTTCGTCATCCTGCTCAACTTCCTGCGCAACCGTACCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTGGCAGGCaa                                                                                                                       <  5:387373/149‑1 (MQ=255)
                             ggCAGCCTGTTCGGCTACTGCCTGCTGTTCGGCATCCTGCTCAACTTCCTGCGCAACCGTACCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTggcaggc                                                                                                                         >  6:382460/1‑147 (MQ=255)
                              ggagacTGTTTCGCTACTGCCTGCTGTTCGGCATACTGCTCAACTTCCTGCGCAACCGTACCGGGCTGTCGCCGGTTACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTGGCAGGCaa                                                                                                                       <  4:302895/143‑1 (MQ=255)
                                aGCCTGTTCGGCTACTGCCTGCTGTTCGGCATCCTGCTCAACTTCCTGCGCAACCGTACCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAGa                                                                                                                     >  3:291943/1‑148 (MQ=255)
                                  ccTGTTCGGCTACTGCCTGCTGTTCGGCATCCTGCTCAACTTCCTGCGCAACCGTACCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAGATc                                                                                                                   >  4:320398/1‑148 (MQ=255)
                                         ggCTACTGCCTGCTGTTCGGCATCCTGCTCAACTTCCTGCGCAACCGTACCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAGATCAACGTGc                                                                                                            <  4:291943/148‑1 (MQ=255)
                                          gCTACTGCCTGCTGTTCGGCATCCTGCTCAACTTCCTGCGCAACCGTACCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAGATCAACGTGca                                                                                                           >  8:218984/1‑148 (MQ=255)
                                            tACTGCCTGCTGTTCGGCATCCTGCTCAACTTCCTGCGCAACCGTACCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAGATCAACGTGCACAt                                                                                                        >  5:337599/1‑149 (MQ=255)
                                            tACTGCCTGCTGTTCGGCATCCTGCTCAACTTCCTGCGCAACCGTACCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCAGGCGCTTGGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAGATCAACGTGcacc                                                                                                         >  4:234095/1‑147 (MQ=255)
                                              ctgcctgcTGTTCGGCATCCTGCTCAACTTCCTGCGCAACCGTACCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAGATCAACGTGCACATcc                                                                                                      >  2:64086/1‑149 (MQ=255)
                                              ctgcctgcTGTTCGGCATCCTGCTCAACTTCCTGCGCAACCGTACCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAGATCAACGTGCACATc                                                                                                       <  3:320398/148‑1 (MQ=255)
                                                 cctgcTGTTCGGCATCCTGCTCAACTTCCTGCGCAACCGTACCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAGATCAACGTGCACATCCTc                                                                                                    <  3:125318/148‑1 (MQ=255)
                                                   tgctgTTCGGCATCCTGCTCAACTTCCTGCGCAACCGTACCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAGATCAACGTGCACATCCTc                                                                                                    >  2:192735/1‑146 (MQ=255)
                                                           ggCATCCTGCTCAACTTCCTGCGCAACCGTACCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAGATCAACGTGcaca                                                                                                         <  1:174282/133‑1 (MQ=255)
                                                           ggCATCCTGCTCAACTTCCTGCGCAACCGTACCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAGATCAACGTGcaca                                                                                                         >  2:174282/1‑133 (MQ=255)
                                                           ggCATCCTGCTCAACTTCCTGCGCAACCGTACCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAGATCAACGTGCACATCCTCGAGAACGTGCt                                                                                         >  3:187257/1‑149 (MQ=255)
                                                            gCATCCTGCTCAACTTCCTGCGCAACCGTACCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAGATCAACGTGCACATCCTCGAGAACGTGCTg                                                                                        <  3:99171/149‑1 (MQ=255)
                                                             cATCCTGCTCAACTTCCTGCGCAACCGTACCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAGATCAACGTGCACATCCTCGAGAACGTGCTgt                                                                                       >  3:387303/1‑149 (MQ=255)
                                                              aTCCTGCTCAACTTCCTGCGCAACCGTACCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTAGGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAGATCAACGTGCACATCCTCGAGAACGTGCTgtt                                                                                      >  8:151468/1‑149 (MQ=255)
                                                               tCCTGCTCAACTTCCTGCGCAACCGTACCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAGATCAACGTGCACATCCTCGAGAACGTGCTgttc                                                                                     >  3:123005/1‑149 (MQ=255)
                                                                 cTGCTCAACTTCCTGCGCAACCGTACCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAGATCAACGTGCACATCCTCGAGAACGTGCTgttcgg                                                                                   <  4:187257/149‑1 (MQ=255)
                                                                 cTGCTCAACTTCCTGCGCAACCGTACCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAGATCAACGTGCACATCCTCGAGAACGTGCTgttcg                                                                                    <  1:227662/148‑1 (MQ=255)
                                                                 cTGCTCAACTTCCTGCGCAACCGTACCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAGATCAACGTGCACATCCTCGAGAACGTGCTgttcg                                                                                    <  7:289988/148‑1 (MQ=255)
                                                                 cTGCTCAACTTCCTGCGCAACCGTACCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTAGGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAGATCAACGTGCACATCCTCGAGAACGTGCTgttcg                                                                                    <  7:151468/148‑1 (MQ=255)
                                                                  tGCTCAACTTCCTGCGCAACCGTACCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAGATCAACGTGCACATCCTCGAGAACGTGCTgttcggt                                                                                  >  6:70910/1‑149 (MQ=255)
                                                                    cTCAACTTCCTGCGCAACCGTACCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAGATCAACGTGCACATCCTCGAGAACGTGCTgttc                                                                                     >  7:153854/1‑144 (MQ=255)
                                                                    cTCAACTTCCTGCGCAACCGTACCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAGATCAACGTGCACATCCTCGAGAACGTGCTgttc                                                                                     <  8:153854/144‑1 (MQ=255)
                                                                     tCAACTTCCTGCGCAACCGTACCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAGATCAACGTGCACATCCTCGAGAACGTGCTgttcggttcg                                                                               >  2:339554/1‑149 (MQ=255)
                                                                       aaCTTCCTGCGCAACCGTACCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAGATCAACGTGCACATCCTCGAGAACGTGCTGttcggttcgg                                                                              <  6:337599/148‑1 (MQ=255)
                                                                        aCTTCCTGCGCAACCGTACCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAGATCAACGTGCACATCCTCGAGAACGTGCTGTTCGGTTCGGTg                                                                            <  1:64086/149‑1 (MQ=255)
                                                                               gcgcAACCGTACCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAGATCAACGTGCACATCCTCGAGAACGTGCTGTTCGGTTCGGTGCTGACCg                                                                     >  5:409990/1‑149 (MQ=255)
                                                                               gcgcAACCGTACCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAGATCAACGTGCACATCCTCGAGAACGTGCTGTTCGGTTCGGTGCTGACCg                                                                     >  5:44210/1‑149 (MQ=255)
                                                                                cgcAACCGTACCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAGATCAACGTGCACATCCTCGAGAACGTGCTGTTCGGTTCGGTGCTGACCGt                                                                    <  6:1345/149‑1 (MQ=255)
                                                                                  cAACCGTACCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAGATCAACGTGCACATCCTCGAGAACGTGCTGTTCGGTTCGGTGCTGACCGTTa                                                                  <  6:44210/149‑1 (MQ=255)
                                                                                   aaCCGTACCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTggcag                                                                                                                           >  4:143876/1‑91 (MQ=255)
                                                                                   aaCCGTACCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTggcag                                                                                                                           <  3:143876/91‑1 (MQ=255)
                                                                                    aCCGTACCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAGATCAACGTGCACATCCTCGAGAACGTGCTGTTCGGTTCGGTGCTGACCGTTAGt                                                                <  3:240766/149‑1 (MQ=255)
                                                                                    aCCGTACCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAGATCAACGTGCACATCCTCGAGAACGTGCTGTTCGGTTCGGTGCTGACCGTTAGt                                                                <  1:136921/149‑1 (MQ=255)
                                                                                     ccGTACCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCCAAGCGCCa                                                                                                                                                >  7:38718/1‑68 (MQ=255)
                                                                                     ccGTACCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCCAAGCGCCa                                                                                                                                                <  8:38718/68‑1 (MQ=255)
                                                                                        tACCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAGATCAACGTGCACATCCTCGAGAACGTGCTGTTCGGTTCGGTGCTGACCGTTAGTGccc                                                            <  7:218984/149‑1 (MQ=255)
                                                                                        tACCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAGATCAACGTGCACATCCTCGAGAACGTGCTGTTCGGTTCGGTGCTGACCGTTAGTGccc                                                            <  5:144251/149‑1 (MQ=255)
                                                                                          ccGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAGATCAACGTGCACATCCTCGAGAACGTGCTGTTCGGTTCGGTGCTGACCGTTAGTGCCCa                                                           <  4:123005/148‑1 (MQ=255)
                                                                                           cGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAGATCAACGTGCACATCCTCGAGAACGTGCTGTTCGGTTCGGTGCTGACCGTTAGTGCCCAg                                                          >  4:365713/1‑148 (MQ=255)
                                                                                            gggCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAGATCAACGTGCACATCCTCGAGAACGTGCTGTTCGGTTCGGTGCTGACCGTTAGTGCCCAGGa                                                        <  5:70910/149‑1 (MQ=255)
                                                                                                  tCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGGTGATGCTGGCAGGCAAGATCAACGTGCACATCCTCGAGAACGTGCTGTTCGGTTCGGTGCTGACCGTTAGTGCCCAGGACTTGc                                                   <  4:348979/148‑1 (MQ=255)
                                                                                                    gCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAGATCAACGTGCACATCCTCGAGAACGTGCTGTTCGGTTCGGTGCTGACCGTTAGTGCCCAGGACTTGCTgg                                                >  6:117608/1‑149 (MQ=255)
                                                                                                     ccGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAGATCAACGTGCACATCCTCGAGAACGTGCTGTTCGGTTCGGTGCTGACCGTTAGTGCCCAg                                                          <  7:267587/138‑1 (MQ=255)
                                                                                                     ccGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAGATCAACGTGCACATCCTCGAGAACGTGCTGTTCGGTTCGGTGCTGACCGTTAGTGCCCAg                                                          >  8:267587/1‑138 (MQ=255)
                                                                                                          tACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAGATCAACGTGCACATCCTCGAGAACGTGCTGTTCGGTTCGGTGCTGACCGTTAGTGCCCAGGACTTGCTGGTACTgg                                          <  5:300388/149‑1 (MQ=255)
                                                                                                             acTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAGATCAACGTGCATATCCTCGAGAACGTGCTGTTCGGTTCGGTGCTGACCGTTAGTGCCCAGGACTTGCTGGTACTGGGCa                                       >  2:70718/1‑149 (MQ=255)
                                                                                                             acTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAGATCAACGTGCACATCCTCGAGAACGTGCTgttc                                                                                     >  1:251457/1‑103 (MQ=255)
                                                                                                             acTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAGATCAACGTGCACATCCTCGAGAACGTGCTgttc                                                                                     <  2:251457/103‑1 (MQ=255)
                                                                                                             acTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAGATCAACGTGCACATCCTCGAGAACGTGCTGTTCGGTTCGGTGCTGACCGTTAGTGCCCAGTACCAGCAAGTCCTGGGCa                                       <  6:54984/149‑1 (MQ=255)
                                                                                                              cTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAGATCAACGTGCACATCCTCGAGAACGTGCTGTTCGGTTCGGTGCTGACCGTTAGTGCCCAGTACCAGCAAGTCCTGGGCac                                      <  1:192735/149‑2 (MQ=255)
                                                                                                                      cGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAGATCAACGTGCACATCCTCGAGAACGTGCTGTTCGGTTCGGTGCTGACCGTTAGTGCCCa                                                           <  5:392336/120‑1 (MQ=255)
                                                                                                                      cGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAGATCAACGTGCACATCCTCGAGAACGTGCTGTTCGGTTCGGTGCTGACCGTTAGTGCCCa                                                           >  6:392336/1‑120 (MQ=255)
                                                                                                                       gtgtTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAGATCAACGTGCATATCCTCGAGAACGTGCTGTTCGGTTCGGTGCTGACCGTTAGTGCCCAGGACTTGCTGGTACTGGGCATCGTCGCggt                             <  1:70718/149‑1 (MQ=255)
                                                                                                                       gtgtTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAGATCAACGTGCACATCCTCGAGAACGTGCTGTTCGGTTCGGTg                                                                            <  7:312604/102‑1 (MQ=255)
                                                                                                                       gtgtTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAGATCAACGTGCACATCCTCGAGAACGTGCTGTTCGGTTCGGTg                                                                            >  8:312604/1‑102 (MQ=255)
                                                                                                                       gtgtTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAGATCAACGTGCACATCCTCGAGAACGTGCTGTTCGGTTCGGTGCTGACCGTTAGTGCCCAGGACTTGCTGGTACTgg                                          >  1:260280/1‑136 (MQ=255)
                                                                                                                       gtgtTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAGATCAACGTGCACATCCTCGAGAACGTGCTGTTCGGTTCGGTGCTGACCGTTAGTGCCCAGGACTTGCTGGTACTgg                                          <  2:260280/136‑1 (MQ=255)
                                                                                                                            ccTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAGATCAACGTGCACATCCTCGAGAACGTGCTGTTCGGTTCGGTGCTGACCGTTAGTGCCCAGGACTTGCTGGTACTGGGCATCGTCGCggtgctg                         <  5:316941/148‑1 (MQ=255)
                                                                                                                            ccTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAGATCAACGTGCACATCCTCGAGAACGTGCTGTTCGGTTCGGTGCTGACCGTTAGTGCCCAGGACTTGCTGGTACTGGGCATCGTCGCggtgctg                         <  5:117608/148‑1 (MQ=255)
                                                                                                                             cTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAGATCAACGTGCACATCCTCGAGAACGTGCTGTTCGGTTCGGTGCTGACCGTTAGTGCCCAGGACTTGCTGGTACTGGGCAtcgtc                                  <  7:198108/138‑1 (MQ=255)
                                                                                                                             cTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAGATCAACGTGCACATCCTCGAGAACGTGCTGTTCGGTTCGGTGCTGACCGTTAGTGCCCAGGACTTGCTGGTACTGGGCAtcgtc                                  >  8:198108/1‑138 (MQ=255)
                                                                                                                             cTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAGATCAACGTGCACATCCTCGAGAACGTGCTGTTCGGTTCGGTGCTGACCGTTAGTGCCCAGGACTTGCTGGTACTGGGCATCGTCGCggtgctggt                       >  4:223425/1‑149 (MQ=255)
                                                                                                                                tCGGTGTCGCTGGCGCTTGGCGCCAGCCTGTTGCTGATGCTGGCAGGCAAGATCAACGTGCACATCCTCGAGAACGTGCTGTTCGGTTCGGTGCTGACCGTTAGTGCCCAGGACTTGCTGGTACTGGGCATCGTCGCggtgctggtgct                    >  1:148119/1‑149 (MQ=255)
                                                                                                                                tCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAGATCAACGTGCACATCCTCGAGAACGTGCTGTTCGGTTCGGTGCTGACCGTTAGTGCCCAGGACTTGCTGGTACTGGGCATCGTCGCggtgctggtg                      <  2:143170/147‑1 (MQ=255)
                                                                                                                                  ggTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAGATCAACGTGCACATCCTCGAGAACGTGCTGTTCGGTTCGTTGCTGACCGTTAGTGCCCAGGACTTGCTGGTACTGGGCATCGTCGCGGTGCTGGTGCtt                   <  5:382460/148‑1 (MQ=255)
                                                                                                                                  ggTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAGATCAACGTGCACATCCTCGAGAACGTGCTGTTCGGTTCGGTGCTGACCGTTAGTGCCCAGGACTTGCTGGTACTGGGCATCGTCGCGGTGCTGGTGCtt                   <  6:409990/148‑1 (MQ=255)
                                                                                                                                   gtgtCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAGATCAACGTGCACATCCTCGAGAACGTGCTGTTCGGTTCGGTGCTGACCGTTAGTGccc                                                            <  6:381009/106‑1 (MQ=255)
                                                                                                                                   gtgtCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAGATCAACGTGCACATCCTCGAGAACGTGCTGTTCGGTTCGGTGCTGACCGTTAGTGccc                                                            >  5:381009/1‑106 (MQ=255)
                                                                                                                                     gtCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAGATCAACGTGCACATCCTCGAGAACGTGCTGTTCGGTTCGGTGCTGACCGTTAGTGCCCAGGACTTGCTGGTACTGGGCATCGTCGCGGTGCTGGTGCTTGc                 >  5:238732/1‑147 (MQ=255)
                                                                                                                                     gtCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAGATCAACGTGCACATCCTCGAGAACGTGCTGTTCGGTTCGGTGCTGACCGTTAGTGCCCAGGACTTGCTGGTACTGGGCATCGTCGCGGTGCTGGTGCTTGc                 <  6:238732/147‑1 (MQ=255)
                                                                                                                                        gCTGGCGCCAAGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAGATCAACGTGCACATCCTCGAGAACGTGCTGTTCGGTTCGGTGCTGACCGTTAGTGCCCAGGACTTGCTGGTACTGGGCATCGTCGCggtgc                           <  3:275226/134‑1 (MQ=255)
                                                                                                                                        gCTGGCGCCAAGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAGATCAACGTGCACATCCTCGAGAACGTGCTGTTCGGTTCGGTGCTGACCGTTAGTGCCCAGGACTTGCTGGTACTGGGCATCGTCGCggtgc                           >  4:275226/1‑134 (MQ=255)
                                                                                                                                        gCTGGCGCCAAGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAGATCAACGTGCACATCCTCGAGAACGTGCTGTTCGGTTCGGTGCTGACCGTTAGTGCCCAGGACTTGCTGGTACTGGGCATCGTCGCGGTGCTGGTGCTTGCCTTgg            >  2:75425/1‑149 (MQ=255)
                                                                                                                                         cTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAGATCAACGTGCACATCCTCGAGAACGTGCTGTTCGGTTCGGTGCTGACCGTTAGTGCCCAGGACTTGCTGGTACTGGGCATCGTCGCGGTGCTGGTGCTTGCCTTGgc           <  1:339554/149‑1 (MQ=255)
                                                                                                                                              gcTTGGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAGATCAACGTGCACATCCTCGAGAACGTGCTGTTCGTTTCGGTGCTGACCGTTAGTGCCCAGGACTTGCTGGTACTGGGCATCGTCGCGGTGCTGGTGCTTGCCTTGGCGCTGc      >  8:180247/1‑149 (MQ=255)
                                                                                                                                              gcTTGGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAGATCAACGTGCACATCCTCGAGAACGTGCTGTTCGGTTCGGTGCTGACCGTTAGTGCCCAGGACTTGCTGGTACTGGGCATCGTCGCGGTGCTGGTGCTTGCCTTGGCGCTGc      >  1:358168/1‑149 (MQ=255)
                                                                                                                                                ttGGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAGATCAACGTGCACATc                                                                                                       <  5:143390/50‑1 (MQ=255)
                                                                                                                                                ttGGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAGATCAACGTGCACATc                                                                                                       >  6:143390/1‑50 (MQ=255)
                                                                                                                                                 tGGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAGATCAACGTGCACATCCTCGAGAACGTGCTGTTCGGTTCGGTGCTGACCGTTAGTGCCCAGGACTTGCTGGTACTGGGCAtcgtc                                  >  7:342348/1‑118 (MQ=255)
                                                                                                                                                 tGGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAGATCAACGTGCACATCCTCGAGAACGTGCTGTTCGGTTCGGTGCTGACCGTTAGTGCCCAGGACTTGCTGGTACTGGGCAtcgtc                                  <  8:342348/118‑1 (MQ=255)
                                                                                                                                                  ggCGCCAGCCTGTTGCTGATGCTGGCAGGCAAGATCAACGTGCACATCCTCGAGAACGTGCTGTTCGGTTCGGTGCTGACCGTTAGTGCCCAGGACTTGCTGGTACTGGGCATCGTCGCGGTGCTGGTGCTTGCCTTGGCGCTGCCg    <  2:148119/147‑1 (MQ=255)
                                                                                                                                                  ggCGCCAGCCTGCTGCTGATGCTGGCAGGCAAGATCAACGTGCACATCCTCGAGAACGTGCTGTTCGGTTCGGTGCTGACCGTTAGTGCCCAGGACTTGCTGGTACTGGGCATCGTCGCGGTGCTGGTGCTTGCCTTGGCGCTGCCGCt  >  2:397314/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
TGCTGGGTGAGCCCTATACCGGCCCTTACGGCAGCCTGTTCGGCTACTGCCTGCTGTTCGGCATCCTGCTCAACTTCCTGCGCAACCGTACCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAGATCAACGTGCACATCCTCGAGAACGTGCTGTTCGGTTCGGTGCTGACCGTTAGTGCCCAGGACTTGCTGGTACTGGGCATCGTCGCGGTGCTGGTGCTTGCCTTGGCGCTGCCGCT  >  NC_002947/4331521‑4331815

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: