Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,027,194 T→A 11.2% D427V (GAC→GTC)  PP_3552 ← two‑component system sensor histidine kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,027,1940TA11.2% 198.9 / 14.1 72D427V (GAC→GTC) PP_3552two‑component system sensor histidine kinase
Reads supporting (aligned to +/- strand):  ref base T (33/30);  new base A (4/4);  total (38/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.74e-01

CTTGCAACACTGCCTCCCCTTATTCCCTGCTGACTGCCAGCCGGCAGCACAGTACCAGACCGCCGTTGTCACCCGCCCGCGCCCACAGGCTGCCACCGAAGCCTTCGATCAGGCTGCGGCTCATGCTCAAGCCCAAGCCCAGGCCGTCGGCCTTGCTGGTGGTAAACGGCGTGAAGATTTCATCCAGCCGCTCGGGTGCCACGCCCGGGCCCTGGTCAGCTACCTCGACCAGGACCCCATCGCCATCGCGCGCCGCGCAGAGCAGGATACGTGACGGCCGCTCCCCCTGCTGCTC  >  NC_002947/4027048‑4027342
                                                                                                                                                  |                                                                                                                                                    
cttgcAACACTGCCTCCCCTTATTCCCTGCTGACTGCCAGCCGGCAGCACAGTACCAGACCGCCGTTGTCACCCGCCCGCGCCCACAGGCTGCCACCGAAGCCTTCGATCAGGCTGCGGCTCATGCTCAAGCCCAAGCCCAGGCCGTCgg                                                                                                                                                   <  1:29418/150‑1 (MQ=255)
 ttgcAACACTGCCTCCCCTTATTCCCTGCTGACTGCCAGCCGGCAGCACAGTACCAGACCGCCGTTGTCACCCGCCCGCGCCCACAGGCTGCCACCGAAGCCTTCGATCAGGCTGCGGCTCATGCTCAAGCCCAAGCCCAGGCCGTCGGc                                                                                                                                                  >  1:509193/1‑150 (MQ=255)
 ttgcAACACTGCCTCCCCTTATTCCCTGCTGACTGCCAGCCGGCAGCACAGTACCAGACCGCCGTTGTCACCCGCCCGCGCCCACAGGCTGCCACCGAAGCCTTCGATCAGGCTGCGGCGCATGCTCAAGCCCAAGCCCAGGCCGTCGGc                                                                                                                                                  >  2:784558/1‑150 (MQ=255)
     cgatcTGCCTCCCCTTATTCCCTGCTGACTGCCAGCCGGCAGCACAGTACCAGACCGCCGTTGTCACCCGCCCGCGCCCACAGGCTGCCACCGAAGCATTCTGTCAGGCTTCGGCTCATGCTCAAGCCCAAGCCCAGGCCGACGGCCtgg                                                                                                                                              <  1:36749/146‑3 (MQ=255)
      acacTGCCTCCCCTTATTCCCTGCTGACTGCCAGCCGGCAGCACAGTACCAGACCGCCGTTGTCACCCGCCCGCGCCCACAGGCTGCCACCGAAGCCTTCGATCAGGCTGCGGCTCATGCTCAAGCCCAAGCCCAGGCCGTCGGCCTTGc                                                                                                                                             <  2:135826/150‑1 (MQ=255)
      acacTGCCTCCCCTTATTCCCTGCTGACTGCCAGCCGGCAGCACAGTACCAGACCGCCGTTGTCACCCGCCCGCGCCCACAGGCTGCCACCGAAGCCTTCGATCAGGCTGCGGCTCATGCTCAAGCCCAAGCCCAGGCCGTCGGCCTTGc                                                                                                                                             >  1:1305318/1‑150 (MQ=255)
         cTGCCTCCCCTTATTCCCTGCTGACTGCCAGCCGGCAGCACAGTACCAGACCGCCGTTGTCACCCGCCCGCGCCCACAGGCTGCCACCGAAGCCTTCGATCAGGCTGCGGCTCATGCTCAAGCCCAAGCCCAGGCCGTCGGCCTTGCtgg                                                                                                                                          <  2:913702/150‑1 (MQ=255)
          tGCCTCCCCTTATTCCCTGCTGACTGCCAGCCGGCAGCACAGTACCAGACCGCCGTTGTCACCCGCCCGCGCCCACAGGCTGCCACCGAAGCCTTCGATCAGGCTGCGGCTCATGCTCAAGCCCAAGCCCAGGCCGACGGCCtgggct                                                                                                                                           >  1:186634/1‑143 (MQ=255)
          tGCCTCCCCTTATTCCCTGCTGACTGCCAGCCGGCAGCACAGTACCAGACCGCCGTTGTCACCCGCCCGCGCCCACAGGCTGCCACCGAAGCCTTCGATCAGGCTGCGGCTCATGCTCAAGCCCAAGCCCAGGCCGACGGCCtgggct                                                                                                                                           <  2:186634/148‑6 (MQ=255)
           gCCTCCCCTTATTCCCTGCTGACTGCCAGCCGGCAGCACAGTACCAGACCGCCGTTGTCACCCGCCCGCGCCCACAGGCTGCCACCGAAGCCTTCGATCAGGCTGCGGCTCATGCTCAAGCCCAAGCCCAGGCCGTCGGCCTTGCtggtg                                                                                                                                        <  1:1133631/150‑1 (MQ=255)
           gCCTCCCCTTATTCCCTGCTGACTGCCAGCCGGCAGCACAGTACCAGACCGCCGTTGTCACCCGCCCGCGCCCACAGGCTGCCACCGAAGCCTTCGATCAGGCTGCGGCTCATGCTCAAGCCCAAGCCCAGGCCGTCGGCCTTGCtggtg                                                                                                                                        <  1:1373285/150‑1 (MQ=255)
           gCCTCCCCTTATTCCCTGCTGACTGCCAGCCGGCAGCACAGTACCAGACCGCCGTTGTCACCCGCCCGCGCCCACAGGCTGCCACCGAAGCCTTCGATCAGGCTGCGGCTCATGCTCAAGCCCAAGCCCAGGCCGTCGGCCTTGCtggtg                                                                                                                                        <  1:182455/150‑1 (MQ=255)
            ccTCCCCTTATTCCCTGCTGACTGCCAGCCGGCAGCACAGTACCAGACCGCCGTTGTCACCCGCCCGCGCCCACAGGCTGCCACCGAAGCCTTCGATCAGGCTTCGGCTCATGCTCAAGCCCAAGCCCACGCCGGCCGCCGTggagttcg                                                                                                                                       >  2:36749/1‑143 (MQ=255)
               ccccTTATTCCCTGCTGACTGCCAGCCAGCCGCCCGGAGCCGGAGCTACTTTGTCGCCCGCCCGCGCCCACAGCATGCCCGCCCCTGCTTTGATCAGGCTGCGGCTCATGCTCAAGCCCAAGCCCAGGCCGTCGGCCTAGCTGGTGGTaa                                                                                                                                    >  1:324664/1‑150 (MQ=255)
                cccTTATTCCCTGCTGACTGCCAGCCGGCAGCACAGTACCAGACCGCCGTTGTCACCCGCCCGCGCCCACAGGCTGCCACCGAAGCCTTCGATCAGGCTGCGGCTCATGCTCAAGCCCAAGCCCAGGCCGTCGGcc                                                                                                                                                 <  2:208127/136‑1 (MQ=255)
                cccTTATTCCCTGCTGACTGCCAGCCGGCAGCACAGTACCAGACCGCCGTTGTCACCCGCCCGCGCCCACAGGCTGCCACCGAAGCCTTCGATCAGGCTGCGGCTCATGCTCAAGCCCAAGCCCAGGCCGTCGGcc                                                                                                                                                 >  1:208127/1‑136 (MQ=255)
                                  tGCCAGCCGGCAGCACAGTACCAGACCGCCGTTGTCACCCGCCCGCGCCCACAGGCTGCCACCGAAGCCTTCGATCAGGCTGCGGCTCATGCTCAAGCCCAAGCCCAGGCCGTCGGCCTTGCTGGTGGTAAACGGCGTGAAGATTTCATc                                                                                                                 >  1:430724/1‑150 (MQ=255)
                                  tGCCAGCCGGCAGCACAGTACCAGACCGCCGTTGTCACCCGCCCGCGCCCACAGGCTGCCACCGAAGCCTTCGATCAGGCTGCGGCTCATGCTCAAGCCCAAGCCCAGGCCGTCGGCCTTGCTGGTGGTAAACGGCGTGAAGATTTCATc                                                                                                                 >  1:1017192/1‑150 (MQ=255)
                                   gCCAGCCGGCAGCACAGTACCAGACCGCCGTTGTCACCCGCCCGCGCCCACAGGCTGCCACCGAAGCCTTCGATCAGGCTGCGGCTCATGCTCAAGCCCAAGCCCAGGCCGTCGGCCTTGCTGGTGGTAAACGGCGTGAAGATTTCATcc                                                                                                                >  2:1252545/1‑150 (MQ=255)
                                    ccAGCCGGCAGCACAGTACCAGACCGCCGTTGTCACCCGCCCGCGCCCACAGGCGGCCACCGAAGCCTTCGATCAGGCTGCGGCTCATGCTCAAGCCCCAGCCCAGGCCGTCGTCCTTACTGGTGGTAAAAGGCGTGAAGATTTCATCCa                                                                                                               >  2:292015/1‑150 (MQ=255)
                                     cAGCCGGCAGCACAGTACCAGACCGCCGTTGTCACCCGCCCGCGCCCACAGGCTGCCACCGAAGCCTTCGATCAGGCTGCGGCTCATGCTCAAGCCCAAGCCCAGGCCGTCGGCCTTGCTGGTGGTAAACGGCGTGAAGATTTCATCCAg                                                                                                              <  2:144914/150‑1 (MQ=255)
                                          ggCAGCACAGTACCAGACCGCCGTTGTCACCCGCCCGCGCCCACAGGCTGCCACCGAAGCCTTCGATCAGGCTGCGGCTCATGCTCAAGCCCAAGCCCAGGCCGTCGGCCTTGCTGGTGGTAAACGGCGTGAAGATTTCATCCAGCCGCt                                                                                                         >  1:441092/1‑150 (MQ=255)
                                               cacaGTACCAGACCGCCGTTGTCACCCGCCCGCGCCCACAGGCTGCCACCGAAGCCTTCGATCAGGCTGCGGCTCATGCTCAAGCCCAAGCCCAGGCCGTCGGCCTTGCTGGTGGTAAACGGCGTGAAGATTTCATCCAGCCGCTCGGGt                                                                                                    >  1:564029/1‑150 (MQ=255)
                                                acaGTACCAGACCGCCGTTGTCACCCGCCCGCGCCCACAGGCTGCCACCGAAGCCTTCGATCAGGCTGCGGCTCATGCTCAAGCCCAAGCCCAGGCCGTCGGCCTTGCTGGTGGTAAACGGCGTGAAGATTTCATCCAGCCGCTCGGGTg                                                                                                   >  2:182795/1‑150 (MQ=255)
                                                 caGTACCAGACCGCCGTTGTCACCCGCCCGCGCCCACAGGCTGCCACCGAAGCCTTCGATCAGGCTGCGGCTCATGCTCAAGCCCAAGCCCAGGCCGTCGGCCTTGCTGGTGGTAAACGGCGTGAAGATTTCATCCAGCCGCTCGGGTGc                                                                                                  <  2:430724/150‑1 (MQ=255)
                                                  aGTACCAGACCGCCGTTGTCACCCGCCCGCGCCCACAGGCTGCCACCGAAGCCTTCGATCAGGCTGCGGCTCATGCTCAAGCCCAAGCCCAGGCCGTCGGCCTTGCTGGTGGTAAACGGCGTGAAGATTTCATCCAGCCGCTCGGGTGcc                                                                                                 <  1:620694/150‑1 (MQ=255)
                                                  aGTACCAGACCGCCGTTGTCACCCGCCCGCGCCCACAGGCTGCCACCGAAGCCTTCGATCAGGCTGCGGCTCATGCTCAAGCCCAAGCCCAGGCCGTCGGCCTTGCTGGTGGTAAACGGCGTGAAGATTTCATCCAGCCGCTCGGGTGcc                                                                                                 <  1:1205122/150‑1 (MQ=255)
                                                    tACCAGACCGCCGTTGTCACCCGCCCGCGCCCACAGGCTGCCACCGAAGCCTTCGATCAGGCTGCGGCTCATGCTCAAGCCCAAGCCCAGGCCGTCGGCCTTGCTGGTGGTAAACGGCGTGAAGATTTCATCCAGCCGCTCGGGTGCCAc                                                                                               <  2:1017192/150‑1 (MQ=255)
                                                       cAGACCGCCGTTGTCACCCGCCCGCGCCCACAGGCTGCCACCGAAGCCTTCGATCAGGCTGCGGCTCATGCTCAAGCCCAAGCCCAGGCCGACGGCCtgggc                                                                                                                                            <  1:170167/102‑5 (MQ=255)
                                                       cAGACCGCCGTTGTCACCCGCCCGCGCCCACAGGCTGCCACCGAAGCCTTCGATCAGGCTGCGGCTCATGCTCAAGCCCAAGCCCAGGCCGACGGCCtgggc                                                                                                                                            >  2:170167/1‑98 (MQ=255)
                                                          aCCGCCGTTGTCACCCGCCCGCGCCCACAGGCTGCCACCGAAGCCTTCGATCAGGCTGCGGCTCATGCTCAAGCCCAAGCCCAGGCCGTCGGCCTTGCTGGTGGTAAACGGCGTGAAGATTTCATCCAGCCGCTCGGGTGCCACGCCCgg                                                                                         <  1:182795/150‑1 (MQ=255)
                                                            cgccgTTGTCACCCGCCCGCGCCCACAGGCTGCCACCGAAGCCTTCGATCAGGCTGCGGCTCATGCTCAAGCCCAAGCCCAGGCCGTCGGCCTTGCTGGTGGTAAACGGCGTGAAGATTTCATCCAGCCGCTCGGGTGCCACGCCCGGGc                                                                                       >  1:908278/1‑150 (MQ=255)
                                                               cgTTGTCACCCGCCCGCGCCCACAGGCTGCCACCGAAGCCTTCGATCAGGCTGCGGCTCATGCTCAAGCCCAAGCCCAGGCCGTCGGCCTTGCTGGTGGTAAACGGCGTGAAGATTTCATCCAGCCGCTCGGGTGCCACGCCCGGGCCCt                                                                                    >  1:343442/1‑150 (MQ=255)
                                                                 ttGTCACCCGCCCGCGCCCACAGGCTGCCACCGAAGCCTTCGATCAGGCTGCGGCTCATGCTCAAGCCCAAAAAGAGGTCGTCGGCCTTCGGGGTTTTTAACGGCGGGAAGATTTCATCCAGCCGCTCGGGTGCCACGCCCGGGCCCTgg                                                                                  <  1:531906/150‑1 (MQ=255)
                                                                   gTCACCCGCCCGCGCCCACAGGCTGCCACCGAAGCCTTCGATCAGGCTGCGGCTTTTGCTCAAGCCCAAGCCCAGGCCGTCGGCCTTGCTGGTGGTAAACGGCGTGAAGATTTCATCCAGCCGCTCGGGTGCCACGCCCGGGCCCTGGTc                                                                                <  2:399636/150‑1 (MQ=255)
                                                                    tCACCCGCCCGCGCCCACAGGCTGCCACCGAAGCCTTCGATCAGGCTGCGGCTCATGCTCAAGCCCAAGCCCAGGCCGTCGGCCTTGCTGGTGGTAAACGGCGTGAAGATTTCATCCAGCCGCTCGGGTGCCACGCCCGGGCCCTGGTCa                                                                               <  2:1251276/150‑1 (MQ=255)
                                                                        ccgcccgcGCCCACAGGCTGCCACCGAAGCCTTCGATCAGGCTGCGGCTCATGCTCAAGCCCAAGCCCAGGCCGTCGGCCTTGCTGGTGGTAAACGGCGTGAAGATTTCATCCAGCCGCTCGGGTGCCACGCCCGGGCCCTGGTCAGCTa                                                                           >  1:1147582/1‑150 (MQ=255)
                                                                        ccgcccgcGCCCACAGGCTGCCACCGAAGCCTTCGATCAGGCTGCGGCTCATGCTCAAGCCCAAGCCCAGGCCGTCGGCCTTGCTGGTGGTAAACGGCGGGAAGATTTCATCCAGCCGCTCGGGTGCCACGCCCGGGCCCTGGTCAGCTa                                                                           >  1:489622/1‑150 (MQ=255)
                                                                         cgcccgcGCCCACAGGCTGCCACCGAAGCCTTCGATCAGGCTGCGGCTCATGCTCAAGCCCAAGCCCAGGCCGTCGGCCTTGCTGGTGGTAAACGGCGTGAAGATTTCATCCAGCCGCTCGGGTGCCACGCCCGGGCCCTGGTCAGCTAc                                                                          >  1:1132917/1‑150 (MQ=255)
                                                                            ccgcGCCCACAGGCTGCCACCGAAGCCTTCGATCAGGCTGCGGCTCATGCTCAAGCCCAAGCCCAGGCCGTCGGCCTTGCTGGTGGTAAACGGCGTGAAGATTTCATCCAGCCGCTCGGGTGCCACGCCCGGGCCCTGGTCAGCTACCTc                                                                       >  1:1277530/1‑150 (MQ=255)
                                                                               cgcCCACAGGCTGCCACCGAAGCCTTCGATCAGGCTGCGGCTCATGCTCAAGCCCAAGCCCAGGCCGTCGGCCTTGCTGGTGGTAAACGGCGTGAAGATTTCATCCAGCCGCTCGGGTGCCACGCCCGGGCCCTGGTCAGCTACCTCGAc                                                                    <  2:1305318/150‑1 (MQ=255)
                                                                                   cacaGGCTGCCACCGAAGCCTTCGATCAGGCTGCGGCTCATGCTCAAGCCCAAGCCCAGGCCGTCGGCCTTGCTGGTGGTAAACGGCGTGAAGATTTCATCCAGCCGCTCGGGTGCCACGCCCGGGCCCTGGTCAGCTACCTCGACCAgg                                                                <  1:1252545/150‑1 (MQ=255)
                                                                                    acaGGCTGCCACCGAAGCCTTCGATCAGGCTGCGGCTCATGCTCAAGCCCAAGCCCTGGCCGTCGGCCTTGCTGGTGGTAAACGGCGTGAAGATTTCATCCAGCCGCTCGGGTGCCACGCCCGGGCCCTGGTCAGCTACCTCGACCAGGa                                                               >  2:923657/1‑150 (MQ=255)
                                                                                    acaGGCTGCCACCGAAGCCTTCGATCAGGCTGCGGCTCATGCTCAAGCCCAAGCCCAGGCCGTCGGCCTTGCtg                                                                                                                                           >  2:456679/1‑74 (MQ=255)
                                                                                    acaGGCTGCCACCGAAGCCTTCGATCAGGCTGCGGCTCATGCTCAAGCCCAAGCCCAGGCCGTCGGCCTTGCtg                                                                                                                                           <  1:456679/74‑1 (MQ=255)
                                                                                     caGGCTGCCACCGAAGCCTTCGATCAGGCTGCGGCTCATGCTCAAGCCCAAGCCCAGGCCGTCGGCCTTGCTGGTGGTAAACGGCGTGAAGATTTCATCCAGCCGCTCGGGTGCCACGCCCGGGCCCTGGTCAGCTACCTCGACCAGGAc                                                              <  2:324664/150‑1 (MQ=255)
                                                                                        gCTGCCACCGAAGCCTTCGATCAGGCTGCGGCTCATGCTCAAGCCCAAGCCCAGGCCGTCGGCCTTGCTGGTGGTAAACGGCGTGAAGATTTCATCCAGCCGCTCGGGTGCCACGCCCGGGCCCTGGTCAGCTACCTCGACCAGGAcccc                                                           >  1:1076988/1‑150 (MQ=255)
                                                                                                 gAAGCCTTCGATCAGGCTGCGGCTCATGCTCAAGCCCAAGCCCAGGCCGTCGGCCTTGCtggtgg                                                                                                                                       >  1:117414/1‑65 (MQ=255)
                                                                                                 gAAGCCTTCGATCAGGCTGCGGCTCATGCTCAAGCCCAAGCCCAGGCCGTCGGCCTTGCtggtgg                                                                                                                                       <  2:117414/65‑1 (MQ=255)
                                                                                                       ttCGATCAGGCTGCGGCTCATGCTCAAGCCCAAGCCCAGGCCGTCGGCCTTGCTGGTGGTAAACGGCGTGAAGATTTCATCCAGCCGCTCGGGTGCCACGCCCGGGCCCTGGTCAGCTACCTCGACCAGGACCCCATCGCCATCgcgcgc                                            >  2:800826/1‑150 (MQ=255)
                                                                                                               ggCTGCGGCTCATGCTCAAGCCCAAGCCCAGGCCGTCGGCCTTGCTGGTGGTAAACGGCGTGAAGATTTCATCCAGCCGCTCGGGTGCCACGCCCGGGCCCTGGTCAGCTACCTCGACCAGGACCCCATCGCCATCGCGCGCCGCGCaga                                    >  2:171436/1‑150 (MQ=255)
                                                                                                                  tGCGGCTCATGCTCAAGCCCAAGCCCAGGCCGTCGGCCTTGCTGGTGGTAAACGGCGTGAAGATTTCATCCAGCCGCTCGGGTGCCACGCCCGGGCCCTGGTCAGCTACCTCGACCAGGACCCCATCGCCATCGCGCGCCGCGCAGAGCa                                 <  2:1122358/150‑1 (MQ=255)
                                                                                                                   gCGGCTCATGCTCAAGCCCAAGCCCAGGCCGTCGGCCTTGCTGGTGGTAAACGGCGTGAAGATTTCATCCAGCCGCTCGGGTGCCACGCCCGGGCCCTGGTCAGCTACCTCGACCAGGACCCCATCGCCATCGCGCGCCGCGCAGAGCAg                                <  2:564029/150‑1 (MQ=255)
                                                                                                                       cTCATGCTCAAGCCCAAGCCCAGGCCGTCGGCCTTGCTGGTGGTAAACGGCGTGAAGATTTCATCCAGCCGCTCGGGTGCCACGCCCGGGCCCTGGTCAGCTACCTCGACCAGGACCCCATCGCCATCGCGCGCCGCGCAGAGCAGGATa                            <  2:1132917/150‑1 (MQ=255)
                                                                                                                        tCATGCTCAAGCCCAAGCCCAGGCCGTCGGCCTTGCTGGTGGTAAACGGCGTGAAGa                                                                                                                        >  1:922080/1‑57 (MQ=255)
                                                                                                                        tCATGCTCAAGCCCAAGCCCAGGCCGTCGGCCTTGCTGGTGGTAAACGGCGTGAAGa                                                                                                                        <  2:922080/57‑1 (MQ=255)
                                                                                                                        tCATGCTCAAGCCCAAGCCCAGGCCGTCGGCCTTGCTGGTGGTAAACGGCGTGAAGAt                                                                                                                       <  1:810003/58‑1 (MQ=255)
                                                                                                                        tCATGCTCAAGCCCAAGCCCAGGCCGTCGGCCTTGCTGGTGGTAAACGGCGTGAAGAt                                                                                                                       >  2:810003/1‑58 (MQ=255)
                                                                                                                            gCTCAAGCCCAAGCCCAGGCCGTCGGCCTTGCTGGTGGTAAACGGCGTGAAGATTTCATCCAGCCGCTCGGGTGCCACGCCCGGGCCCTGGTCAGCTACCTCGACCAGGACCCCATCGCCATCGCGCGCCGCGCAGAGCAGGATACGTGa                       <  2:1147582/150‑1 (MQ=255)
                                                                                                                              tCAAGCCCAAGCCCAGGCCGTCGGCCTTGCTGGTGGTAAACGGCGTGAAGATTTCATCCAGCCGCTCGGGTGCCACGCCCGGGCCCTGGTCAGCTACCTCGACCAGGa                                                               >  2:430111/1‑108 (MQ=255)
                                                                                                                              tCAAGCCCAAGCCCAGGCCGTCGGCCTTGCTGGTGGTAAACGGCGTGAAGATTTCATCCAGCCGCTCGGGTGCCACGCCCGGGCCCTGGTCAGCTACCTCGACCAGGa                                                               <  1:430111/108‑1 (MQ=255)
                                                                                                                              tCAAGCCCAAGCCCAGGCCGTCGGCCTTGCTGGTGGTAAACGGCGTGAAGATTTCATCCAGCCGCTCGGGTGCCACGCCCGGGCCCTGGTCAGCTACCTCGACCAGGACCCCATCGCCATCGCGCGCCGCGCAGAGCAGGATACGTGACg                     >  1:781024/1‑150 (MQ=255)
                                                                                                                              tCAAGCCCAAGCCCAGGCCGTCGGCCTTGCTGGTGGTAAACGGCGTGAAGATTTCATCCAGCCGCTCGGGTGCCACGCCCGGGCCCTGGTCAGCTACATCGACCAGGACCCCATCGCCATCGCGCGCCGCGCAGAGCAGGATACGTGACg                     >  1:983068/1‑150 (MQ=255)
                                                                                                                                     caagcccaGGCCGTCGGCCTTGCTGGTGGTAAACGGCGTGAAGATTTCATCCAGCCGCTCGGGTGCCACGCCCGGGCCCTGGTCAGCTACCTCGACCAGGACCCCATCGCCATCGCGCGCCGCGCAGAGCAGGATACGTGACGGCCGCTc              >  2:1153977/1‑150 (MQ=255)
                                                                                                                                       cagcaaGGCCGACGGCCTTGCTGGTGGTAAACGGCGTGAAGATTTCATCCAGCCGCTCGGGTGCCACGCCCGGGCCCTGGTCAGCTACCTCGACCAGGACCCCATCGCCATCGCGCGCCGCGCAGAGCAGGATACGTGACGGCCGCTccc            >  2:1019759/6‑150 (MQ=255)
                                                                                                                                         gcaaGGCCGACGGCCTTGCTGGTGGTAAACGGCGTGAAGATTTCATCCAGCCGCTCGGGTGCCACGCCCGGGCCCTGGTCAGCTACCTCGACCAGGACCCCATCGCCATCGCGCGCCGCGCagag                                   <  2:832784/122‑1 (MQ=255)
                                                                                                                                         gcaaGGCCGACGGCCTTGCTGGTGGTAAACGGCGTGAAGATTTCATCCAGCCGCTCGGGTGCCACGCCCGGGCCCTGGTCAGCTACCTCGACCAGGACCCCATCGCCATCGCGCGCCGCGCagag                                   >  1:832784/4‑125 (MQ=255)
                                                                                                                                          ccaGGCCGTCGGCCTTGCTGGTGGTAAACGGCGTGAAGATTTCATCCAGCCGCTCGGGTGCCACGCCCGGGCCCTGGTCAGCTACCTCGACCAGGACCCCATCGCCATCGCGCGCCGCGCAGAGCAGGATACGTGACGGCCGCTCCCCct         >  2:225931/1‑150 (MQ=255)
                                                                                                                                            aGGCCGTCGGCCTTGCTGGTGGTAAACGGCGTGAAGATTTCATCCAGCCGCTCGGGTGCCACGCCCGGGCCCTGGTCAGCTACCTCGACCAGGACCCCATCGCCATCGCGCGCCGCGCAGAGCAGGATACGTGACGGCCGCTCCCCctgc       <  1:955846/150‑1 (MQ=255)
                                                                                                                                             ggCCGTCGGCCTTGCTGGTGGTAAACGGCGTGAAGATTTCATCCAGCCGCTCGGGTGCCACGCCCGGGCCCTGGTCAGCTACCTCGACCAGGACCCCATCGCCATCGCGCGCCGCGCAGAGCAGGATACGTGACGGCCGCTCCCCctgct      >  2:1163330/1‑150 (MQ=255)
                                                                                                                                               ccGTCGGCCTTGCTGGTGGTAAACGGCGTGAAGATTTCATCCAGCCGCTCGGGTGCCACGCCCGGGCCCTGGTCAGCTACCTCGACCAGGACCCCATCGCCATCGCGCGCCGCGCAGAGCAGGATACGTGACGGCCGCTCCCCCtgctgc    <  2:441092/150‑1 (MQ=255)
                                                                                                                                                 gTCGGCCTTGCTGGTGGTAAACGGCGTGAAGATTTCATCCAGCCGCTCGGGTGCCACGCCCGGGCCCTGGTCAGCTACCTCGACCAGGACCCCATCGCCATCGCGCGCCGCGCAGAGCAGGATACGTGACGGCCGCTCCCCCTGCTGCTc  >  1:1233371/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
CTTGCAACACTGCCTCCCCTTATTCCCTGCTGACTGCCAGCCGGCAGCACAGTACCAGACCGCCGTTGTCACCCGCCCGCGCCCACAGGCTGCCACCGAAGCCTTCGATCAGGCTGCGGCTCATGCTCAAGCCCAAGCCCAGGCCGTCGGCCTTGCTGGTGGTAAACGGCGTGAAGATTTCATCCAGCCGCTCGGGTGCCACGCCCGGGCCCTGGTCAGCTACCTCGACCAGGACCCCATCGCCATCGCGCGCCGCGCAGAGCAGGATACGTGACGGCCGCTCCCCCTGCTGCTC  >  NC_002947/4027048‑4027342

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