Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,987,993 G→C 6.5% V153L (GTC→CTC) ‡ PP_1779 → lipopolysaccharide ABC transporter ATP‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,987,9930GC 6.5% 240.8 / 6.3 77V153L (GTC→CTC) ‡PP_1779lipopolysaccharide ABC transporter ATP‑binding protein
Reads supporting (aligned to +/- strand):  ref base G (35/37);  new base C (3/2);  total (38/39)
Fisher's exact test for biased strand distribution p-value = 6.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.58e-01

CGCCAGAACGTGTACCACGCTGCGGGTCTGATGGGTTTCGATGTACAGCAGATCAATGATGTCATTGATGACATTGAAGCCTTTGCCGAGATTGGCGAGTACTTTGATGAGGCAGTGCGTACCTATTCCAGCGGTATGCAAATGCGCGTCGCATTTGCCGTGGCGACTGCGTTCCGTCCAGAAATCCTTATCGTCGATGAGGCGCTGTCAGTAGGGGACAGCTATTTCCAGCACAAGAGTTTTGAGCGCATCCGTGAGTTCCAGCGGCTGGGAACATCGCTTTTGATAGTGT  >  NC_002947/1987846‑1988137
                                                                                                                                                   |                                                                                                                                                
cgccTGAACGTGTACCACGCTGCGGGTCTGATGGGTTTCGATGTACAGCAGATCAATGATGTCATTGATGACATTGAAGCCTTTGCCGAGATTGGCGAGTTCTTTGATGAGTCAGTGCGTACCTATTCCAGCGTTATGCAAATGcgcg                                                                                                                                                  <  3:401255/148‑1 (MQ=255)
cgccAGAACGTGTACAACGCTGCGGGTCTGATGGGTTTCGATGTACCGCAGATCAATGATGTCATTGATGACATTGAAGCCTTTGCCGAGATTGGCGAGTACTTTGATGAGGCAGTGCGTACCTATTCCAGCGGTATGCAAATGCGCGt                                                                                                                                                 >  2:263935/1‑149 (MQ=255)
      aaCGTGTACCACGCTGCGGGTCTGATGGGTTTCGATGTACAGCAGATCAATGATGTCATTGATGACATTGAAGCCTTTGCCGAGATTGGCGAGTACTTTGATGAGGCAGTGCGTACCTATTCCAGCGGTATGCAAATGCGCGTCGCAtt                                                                                                                                           <  8:252158/149‑1 (MQ=255)
          tgtACCACGCTGCGGGTCTGATGGGTTTCGATGTACAGCAGATCAATGATGTCATTGATGACATTGAAGCCTTTGCCGAGATTGGCGAGTACTTTGATGAGGCAGTGCGTACCTATTCCAGCGGTATGCAAATGCGACGCGCATTTGca                                                                                                                                       <  3:240512/149‑2 (MQ=255)
            tACCACGCTGCGGGTCTGATGGGTTTCGATGTACAGCAGATCAATGATGTCATTGATGACATTGAAGCCTTTGCCGAGATTGGCGAGTACTTTGATGAGGCAGTGCGTACCTATTCCAGCGGTATGCAAATGCGCGTCGCATTTGCCGt                                                                                                                                     >  5:327816/1‑149 (MQ=255)
                    tGCGGGTCTGATGGGTTTCGATGTACAGCAGATCAATGATGTCATTGATGACATTGAAGCCTTTGCCGAGATTGGCGAGTACTTTGATGAGGCAGTGCGTACCTATTCCAGCGGTATGCAAATGCGCGTCGCATTTGCCGTGGCGACTg                                                                                                                             <  4:56948/149‑1 (MQ=255)
                               tGGGTTTCGATGTACAGCAGATCAATGATGTCATTGATGACATTGAAGCCTTTGCCGAGATTGGCGAGTACTTTGATGAGGCAGTGCGTACCTATTCCAGCGGTATGCAAATGCGCGTCGCATTTGCCGTGGCGACTGCTTTCCGTcc                                                                                                                   <  6:38738/148‑1 (MQ=255)
                                gggTTTCGATGTACAGCAGATCAATGATGTCATTGATGACATTGAAGCCTTTGCCGAGATTGGCGAGTACTTTGATGAGGCAGTGCGTACCTATTCCAGCGGTATGCAAATGCGCGTCGCATTTGCCGTGGCGACTGCGTTCCGTc                                                                                                                    >  7:126161/1‑146 (MQ=255)
                                  gTTTCGATGTACAGCAGATCAATGATGTCATTGATGACATTGAAGCCTTTGCCGAGATTGGCGAGTACTTTGATGAGGCAGTGCGTACCTATTCCAGCGGTATGCAAATGCGCGTCGCATTTGCCGTGGCGACTGCGTTCCGTCCAg                                                                                                                 <  8:116924/147‑1 (MQ=255)
                                  gTTTCGATGTACAGCAGATCAATGATGTCATTGATGACATTGAAGCCTTTGCCGAGATTGGCGAGTACTTTGATGAGGCAGTGCGTACCTATTCCAGCGGTATGCAAATGCGCGTCGCATTTGCCGTGGCGACTGCGTTCCGTCCAg                                                                                                                 >  7:116924/1‑147 (MQ=255)
                                   tttCGATGTACAGCAGATCAATGATGTCATTGATGACATTGAAGCCTTTGCCGAGATTGGCGAGTACTTTGATGAGGCAGTGCGTACCTATTCCAGCGGTATGCAAATGCGCGTCGCATTTGCCGTGGCGACTGCGTTCCGTCCAGaaa                                                                                                              >  1:29988/1‑149 (MQ=255)
                                   tttCGATGTACAGCAGATCAATGATGTCATTGATGACATTGAAGCCTTTGCCGAGATTGGCGAGTACTTTGATGAGGCAGTGCGTACCTATTCCAGCGGTATGCAAATGCGCGTCGCATTTGCCGTGGCGACTGCGTTCCGTCCAGaa                                                                                                               >  6:58003/1‑148 (MQ=255)
                                    ttCGATGTACAGCAGATCAATGATGTCATTGATGACATTGAAGCCTTTGCCGAGATTGGCGAGTACTTTGATGAGGCAGTGCGTACCTATTCCAGCGGTATGCAAATGCGCGTCGCATTTGCCGTGGCGACTGCGTTCCGTCCAGAAAt                                                                                                             >  5:29491/1‑149 (MQ=255)
                                       gATGTACAGCAGATCAATGATGTCATTGATGACATTGAAGCCTTTGCCGAGATTGGCGAGTACTTTGATGAGGCAGTGCGTACCTATTCCAGCGGTATGCAAATGCGCGTCGCATTTGCCGTGGCGACTGCGTTCCGTCCAGAAATcc                                                                                                           <  5:58003/148‑1 (MQ=255)
                                       gATGTACAGCAGATCAATGATGTCATTGATGACATTGAAGCCTTTGCCGAGATTGGCGAGTACTTTGATGAGGCAGTGCGTACCTATTCCAGCGGTATGCAAATGCGCGTCGCATTTGCCGTGGCGACTGCGTTCCGTCCAGAAATcc                                                                                                           <  8:204190/148‑1 (MQ=255)
                                       gATGTACAGCAGATCAATGATGTCATTGATGACATTGAAGCCTTTGCCGAGATTGGCGAGTACTTTGATGAGGCAGTGCGTACCTATTCCAGCGGTATGCAAATGCGCGTCGCATTTGCCGTGGCGACTGCGTTCCGTCCAGAAATCCt                                                                                                          >  3:345755/1‑149 (MQ=255)
                                         tGTACAGCAGATCAATGATGTCATTGATGACATTGAAGCCTTTGCCGAGATTGGCGAGTACTTTGATGAGGCAGTGCGTACCTATTCCAGCGGTATGCAAATGCGCGTCGCATTTGCCGTGGCGACTGCGTTCCGTCCAGAAATCCTTa                                                                                                        <  1:53692/149‑1 (MQ=255)
                                         tGTACAGCAGATCAATGATGTCATTGATGACATTGAAGCCTTTGCCGAGATTGGCGAGTACTTTGATGAGGCAGTGCGTACCTATTCCAGCGGTATGCAAATGCGCGTCGCATTTGCCGTGGCGACTGCGTTCCGTCCAGAAATCCTTa                                                                                                        <  2:62098/149‑1 (MQ=255)
                                            aCAGCAGATCAATGATGTCATTGATGACATTGAAGCCTTTGCCGAGATTGGCGAGTACTTTGATGAGGCAGTGCGTACCTATTCCAGCGGTATGCAAATGCGCGTCGCATTTGCCGCGGCGACTGCGTTCCGTCCAGAAATCCTTAtcg                                                                                                     >  6:119820/1‑149 (MQ=255)
                                               gcagATCAATGATGTCATTGATGACATTGAAGCCTTTGCCGAGATTGGCGAGTACTTTGATGAGGCAGTGCGTACCTATTCCAGCGGTATGCAAATGCGCGTCGCATTTGCCGTGGCGACTGCGTTCCGTCCAGAAATCCTTAtcgtc                                                                                                   >  8:278505/1‑148 (MQ=255)
                                                cagATCAATGATGTCATTGATGACATTGAAGCCTTTGCCGAGATTGGCGAGTACTTTGATGAGGCAGTGCGTACCTATTCCAGCGGTATGCAAATGCGCGTCGCATTTGCCGTGGCGACTGCGTTCCGTCCAGAAATCCTTATCGTCGa                                                                                                 >  1:396899/1‑149 (MQ=255)
                                                   aTCAATGATGTCATTGATGACATTGAAGCCTTTGCCGAGATTGGCGAGTACTTTGATGAGGCAGTGCGTACCTATTCCAGCGGTATGCAAATGCGCGTCGCATTTGCCGTGGCGACTGCGTTCCGTCCAGAAATCCTTATCGTCGAt                                                                                                >  7:338395/1‑147 (MQ=255)
                                                   aTCAATGATGTCATTGATGACATTGAAGCCTTTGCCGAGATTGGCGAGTACTTTGATGAGGCAGTGCGTACCTATTCCAGCGGTATGCAAATGCGCGTCGCATTTGCCGTGGCGACTGCGTTCCGTCCAGAAATCCTTATCGTCGATGa                                                                                              <  6:81865/149‑1 (MQ=255)
                                                           tgTCATTGATGACATTGAAGCCTTTGCCGAGATTGGCGAGTACTTTGATGAGGCAGTGCGTACCTATTCCAGCGGTATGCAAATGCGCGTCGCATTTGCCGTGGCGACTGCGTTCCGTCCAGAAATCCTTATCGTCGATGAGGCGCTGt                                                                                      <  6:29491/149‑1 (MQ=255)
                                                           tgTCATTGATGACATTGAAGCCTTTGCCGAGATTGGCGAGTACTTTGATGAGGCAGTGCGTACCTATTCCAGCGGTATGCAAAAGCGCGTCGCATTTGCCGTGGCGACTGCGTTCCGTCCAGAAATCCTTATCGTCGATGAGGCGCTGt                                                                                      <  3:9825/149‑1 (MQ=255)
                                                                    tgaCATTGAAGCCTTTGCCGAGATTGGCGAGTACTTTGATGAGGCAGTGCGTACCTATTCCAGCGGTATGCAAATGCGCGTCGCATTTGCCGTGGCGACTGCGTTCCGTCCAGAAATCCTTATCGTCGATGAGGCGCTGTCAGTAggg                                                                              >  2:21284/1‑148 (MQ=255)
                                                                    tgaCATTGAAGCCTTTGCCGAGATTGGCGAGTACTTTGATGAGGCAGTGCGTACCTATTCCAGCGGTATGCAAATGCGCGTCGCATTTGCCGTGGCGACTGCGTTCCGTCCAGAAATCCTTATCGTCGATGAGGCGCTGTCAGTAggg                                                                              <  4:117273/148‑1 (MQ=255)
                                                                    tgaCATTGAAGCCTTTGCCGAGATTGGCGAGTACTTTGATGAGGCAGTGCGTACCTATTCCAGCGGTATGCAAATGCGCGTCGCATTTGCCGTGGCGACTGCGTTCCGTCCAGAAATCCTTATCGTCGATGAGGCGCTGTCAGTAggg                                                                              <  4:3287/148‑1 (MQ=255)
                                                                    tgaCATTGAAGCCTTTGCCGAGATTGGCGAGTACTTTGATGAGGCAGTGCGTACCTATTCCAGCAGTATGCAAATGCGCGTCGCATTTGCCGTGGCGACTGCGTTCCGTCCAGAAATCCTTATCGTCGATGAGGCGCTGTCAGTAggg                                                                              >  1:349341/1‑148 (MQ=255)
                                                                     gaCATTGAAGCCTTTGCCGAGATTGGCGAGTACTTTGATGAGGCAGTGCGTACCTATTCCAGCGGTATGCAAATGCGCGTCGCATTTGCCGTGGCGACTGCGTTCCGTCCAGAAATCCTTATCGTCGATGAGGCGCTGTCAGTAGGGGa                                                                            >  4:158255/1‑149 (MQ=255)
                                                                      aCATTGAAGCCTTTGCCGAGATTGGCGAGTACTTTGATGAGGCAGTGCGTACCTATTCCAGCGGTATGCAAATGCGCGTCGCATTTGCCGTGGCGACTGCGTTCCGTCCAGAAATCCTTATCGTCGATGAGGCGCTGTCAGTAGGGGa                                                                            <  8:126161/148‑1 (MQ=255)
                                                                      aCATTGAAGCCTTTGCCGAGATTGGCGAGTACTTTGATGAGGCAGTGCGTACCTATTCCAGCGGTATGCAAATGCGCGTCGCATTTGCCGTGGCGACTGCGTTCCGTCCAGAAATCCTTATCGTCGATGAGGCGCTGTCAGTAGGGGAc                                                                           <  8:338395/149‑1 (MQ=255)
                                                                       cATTGAAGCCTTTGCCGAGATTGGCGAGTACTTTGATGAGGCAGTGCGTACCTATTCCAGCGGTATGCAAATGCGCGTCGCATTTGCCGTGGCGACTGCGTTCCGTCCAGAAATCCTTATCGTCGATGAGGCGCTGTCAGTAGGGGACa                                                                          <  6:327816/149‑1 (MQ=255)
                                                                       cATTGAAGCCTTTGCCGAGATTGGCGAGTACTTTGATGAGGCAGTGCGTACCTATTCCAGCGGTATGCAAATGCGCGTCGCATTTGCCGTGGCGACTGCGTTCCGTCCAGAAATCCTTATCGTCGATGAGGCACTGTCAGTAGGGGa                                                                            >  5:330724/1‑147 (MQ=255)
                                                                        aTTGAAGCCTTTGCCGAGATTGGCGAGTACTTTGATGAGGCAGTGCGTACCTATTCCAGCGGTATGCAAATGCGACGCGCATTTGcatac                                                                                                                                    <  3:3577/90‑5 (MQ=38)
                                                                        aTTGAAGCCTTTGCCGAGATTGGCGAGTACTTTGATGAGGCAGTGCGTACCTATTCCAGCGGTATGCAAATGCGACGCGCATTTGcatac                                                                                                                                    >  4:3577/1‑86 (MQ=38)
                                                                         ttGAAGCCTTTGCCGAGATTGGCGAGTACTTTGATGAGGCAGTGCGTACCTATTCCAGCGGTATGCAAATGCGCGTCGCATTTGCCGTGGCGACTGCGTTCCGTCCAGAAATCCTTATCGTCGATGAGGCGCTGTCAGTAGGGGACAGc                                                                        <  2:29988/149‑1 (MQ=255)
                                                                            aaGCCTTTGCCGAGATTGGCGAGTACTTTGATGAGGCAGTGCGTACCTATTCCAGCGGTATGCAAATGCGCGTCGCATTTGCCGTGGCGACTGCGTTCCGTCCAGAAATCCTTATCGTCGATGAGGCGCTGTCAGTAGGGGACAGCTAt                                                                     >  1:154532/1‑149 (MQ=255)
                                                                            aaGCCTTTGCCGAGATTGGCGAGTACTTTGATGAGGCAGTGCGTACCTATTCCAGCGGTATGCAAATGCGCGTCGCATTTGCCGTGGCGACTGCGTTCCGTCCAGAAATCCTTATCGTCGATGAGGCGCTGTCAGTAGGGGACAGCTAt                                                                     >  2:71657/1‑149 (MQ=255)
                                                                             aGCCTTTGCCGAGATTGGCGAGTACTTTGATGAGGCAGTGCGTACCTATTCCAGCGGTATGCAAATGCGCGTCGCATTTGCCGTGGCGACTGCGTTCCGTCCAGAAATCCTTATCGTCGATGAGGCGCTGTCAGTAGGGGACAGCTAtt                                                                    <  1:437212/149‑1 (MQ=255)
                                                                              gCCTTTGCCGAGATTGGCGAGTACTTTGATGAGGCAGTGCGTACCTATTCCAGCGGTATGCAAATGCGCGTCGCATTTGCCGTGGCGACTGCGTTCCGTCCAGAAATCCTTATCGTCGATGAGGCGCTGTCAGTAGGGGACAGCTAttt                                                                   >  1:115838/1‑149 (MQ=255)
                                                                               ccTTTGCCGAGATTGGCGAGTACTTTGATGAGGCAGTGCGTACCTATTCCAGCGGTATGCAAATGCGCGTCGCATTTGCCGTGGCGACTGCTTTCCGTCCAGAAATCCTTATCGTCGATGAGGCGCTGTCAGTAGGGGACAGCTATTTc                                                                  <  5:285814/149‑1 (MQ=255)
                                                                               ccTTTGCCGAGATTGGCGAGTACTTTGATGAGGCAGTGCGTACCTATTCCAGCGGTATGCAAATGCGCGTCGCATTTGCCGTGGCGACTGCGTTCCGTCCAGAAATCCTTATCGTCGATGAGGCGCTGTCAGTAGGGGACAGCTATTTc                                                                  <  8:251261/149‑1 (MQ=255)
                                                                                cTTTGCCGAGATTGGCGAGTACTTTGATGAGGCAGTGCGTACCTATTCCAGCGGTATGCAAATGCGCGTCGCATTTGCCGTGGCGACTGCGTTCCGTCCAGAAATCCTTATCGTCGATGAGGCGCTGTCAGTAGGGGACAGCTATTTc                                                                  >  4:21145/1‑148 (MQ=255)
                                                                                    gCCGAGATTGGCGAGTACTTTGATGAGGCAGTGCGTACCTATTCCAGCGGTATGCAAATGCGCGTCGCATTTGCCGTGGCGACTGCGTTCCGTCCAGAAATCCTTATCGTCGATGAGGCGCTGTCAGTAGGGGACAGCTATTTCCAGc                                                              <  2:396899/148‑1 (MQ=255)
                                                                                      cGAGATTGGCGAGTACTTTGATGAGGCAGTGCGTACCTATTCCAGCGGTATGCAAATGCGCGTCGCATTTGCCGTGGCGACTGCGTTCCGTCCAGAAATCCTTATCGTCGATGAGGCGCTGTCa                                                                                    <  1:167386/124‑1 (MQ=255)
                                                                                      cGAGATTGGCGAGTACTTTGATGAGGCAGTGCGTACCTATTCCAGCGGTATGCAAATGCGCGTCGCATTTGCCGTGGCGACTGCGTTCCGTCCAGAAATCCTTATCGTCGATGAGGCGCTGTCa                                                                                    >  2:167386/1‑124 (MQ=255)
                                                                                      cGAGATTGGCGAGTACTTTGATGAGGCAGTGCGTACCTATTCCAGCGGTATGCAAATGCGCGTCGCATTTGCCGTGGCGACTGCGTTCCGTCCAGAAATCCTTATCGTCGATGAGGCGCTGTCAGTAGGGGACAGCTATTTCCAGcac                                                            >  1:67374/1‑148 (MQ=255)
                                                                                            tGGCGAGTACTTTGATGAGGCAGTGCGTACCTATTCCAGCGGTATGCAAATGCGCGTCGCATTTGCCGTGGCGACTGCGTTCCGTCCAGAAATCCTTATCGTCGATGAGGCGCTGTCAGTAGGGGACAGCTATTTCCAGCACAAGAGt                                                      <  2:349341/148‑1 (MQ=255)
                                                                                              gCGAGTACTTTGATGAGGCAGTGCGTACCTATTCCAGCGGTATGCAAATGCGCGTCGCATTTGCCGTGGCGACTGCGTTCCGTCCAGAAATCCTTATCGTCGATGAGGCGCTGTCAGTAGGGGACAGCTATTTCCAGCACAAGAGtttt                                                   <  4:345755/149‑1 (MQ=255)
                                                                                                gAGTACTTTGATGAGGCAGTGCGTACCTATTCCAGCGGTATGCAAATGCGCGTCGCATTTGCCGTGGCGACTGCGTTCCGTCCAGAAATCCTTATCGTCGATGAGGCGCTGTCAGTAGGGGACAGCTATTTCCAGCACAAGAGtttt                                                   >  6:42255/1‑147 (MQ=255)
                                                                                                gAGTACTTTGATGAGGCAGTGCGTACCTATTCCAGCGGTATGCAAATGCGCGTCGCATTTGCCGTGGCGACTGCGTTCCGTCCAGAAATCCTTATCGTCGATGAGGCGCTGTCAGTAGGGGACAGCTATTTCCAGCACAAGAGTTTTGa                                                 <  5:83423/149‑1 (MQ=255)
                                                                                                gAGTACTTTGATGAGGCAGTGCGTACCTATTCCAGCGGTATGCAAATGCGCGTCGCATTTGCCGTGGCGACTGCGTTCCGTCCAGAAATCCTTATCGTCGATGAGGCGCTGTCAGTAGGGGACAGCTATTTCCAGCACAAGAGTTTTGa                                                 <  7:278505/149‑1 (MQ=255)
                                                                                                  gTACTTTGATGAGGCAGTGCGTACCTATTCCAGCGGTATGCAAATGCGCGTCGCATTTGCCGTGGCGACTGCGTTCCGTCCAGAAATCCTTATCGTCGATGAGGCGCTGTCAGTAGGGGACAGCTATTTCCAGCACAAGAGTTTTGAgc                                               <  5:42255/149‑1 (MQ=255)
                                                                                                   tACTTTGATGAGGCAGTGCGTACCTATTCCAGCGGTATGCAAATGCGCGTCGCATTTGCCGTGGCGACTGCGTTCCGTCCAGAAATCCTTATCGTCGATGAGGCGCTGTCAGTAGGGGACAGCTATTTCCAGCACAAGAGTTTTGAgcg                                              >  7:101746/1‑149 (MQ=255)
                                                                                                     cTTTGATGAGGCAGTGCGTACCTATTCCAGCGGTATGCAAATGCGCGTCGCATTTGCCGTGGCGACTGCGTTCCGTCCAGAAATCCTTATCGTCGATGAGGCGCTGTCAGTAGGGGACAGCTATTTCCAGCACAAGAGTTTTGAGCGCa                                            >  8:84707/1‑149 (MQ=255)
                                                                                                      tttGATGAGGCAGTGCGTACCTATTCCAGCGGTATGCAAATGCGCGTCGCATTTGCCGTGGCGACTGCGTTCCGTCCAGAAATCCTTATCGTCGATGAGGCGCTGTCAGTAGGGGACAGCTATTTCCAGCACAAGAGTTTTGAGCGCAt                                           >  8:293914/1‑149 (MQ=255)
                                                                                                        tgatgaGGCAGTGCGTACCTATTCCAGCGGTATGCAAATGCGCGTCGCATTTGCCGTGGCGACTGCGTTCCGTCCAGAAATCCTTATCGTCGATGAGGCGCTGTCAGTAGGGGACAGCTATTTCCAGCACAAGAGTTTTGAGCGCAt                                           <  2:115838/147‑1 (MQ=255)
                                                                                                          atgaGGCAGTGCGTACCTATTCCAGCGGTATGCAAATGCGCGTCGCATTTGCCGTGGCGACTGCGTTCCGTCCAGAAATCCTTATCGTCGATGAGGCGCTGTCAGTAGGGGACAGCTATTTCCAGCACAAGAGTTTTGAGCGCATCCGt                                       >  6:278439/1‑149 (MQ=255)
                                                                                                           tgaGGCAGTGCGTACCTATTCCAGCGGTATGCAAATGCGCGTCGCATTTGCCGTGGCGACTGCGTTCCGTCCAGAAATCCTTATCGTCGATGAGGCGCTGTCAGTAGGGGACAGCTATTTCCAGCACAAGAGTTTTGAGCGCATCCGt                                       <  7:84707/148‑1 (MQ=255)
                                                                                                            gaGGCAGTGCGTACCTATTCCAGCGGTATGCAAATGCGCGTCGCATTTGCCGTGGCGACTGCGTTCCGTCCAGAAATCCTTATCGTCGATGAGGCGCTGTCAGTAGGGGACAGCTATTTCCAGCACAAGAGTTTTGAg                                                <  6:166923/138‑1 (MQ=255)
                                                                                                            gaGGCAGTGCGTACCTATTCCAGCGGTATGCAAATGCGCGTCGCATTTGCCGTGGCGACTGCGTTCCGTCCAGAAATCCTTATCGTCGATGAGGCGCTGTCAGTAGGGGACAGCTATTTCCAGCACAAGAGTTTTGAg                                                >  5:166923/1‑138 (MQ=255)
                                                                                                                cAGTGCGTACCTATTCCAGCGGTATGCAAATGCGCGTCGCATTTGCCGTGGCGACTGCGTTCCGTCCAGAAATCCTTATCGTCGATGAGGCGCTGTCAGTAGGGGACAGCTATTTCCAGCACAAGAGTTTTGAGCGCATCCGTGAGTTc                                 <  2:155488/149‑1 (MQ=255)
                                                                                                                         ccTATTCCAGCGGTATGCAAATGCGCGTCGCATTTGCCGTGGCGACTGCGTTCCGTCCAGAAATCCTTATCGTCGATGAGGCGCTGTCAGTAGGGGACAGCTATTTCCAGCACAAGAGTTTTGAGCGCATCCGTGAGTTCCAGCGGCt                         <  1:21284/148‑1 (MQ=255)
                                                                                                                          cTATTCCAGCGGTATGCAAATGCGCGTCGCATTTGCCGTGGCGACTGCGTTCCGTCCAGAAATCGTTATCGTCGATGAGGCGCTGTCAGTAGGGGAGAGCTATTTCCAGCACAAGAGTTTTGAGCGCATCAGAGAGTTCAAGCGGCTgg                       >  8:101107/1‑149 (MQ=255)
                                                                                                                          cTATTCCAGCGGTATGCAAATGCGCGTCGCATTTGCCGTGGCGACTGCGTTCCGTCCAGAAATCCTTATCGTCGATGAGGCGCTGTCAGTAGGGGACAGCTATTTCCAGCACAAGAGTTTTGAGCGCATCCGTGAGTTCCAGCGGCTgg                       >  5:177110/1‑149 (MQ=255)
                                                                                                                          cTATTCCAGCGGTATGCAAATGCGCGTCGCATTTGCCGTGGCGACTGCGTTCCGTCCAGAAATCCTTATCGTCGATGAGGCGCTGTCAGTAGGGGACAGCTATTTCCAGCACAAGAGTTTTGAGCGCATCCGTGAGTTCCAGCGGCTgg                       >  8:236805/1‑149 (MQ=255)
                                                                                                                          cTATTCCAGCGGTATGCAAATGCGCGTCGCATTTGCCGTGGCGACTGCGTTCCGTCCAGAAATCCTTATCGTCGATGAGGCGCTGTCAGTAGGGGACAGCTATTTACAGCACAAGAGTTTTGAGCGCATCCGTGAGTTCCAGCGGCTgg                       >  3:362406/1‑149 (MQ=255)
                                                                                                                             ttCCAGCGGTATGCAAATGCGCGTCGCATTTGCCGTGGCGACTGCGTTCCGTCCAGAAATCCTTATCGTCGATGAGGCGCTGTCAGTAGGGGACAGCTATTTCCAGCACAAGAGTTTTGAGCGCATCCGTGAGTTCCAGCGGCTGGGa                     <  2:154532/148‑1 (MQ=255)
                                                                                                                               ccAGCGGTATGCAAATGCGCGTCGCATTTGCCGTGGCGACTGCGTTCCGTCCAGAAATCCTTATCGTCGATGAGGCGCTGTCAGTAGGGGACAGCTATTTCCAGCACAAGAGTTTTGAGCGCATCCGTGAGTTCCAGCGGCTGGGAAc                   <  5:278439/148‑1 (MQ=255)
                                                                                                                                 aGCGGTATGCAAATGCGCGTCGCATTTGCCGTGGCGACTGCGTTCCGTCCAGAAATCCTTATCGTCGATGAGGCGCTGTCAGTAGGGGACAGCTATTTCCAGCACAAGAGTTTTGAGCGCATCCGTGAGTTCCAGCGGCTGGGAACATc                <  7:236805/149‑1 (MQ=255)
                                                                                                                                        tGCAAATGCGCGTCGCATTTGCCGTGGCGACTGCGTTCCGTCCAGAAATCCTTATCGTCGATGAGGCGCTGTCAGTAGGGGACAGCTATTTCCAGCACAAGAGTTTTGAGCGCATCCGTGAGTTCCAGCGGCTGGGAACATCGCtttt          <  6:325500/148‑1 (MQ=255)
                                                                                                                                        ggCAAATGCGACGCGCATTTGCCGTGGCGACTGCGTTCCGTCCAGAAATCCTTATCGTCGATGAGGCGCTGTCAGTAGGGGACAGCTATTTCCAGCACAAGAGTTTTGAGCGCATCCGTGAGTTCCAGCGGCTGGGAACATCGCtttt          >  1:389238/2‑148 (MQ=255)
                                                                                                                                          cAAATGCGACGCGCATTTGCCGTGGCGACTGCGTTCCGTCCAGAAATCCTTATCGTCGATGAGGCGCTGTCAGTAGGGGACAGCTATTTCCAGCACAAGAGTTTTGAGCGCATCCGTGAGTTCCAGCGGCTGGGAACATCGCTTTTGa        >  7:40396/1‑148 (MQ=255)
                                                                                                                                           aaaTGCGCGTCGCATTTGCCGTGGCGACTGCGTTCCGTCCAGAAATCCTTATCGTCGATGAGGCGCTGTCAGTAGGGGACAGCTATTTCCAGCACAAGAGTTTTGAGCGCATCCGTGAGTTCCAGCGGCTGGGAACATCGCTTTTGATa      >  8:173472/1‑149 (MQ=255)
                                                                                                                                           aaaTGCGCGTCGCATTTGCCGTGGAGACTGCGTTCCGTCCAGAAATACTTATCGTCGATGAGGCGCTGTCAGTTGGGGACAGCTATTTCCAGAACAAGAGTTTTGAGCGCATCCGTGAGTTCCAGCGGCTGGGAACATCGCTTTTGa        >  6:400740/1‑147 (MQ=255)
                                                                                                                                           aaaTGCGCGTAGCATTTGCCGTGGCGACTGCGTTCAGTCCAGAAAACCTTATCGTCGATGAGGCGCTGTCAGTAGGGGACAGCTATTTCCAGCACAAGAGATTAGAGCGAAAACGTCAGTTCCAGCGGCTGGGAACATCGCTTTTGa        >  3:96192/1‑147 (MQ=255)
                                                                                                                                             aTGCGCGTCGCATTTGCCGTGGCGACTGCGTTCCGTCCAGAAATCCTTATCGTCGATGAGGCGCTGTCAGTAGGGGACAGCTATTTCCAGCACAAGAGTTTTGAGCGCATCCGTGAGTTCCAGCGGCTGGGAACATCGCTTTTGATAgt    <  3:21145/149‑1 (MQ=255)
                                                                                                                                                cgcgTCGCATTTGCCGTGGCGACTGCGTTCCGTCCAGAAATCCTTATCGTCGATGAGGCGCTGTCAGTAGGGGACAGCTATTTCCAGCACAAGAGTTTTGAGCGCATCCGTGAGTTCCAGCGGCTGGGAACATCGCTTTTGATAgtgt  <  7:173472/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
CGCCAGAACGTGTACCACGCTGCGGGTCTGATGGGTTTCGATGTACAGCAGATCAATGATGTCATTGATGACATTGAAGCCTTTGCCGAGATTGGCGAGTACTTTGATGAGGCAGTGCGTACCTATTCCAGCGGTATGCAAATGCGCGTCGCATTTGCCGTGGCGACTGCGTTCCGTCCAGAAATCCTTATCGTCGATGAGGCGCTGTCAGTAGGGGACAGCTATTTCCAGCACAAGAGTTTTGAGCGCATCCGTGAGTTCCAGCGGCTGGGAACATCGCTTTTGATAGTGT  >  NC_002947/1987846‑1988137

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