Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,698,157 G→T 7.9% R264L (CGT→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,1570GT 7.9% 234.3 / 8.3 76R264L (CGT→CTT) hslUprotease HslVU ATPase subunit
Reads supporting (aligned to +/- strand):  ref base G (34/36);  new base T (3/3);  total (37/39)
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.98e-01

GCTGAAGGTGAAAGAAGCGCTGAAGATGGTTCGCGATGAAGAAGCGAGCCGCCTGGTCAACGAGGAAGAGCTCAAGGCCAAGGCCCTGGAAGCAGTCGAACAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGACCTGCTGCCGCTGATCGAAGGCTGCACCGTCAACACCAAGCTGGGCATGGTCAAGACCGACCACATTCTGTTCATTGCCTCGGGTGCATTCCACCT  >  NC_002947/5698011‑5698301
                                                                                                                                                  |                                                                                                                                                
gCTGAAGGTGAAAGAAGCGCTGAAGATGGTTCGCGATGAAGAAGCGAGCCGCCTGGTCAACGAGGAAGAGCTCAAGGCCAAGGCCCTGGAAGCAGTCGAACAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTg                                                                                                                                                <  4:162400/149‑1 (MQ=255)
gCTGAAGGTGAAAGAAGCGCTGAAGATGGTTCGCGATGAAGAAGCGAGCCGCCTGGTCAACGAGGAAGAGCTCAAGGCCAAGGCCCTGGAAGCAGTCGAACAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTg                                                                                                                                                >  2:331994/1‑149 (MQ=255)
     aGGTGAAAGAAGCGCTGAAGATGGTTCGCGATGAAGAAGCGAGCCGCCTGGTCAACGAGGAAGAGCTCAAGGCCAAGGCCCTGGAAGCAGTCGAACAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAAc                                                                                                                                           <  2:220279/149‑1 (MQ=255)
     aGGTGAAAGAAGCGCTGAAGATGGTTCGCGATGAAGAAGCGAGCCGCCTGGTCAACGAGGAAGAGCTCAAGGCCAAGGCCCTGGAAGCAGTCGAACAGCACGATATCGTGTTCATCGACGAAATCGACAACGTTGCCACGCTTGGCAAc                                                                                                                                           <  1:78569/149‑1 (MQ=255)
         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)
                              tCGCGATGAAGAAGCGAGCCGCCTGGTCAACGAGGAAGAGCTCAAGGCCAAGGCCCTGGAAGCAGTCGAACAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTg                                                                                                                                                <  6:182093/119‑1 (MQ=255)
                                 cgATGAAGAAGCGAGCCGCCTGGTCAACGAGGAAGAGCTCAAGGCCAAGGCCCTGGAAGCAGTCGAACAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAg                                                                                                               >  3:51896/1‑149 (MQ=255)
                                 cgATGAAGAAGCGAGCCGCCTGGTCAACGAGGAAGAGCTCAAGGCCAAGGCCCTGGAAGCAGTCGAACAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAg                                                                                                               <  1:196375/149‑1 (MQ=255)
                                  gATGAAGAAGCGAGCCGCCTGGTCAACGAGGAAGAGCTCAAGGCCAAGGCCCTGGAAGCAGTCGAACAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGt                                                                                                                                                 <  8:14820/114‑1 (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                                                                                                                           >  8:240354/1‑86 (MQ=255)
                                                                                    ccTGGAAGCAGTCGAACAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACg                                                                                                                           <  7:240354/86‑1 (MQ=255)
                                                                                          aGCAGTCGAACAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGACCTGCTGCCGCTGATCGAAGGCTGCACCGTCAACACCAAGCTgg                                                      >  4:88919/1‑149 (MQ=255)
                                                                                             aGTCGAACAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGACCTGCTGCCGCTGa                                                                                    >  4:332913/1‑116 (MQ=255)
                                                                                             aGTCGAACAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGACCTGCTGCCGCTGa                                                                                    <  3:332913/116‑1 (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                                                                >  1:232590/1‑125 (MQ=255)
                                                                                                        aCGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGACCTGCTGCCGCTGATCGAAGGCTGCACCGTCAac                                                                <  2:232590/125‑1 (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                      <  8:185906/130‑1 (MQ=255)
                                                                                                                                             cAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGACCTGCTGCCGCTGATCGAAGGCTGCACCGTCAACACCAAGCTGGGCATGGTCAAGACCGACCACATTCTGTTCAtt                      >  7:185906/1‑130 (MQ=255)
                                                                                                                                              aaGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGAACAGCGCGACCTGCTGCCGCTGATCGAAGGCTGCACCGTCAACACCAAGCTGGGCATGGTCAAGACCGACCACATTCTGTTCATTGCCTCGGGTGCATTCCACCt  >  1:243843/1‑149 (MQ=255)
                                                                                                                                                gCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGACCTGCTGCCGCTGATCGAAGGCTGCACCGTCAACACCAAGCTGGGCATGGTCAAGACCGACCACATTCTGTTCATTGCCTCGGGTGCATTCCACCt  >  8:228124/1‑147 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
GCTGAAGGTGAAAGAAGCGCTGAAGATGGTTCGCGATGAAGAAGCGAGCCGCCTGGTCAACGAGGAAGAGCTCAAGGCCAAGGCCCTGGAAGCAGTCGAACAGCACGGTATCGTGTTCATCGACGAAATCGACAAGGTTGCCAAGCGTGGCAACGTTGGCGGCGCCGACGTGTCCCGTGAAGGCGTACAGCGCGACCTGCTGCCGCTGATCGAAGGCTGCACCGTCAACACCAAGCTGGGCATGGTCAAGACCGACCACATTCTGTTCATTGCCTCGGGTGCATTCCACCT  >  NC_002947/5698011‑5698301

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: