Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,818,534 A→G 5.7% N212D (AAC→GAC)  PP_1622 → M23B subfamily metallopeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,818,5340AG 5.7% 289.4 / 4.0 95N212D (AAC→GAC) PP_1622M23B subfamily metallopeptidase
Reads supporting (aligned to +/- strand):  ref base A (40/49);  new base G (2/4);  total (42/53)
Fisher's exact test for biased strand distribution p-value = 6.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATAAAGGCATTGATATCGCCGGTGATTTGGGACAGCCTGTTTTTGCTGCGTCTGATGGTGCGGTGGTTTATGCCGGGAGTGGTTTGAGGGGCTACGGCGAGCTGATCATCATCAAGCACAGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCGATTGCTGAAATGGGGTCTACGGGCACAGATCGGGTGAAGCTGCATTTCGAGATTCGCCGCCAGGGCAAACCCGTCGATCCACTCCAGTTCCTGC  >  NC_002947/1818385‑1818678
                                                                                                                                                     |                                                                                                                                                 
aTAAAGGCATTGATATCGCCGGTGATTTGGGACAGCCTGTTTTTGCTGCGTCTGATGGTGCGGTGGTTTATGCCGGGAGTGGTTTGAGGGGCTACGGCGAGCTGATCATCATCAAGCACAGCGACACCTACGTCAGTGCCTACGGTCATa                                                                                                                                                   <  1:1895514/150‑1 (MQ=255)
      gCATTGATATCGCCGGTGATTTGGGACAGCCTGTTTTTGCTGCGTCTGATGGTGCGGTGGTTTATGCCGGGAGTGGTTTGAGGGGCTACGGCGAGCTGATCATCATCAAGCACAGCGACACCTACGTCAGTGCCTACGGTCATAACCGCa                                                                                                                                             >  2:340639/1‑150 (MQ=255)
        aTTGATATCGCCGGTGATTTGGGACAGCCTGTTTTTGCTGCGTCTGATGGTGCGGTGGTTTATGCCGGGAGTGGTTTGAGGGGCTACGGCGAGCTGATCATCATCAAGCACAGCGACACCTACGTCAGTGCCTGCGGTTATGACCGTAgg                                                                                                                                           <  2:1652331/150‑1 (MQ=255)
         ttGATATCGCCGGTGATTTGGGACAGCCTGTTTTTGCTGCGTCTGATGGTGCGGTGGTTTATGCCGGGAGTGGTTTGAGGGGCTACGGCGAGCTGATCATCATCAAGCACAGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGc                                                                                                                                          <  2:1183536/150‑1 (MQ=255)
            atatCGCCGGTGATTTGGGACAGCCTGTTTTTGCTGCGTCTGATGGTGCGGTGGTTTATGCCGGGAGTGGTTTGAGGGGCTACGGCGAGCTGATCATCATCAAGCACAGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGt                                                                                                                                       >  2:232935/1‑150 (MQ=255)
                  ccGGTGATTTGGGACAGCCTGTTTTTGCTGCGTCTGATGGTGCGGTGGTTTATGCCGGGAGTGGTTTGAGGGGCTACGGCGAGCTGATCATCATCAAGCACAGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTc                                                                                                                                 >  2:1583548/1‑150 (MQ=255)
                  ccGGTGATTTGGGACAGCCTGTTTTTGCTGCGTCTGATGGTGCGGTGGTTTATGCCGGGAGTGGTTTGAGGGGCTACGGCGAGCTGATCATCATCAAGCACAGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTc                                                                                                                                 <  1:1224857/150‑1 (MQ=255)
                   cGGTGATTTGGGACAGCCTGTTTTTGCTGCGTCTGATGGTGCGGTGGTTTATGCCGGGAGTGGTTTGAGGGGCTACGGCGAGCTGATCATCATCAAGCACAGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCg                                                                                                                                >  1:147409/1‑150 (MQ=255)
                    ggTGATTTGGGACAGCCTGTTTTTGCTGCGTCTGATGGTGCGGTGGTTTATGCCGGGAGTGGTTTGAGGGGCTACGGCGAGCTGATCATCATCAAGCACAGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCgg                                                                                                                               <  2:487117/150‑1 (MQ=255)
                    ggTGATTTGGGACAGCCTGTTTTTGCTGCGTCTGATGGTGCGGTGGTTTATGCCGGGAGTGGTTTGAGGGGCTACGGCGAGCTGATCATCATCAAGCACAGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCgg                                                                                                                               <  1:1774363/150‑1 (MQ=255)
                     gTGATTTGGGACAGCCTGTTTTTGCTGCGTCTGATGGTGCGGTGGTTTATGCCGGGAGTGGTTTGAGGGGCTACGGCGAGCTGATCATCATCAAGCACAGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCggg                                                                                                                              <  1:1082850/150‑1 (MQ=255)
                      tGATTTGGGACAGCCTGTTTTTGCTGCGTCTGATGGTGCGGTGGTTTATGCCGGGAGTGGTTTGAGGGGCTACGGCGAGCTGATCATCATCAAGCACAGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGa                                                                                                                             <  2:147409/150‑1 (MQ=255)
                      tGATTTGGGACAGCCTGTTTTTGCTGCGTCTGATGGTGCGGTGGTTTATGCCGGGAGTGGTTTGAGGGGCTACGGCGAGCTGATCATCATCAAGCACAGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGa                                                                                                                             <  2:67790/150‑1 (MQ=255)
                       gATTTGGGACAGCCTGTTTTTGCTGCGTCTGATGGTGCGGTGGTTTATGCCGGGAGTGGTTTGAGGGGCTACGGCGAGCTGATCATCATCAAGCACAGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAg                                                                                                                            <  2:11682/150‑1 (MQ=255)
                        aTTTGGGACAGCCTGTTTTTGCTGCGTCTGATGGTGCGGTGGTTTATGCCGGGAGTGGTTTGAGGGGCTACGGCTAGCTGATCATCATCAAGCACAGCGACACCTACGTCAGTGCTCACGGTCATAACCGCAGGCTGTTGGTTCGGGAgg                                                                                                                           <  1:1590336/150‑1 (MQ=255)
                          ttGGGACAGCCTGTTTTTGCTGCGTCTGATGGTGCGGTGGTTTATGCCGGGAGTGGTTTGAGGGGCTACGGCGAGCTGATCATCATCAAGCACAGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAgggg                                                                                                                         <  2:1527763/150‑1 (MQ=255)
                          ttGGGACAGCCTGTTTTTGCTGCGTCTGATGGTGCGGTGGTTTATGCCGGGAGTGGTTTGAGGGGCTACGGCGAGCTGATCATCATCAAGCACAGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAgggg                                                                                                                         >  2:1302022/1‑150 (MQ=255)
                                  gCCTGTTTTTGCTGCGTCTGATGGTGCGGTGGTTTATGCCGGGAGTGGTTTGAGGGGCTACGGCGAGCTGATCATCATCAAGCACAGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGcagcagtt                                                                                                                <  1:1013120/150‑3 (MQ=255)
                                   ccTGTTTTTGCTGCGTCTGATGGTGCGGTGGTTTATGCCGGGAGTGGTTTGAGGGGCTACGGCGAGCTGATCATCATCAAGCACAGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTc                                                                                                               >  1:1292193/1‑150 (MQ=255)
                                      gTTTTTGCTGCGTCTGATGGTGCGGTGGTTTATGCCGGGAGTGGTTTGAGGGGCTACGGCGAGCTGATCATCATCAAGCACAGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAg                                                                                                            <  1:340639/150‑1 (MQ=255)
                                      gTTTTTGCTGCGTCTGATGGTGCGGTGGTTTATGCCGGGAGTGGTTTGAGGGGCTACGGCGAGCTGATCATCATCAAGCACAGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAg                                                                                                            <  1:232935/150‑1 (MQ=255)
                                          ttGCTGCGTCTGATGGTGCGGTGGTTCATGCCGGGAGTGGTTTGAGGGGCTACGGCGAGCTGATCATCATCAAGCACAGCGACACCTACGTCAGTGCATACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAggcag                                                                                                        <  2:310615/150‑1 (MQ=255)
                                           tgctgcGTCTGATGGTGCGGTGGTTTATGCCGGGAGTGGTTTGAGGGGCTACGGCGAGCTGATCATCATCAAGCACAGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAggcagg                                                                                                       <  2:1231016/150‑1 (MQ=255)
                                            gctgcGTCTGATGGTGCGGTGGTTTATGCCGGGAGTGGTTTGAGGGGCTACGGCGAGCTGATCATCATCAAGCACAGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAggcaggg                                                                                                      <  1:1302022/150‑1 (MQ=255)
                                              tgcGTCTGATGGTGCGGTGGTTTATGCCGGGAGTGGTTTGAGGGGCTACGGCGAGCTGATCATCATCAAGCACAGCGACACCTACGTCAGTGCCTGCGGTTATGACCgta                                                                                                                                             >  1:369785/1‑108 (MQ=255)
                                              tgcGTCTGATGGTGCGGTGGTTTATGCCGGGAGTGGTTTGAGGGGCTACGGCGAGCTGATCATCATCAAGCACAGCGACACCTACGTCAGTGCCTGCGGTTATGACCgta                                                                                                                                             <  2:369785/110‑3 (MQ=255)
                                              tgcGTCTGATGGTGCGGTGGTTTATGCCGGGAGTGGTTTGAGGGGCTACGGCGAGCTGATCATCATCAAGCACAGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAggcagggca                                                                                                    <  2:617588/150‑1 (MQ=255)
                                                cGTCTGATGGTGCGGTGGTTTATGCCGGGAGTGGTTTGAGGGGCTACGGCGAGCTGATCATCATCAAGCACAGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGt                                                                                                  <  2:890626/150‑1 (MQ=255)
                                                 gTCTGATGGTGCGGTGGTTTATGCCGGGAGTGGTTTGAGGGGCTACGGCGAGCTGATCATCATCAAGCACAGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTc                                                                                                 <  2:1725715/150‑1 (MQ=255)
                                                  tCTGATGGTGCGGTGGTTTATGCCGGGAGTGGTTTGAGGGGCTACGGCGAGCTGATCATCATCAAGCACAGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCg                                                                                                >  1:1713618/1‑150 (MQ=255)
                                                    tGATGGTGCGGTGGCTTATGCCGGGAGTGGTTTGAGGGGCAAGGGCGAGCTGATCAAAGTGAATCACAGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCGAt                                                                                              >  2:1005905/1‑150 (MQ=255)
                                                     gATGGTGCGGTGGTTTATGCCGGGAGTGGTTTGAGGGGCTACGGCGAGCTGATCATCATCAAGCACAGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCGAtt                                                                                             <  2:1292193/150‑1 (MQ=255)
                                                      aTGGTGCGGTGGTTTATGCCGGGAGTGGTTTGAGGGGCTACGGCGAGCTGATCATCATCAAGCACAGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAggg                                                                                                                          >  1:1250665/1‑121 (MQ=255)
                                                      aTGGTGCGGTGGTTTATGCCGGGAGTGGTTTGAGGGGCTACGGCGAGCTGATCATCATCAAGCACAGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAggg                                                                                                                          <  2:1250665/121‑1 (MQ=255)
                                                       tGGTGCGGTGGTTTATGCCGGGAGTGGTTTGAGGGGCTACGGCGAGCTGATCATCATCAAGCACAGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCGATTGc                                                                                           <  1:1513381/150‑1 (MQ=255)
                                                             ggtggtTTATGCCGGGAGAGGATTGAGGGGCTACGGCGAGCTGATCATCATCAAGCACAGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCGATTGCTGAAAt                                                                                     >  2:1664945/1‑150 (MQ=255)
                                                              gtggtTTATGCCGGGAGTGGTTTGAGGGGCTACGGCGAGCTGATCATCATCAAGCACAGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCGATTGCTGANATg                                                                                    <  2:1188394/150‑1 (MQ=255)
                                                                   ttATGCCGGGAGTGGTTTGAGGGGCTACGGCGAGCTGATCATCATCAAGCACAGCGACACCTACGTCAGTGCCTGCGGTTATGACCGTAgg                                                                                                                                           <  1:161942/91‑1 (MQ=255)
                                                                   ttATGCCGGGAGTGGTTTGAGGGGCTACGGCGAGCTGATCATCATCAAGCACAGCGACACCTACGTCAGTGCCTGCGGTTATGACCGTAgg                                                                                                                                           >  2:161942/1‑91 (MQ=255)
                                                                    tATGCCGGGAGTGGTTTGAGGGGCTACGGCGAGCTGATCATCATCAAGCACAGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCGATTGCTGAAATGGGGTCt                                                                              <  1:467284/150‑1 (MQ=255)
                                                                     aTGCCGGGAGTGGTTTGAGGGGCTACGGCGAGCTGATCATCATCAAGCACAGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCGATTGCTGAAATGGGGTCTa                                                                             <  1:629757/150‑1 (MQ=255)
                                                                        ccGGGAGTGGTTTGAGGGGCTACGGCGAGCTGATCATCATCAAGCACAGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAg                                                                                                                                            <  2:1045375/85‑1 (MQ=255)
                                                                        ccGGGAGTGGTTTGAGGGGCTACGGCGAGCTGATCATCATCAAGCACAGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAg                                                                                                                                            >  1:1045375/1‑85 (MQ=255)
                                                                        ccGGGAGTGGTTTGAGGGGCTACGGCGAGCTGATCATCATCAAGCACAGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCGATTGCTGAAATGGGGTCTACgg                                                                          <  1:914369/150‑1 (MQ=255)
                                                                        ccGGGAGTGGTTTGAGGGGCTACGGCGAGCTGATCATCATCAAGCACAGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCGATTGCTGAAATGGGGTCTACgg                                                                          <  2:14725/150‑1 (MQ=255)
                                                                         cGGGAGTGGTTTGAGGGGCTACGGCGAGCTGATCATCATCAAGCACAGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCGATTGCTGAAATgg                                                                                   >  1:609742/1‑140 (MQ=255)
                                                                         cGGGAGTGGTTTGAGGGGCTACGGCGAGCTGATCATCATCAAGCACAGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCGATTGCTGAAATgg                                                                                   <  2:609742/140‑1 (MQ=255)
                                                                              gTGGTTTGAGGGGCTACGGCGAGCTGATCATCATCAAGCACAGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCCGCATGTCAAGGCAGGGCAGTCGATTGCTGAAATGGGGTCTACGGGCACAg                                                                    >  1:115499/1‑150 (MQ=255)
                                                                                 gTTTGAGGGGCTACGGCGAGCTGATCATCATCAAGCACAGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCGAGGCAGGGCAGTCGATTGCTGAAATGGGGTCTACGGGCACAGATc                                                                 >  2:1792136/1‑150 (MQ=255)
                                                                                 gTTTGAGGGGCTACGACGAGCTGATCATCCTCAAGCACCGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAAGCCGGGGGTTCGGGAGCGGCAGCAGCTCACGGCAGGGCCGTCCATTGCTGAACTGGGGTCTACCGGCACAGacc                                                                 >  2:630125/1‑148 (MQ=255)
                                                                                  tttGAGGGGCTACGGCGAGCTGATCATCATCAAGCACAGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCTATTGCTGAAATGGGGTCTACGGGCACAGATCg                                                                >  1:1844684/1‑150 (MQ=255)
                                                                                  tttGAGGGGCTACGGCGAGCTGATCATCATCAAGCACAGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCGATTGCTGAAATGGGGTCTACGGGCACAGATCg                                                                >  2:1009241/1‑150 (MQ=255)
                                                                                   ttGAGGGGCTACGGCGAGCTGATCATCATCAAGCACAGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCGATTGCTGAAATGGGGTCTACGGGCACAGATCgg                                                               >  2:179957/1‑150 (MQ=255)
                                                                                      aGGGGCTACGGCGAGCTGATCATCATCAAGCACAGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTGGGGCGGGGAGTGCCGGTTATGCAGGGCAGTCGATTGCTGAAATGGGGTCTACGGGCACAGATCGtgtg                                                            >  1:697343/1‑150 (MQ=255)
                                                                                          gCTACGGCGAGCTGATCATCATCAAGCACAGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAggg                                                                                                                          <  1:1300807/85‑1 (MQ=255)
                                                                                          gCTACGGCGAGCTGATCATCATCAAGCACAGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAggg                                                                                                                          >  2:1300807/1‑85 (MQ=255)
                                                                                           cTACGGCGAGCTGATCATCATCAAGCACAGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCGATTGCTGAAATGGGGTCTACGGGCACAGATCGGGTGAAGCt                                                       >  2:92445/1‑150 (MQ=255)
                                                                                           cTACGGCGAGCTGATCATCATCAAGCACAGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCGATTGCTGAAATGGGGTCTACGGGCACAGATCGGGTGAAGCt                                                       >  1:907434/1‑150 (MQ=255)
                                                                                             aCGGCGAGCTGATCATCATCAAGCACAGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCGATTGCTGAAATGGGGTCTACGGGCACAGATCGGGTGAAGCTGc                                                     >  1:3179/1‑150 (MQ=255)
                                                                                                 cGAGCTGATCATCATCAAGCACAGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCGATTGCTGAAATGGGGTCTACGGGCACAGATCGGGTGAAGCTGCAttt                                                 <  2:197912/150‑1 (MQ=255)
                                                                                                        atcatcatcAAGCACAGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCGATTGCTGAAATGGGGTCTACGGGCACAGATCGGGTGAAGCTGCATTTCgaga                                            <  2:1420710/148‑1 (MQ=255)
                                                                                                        atcatcatcAAGCACAGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCGATTGCTGAAATGGGGTCTACGGGCACAGATCGGGTGAAGCTGCATTTCgaga                                            >  1:1420710/1‑148 (MQ=255)
                                                                                                        atcatcatcAAGCACAGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCGATTGCTGAAATGGGGTCTACGGGCACAGATCGGGTGAAGCTGCATTTCGAGAtt                                          >  1:1699381/1‑150 (MQ=255)
                                                                                                         tcatcatcaAGCACAGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCGATTGCTGAAATGGGGTCTACGGGCACAGATCGGGTGAAGCTGCATTTCGAGATTc                                         <  2:819789/150‑1 (MQ=255)
                                                                                                            tcatcaAGCACAGCGACACCTACGTCAGTGCCTAAGGCAATAACCACAAGCGGTTGGTTCTGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCGATTGCTGAAATGGGGTCTACGGGCACAGATCGGGTGAAGCTGCATTTCGAGATTcgcc                                      >  1:1775505/1‑150 (MQ=255)
                                                                                                               tcaAGCACAGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTg                                                                                                                                     >  2:344451/1‑53 (MQ=255)
                                                                                                               tcaAGCACAGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTg                                                                                                                                     <  1:344451/53‑1 (MQ=255)
                                                                                                               tcaAGCACAGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCGATTGCTGAAATGGGGTCTACGGGCACAGATCGGGTGAAGCTGCATTTCGAGATTcg                                        >  2:85757/1‑145 (MQ=255)
                                                                                                               tcaAGCACAGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCGATTGCTGAAATGGGGTCTACGGGCACAGATCGGGTGAAGCTGCATTTCGAGATTcg                                        <  1:85757/145‑1 (MQ=255)
                                                                                                               tcaAGCACAGAGACACCTACATCAGTAACTACGGTCATAAAAGCAGTCTGTAGTGGCAGGGGTGCGAGCTTATCTTTGGATGGCAGTCTAGTTCTGAAATGGGGTCTACGGGCACAGATCGGGTGAAGCTGCATTTCGAGATTCgccgcc                                   >  2:1163221/1‑150 (MQ=255)
                                                                                                                caAGCACAGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCGATTGCTGAAATGGGGTCTACGGGCACAGATCGGGTGAAGCTGCATTTCGAGATTCGCCGCCa                                  <  1:1047720/150‑1 (MQ=255)
                                                                                                                 aaGCACAGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCGATTGCTGAAATGGGGTCTACGGGCACAGATCGGGTGAAGCTGCATTTCGAGATTCGCCGCCAg                                 >  2:1921714/1‑150 (MQ=255)
                                                                                                                   tcacaGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCGATTGCTGAAATGGGGTCTACGGGCACAGATCGGGTGAAGCTGCATTTCGAGATTCGCCGCCAggg                               <  1:240507/149‑1 (MQ=255)
                                                                                                                   gCACAGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCGATTGCTGAAATGGGGTCTACGGGCACAGATCGGGTGAAGCTGCATTTCGAGATTCGCCGCCAggg                               <  1:92445/150‑1 (MQ=255)
                                                                                                                   gCACAGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCGATTGCTGAAATGGGGTCTACGGGCACAGATCGGGTGAAGCTGCATTTCGAGATTCGCCGCCAggg                               <  2:1699381/150‑1 (MQ=255)
                                                                                                                     acaGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCGATTGCTGAAATGGGGTCTACGGGCACAGATCGGGTGACGCTTCTTTTCGAGATTCGCCGCCAGGGCa                             >  2:240507/1‑150 (MQ=255)
                                                                                                                      caGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCGATTGCTGAAATGGGGTCTACGGGCACAGATCGGGTGAAGCTGCATTTCGAGATTCGCCGCCAGGGCaa                            <  1:1321795/150‑1 (MQ=255)
                                                                                                                               cTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAg                                                                                                            <  1:1690153/61‑1 (MQ=255)
                                                                                                                               cTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAg                                                                                                            >  2:1690153/1‑61 (MQ=255)
                                                                                                                               cTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCGATTGCTGAAATGGGGTCTACGGGCACAGATCGGGTGAAGCTGCATTTCGAGATTCGCCGCCAGGGCAAACCCGTCGa                   <  2:900488/150‑1 (MQ=255)
                                                                                                                                    tCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCGATTGCTGAAATGGGGTCTACGGGCACAGATCGGGTGAAGCTGCATTTCGAGATTCGCCGCCAGGGCAAAcc                         >  1:237267/1‑139 (MQ=255)
                                                                                                                                    tCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCGATTGCTGAAATGGGGTCTACGGGCACAGATCGGGTGAAGCTGCATTTCGAGATTCGCCGCCAGGGCAAAcc                         <  2:237267/139‑1 (MQ=255)
                                                                                                                                    tCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCGATTGCTGAAATGGGGTCTACGGGCACAGATCGGGTGAAGCTGCATTTCGAGATTCGCCGCCAGGGCAAACCCGTCGATCCAc              >  1:487152/1‑150 (MQ=255)
                                                                                                                                    tCAGGGCCTACGGTCATGACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCGATTGCTGAAATGGGGTCTACGGGCACAGATCGGGTGAAGCTGCATTTCGAGATTCGCCGCCAGGTCAAACCCGTCGATCCAc              <  2:645351/150‑1 (MQ=255)
                                                                                                                                       gTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAggcag                                                                                                        <  1:622642/57‑1 (MQ=255)
                                                                                                                                       gTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAggcag                                                                                                        >  2:622642/1‑57 (MQ=255)
                                                                                                                                        tGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCTATTGCTGAAATGGGGTCTACGGGCACAGATCGGGTGAAGCTGCATTTCGAGATTCGCCGCCAGGGCAAACCCGTCGATCCACTCCa          <  2:1844684/150‑1 (MQ=255)
                                                                                                                                        tGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCGATTGCTGAAATGGGGTCTACGGGCACAGATCGGGTGAAGCTGCATTTCga                                              <  2:412278/114‑1 (MQ=255)
                                                                                                                                        tGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCGATTGCTGAAATGGGGTCTACGGGCACAGATCGGGTGAAGCTGCATTTCga                                              >  1:412278/1‑114 (MQ=255)
                                                                                                                                           cTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCGATTGCTGAAATGGGGTCTACGCGCACAGATCGGGTGAAGCTACATTTCGAGATTCGCCGCCAGGGCAAACCCGTCGATCCACTCCAGtt       >  2:124428/1‑150 (MQ=255)
                                                                                                                                            tACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCGATTg                                                                                            >  2:449588/1‑64 (MQ=255)
                                                                                                                                            tACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCGATTg                                                                                            <  1:449588/64‑1 (MQ=255)
                                                                                                                                            tACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGGCGATTGCTGAAATGGGGTCTACGGGCACAGATCGGGTGAAGCTGCATTTCGAGATTCGCCGCCAGGGCAAACCCGTCGATCCACTCCAGTTc      >  1:516911/1‑150 (MQ=255)
                                                                                                                                                gTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCGATTGCTGAAATGGGGTCTACGGGCACAGATCGGGTGAAGCTGCATTTCGAGATTcgcc                                      >  1:526558/1‑114 (MQ=255)
                                                                                                                                                gTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCGATTGCTGAAATGGGGTCTACGGGCACAGATCGGGTGAAGCTGCATTTCGAGATTcgcc                                      <  2:526558/114‑1 (MQ=255)
                                                                                                                                                gTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCGATTGCTGAAATGGGGTCTACGGGCACAGATCGGGTGAAGCTGCATTTCGAGATTCGCCGCCAGGGCAAACCCGTCGATCCACTCCAGTTCCTGc  <  1:179957/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                 
ATAAAGGCATTGATATCGCCGGTGATTTGGGACAGCCTGTTTTTGCTGCGTCTGATGGTGCGGTGGTTTATGCCGGGAGTGGTTTGAGGGGCTACGGCGAGCTGATCATCATCAAGCACAGCGACACCTACGTCAGTGCCTACGGTCATAACCGCAGGCTGTTGGTTCGGGAGGGGCAGCAGGTCAAGGCAGGGCAGTCGATTGCTGAAATGGGGTCTACGGGCACAGATCGGGTGAAGCTGCATTTCGAGATTCGCCGCCAGGGCAAACCCGTCGATCCACTCCAGTTCCTGC  >  NC_002947/1818385‑1818678

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: