Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,331,665 T→C 5.5% L122P (CTT→CCT) ‡ 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,6650TC 5.5% 271.7 / 4.5 91L122P (CTT→CCT) ‡PP_3803cation ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base T (39/47);  new base C (3/2);  total (42/49)
Fisher's exact test for biased strand distribution p-value = 6.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.50e-01

TGCTGGGTGAGCCCTATACCGGCCCTTACGGCAGCCTGTTCGGCTACTGCCTGCTGTTCGGCATCCTGCTCAACTTCCTGCGCAACCGTACCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAGATCAACGTGCACATCCTCGAGAACGTGCTGTTCGGTTCGGTGCTGACCGTTAGTGCCCAGGACTTGCTGGTACTGGGCATCGTCGCGGTGCTGGTGCTTGCCTTGGCGCTGC  >  NC_002947/4331521‑4331811
                                                                                                                                                |                                                                                                                                                  
tgctgGGTGAGCCCTATACCGGCCCTTACGGCAGCCTGTTCGGCTACTGCCTGCTGTTCGGCATCCTGCTCAACTTCCTGCGCAACCGTACCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCCGTGTCGCTGGCGCTTGgc                                                                                                                                                >  3:348979/1‑149 (MQ=255)
      gTGAGCCCCATACCGGCCCTTACGGCAGCCTGTTCGGCTACTGCCTGCTGTTCGGCATCCTGCTCAACTTCCTGCGCAACCGTACCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGc                                                                                                                                          <  5:351969/149‑1 (MQ=255)
           cccTATACCGGCCCTTACGGCAGCCTGTTCGGCTACTGCCTGCTGTTCGGCATCCTGCTCAACTTCCTGCGCAACCGTACCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCctgct                                                                                                                                     >  5:81993/1‑149 (MQ=255)
           cccTATACCGGCCCTTACGGCAGCCTGTTCGGCTACTGCCTGCTGTTCGGCATCCTGCTCAACTTCCTGCGCAACCGTACCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCctgct                                                                                                                                     >  6:76499/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                                                                                                     >  2:174282/1‑133 (MQ=255)
                                                           ggCATCCTGCTCAACTTCCTGCGCAACCGTACCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAGATCAACGTGcaca                                                                                                     <  1:174282/133‑1 (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                                                                                <  7:289988/148‑1 (MQ=255)
                                                                 cTGCTCAACTTCCTGCGCAACCGTACCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAGATCAACGTGCACATCCTCGAGAACGTGCTgttcg                                                                                <  1:227662/148‑1 (MQ=255)
                                                                 cTGCTCAACTTCCTGCGCAACCGTACCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTAGGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAGATCAACGTGCACATCCTCGAGAACGTGCTgttcg                                                                                <  7:151468/148‑1 (MQ=255)
                                                                  tGCTCAACTTCCTGCGCAACCGTACCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAGATCAACGTGCACATCCTCGAGAACGTGCTgttcggt                                                                              >  6:70910/1‑149 (MQ=255)
                                                                    cTCAACTTCCTGCGCAACCGTACCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAGATCAACGTGCACATCCTCGAGAACGTGCTgttc                                                                                 <  8:153854/144‑1 (MQ=255)
                                                                    cTCAACTTCCTGCGCAACCGTACCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAGATCAACGTGCACATCCTCGAGAACGTGCTgttc                                                                                 >  7:153854/1‑144 (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:117608/148‑1 (MQ=255)
                                                                                                                            ccTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAGATCAACGTGCACATCCTCGAGAACGTGCTGTTCGGTTCGGTGCTGACCGTTAGTGCCCAGGACTTGCTGGTACTGGGCATCGTCGCggtgctg                     <  5:316941/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             <  6:238732/147‑1 (MQ=255)
                                                                                                                                     gtCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAGATCAACGTGCACATCCTCGAGAACGTGCTGTTCGGTTCGGTGCTGACCGTTAGTGCCCAGGACTTGCTGGTACTGGGCATCGTCGCGGTGCTGGTGCTTGc             >  5:238732/1‑147 (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                                                                                                   >  6:143390/1‑50 (MQ=255)
                                                                                                                                                ttGGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAGATCAACGTGCACATc                                                                                                   <  5:143390/50‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                  
TGCTGGGTGAGCCCTATACCGGCCCTTACGGCAGCCTGTTCGGCTACTGCCTGCTGTTCGGCATCCTGCTCAACTTCCTGCGCAACCGTACCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAGATCAACGTGCACATCCTCGAGAACGTGCTGTTCGGTTCGGTGCTGACCGTTAGTGCCCAGGACTTGCTGGTACTGGGCATCGTCGCGGTGCTGGTGCTTGCCTTGGCGCTGC  >  NC_002947/4331521‑4331811

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: