Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,698,165 G→C 5.5% V267L (GTT→CTT)  hslU → protease HslVU ATPase subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,698,1650GC 5.5% 235.9 / 4.0 73V267L (GTT→CTT) hslUprotease HslVU ATPase subunit
Reads supporting (aligned to +/- strand):  ref base G (36/33);  new base C (2/2);  total (38/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.65e-01

GAAAGAAGCGCTGAAGATGGTTCGCGATGAAGAAGCGAGCCGCCTGGTCAACGAGGAAGAGCTCAAGGCCAAGGCCCTGGAAGCAGTCGAACAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGACCTGCTGCCGCTGATCGAAGGCTGCACCGTCAACACCAAGCTGGGCATGGTCAAGACCGACCACATTCTGTTCATTGCCTCGGGTGCATTCCACCTGAGCAAGCCGAG  >  NC_002947/5698020‑5698313
                                                                                                                                                 |                                                                                                                                                    
gAAAGAAGCGCTGAAGATGGTTCGCGATGAAGAAGCGAGCCGCCTGGTCAACGAGGAAGAGCTCAAGGCCAAGGCCCTGGAAGCAGTCGAACAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCACGCTTGGCAACCTTg                                                                                                                                                   >  4:18597/1‑149 (MQ=255)
gAAAGAAGCGCTGAAGATGGTTCGCGATGAAGAAGCGAGCCGCCTGGTCAACGAGGAAGAGCTCAAGGCCAAGGCCCTGGAAGCAGTCGAACAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTg                                                                                                                                                   >  3:326091/1‑149 (MQ=255)
 acagaagCGCTGAAGATGGAGCGCGATGAAGAAGCGAGCCGCCTGGTCAACGAGGAAGAGCTCAAGGCCAAGGCCCTGGAAGCAGTCGAACAGCACGGTATCGTGTTCATCGACGAAATCGCCAACGTTGCACGCTTGGCAACCTTgtc                                                                                                                                                 <  8:223423/147‑3 (MQ=255)
 aaaGAAGCGCTGAAGATGGTTCGCGATGAAGAAGCGAGCCGCCTGGTCAACGAGGAAGAGCTCAAGGCCAAGGCCCTGGAAGCAGTCGAACAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCACGCTTGGCAACCTTg                                                                                                                                                   <  3:18597/148‑1 (MQ=255)
  aagaagCGCTGAAGATGGTTCGCGATGAAGAAGCGAGCCGCCTGGTCAACGAGGAAGAGCTCAAGGCCAAGGCCCTGGAAGCAGTCGAACAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTgg                                                                                                                                                  <  1:331994/148‑1 (MQ=255)
   agaagCGCTGAAGATGGTTCGCGATGAAGAAGCGAGCCGCCTGGTCAACGAGGAAGAGCTCAAGGCCAAGGCCCTGGAAGCAGTCGAACAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGtt                                                                                                                                                    >  6:295621/1‑145 (MQ=255)
   agaagCGCTGAAGATGGTTCGCGATGAAGAAGCGAGCCGCCTGGTCAACGAGGAAGAGCTCAAGGCCAAGGCCCTGGAAGCAGTCGAACAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGtt                                                                                                                                                    <  5:295621/145‑1 (MQ=255)
        cgcTGAAGATGGTTCGCGATGAAGAAGCGAGCCGCCTGGTCAACGAGGAAGAGCTCAAGGCCAAGGCCCTGGAAGCAGTCGAACAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGcc                                                                                                                                           <  7:243194/149‑1 (MQ=255)
            gAAGATGGTTCGCGATGAAGAAGCGAGCCGCCTGGTCAACGAGGAAGAGCTCAAGGCCAAGGCCCTGGAAGCAGTCGAACAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACg                                                                                                                                       >  7:315587/1‑149 (MQ=255)
                        cgATGAAGAAGCGAGCCGCCTGGTCAACGAGGAAGAGCTCAAGGCCAAGGCCCTGGAAGCAGTCGAACAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAg                                                                                                                           <  1:196375/149‑1 (MQ=255)
                        cgATGAAGAAGCGAGCCGCCTGGTCAACGAGGAAGAGCTCAAGGCCAAGGCCCTGGAAGCAGTCGAACAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAg                                                                                                                           >  3:51896/1‑149 (MQ=255)
                         gATGAAGAAGCGAGCCGCCTGGGCAACGAGGAAGAGCGCAATGCCAAGGCCCTGGAAGCAGTCGAACACCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAgg                                                                                                                          <  4:276277/149‑1 (MQ=255)
                             aagaagCGAGCCGCCTGGTCAACGAGGAAGAGCTCAAGGCCAAGGCCCTGGAAGCAGTCGAACAGCGCGGTATCGTGTTCATCGACGAAATCGCCAACGTTGCCACGCTTGGCAACCTTgtcga                                                                                                                                               >  8:311495/1‑123 (MQ=38)
                                 agCGAGCCGCCTGGTCAACGAGGAAGAGCTCAAGGCCAAGGCCCTGGAAGCAGTCGAACAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAgc                                                                                                                  <  4:326091/149‑1 (MQ=255)
                                  gCGAGCCGCCTGGTCAACGAGGAAGAGCTCAAGGCCAAGGCCCTGGAAGCAGTCGAACAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAgcg                                                                                                                 <  6:257604/149‑1 (MQ=255)
                                    gAGCCGCCTGGTCAACGAGGAAGAGCTCAAGGCCAAGGCCCTGGAAGCAGTCGAACAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGcgcg                                                                                                               >  4:308037/1‑149 (MQ=255)
                                      gccgccTGGTCAACGAGGAAGAGCTCAAGGCCAAGGCCCTGGAAGCAGTCGAACAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGAc                                                                                                             <  8:161320/149‑1 (MQ=255)
                                         gccTGGTCAACGAGGAAGAGCTCAAGGCCAAGGCCCTGGAAGCAGTCGAACAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGACctg                                                                                                          <  4:51896/149‑1 (MQ=255)
                                         gccTGGTCAACGAGGAAGAGCTCAAGGCCAAGGCCCTGGAAGCAGTCGAACAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGACctg                                                                                                          <  6:344635/149‑1 (MQ=255)
                                            tGGTCAACGAGGAAGAGCTCAAGGCCAAGGCCCTGGAAGCAGTCGAACAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGCCGCCGACGTGTCCCGTGAAGGCGTACAGCGCCACctgctg                                                                                                       >  3:126025/1‑149 (MQ=255)
                                                 aaCGAGGAAGAGCTCAAGGCCAAGGCCCTGGAAGCAGTCGAACAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACa                                                                                                                    <  1:153712/131‑1 (MQ=255)
                                                 aaCGAGGAAGAGCTCAAGGCCAAGGCCCTGGAAGCAGTCGAACAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACa                                                                                                                    >  2:153712/1‑131 (MQ=255)
                                                 aaCGAGGAAGAGCTCAAGGCCAAGGCCCTGGAAGCAGTCGAACAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGACCTGCTGCCGCt                                                                                                  >  2:122870/1‑149 (MQ=255)
                                                 aaCGAGGAAGAGCTCAAGGCCAAGGCCCTGGAAGCAGTCGAACAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGACCTGCTGCCGCt                                                                                                  >  4:151876/1‑149 (MQ=255)
                                                        aaGAGCTCAAGGCCAAGGCCCTGGAAGCAGTCGAACAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGACCTGCTGCCGCTGATCGaa                                                                                           <  3:151876/149‑1 (MQ=255)
                                                        aaGAGCTCAAGGCCAAGGCCCTGGAAGCAGTCGAACAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGACCTGCTGCCGCTGATCGaa                                                                                           <  8:101142/149‑1 (MQ=255)
                                                         agagCTCAAGGCCAAGGCCCTGGAAGCAGTCGAACAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAg                                                                                                                           <  3:245092/116‑1 (MQ=255)
                                                         agagCTCAAGGCCAAGGCCCTGGAAGCAGTCGAACAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAg                                                                                                                           >  4:245092/1‑116 (MQ=255)
                                                         agagCTCAAGGCCAAGGCCCTGGAAGCAGTCGAACAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGACCTGCTGCCGCTGATCGAAg                                                                                          >  6:144900/1‑149 (MQ=255)
                                                          gagCTCAAGGCCAAGGCCCTGGAAGCAGTCGAACAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGACCTGCTGCCGCTGATCGAAgg                                                                                         >  3:223346/1‑149 (MQ=255)
                                                          gagCTCAAGGCCAAGGCCCTGGAAGCAGTCGAACAGCACGGTATCCTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTAGCCGCCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGACCTGCTGCAGCTGATCGAAgg                                                                                         >  5:291324/1‑149 (MQ=255)
                                                           agCTCAAGGCCAAGGCCCTGGAAGCAGTCGAACAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGACCTGCTGCCGCTGATCGAAGGc                                                                                        <  8:315587/149‑1 (MQ=255)
                                                             cTCAAGGCCAAGGCCCTGGAAGCAGTCGAACAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTc                                                                                                                                   <  5:319305/104‑1 (MQ=255)
                                                             cTCAAGGCCAAGGCCCTGGAAGCAGTCGAACAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTc                                                                                                                                   >  6:319305/1‑104 (MQ=255)
                                                                   gccaaggccCTGGAAGCAGTCGAACAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGACCTGCTGCCGCTGATCGAAGGCTGCACCGt                                                                                <  5:144900/149‑1 (MQ=255)
                                                                    ccaaggccCTGGAAGCAGTCGAACAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGACCTGCTGCCGCTGATCGAAGGCTGCACCGTc                                                                               <  3:140398/149‑1 (MQ=255)
                                                                           ccTGGAAGCAGTCGAACAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACg                                                                                                                                       <  7:240354/86‑1 (MQ=255)
                                                                           ccTGGAAGCAGTCGAACAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACg                                                                                                                                       >  8:240354/1‑86 (MQ=255)
                                                                                 aGCAGTCGAACAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGACCTGCTGCCGCTGATCGAAGGCTGCACCGTCAACACCAAGCTgg                                                                  >  4:88919/1‑149 (MQ=255)
                                                                                    aGTCGAACAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGACCTGCTGCCGCTGa                                                                                                <  3:332913/116‑1 (MQ=255)
                                                                                    aGTCGAACAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGACCTGCTGCCGCTGa                                                                                                >  4:332913/1‑116 (MQ=255)
                                                                                     gTCGAACAGCACGGTATCGTGTTCATCGTCGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGACCTGCTGcc                                                                                                     >  1:194955/1‑110 (MQ=255)
                                                                                     gTCGAACAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGACCTGCTGcc                                                                                                     <  2:194955/110‑1 (MQ=255)
                                                                                          aCAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGACCTGCTGCCGCTGATCGAAGGCTGCACCGTCAACACCAAGCTGGGCATGGTc                                                          >  4:311420/1‑148 (MQ=255)
                                                                                             ccgatctATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGACCTGCTGCCGCTGATCGAAGGCTGCTCCGTCAACACCAAGCTGGGCATGGTCAAga                                                      <  1:32758/143‑1 (MQ=255)
                                                                                               aCGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGACCTGCTGCCGCTGATCGAAGGCTGCACCGTCAac                                                                            <  2:232590/125‑1 (MQ=255)
                                                                                               aCGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGACCTGCTGCCGCTGATCGAAGGCTGCACCGTCAac                                                                            >  1:232590/1‑125 (MQ=255)
                                                                                                cGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGACCTGCTGCCGCTGATCGAAGGCTGCACCGTCAACACCAAGCTGGGCATGGTCAAgaccg                                                   <  1:122870/149‑1 (MQ=255)
                                                                                                cGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGACCTGCTGCCGCTGATCGAAGGCTGCACCGTCAACACCAAGCTGGGCATGGTCAAgacc                                                    >  1:340958/1‑148 (MQ=255)
                                                                                                   tATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGACCTGCTGCCGCTGATCGAAGGCTGCACCGTCAACACCAAGCTGGGCATGGTCAAgaccgacc                                                <  3:88919/149‑1 (MQ=255)
                                                                                                     tCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGACCTGCTGCCGCTGATCGAAGGCCGCACCGTCAACACCAAGCTGGGCATGGTCAAga                                                      <  6:290010/141‑1 (MQ=255)
                                                                                                     tCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGACCTGCTGCCGCTGATCGAAGGCAGCACCGTCAACACCAAGCTGGGCATGGTCAAga                                                      >  5:290010/1‑141 (MQ=255)
                                                                                                     tCGTGTTAATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGACCTGCTGACGCTGATCGAAGGCTGCACCGTCAACACCAAGCTGGGCATGGTCaaaa                                                      >  2:32758/1‑139 (MQ=255)
                                                                                                        tgtTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGACCTGCTGCCGCTGATAGAAGGCTGCACCGTCAACACCAAGCTGGGCATGGTCAAGACCGACCACAtt                                           >  7:143727/1‑149 (MQ=255)
                                                                                                         gtTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGACCTGCTGCCGCTGATCGAAGGCTGCACCGTCAACACCAAGCTGGGCATGGTCAAGACCGACCACATTc                                          >  7:34820/1‑149 (MQ=255)
                                                                                                            cATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGACCTGCTGCCGCTGATCGAAGGCTGCACCGTCAACACCAAGCTGGGCATGGTCAAGACCGACCACATTCTGt                                       >  7:183401/1‑149 (MQ=255)
                                                                                                             aTCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGACCTGCTGCCGCTGATCGAAGGCTGCACCGTCAACACCAAGCTGGGCATGGTCAAGACCGACCACATTCTGtt                                      >  3:196775/1‑149 (MQ=255)
                                                                                                                 acgaAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGACCTGCTGCCGCTGATCGAAGGCTGCACCGTCAACACCAAGCTGGGCATGGTCAAGACCGACCACATTCTGTTCAtt                                  <  3:311420/149‑1 (MQ=255)
                                                                                                                   gaAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGACCTGCTGCCGCTGATCGAAGGCTGCACCGTCAACACCAAGCTGGGCATGGTCAAGACCGACCACATTCTGTTCATTGc                                >  7:69448/1‑149 (MQ=255)
                                                                                                                      aTCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGACCTGCTGCCGCTGATCGAAGGCTGCACCGTCAACACCAAGCTGGGCATGGTCAAGACCGACCACATTCTGTTCATTGCCTc                             <  3:308037/149‑1 (MQ=255)
                                                                                                                         gACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGACCTGCTGCCGCTGATCGAAGGCTGCACCGTCAACACCAAGCTGGGCATGGTCAAGACCGACCACATTCTGTTCATTGCCTCggg                          >  2:287505/1‑149 (MQ=255)
                                                                                                                          aCAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGACCTGCTGCCGCTGATCGAAGGCTGCACCGTCAACACCAAGCTGGGCATGGTCAAGACCGACCACATTCTGTTCATTGCCTCGGGt                         <  8:183401/149‑1 (MQ=255)
                                                                                                                             acgTTGCCACGCTTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGACCTGCTGCCGCTGATCGAAGGCTGCACCGTCAACACCAAGCTGGGCATGGTCAAGACCGACCACATTCTGTTCATTGCCTCGGGGGCa                      >  4:195665/3‑149 (MQ=255)
                                                                                                                               gTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGACCTGCTGCCGCTGATCGAAGGCTGCACCGTCAACACCAAGCTGGGCATGGTCAAGACCGACCACATTCTGTTCATTGCCTCGGGTGCAtt                    >  6:97323/1‑149 (MQ=255)
                                                                                                                                   ccAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGACCTGCTGCCGCTGATCGAAGGCTGCACCGTCAACACCAAGCTGGGCATGGTCAAGACCGACCACATTCTGTTCATTGCCTCGGGTGCATTCCa                 <  4:196775/148‑1 (MQ=255)
                                                                                                                                   ccAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGACCTGCTGCCGCTGATCGAAGGCTGCACCGTCAACACCAAGCTGGGCATGGTCAAGACCGACCACATTCTGTTCATTGCCTCGGGTGCATTCCAc                <  2:340958/149‑1 (MQ=255)
                                                                                                                                    cAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGACCTGCTGCCGCTGATCGAAGGCTGCACCGTCAACACCAAGCTGGGCATGGTCAAGACCGACCACATTCTGTTCAtt                                  >  7:185906/1‑130 (MQ=255)
                                                                                                                                    cAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGACCTGCTGCCGCTGATCGAAGGCTGCACCGTCAACACCAAGCTGGGCATGGTCAAGACCGACCACATTCTGTTCAtt                                  <  8:185906/130‑1 (MQ=255)
                                                                                                                                     aaGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGAACAGCGCGACCTGCTGCCGCTGATCGAAGGCTGCACCGTCAACACCAAGCTGGGCATGGTCAAGACCGACCACATTCTGTTCATTGCCTCGGGTGCATTCCACCt              >  1:243843/1‑149 (MQ=255)
                                                                                                                                       gCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGACCTGCTGCCGCTGATCGAAGGCTGCACCGTCAACACCAAGCTGGGCATGGTCAAGACCGACCACATTCTGTTCATTGCCTCGGGTGCATTCCACCt              >  8:228124/1‑147 (MQ=255)
                                                                                                                                             cAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGACCTGCTGCCGCTGATCGAAGGCTTCACCGTCAACACCAAGCTGGGCATGGTCAAGACCGACCACAtt                                           <  1:242721/112‑1 (MQ=255)
                                                                                                                                             cAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCCCGACCTGCTGCCGCTGATCGAACGCTGCACCGTCAACACCAACCTGGGCATGGTCAAAACCGACCACAtt                                           >  2:242721/1‑112 (MQ=255)
                                                                                                                                              aaCGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGACCTGCTGCCGCTGATCGAAGGCTGCACCGTCAACACCAAGCTGGGCATGGTCAAGACCGACCACATTCTGTTCATTGCCTCGGGTGCATTCCACCTGAGCAAGcc     >  4:234085/1‑149 (MQ=255)
                                                                                                                                                 gTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGACCTGCTGCCGCTGATCGAAGGCTGCACCGTCAACACCAAGCTGGGCATGGTCAAGACCGACCACATTCTGTTCATTGCCTCGGGTGCATTCCACCTGAGCAAGCCGAg  <  5:97323/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
GAAAGAAGCGCTGAAGATGGTTCGCGATGAAGAAGCGAGCCGCCTGGTCAACGAGGAAGAGCTCAAGGCCAAGGCCCTGGAAGCAGTCGAACAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGACCTGCTGCCGCTGATCGAAGGCTGCACCGTCAACACCAAGCTGGGCATGGTCAAGACCGACCACATTCTGTTCATTGCCTCGGGTGCATTCCACCTGAGCAAGCCGAG  >  NC_002947/5698020‑5698313

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: