Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 461,033 C→G 5.4% V479L (GTG→CTG)  pqqF ← coenzyme PQQ synthesis protein F

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947461,0330CG 5.4% 355.2 / 8.8 111V479L (GTG→CTG) pqqFcoenzyme PQQ synthesis protein F
Reads supporting (aligned to +/- strand):  ref base C (51/54);  new base G (4/2);  total (55/56)
Fisher's exact test for biased strand distribution p-value = 4.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.20e-01

Only 100 of 115 total aligned reads displayed.

GGGCAGGCTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTCACAACGTTCCTGCAGTGGCGCCAGGGCCTGTTCCAGCACACGCTGTAACGGCTGGCGCATAGGTGAAGG  >  NC_002947/460887‑461180
                                                                                                                                                   |                                                                                                                                                    
gggCAGGCTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCActgc                                                                                                                                                   >  2:632373/1‑150 (MQ=255)
      gCTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGCTCCAGGTCCCTGCGTTTTTCCACCATGTCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTgg                                                                                                                                             <  1:1694657/150‑1 (MQ=255)
       cTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGc                                                                                                                                            <  2:2234028/150‑1 (MQ=255)
        tACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCa                                                                                                                                           <  1:1209130/150‑1 (MQ=255)
        tACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCATGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCa                                                                                                                                           >  1:266386/1‑150 (MQ=255)
         aCCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAg                                                                                                                                          <  1:594942/150‑1 (MQ=255)
         aCCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCAGCAGCAGTGCTGCGGGCAg                                                                                                                                          <  1:201544/150‑1 (MQ=255)
           cccGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGc                                                                                                                                        >  1:2168994/1‑150 (MQ=255)
            ccGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGcc                                                                                                                                       >  2:1195752/1‑150 (MQ=255)
            ccGTAAGGCGTCGGGCAGTTGCGTGAGCAGGGCCCTGATGGGGCTCAGGGCCGGGGTCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTAGCTTGAGCGTCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGcc                                                                                                                                       >  2:146818/1‑150 (MQ=255)
                   gcgTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAg                                                                                                                                <  1:1022565/150‑1 (MQ=255)
                   gcgTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAg                                                                                                                                <  2:1562068/150‑1 (MQ=255)
                    cgTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGt                                                                                                                               <  1:1970456/150‑1 (MQ=255)
                        tGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCaa                                                                                                                           <  2:423628/150‑1 (MQ=255)
                            aggtGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCGGGGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTgg                                                                                                                       <  2:929578/147‑2 (MQ=255)
                            aGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTTGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGc                                                                                                                       <  1:2028390/150‑1 (MQ=255)
                            aGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGc                                                                                                                       <  1:236314/150‑1 (MQ=255)
                              ttGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCa                                                                                                                     <  2:1785476/150‑1 (MQ=255)
                              ttGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTCCCACTGCCCGCAGCActgct                                                                                                                                                  <  2:2197638/121‑1 (MQ=255)
                              ttGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTCCCACTGCCCGCAGCActgct                                                                                                                                                  >  1:2197638/1‑121 (MQ=255)
                               tGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAg                                                                                                                    <  1:1624469/150‑1 (MQ=255)
                                  ttGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTAt                                                                                                                 >  1:1368654/1‑150 (MQ=255)
                                      gCAGGGCCCTGAAGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACAGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTccc                                                                                                             <  1:2229909/150‑1 (MQ=255)
                                         gggCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCAAGAGTCTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCgg                                                                                                          >  1:1227650/1‑150 (MQ=255)
                                            cccTGATGGGGATCAGGGCCGGGGGAAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCAc                                                                                                       >  1:427521/1‑150 (MQ=255)
                                              cTGATGGGGATCAGGGCCGGGGGGAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCActg                                                                                                    <  2:38862/150‑1 (MQ=255)
                                              cTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCGGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCActg                                                                                                    >  1:141607/1‑150 (MQ=255)
                                               tGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCActgc                                                                                                   >  1:104321/1‑150 (MQ=255)
                                                  tGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGt                                                                                                >  2:2244134/1‑150 (MQ=255)
                                                   gaggATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCACTGCAACTGCCCGTATTCCCCCGCACTTCTGtc                                                                                               >  2:138555/3‑149 (MQ=255)
                                                    gggATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTAc                                                                                              >  2:1172169/1‑150 (MQ=255)
                                                        tCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCa                                                                                                                     <  2:156668/124‑1 (MQ=255)
                                                        tCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCa                                                                                                                     >  1:156668/1‑124 (MQ=255)
                                                        tCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTACTGCa                                                                                          <  1:1681180/150‑1 (MQ=255)
                                                                   gggCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTACTGCAGTTGCACGGa                                                                                <  1:2104930/150‑1 (MQ=255)
                                                                    ggCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGCTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTACTGCAGTTGCACGGAg                                                                               <  1:1403022/150‑1 (MQ=255)
                                                                                cAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTCACaa                                                                   <  2:266386/150‑1 (MQ=255)
                                                                                cAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTCACaa                                                                   <  1:1156372/150‑1 (MQ=255)
                                                                                   ggCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCAct                                                                                                      <  2:562432/113‑1 (MQ=255)
                                                                                   ggCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCAct                                                                                                      >  1:562432/1‑113 (MQ=255)
                                                                                       gccagcAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTACCGCAGTTGCACGGAGGCGCGGTCACAACGTTCCt                                                            <  2:141607/150‑1 (MQ=255)
                                                                                        ccagcAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTCACAACGTTCCTg                                                           >  1:1873332/1‑150 (MQ=255)
                                                                                        ccagcAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTCACAACGTTCCTg                                                           >  2:518130/1‑150 (MQ=255)
                                                                                        ccagcAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTCACAACGTTCCTg                                                           <  1:518130/150‑1 (MQ=255)
                                                                                           gcAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTCACAACGTTCCt                                                            >  1:469346/1‑146 (MQ=255)
                                                                                           gcAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTCACAACGTTCCt                                                            <  2:469346/146‑1 (MQ=255)
                                                                                              aCTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTACTACAGTTGCACGGAGGCGCGGTCACAACGTTCCTGCAGTgg                                                     >  1:255134/1‑150 (MQ=255)
                                                                                               cTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCTCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTCACAACGTTCCTGCAGTGgc                                                    <  1:2071874/150‑1 (MQ=255)
                                                                                                tGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTCACAACGTTCCTGCAGTgg                                                     <  2:681999/148‑1 (MQ=255)
                                                                                                  gCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTCACAACGTTCCTGCAGTGGCGcc                                                 >  2:1981850/1‑150 (MQ=255)
                                                                                                      gTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTCACAACGTTCCTGCAGTGGCGCCAggg                                             >  2:1141209/1‑150 (MQ=255)
                                                                                                       tGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTCACAACGTTCCTGCAGTGGCGCCAGGGc                                            <  1:316361/150‑1 (MQ=255)
                                                                                                       tGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTCACAACGTTCCTGCAGTGGCGCCAGGGc                                            <  1:1424015/150‑1 (MQ=255)
                                                                                                         gCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTCACAACGTTCCTGCAGTGGCGCCAGGGCCt                                          <  1:947949/150‑1 (MQ=255)
                                                                                                         gCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCACTAACTACGTTCCTGCAGTGGCGCCAGGGCCt                                          <  1:652380/150‑1 (MQ=255)
                                                                                                           ttGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTCACAACGTTCCTGCAGTGGCGCCAGGGCCTGt                                        >  1:789420/1‑150 (MQ=255)
                                                                                                           ttGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTCACAACGTTCCTGCAGTGGCGCCAGGGCCTGt                                        >  1:908239/1‑150 (MQ=255)
                                                                                                           ttGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTCACAACGTTCCTGCAGTGGCGCCAGGGCCTGt                                        >  2:116861/1‑150 (MQ=255)
                                                                                                             gAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTCACAACGTTCCTGCAGTGGCGcc                                                 <  2:2054748/139‑1 (MQ=255)
                                                                                                             gAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTCACAACGTTCCTGCAGTGGCGcc                                                 >  1:2054748/1‑139 (MQ=255)
                                                                                                               gCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTCACAACGTTCCTGCAGTGGCGCCAGGGCCTGTTCCa                                    <  2:1873332/150‑1 (MQ=255)
                                                                                                                 aTCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTCACAACGTTCCTGCAGTGGCGCCAGGGCCTGTTCCAGc                                  <  1:1195752/150‑1 (MQ=255)
                                                                                                                 aTCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTCACAACGTTCCTGCAGTGGCGCCAGGGCCTGTTCCAGc                                  <  1:116861/150‑1 (MQ=255)
                                                                                                                  tCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTCACAACGTTCCTGCAGTGGCGCCAg                                               <  2:912582/136‑1 (MQ=255)
                                                                                                                  tCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTCACAACGTTCCTGCAGTGGCGCCAg                                               >  1:912582/1‑136 (MQ=255)
                                                                                                                   caaggccaaggcTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTCACAACGTTCCTGCAGTGGCGCCAGGGCCTGTTCCAGcac                                >  1:1544589/1‑150 (MQ=255)
                                                                                                                    aaggccaaggcTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTCACAACGTTCCTGCAGTGGCGCCAGGGCCTGTTCCAGcaca                               >  1:1629206/1‑150 (MQ=255)
                                                                                                                     aggccaagtcTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTCACAACGTTCCTGCAGTGGCGCCAGGGCCTGTTCCAGCacac                              >  1:98313/1‑150 (MQ=255)
                                                                                                                     aggccaaggcTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTCACAACGTTCCTGCAGTGGCGCCAGGGCCTGTTCCAGCacac                              >  2:1829822/1‑150 (MQ=255)
                                                                                                                       gccaaggcTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTCACAACGTTCCTg                                                           <  2:1161561/119‑1 (MQ=255)
                                                                                                                       gccaaggcTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTCACAACGTTCCTg                                                           >  1:1161561/1‑119 (MQ=255)
                                                                                                                        ccaaggcTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTACTGCAGTTGCACgaag                                                                               >  1:172305/1‑95 (MQ=255)
                                                                                                                        ccaaggcTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTCACAACGTTCCTGCAGTGGCGCCAGGGCCTGTTCCGGCACACGCt                           >  1:394179/1‑150 (MQ=255)
                                                                                                                        ccaaggcTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTCACAACGTTCCTGCAGTGGCGCCAGGGCCTGTTCCAGCACACGCt                           >  2:1206042/1‑150 (MQ=255)
                                                                                                                        ccaaggcTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTACTGCAGTTGCACGAAg                                                                               <  2:172305/98‑1 (MQ=255)
                                                                                                                          aaggcTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTCACAACGTTCCTGCAg                                                        >  2:472862/1‑119 (MQ=255)
                                                                                                                           aggcTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTCACAACGTTCCTGCAGTGGCGCCAGGGCCTGTTCCAGCACACGCTGTa                        <  1:1141209/150‑1 (MQ=255)
                                                                                                                            ggcTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTccc                                                                                                             <  2:2057824/65‑1 (MQ=255)
                                                                                                                            ggcTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTccc                                                                                                             >  1:2057824/1‑65 (MQ=255)
                                                                                                                            ggcTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTCACAAc                                                                  >  1:2209210/1‑107 (MQ=255)
                                                                                                                            ggcTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTCACAAc                                                                  <  2:2209210/107‑1 (MQ=255)
                                                                                                                                tGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCACTCGCTGACTGCAGTTGCACGGAGGCGCGGTCACAACGTTCCTGCAGTGGCGCCAGGGCCTGTTCCAGCACACGCTGTAACGGc                   <  1:632373/150‑1 (MQ=255)
                                                                                                                                  ttCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTCACAACGTTCCTGCAGTGGCGCCAGGGCCTGTTCCAGCACACGCTGTAACGGCTg                 >  1:1136284/1‑150 (MQ=255)
                                                                                                                                   tCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTCACAACGTTCCTGCAGTGGCGCCAGGGCCTGTTCCAGCACACGCTGTAACGGCTgg                >  1:1457265/1‑150 (MQ=255)
                                                                                                                                   tCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTCACAACGTTCCTGCAGTGGCGCCAGGGCCTGTTCCAGCACACGCTGTAACGGCTgg                <  1:335864/150‑1 (MQ=255)
                                                                                                                                     cACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTACTGCAGTTGCAc                                                                                   <  1:1760000/81‑1 (MQ=255)
                                                                                                                                        tGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTCACAACGTTCCTGCAGTGGCGCCAGGGCCTGTTc                                      <  1:578140/123‑1 (MQ=255)
                                                                                                                                        tGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTCACAACGTTCCTGCAGTGGCGCCAGGGCCTGTTCCAGCACACGCTGTAACGGCTGGCGCAt           >  2:7293/1‑150 (MQ=255)
                                                                                                                                        tGCCAGCAGCAGTGCTGCGGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTCACAACGTTCCTGCAGTGGCGCCAGGGCCTGTTCCAGCACACGCTGTAACGGCTGGCGCAt           >  1:1300533/1‑150 (MQ=255)
                                                                                                                                          cccGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTCACAACGTTCCTGCAGTGGCGCCAGGGCCTGTTCCAGCACACGCTGTAACGGCTGGCGCATAg         >  1:492491/1‑150 (MQ=255)
                                                                                                                                          cccGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTCACAACGTTCCTGCAGTGGCGCCAGGGCCTGTTCCAGCACACGCTGTAACGGCTGGCGCATAg         >  1:1473093/1‑150 (MQ=255)
                                                                                                                                          cccGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTCACAACGTTCCTGCAGTGGCGCCAGGGCCTGTTCCAGCACACGCTGTAACGGCTGGCGCATAg         >  1:2276611/1‑150 (MQ=255)
                                                                                                                                          ccagcagcaGTGCTGCGGGCAGTCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTCACAACGTTCCTGCa                                                         >  1:698909/4‑102 (MQ=255)
                                                                                                                                          ccagcagcaGTGCTGCGGGCAGTCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTCACAACGTTCCTGCa                                                         <  2:698909/99‑1 (MQ=255)
                                                                                                                                           cagcagcaGTGCTGCGGGCAGTCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTCACAACGTTCCTGCAGTGGCGCCAGGGCCTGTTCCAGCACACGCTGTAACGGCTGGCGCATAgg        >  2:2310496/3‑150 (MQ=255)
                                                                                                                                            cGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTCACAACGTTCCTGCAGTGGCGCCAGGGCCTGTTCCAGCACACGCTGTAACGGCTGGCGCATAGGt       <  2:427521/150‑1 (MQ=255)
                                                                                                                                             gcagcaCTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTACt                                                                                             <  2:1414610/63‑1 (MQ=255)
                                                                                                                                                ggaCTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTCACAACGTTCCTGCAGTGGCGCCAGGGCCTGTTCCAGCACACGCTGTAACGGCTGGCGCATAGGTGAAg   <  2:227190/148‑1 (MQ=255)
                                                                                                                                                 caCTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTCACAACGTTCCTGCAGTGGCGCCAGGGCCTGTTCCAGCACACGCTGTAACGGCTGGCGCATAGGTGAAgg  >  2:1264296/1‑150 (MQ=255)
                                                                                                                                                  aCTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTCACAACGTTCCTGCAGTGgc                                                    >  1:1115654/1‑99 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GGGCAGGCTACCCCGTACCGCGTCTGGCAGTTGCTTGAGCAGGGCCCTGATGGGGATCAGGGCCGGGGGCAGCGCCGTGCAGGGCAGCCAGCAACTGGCAGGTGGCTTGAGCATCAAGGCCAAGGCTTGTTCCACTGCCCGCAGCACTGCTGCTGGCAGGCCGGCGCAGTGCAACTGCCAGTATTCCCCGGCACTGCTGTACTGCAGTTGCACGGAGGCGCGGTCACAACGTTCCTGCAGTGGCGCCAGGGCCTGTTCCAGCACACGCTGTAACGGCTGGCGCATAGGTGAAGG  >  NC_002947/460887‑461180

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