Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_008463 2,158,823 G→C Y90* (TAC→TAG dacB ← D‑alanyl‑D‑alanine carboxypeptidase/D‑alanyl‑D‑alanine‑endopeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0084632,158,8230GC98.6% 1082.1 / NA 290Y90* (TAC→TAGdacBD‑alanyl‑D‑alanine carboxypeptidase/D‑alanyl‑D‑alanine‑endopeptidase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base C (120/166);  minor base A (3/1);  total (123/167)
Fisher's exact test for biased strand distribution p-value = 3.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 290 total aligned reads displayed.

GGTCACGCCATTGGCGCGCAGGTCACGCATCAACAGCCAGAGCTTCTCCATGTTCAGTTTCGGATCGCCGCCGCCCTTGAGGTAGAGGTTGCCGTTGAGCACTCCGTTCTTGAGCTGGCCGTCGGTGTAGAACTCGGTCTTCCACTGGTAGGTCGGGCCGAGCATTTCCAGGGCGGCATAGGTGGTGAACAGCTTCATGGTCGAGGCCGGGTTCACCGATACGTCGGCGTTGTAGTAGGTGGGGTTGCCCGGGCCGTCGAGCGGAATCAGCACCAGCGACAGATCGTTGCCGGTC  >  NC_008463/2158676‑2158970
                                                                                                                                                   |                                                                                                                                                   
ggTCACGCCATTGGCGCGCAGGTCACGCATCAACAGCCAGAGCTTCTCCATGTTCAGTTTCGGATCGCCGCCGCCCTTGAGGTAGAGGTTGCCGTTGAGCACTCCGTTCTTGAGCTGGCCGTCGGTGTAGAACTCGGTCTTCCACTGCTAg                                                                                                                                                  >  2:4310621/1‑151 (MQ=255)
  tCACGCCATTGGCGCGCAGGTCACGCATCAACAGCCAGAGCTTCTCCATGTTCAGTTTCGGATCGCCGCCGCCCTTGAGGTAGAGGTTGCCGTTGAGCACTCCGTACTTGAGCTGGCCGTCGGTGTAGAACTCGGTCTTCCACTGCTAGGt                                                                                                                                                <  2:3166620/151‑1 (MQ=255)
    aCGCCATTGGCGCGCAGGTCACGCATCAACAGCCAGAGCTTCTCCATGTTCAGTTTCGGATCGCCGCCGCCCTTGAGGTAGAGGTTGCCGTTGAGCACTCCGTTCTTGAGCTGGCCGTCGGTGTAGAACTCGGTCTTCCACTGCTAGGTCg                                                                                                                                              <  1:6707597/151‑1 (MQ=255)
    aCGCCATTGGCGCGCAGGTCACGCATCAACAGCCAGAGCTTCTCCATGTTCAGTTTCGGATCGCCGCCGCCCTTGAGGTAGAGGTTGCCGTTGAGCACTCCGTTCTTGAGCTGGCCGTCGGTGTAGAACTCGGTCTTCCACTGCTAGGTCg                                                                                                                                              >  1:6029954/1‑151 (MQ=255)
       cccttGGCGCGCAGGTCACGCATCAACAGCCAGAGCTTCTCCATGTTCAGTTTCGGATCGCCGCCGCCCTTGAGGTAGAGGTTGCCGTTGAGCACTCCGTTCTTGAGCTGGCCGTCGGTGTAGAACTCGGTCTTCCACTGCTAGGTCGGGc                                                                                                                                           <  2:8365086/148‑1 (MQ=255)
       ccATTGGCGCGCAGGTCACGCATCAACAGCCAGAGCTTCTCCATGTTCAGTTTCGGATCGCCGCCGCCCTTGAGGTAGAGGTTGCCGTTGAGCACTCCGTTCTTGAGCTGGCCGTCGGTGTAGAACTCGGTCTTCCACTGCTAGTTCGGGc                                                                                                                                           >  2:11064687/1‑151 (MQ=255)
       ccATTGGCGCGCAGGTCACGCATCAACAGCCAGAGCTTCTCCATGTTCAGTTTCGGATCGCCGCCGCCCTTGAGGTAGAGGTTGCCGTTGAGCACTCCGTTCTTGAGCTGGCCGTCGGTGTAGAACTCGGTCTTCCACTGCTAGGTCGGGc                                                                                                                                           <  1:1108673/151‑1 (MQ=255)
       ccATTGGCGCGCAGGTCACGCATCAACAGCCAGAGCTTCTCCATGTTCAGTTTCGGATCGCCGCCGCCCTTGAGGTAGAGGTTGCCGTTGAGCACTCCGTTCTTGAGCTGGCCGTCGGTGTAGAACTCGGTCTTCCACTGCTAGGTCGGGc                                                                                                                                           >  2:1688292/1‑151 (MQ=255)
        cATTGGCGCGCAGGTCACGCATCAACAGCCAGAGCTTCTCCATGTTCAGTTTCGGATCGCCGCCGCCCTTGAGGTAGAGGTTGCCGTTGAGCACTCCGTTCTTGAGCTGGCCGTCGGTGTAGAACTCGGTCTTCCACTGCTAGGTCGGGcc                                                                                                                                          >  2:192281/1‑151 (MQ=255)
         aTTGGCGCGCAGGTCACGCATCAACAGCCAGAGCTTCTCCATGTTCAGTTTCGGATCGCCGCCGCCCTTGAGGTAGAGGTTGCCGTTGAGCACTCCGTTCTTGAGCTGGCCGTCGGTGTAGAACTCGGTCTTCCACTGCTAGGTCGGGCCg                                                                                                                                         <  1:7527904/151‑1 (MQ=255)
          ttGGCGCGCAGGTCACGCATCAACAGCCAGAGCTTCTCCATGTTCAGTTTCGGATCGCCGCCGCCCTTGAGGTAGAGGTTGCCGTTGAGCACTCCGTTCTTGAGCTGGCCGTCGGTGTAGAACTCGGTCTTCCACTGCTAGGTCGGGCCGa                                                                                                                                        <  2:10895839/151‑1 (MQ=255)
          ttGGCGCGCAGGTCACGCATCAACAGCCAGAGCTTCTCCATGTTCAGTTTCGGATCGCCGCCGCCCTTGAGGTAGAGGTTGCCGTTGAGCACTCCGTTCTTGAGCTGGCCGTCGGTGTAGAACTCGGTCTTCCACTGCTAGGTCGGGCCGa                                                                                                                                        >  1:4715733/1‑151 (MQ=255)
           tGGCGCGCAGGTCACGCATCAACAGCCAGAGCTTCTCCATGTTCAGTTTCGGATCGCCGCCGCCCTTGAGGTAGAGGTTGCCGTTGAGCACTCCGTTCTTGAGCTGGCCGTCGGTGTAGAACTCGGTCTTCCACTGCTAGGTCGGGCCGAg                                                                                                                                       <  2:4504959/151‑1 (MQ=255)
           tGGCGCGCAGGTCACGCATCAACAGCCAGAGCTTCTCCATGTTCAGTTTCGGATCGCCGCCGCCCTTGAGGTAGAGGTTGCCGTTGAGCACTCCGTTCTTGAGCTGGCCGTCGGTGTAGAACTCGGTCTTCCACTGCTAGGTCGGGCCGAg                                                                                                                                       <  1:1051411/151‑1 (MQ=255)
           tGGCGCGCAGGTCACGCATCAACAGCCAGAGCTTCTCCATGTTCAGTTTCGGATCGCCGCCGCCCTTGAGGTAGAGGTTGCCGTTGAGCACTCCGTTCTTGAGCTGGCCGTCGGTGTAGAACTCGGTCTTCCACTGCTAGGTCGGGCCGAg                                                                                                                                       >  1:330567/1‑151 (MQ=255)
            ggCGCGCAGGTCACGCATCAACAGCCAGAGCTTCTCCATGTTCAGTTTCGGATCGCCGCCGCCCTTGAGGTAGAGGTTGCCGTTGAGCACTCCGTTCTTGAGCTGGCCGTCGGTGTAGAACTCGGTCTTCCACTGCTAGGTCGGGCCGAGc                                                                                                                                      >  1:7851415/1‑151 (MQ=255)
              cgcgcAGGTCACGCATCAACAGCCAGAGCTTCTCCATGTTCAGTTTCGGATCGCCGCCGCCCTTGAGGTAGAGGTTGCCGTTGAGCACTCCGTTCTTGAGCTGGCCGTCGGTGTAGAACTCGGTCTTCCACTGCTAGGTCGGGCCGAGCAt                                                                                                                                    >  1:5114142/1‑151 (MQ=255)
              cgcgcAGGTCACGCATCAACAGCCAGAGCTTCTCCATGTTCAGTTTCGGATCGCCGCCGCCCTTGAGGTAGAGGTTGCCGTTGAGCACTCCGTTCTTGAGCTGGCCGTCGGTGTAGAACTCGGTCTTCCACTGCTAGGTCGGGCCGAGCAt                                                                                                                                    <  2:1369230/151‑1 (MQ=255)
                 gcAGGTCACGCATCAACAGCCAGAGCTTCTCCATGTTCAGTTTCGGATCGCCGCCGCCCTTGAGGTAGAGGTTGCCGTTGAGCACTCCGTTCTTGAGCTGGCCGTCGGTGTAGAACTCGGTCTTCCACTGCTAGGTCGGGCCGAGCATTTc                                                                                                                                 <  2:5222212/151‑1 (MQ=255)
                   aGGTCACGCATCAACAGCCAGAGCTTCTCCATGTTCAGTTTCGGATCGCCGCCGCCCTTGAGGTAGAGGTTGCCGTTGAGCACTCCGTTCTTGAGCTGGCCGTCGGTGTAGAACTCGGTCTTCCACTGCTAGGTCGGGCCGAGCATTTCCa                                                                                                                               <  2:7604723/151‑1 (MQ=255)
                     gTCACGCATCAACAGCCAGAGCTTCTCCATGTTCTGTTTCGGATCGCCGCCGCCCTTGAGGTAGAGGTTGCCGTTGAGCACTCCGTTCTTGAGCTGGCCGTCTGTGTAGAACTCGGTCTTCCACTGCTAGGTCGGGCCGAGCATTTCCAgg                                                                                                                             <  1:9874330/151‑1 (MQ=255)
                          gCATCAACAGCCAGAGCTTCTCCATGTTCAGTTTCGGATCGCCGCCGCCCTTGAGGTAGAGGTTGCCGTTGAGCACTCCGTTCTTGAGCTGGCCGTCGGTGTAGAACTCGGTCTTCCACTGCTAGGTCGGGCCGAGCATTTCCAGggcggc                                                                                                                        >  1:9562888/1‑151 (MQ=255)
                          gCATCAACAGCCAGAGCTTCTCCATGTTCAGTTTCGGATCGCCGCCGCCCTTGAGGTAGAGGTTGCCGTTGAGCACTCCGTTCTTGAGCTGGCCGTCGGTGTAGAACTCGGTCTTCCACTGCTAGGTCGGGCCGAGCATTTCCAGggcggc                                                                                                                        >  1:3617159/1‑151 (MQ=255)
                            aTCAACAGCCAGAGCTTCTCCATGTTCAGTTTCGGATCGCCGCCGCCCTTGAGGTAGAGGTTGCCGTTGAGCACTCCGTTCTTGAGCTGGCCGTCGGTGTAGAACTCGGTCTTCCACTGCTAGGTCGGGCCGAGCATTTCCAGGGCGGCAt                                                                                                                      >  1:9631253/1‑151 (MQ=255)
                                 cAGCCAGAGCTTCTCCATGTTCAGTTTCGGATCGCCGCCGCCCTTGAGGTAGAGGTTGCCGTTGAGCACTCCGTTCTTGAGCTGGCCGTCGGTGTAGAACTCGGTCTTCCACTGCTAGGTCGGGCCGAGCATTTCCAGGGCGGCATAggtg                                                                                                                 >  1:8119757/1‑151 (MQ=255)
                                 cAGCCAGAGAATCTCAATGGTCAGTTTCCGAGCGCCGCCGCCCTTGAGGTAGTGGTTGCCGTCGAGCACTCCGTTCTTGCGCGGGCCGTCGGTGTAGAACTCGGTCTGCCACTGCTAGGTCGGGCCGAGCATTTCCAGGGCGGAATCggtg                                                                                                                 <  2:10234063/151‑1 (MQ=255)
                                   gCCAGAGCTTCTCCATGTTCAGTTTCGGATCGCCGCCGCCCTTGAGGTAGAGGTTGCCGTTGAGCACTCCGTTCTTGAGCTGGCCGTCGGTGTAGAACTCGGTCTTCCACTGCTAGGTCGGGCCGAGCATTTCCAGGGCGGCATAggtggt                                                                                                               <  1:9034911/151‑1 (MQ=255)
                                   cccAGTGCTTATCGTTGTTCAGAGTCGCTGCTCCGCCGCCCTTGCGGGCGTCTTCGCGCTTGTGCACTCCGGTGTTGAGCGGTCGGTCGGTGTAGACGTCGGTCTGCCAATGCTAGGTCTGGCCGAGCATTTCCAGGGCGGCATAggtggt                                                                                                               <  2:59801/150‑1 (MQ=255)
                                    ccAGAGCTTCTCCATGTTCAGTTTCGGATCGCCGCCGCCCTTGAGGTAGAGGTTGCCGTTGAGCACTCCGTTCTTGAGCTGGCCGTCGGTGTAGAACTCGGTCTTCCACTGCTAGGTCGGGCCGAGCATTTCCAGGGCGGCATAGgtggtg                                                                                                              <  2:8093428/151‑1 (MQ=255)
                                     cAGAGCTTCTCCATGTTCAGTTTCGGATCGCCGCCGCCCTTGAGGTAGAGGTTGCCGTTGAGCACTCCGTTCTTGAGCTGGCCGTCGGTGTAGAACTCGGTCTTCCACTGCTAGGTCGGGCCGAGCATTTCCAGGGCGGCATAGGTGGTGa                                                                                                             >  1:8578992/1‑151 (MQ=255)
                                     cAGAGCTTCTCCATGTTCAGTTTCGGATCGCCGCCGCCCTTGAGGTAGAGGTTGCCGTTGAGCACTCCGTTCTTGAGCTGGCCGTCGGTGTAGAACTCGGTCTTCCACTGCTAGGTCGGGCCGAGCATTTCCAGGGCGGCATAGGTGGTGa                                                                                                             <  2:8450583/151‑1 (MQ=255)
                                     cAGAGCTTCTCCATGTTCAGTTTCGGATCGCCGCCGCCCTTGAGGTAGAGGTTGCCGTTGAGCACTCCGTTCTTGAGCTGGCCGTCGGTGTAGAACTCGGTCTTCCACTGCTAGGTCGGGCCGAGCATTTCCAGGGCGGCATAGGTGGTGa                                                                                                             >  1:4283479/1‑151 (MQ=255)
                                      agagCTTCTCCATGTTCAGTTTCGGATCGCCGCCGCCCTTGAGGTAGAGGTTGCCGTTGAGCACTCCGTTCTTGAGCTGGCCGTCGGTGTAGAACTCGGTCTTCCACTGCTAGGTCGTGCCGAGCATTTCCAGGGCGGCATAGGTGGTGaa                                                                                                            >  2:5365123/1‑151 (MQ=255)
                                      agagCTTCTCCATGTTCAGTTTCGGATCGCCGCCGCCCTTGAGGTAGAGGTTGCCGTTGAGCACTCCGTTCTTGAGCTGGCCGTCGGTGTAGAACTCGGTCTTCCACTGCTAGGTCGGGCCGAGCATTTCCAGGGCGGCATAGGTGGTGaa                                                                                                            <  1:9163106/151‑1 (MQ=255)
                                      agagCTTCTCCATGTTCAGTTTCGGATCGCCGCCGCCCTTGAGGTAGAGGTTGCCGTTGAGCACTCCGTTCTTGAGCTGGCCGTCGGTGTAGAACTCGGTCTTCCACTGCTAGGTCGGGCCGAGCATTTCCAGGGCGGCATAGGTGGTGaa                                                                                                            <  2:1563505/151‑1 (MQ=255)
                                         gCTTCTCCATGTTCAGTTTCGGATCGCCGCCGCCCTTGAGGTAGAGGTTGCCGTTGAGCACTCCGTTCTTGAGCTGGCCGTCGGTGTAGAACTCGGTCTTCCACTGCTAGGTCGGGCCGAGCATTTCCAGGGCGGCATAGGTGGTGAACAg                                                                                                         >  1:1677172/1‑151 (MQ=255)
                                             ctcCATGTTCAGTTTCGGATCGCCGCCGCCCTTGAGGTAGAGGTTGCCGTTGAGCACTCCGTTCTTGAGCTGGCCGTCGGTGTAGAACTCGGTCTTCCACTGCTAGGTCGGGCCGAGCATTTCCAGGGCGGCATAGGTGGTGAACAGCTTc                                                                                                     <  1:1433005/151‑1 (MQ=255)
                                              tcCATGTTCAGTTTCGGATCGCCGCCGCCCTTGAGGTAGAGGTTGCCGTTGAGCACTCCGTTCTTGAGCTGGCCGTCGGTGTAGAACTCGGTCTTCCACTGCTAGGTCGGGCCGAGCATTTCCAGGGCGGCATAGGTGGTGAACAGCTTCa                                                                                                    <  1:4121622/151‑1 (MQ=255)
                                              tcCATGTTCAGTTTCGGATCGCCGCCGCCCTTGAGGTAGAGGTTGCCGTTGAGCACTCCGTTCTTGAGCTGGCCGTCGGTGTAGAACTCGGTCTTCCACTGCTAGGTCGGGCCGAGCATTTCCAGGGCGGCATAGGTGGTGAACAGCTTCa                                                                                                    >  1:1009765/1‑151 (MQ=255)
                                                 aTGTTCAGTTTCGGATCGCCGCCGCCCTTGAGGTAGAGGTTGCCGTTGAGCACTCCGTTCTTGAGCTGGCCGTCGGTGTAGAACTCGGTCTTCCACTGCTAGGTCGGGCCGAGCATTTCCAGGGCGTCATAGGTGGTGAACAGCCTCATgg                                                                                                 >  1:5371819/1‑151 (MQ=255)
                                                 aTGTTCAGTTTCGGATCGCCGCCGCCCTTGAGGTAGAGGTTGCCGTTGAGCACTCCGTTCTTGAGCTGGCCGTCGGTGTAGAACTCGGTCTTCCACTGCTAGGTCGGGCCGAGCATTTCCAGGGCGGCATAGGTGGTGAACAGCTTCATgg                                                                                                 >  1:10971690/1‑151 (MQ=255)
                                                  tGTTCAGTTTCGGATCGCCGCCGCCCTTGAGGTAGAGGTTGCCGTTGAGCACTCCGTTCTTGAGCTGGCCGTCGGTGTAGAACTCGGTCTTCCACTGCTAGGTCGGGCCGAGCATTTCCAGGGCGGCATAGGTGGTGAACAGCTTCATGGt                                                                                                <  1:2461197/151‑1 (MQ=255)
                                                       aGTTTCGGATCGCCGCCGCCCTTGAGGTAGAGGTTGCCGTTGAGCACTCCGTTCTTGAGCTGGCCGTCGGTGTAGAACTCGGTCTTCCACTGCTAGGTCGGGCCGAGCATTTCCAGGGCGGCATAGGTGGTGAACAGCTTCATGGTCGAgg                                                                                           >  1:6859050/1‑151 (MQ=255)
                                                       aGTTTCGGATCGCCGCCGCCCTTGAGGTAGAGGTTGCCGTTGAGCACTCCGTTCTTGAGCTGGCCGTCGGTGTAGAACTCGGTCTTCCACTGCTAGGTCGGGCCGAGCATTTCCAGGGCGGCATAGGTGGTGAACAGCTTCATGGTCGAgg                                                                                           >  1:9914207/1‑151 (MQ=255)
                                                       aGTTTCGGATCGCCGCCGCCCTTGAGGTAGAGGTTGCCGTTGAGCACTCCGTTCTTGAGCTGGCCGTCGGTGTAGAACTCGGTCTTCCACTGCTAGGTCGGGCCGAGCATTTCCAGGGCGGCATAGGTGGTGAACAGCTTCATGGTCGAgg                                                                                           <  1:10198888/151‑1 (MQ=255)
                                                       aGTTTCGGATCGCCGCCGCCCTTGAGGTAGAGGTTGCCGTTGAGCACTCCGTTCTTGAGCTGGCCGTCGGTGTAGAACTCGGTCTTCCACTGCTAGGTCGGGCCGAGCATTTCCAGGGCGGCATAGGTGGTGAACAGCTTCATGGTCGAgg                                                                                           <  2:6107427/151‑1 (MQ=255)
                                                        gTTTCGGATCGCCGCCGCCCTTGAGGTAGAGGTTGCCGTTGAGCACTCCGTTCTTGAGCTGGCCGTCGGTGTAGAACTCGGTCTTCCACTGCTAGGTCGGGCCGAGCATTTCCAGGGCGGCATAGGTGGTGAACAGCTTCATGGTCGAGGc                                                                                          >  1:9799764/1‑151 (MQ=255)
                                                         tttCGGATCGCCGCCGCCCTTGAGGTAGAGGTTGCCGTTGAGCACTCCGTTCTTGAGCTGGCCGTCGGTGTAGAACTCGGTCTTCCACTGCTAGGTCGGGCCGAGCATTTCCAGGGCGGCATAGGTGGTGAACAGCTTCATGGTCGAGGcc                                                                                         >  2:3796224/1‑151 (MQ=255)
                                                          ttCGGATCGCCGCCGCCCTTGAGGTAGAGGTTGCCGTTGAGCACTCCGTTCTTGAGCTGGCCGTCGGTGTAGAACTCGGTCTTCCACTGCTAGGTCGGGCCGAGCATTTCCAGGGCGGCATAGGTGGTGAACAGCTTCATGGTCGAGGCCg                                                                                        <  2:2699851/151‑1 (MQ=255)
                                                           tCGGATCGCCGCCGCCCTTGAGGTAGAGGTTGCCGTTGAGCACTCCGTTCTTGAGCTGGCCGTCGGTGTAGAACTCGGTCTTCCACTGCTAGGTCGGGCCGAGCATTTCCAGGGCGGCATAGGTGGTGAACAGCTTCATGGTCGAGGCCgg                                                                                       >  1:4921846/1‑151 (MQ=255)
                                                               aTCGCCGCCGCCCTTGAGGTAGAGGTTGCCGTTGAGCACTCCGTTCTTGAGCTGGCCGTCGGTGTAGAACTCGGTCTTCCACTGCTAGGTCGGGCCGAGCATTTCCAGGGCGGCATAGGTGGTGAACAGCTTCATGGTCGAGGCCGGGTTc                                                                                   <  2:977254/151‑1 (MQ=255)
                                                                tCGCCGCCGCCCTTGAGGTAGAGGTTGCCGTTGAGCACTCCGTTCTTGAGCTGGCCGTCGGTGTAGAACTCGGTCTTCCACTGCTAGGTCGGGCCGAGCATTTCCAGGGCGGCATAGGTGGTGAACAGCTTCATGGTCGAGGCCGGGTTCa                                                                                  <  2:5742608/151‑1 (MQ=255)
                                                                   ccgccgccCTTGAGGTAGAGGTTGCCGTTGAGCACTCCGTTCTTGAGCTGGCCGTCGGTGTAGAACTCGGTCTTCCACTGCTAGGTCGGGCCGAGCATTTCCAGGGCGGCATAGGTGGTGAACAGCTTCATGGTCGAGGCCGGGTTCACCg                                                                               <  1:8820995/151‑1 (MQ=255)
                                                                     gccgccCTTGAGGTAGAGGTTGCCGTTGAGCACTCCGTTCTTGAGCTGGCCGTCGGTGTAGAACTCGGTCTTCCACTGCTAGGTCGGGCCGAGCATTTCCAGGGCGGCATAGGTGGTGAACAGCTTCATGGTCGAGGCCGGGTTCACCGAt                                                                             <  2:6684003/151‑1 (MQ=255)
                                                                        gccCTTGAGGTAGAGGTTGCCGTTGAGCACTCCGTTCTTGAGCTGGCCGTCGGTGTAGAACTCGGTCTTCCACTGCTAGGTCGGGCCGAGCATTTCCAGGGCGGCATAGGTGGTGAACAGCTTCATGGTCGAGGCCGGGTTCACCGATACg                                                                          >  2:6469209/1‑151 (MQ=255)
                                                                        gccCTTGAGGTAGAGGTTGCCGTTGAGCACTCCGTTCTTGAGCTGGCCGTCGGTGTAGAACTCGGTCTTCCACTGCTAGGTCGGGCCGAGCATTTCCAGGGCGGCATAGGTGGTGAACAGCTTCATGGTCGAGGCCGGGTTCACCGATACg                                                                          <  1:2125221/151‑1 (MQ=255)
                                                                         cccTTGAGGTAGAGGTTGCCGTTGAGCACTCCGTTCTTGAGCTGGTCGTCGGTGTAGAACTCGGTGGTCCCCGTATCATTAAAAAAGAGAATTTACAGGGCGGCGTCGGGGGTGAAGAGCTTCCTGGGGGAGGCGGGGTGCACCgtgaggt                                                                         >  2:9467027/1‑145 (MQ=255)
                                                                          ccTTGAGGTAGAGGTTGCCGTTGAGCACTCCGTTCTTGAGCTGGCCGTCGGTGTAGAACTCGGTCTTCCACTGCTAGGTCGGGCCGAGCATTTCCAGGGCGGCATAGGTGGTGAACAGCTTCATGGTCGAGGCCGGGTTCACCGATACGTc                                                                        <  2:3009369/151‑1 (MQ=255)
                                                                           cTTGAGGTAGAGGTTGCCGTTGAGCACTCCGTTCTTGAGCTGGCCGTCGGTGTAGAACTCGGTCTTCCACGGCTAGGTCTGGCCGAGCATTTCCAGGGCGGCATAGGTGGTGAACAGCTTCTTGGTCGAGGCCGGGTTCTCCGATACGTCg                                                                       >  2:7344001/1‑151 (MQ=255)
                                                                            ttGAGGTAGAGGTTGCCGTTGAGCACTCCGTTCTTGAGCTGGCCGTCGGTGTAGAACTCGGTCTTCCACTGCTAGGTCGGGCCGAGCATTTCCAGGGCGGCATAGGTGGTGAACAGCTTCATGGTCGAGGCCGGGTTCACCGATACGTCgg                                                                      <  1:7645460/151‑1 (MQ=255)
                                                                            ttGAGGTAGAGGTTGCCGTTGAGCACTCCGTTCTTGAGCTGGCCGTCGGTGTAGAACTCGGTCTTCCACTGCTAGGTCGGGCCGAGCATTTCCAGGGCGGCATAGGTGGTGAACAGCTTCATGGTCGAGGCCGGGTTCACCGATACGTCgg                                                                      >  1:4245444/1‑151 (MQ=255)
                                                                                ggTAGAGGTTGCCGTTGAGCACTCCGTTCTTGAGCTGGCCGTCGGTGTAGAACTCGGTCTTCCACTGCTAGGTCGGGCCGAGCATTTCCAGGGCGGCATAGGTGGTGAACAGCTTCATGGTCGAGGCCGGGTTCACCGATACGTCGGCGtt                                                                  <  1:2832274/151‑1 (MQ=255)
                                                                                       gTTGCCGTTGAGCACTCCGTTCTTGAGCTGGCCGTCGGTGTAGAACTCGGTCTTCCACTGCTAGGTCGGGCCGAGCATTTCCAGGGCGGCATAGGTGGTGAACAGCTTCATGGTCGAGGCCGGGTTCACCGATACGTCGGCGTTGtagtag                                                           >  1:7596403/1‑151 (MQ=255)
                                                                                       gTTGCCGTTGAGCACTCCGTTCTTGAGCTGGCCGTCGGTGTAGAACTCGGTCTTCCACTGCTAGGTCGGGCCGAGCATTTCCAGGGCGGCATAGGTGGTGAACAGCTTCATGGTCGAGGCCGGGTTCACCGATACGTCGGCGTTGtagtag                                                           <  2:9609241/151‑1 (MQ=255)
                                                                                          gCCGTTGAGCACTCCGTTCTTGAGCTGGCCGTCGGTGTAGAACTCGGTCTTCCACTGCTAGGTCGGGCCGAGCATTTCCAGGGCGGCATAGGTGGTGAACAGCTTCATGGTCGAGGCCGGGTTCACCGATACGTCGGCGTTGTAGTAGGTg                                                        <  1:4446137/151‑1 (MQ=255)
                                                                                           ccGTTGAGCACTCCGTTCTTGAGCTGGCCGTCGGTGTAGAACTCGGTCTTCCACTGCTAGGTCGGGCCGAGCATTTCCAGGGCGGCATAGGTGGTGAACAGCTTCATGGTCGAGGCCGGGTTCACCGATACGTCGGCGTTGTAGTAGGTgg                                                       >  1:569546/1‑151 (MQ=255)
                                                                                           ccGTTGAGCACTCCGTTCTTGAGCTGGCCGTCGGTGTAGAACTCGGTCTTCCACTGCTAGGTCGGGCCGAGCATTTCCAGGGCGGCATAGGTGGTGAACAGCTTCATGGTCGAGGCCGGGTTCACCGATACGTCGGCGTTGTAGTAGGTgg                                                       >  1:569408/1‑151 (MQ=255)
                                                                                            cGTTGAGCACTCCGTTCTTGAGCTGGCCGTCGGTGTAGAACTCGGTCTTCCACTGCTAGGTCGGGCCGAGCATTTCCCGGGCGGCATAGGTGGTGAACAGCTTCATGGTCGAGGCCGGGTTCACCGATACGTCGGCGTTGTAGTAGGTggg                                                      <  2:4166163/151‑1 (MQ=255)
                                                                                            cGTTGAGCACTCCGTTCTTGAGCTGGCCGTCGGTGTAGAACTCGGTCTTCCACTGCTAGGTCGGGCCGAGCATTTCCAGGGCGGCATAGGTGGTGAACAGCTTCATGGTCGAGGCCGGGTTCACCGATACGTCGGCGGTGTAGTAGGTggg                                                      >  1:10631002/1‑151 (MQ=255)
                                                                                             gTTGAGCACTCCGTTCTTGAGCTGGCCGTCGGTGTAGAACTCGGTCTTCCACTGCTAGGTCGGGCCGAGCATTTCCAGGGCGGCGTAGGTGGTGAACAGCTTCAGGGTCGAGGCCGGGTTCTCCGGGACGGCGGGGGGGGGGTAGGTgggg                                                     >  2:9278591/1‑151 (MQ=255)
                                                                                             gTTGAGCACTCCGTTCTTGAGCTGGCCGTCGGTGTAGAACTCGGTCTTCCACTGCTAGGTCGGGCCGAGCATTTCCAGGGCGGCATAGGTGGTGAACAGCTTCATGGTCGAGGCCGGGTTCACCGATACGTCGGCGTTGTAGTAGGTgggg                                                     <  2:10445942/151‑1 (MQ=255)
                                                                                              ttGAGCACTCCGTTCTTGAGCTGGCCGTCGGTGTAGAACTCGGTCTTCCACTGCTAGGTCGGGCCGAGCATTTCCAGGGCGGCATAGGTGGTGAACAGCTTCATGGTCGAGGCCGGGTTCACCGATACGTCGGCGTTGTAGTAGGTGGGGt                                                    <  1:1735148/151‑1 (MQ=255)
                                                                                               tGAGCACTCCGTTCTTGATCTGGCCGTCGGTGTAGTACTCGGTCTTCCACTGCTAGGTCGGGCCGAGCATTTCCAGGGCGGCATAGGTGGTGAACAGCTTCATTGTCGAGGCCGTGTTCACCGATACGTCGGCTTTGTAGTAGGTGGGGtt                                                   <  1:4854881/151‑1 (MQ=255)
                                                                                                 aGCACTCCGTTCTTGAGCTGGCCGTCGGTGTAGAACTCGGTCTTCCACTGCTAGGTCGGGCCGAGCATTTCCAGGGCGGCATAGGTGGTGAACAGCTTCATGGTCGAGGCCGGGTTCACCGATACGTCGGCGTTGTAGTAGGTGGGGTTGc                                                 <  2:3328869/151‑1 (MQ=255)
                                                                                                   cACTCCGTTCTTGAGCTGGCCGTCGGTGTAGAACTCGGTCTTCCACTGCTAGGTCGGGCCGAGCATTTCCAGGGCGGCATAGGTGGTGAACAGCTTCATGGTCGAGGCCGGGTTCACCGATACGTCGGCGTTGTAGTAGGTGGGGTTGccc                                               >  1:9336535/1‑151 (MQ=255)
                                                                                                      tCCGTTCTTGAGCTGGCCGTCGGTGTAGAACTCGGTCTTCCACTGCTAGGTCGGGCCGAGCATTTCCAGGGCGGCATAGGTGGTGAACAGCTTCATGGTCGAGGCCGGGTTCACCGATACGTCGGCGTTGTAGTAGGTGGGGTTGCCCggg                                            <  1:5761905/151‑1 (MQ=255)
                                                                                                         gTTCTTGAGCTGGCCGTCGGTGTAGAACTCGGTCTTCCACTGCTAGGTCGGGCCGAGCATTTCCAGGGCGGCATTGGTGGTGAACAGCTTCATGGTCGAGGCCGGGTTCACCGATACGTCGGCGTTGTAGTAGGTGGGGTTGCCCGGGCCg                                         >  1:4573757/1‑151 (MQ=255)
                                                                                                         gTTCTTGAGCTGGCCGTCGGTGTAGAACTCGGTCTTCCACTGCTAGGTCGGGCCGAGCATTTCCAGGGCGGCATAGGTGGTGAACAGCTTCATGGTCGAGGCCGGGTTCACCGATACGTCGGCGTTGTAGTAGGTGGGGTTGCCCGGGCCg                                         <  2:5983790/151‑1 (MQ=255)
                                                                                                           tCTTGAGCTGGCCGTCGGTGTAGAACTCGGTCTTCCACTGCTAGGTCGGGCCGAGCATTTCCAGGGCGGCATAGGTGGTGAACAGCTTCATGGTCGAGGCCGGGTTCACCGATACGTCGGCGTTGTAGTAGGTGGGGTTGCCCGGGCCGTc                                       >  1:7169235/1‑151 (MQ=255)
                                                                                                           tCTTGAGCTGGCCGTCGGTGTAGAACTCGGTCTTCCACTGCTAGGTCGGGCCGAGCATTTCCAGGGCGGCATAGGTGGTGAACAGCTTCATGGTCGAGGCCGGGTTCACCGATACGTCGGCGTTGTAGTAGGTGGGGTTGCCCGGGCCGTc                                       <  2:7013473/151‑1 (MQ=255)
                                                                                                           tCTTGAGCTGGCCGTCGGTGTAGAACTCGGTCTTCCACTGCTAGGTCGGGCCGAGCATTTCCAGGGCGGCATAGGTGGTGAACAGCTTCATGGTCGAGGCCGGGTTCACCGATACGTCGGCGTTGTAGTAGGTGGGGTTGCCCGGGCCGTc                                       <  2:10706861/151‑1 (MQ=255)
                                                                                                               gAGCTGGCCGTCGGTGTAGAACTCGGTCTTCCACTGCTAGGTCGGGCCGAGCATTTCCAGGGCGGCATAGGTGGTGAACAGCTTCATGGTCGAGGCCGGGTTCACCGATACGTCGGCGTTGTAGTAGGTGGGGTTGCCCGGGCCGTCGAGc                                   <  2:5641737/151‑1 (MQ=255)
                                                                                                                  cTGGCCGTCGGTGTATAACTCGTTCTTCCACTGCTAGGTCGGGCCGAGCATTTCCAGGGCGGCATAGGTGGTGAACAGCTTCATGGTCGAGGCCGGGTTCACCGATACGTCGTCTTTGTAGTAGGTGGGGTTGACCGGGACGTCGAGCGGa                                >  1:2729654/1‑151 (MQ=255)
                                                                                                                    ggCCGTCGGTGTAGAACTCGGTCTTCCACTGCTAGGTCGGGCCGAGCATTTCCAGGGCGGCATAGGTGGTGAACAGCTTCATGGTCGAGGCCGGGTTCACCGATACGTCGGCGTTGTAGTAGGTGGGGTTGCCCGGGCCGTCGAGCGGAAt                              <  2:7790364/151‑1 (MQ=255)
                                                                                                                     gCCGTCGGTGTAGAACTCGGTCTTCCACTGCTAGGTCGGGCCGAGCATTTCCAGGGCGGCATAGGTGGTGAACAGCTTCATGGTCGAGGCCGGGTTCACCGATACGTCGGCGTTGTAGTAGGTGGGGTTGCCCGGGCCGTCGAGCGGAATc                             <  1:2280959/151‑1 (MQ=255)
                                                                                                                      ccGTCGGTGTAGAACTCGGTCTTCCACTGCTAGGTCGGGCCGAGCATTTCCAGGGCGGCATAGGTGGTGAACAGCTTCATGGTCGAGGCCGGGTTCACCGATACGTCGGCGTTGTAGTAGGTGGGGTTGCCCGGGCCGTCGAGCGGAATCa                            >  1:1998280/1‑151 (MQ=255)
                                                                                                                        gTCGGTGTAGAACTCGGTCTTCCACTGCTAGGTCGGGCCGAGCATTTCCAGGGCGGCATAGGTGGTGAACAGCTTCATGGTCGAGGCCGGGTTCACCGATACGTCGGCGGTGTAGTAGGTGGGGTTGCCCGGGCCGTCGAGCGGAATCAGc                          >  2:10582880/1‑151 (MQ=255)
                                                                                                                        gTCGGTGTAGAACTCGGTCTTCCACTGCTAGGTAGGGCCGAGCATTTCCAGGGCGGCATGGGTGGTGAACAGCTTCGTGGTCGAGGCCGGGTTCACCGATTCGTCGGCGGTGTAGGAGGGGGGGTTGCCCGGGCCGTCGAGCGGAATCTGc                          >  2:908174/1‑151 (MQ=255)
                                                                                                                               tAGAACTCGGTCTTCCACTGCTAGGTCGGGCCGAGCATTTCCAGGGCGGCATAGGTGGTGAACAGCTTCATGGTCGAGGCCGGGTTCACCGATACGTCGGCGTTGTAGTAGGTGGGGTTGCCCGGGCCGTCGAGCGGAATCAGCACCAGCg                   <  2:4938087/151‑1 (MQ=255)
                                                                                                                                  aaCTCGGTCTTCCACTGCTAGGTCGGGCCGAGCATTTCCAGGGCGGCATAGGTGGTGCACAGCTTCATGGTCGAGGCCGGGTTCACCGATACGTCGGCGTTGTAGTAGGTGGGGTTGCCCGGGCCGTCGAGCGGAATCAGCACCAGCGACa                <  2:2519760/151‑1 (MQ=255)
                                                                                                                                    cTCGGTCTTCCACTGCTAGGTCGGGCCGAGCATTTCCAGGGCGGCATAGGTGGTGAACAGCTTCATGGTCGAGGCCGGGTTCACCGATACGTCGGCGTTGTAGTAGGTGGGGTTGCCCGGGCCGTCGAGCGGAATCAGCACCAGCGACAGa              >  1:5228394/1‑151 (MQ=255)
                                                                                                                                     tCGGTCTTCCACTGCTAGGTCGGGCCGAGCATTTCCAGGGCGGCATAGGTGGTGAACAGCTTCATGGTCGAGGCCGGGTTCACCGATACGTCGGCGTTGTAGTAGGTGGGGTTGCCCGGGCCGTCGAGCGGAATCAGCACCAGCGACAGAt             >  1:8788275/1‑151 (MQ=255)
                                                                                                                                      cGGTCTTCCACTGCTAGGTCGGGCCGAGCATTTCCAGGGCGGCATAGGTGGTGAACAGCTTCATGGTCGAGGCCGGGTTCACCGATACGTCGGCGTTGTAGTAGGTGGGGTTGCCCGGGCCGTCGAGCGGAATCAGCACCAGCGACAGATc            <  2:8970925/151‑1 (MQ=255)
                                                                                                                                         tCTTCCACTGCTAGGTCGGGCCGAGCTTTTCCAGGGCGGCATAGGTGGTGAACAGCTTCATGGTCGAGGCCGGGTTCACCGATACGTCGGCGTTGTAGTAGGTGGGGTTGCCCGGGCCGTCGAGCGGAATCAGCACCAGCGACAGATCGtt         <  1:6418352/151‑1 (MQ=255)
                                                                                                                                         tCTTCCACTGCTAGGTCGGGCCGAGCATTTCCAGGGCGGCATAGGTGGTGAACAGCTTCATGGTCGAGGCCGGGTTCACCGATACGTCGGCGTGGTAGTAGGTGGGGTTGCCCGGGCCGTCGAGCGGAATCAGCACCAGCGACAGATCGtt         >  1:6042608/1‑151 (MQ=255)
                                                                                                                                             ccACTGCTAGGTCGGGCCGAGCATTTCCAGGGCGGCATAGGTGGTGAACAGCTTCATGGTCGAGGCCGGGTTCACCGATACGTCGGCGTTGTAGTAGGTGGGGTTGCCCGGGCCGTCGAGCGGAATCAGCACCAGCGACAGATCGTTGCCg     >  1:30285/1‑151 (MQ=255)
                                                                                                                                              cACTGCTAGGTCGGGCCGAGCATTTCCAGGGCGGCATAGGTGGTGAACAGCTTCATGGTCGAGGCCGGGTTCACCGATACGTCGGCGTTGTAGTAGGTGGGGTTGCCCGGGCCGTCGAGCGGAATCAGCACCAGCGACAGATCGTTGCCgg    >  2:1783067/1‑151 (MQ=255)
                                                                                                                                              cACTGCTAGGTCGGGCCGAGCATTTCCAGGGCGGCATAGGTGGTGAACAGCTTCATGGTCGAGGCCGGGTTCACCGATACGTCGGCGTTGTAGTAGGTGGGGTTGCCCGGGCCGTCGAGCGGAATCAGCACCAGCGACAGATCGGTGCCgg    >  2:1794242/1‑151 (MQ=255)
                                                                                                                                               aCTGCTAGGTCGGGCCGAGCATTTCCAGGGCGGCATAGGTGGTGAACAGCTTCATGGTCGAGGCCGGGTTCACCGATACGTCGGCGTTGTAGTAGGTGGGGTTGCCCGGGCCGTCGAGCGGAATCAGCACCAGCGACAGATCGTTGCCGGt   >  1:4759383/1‑151 (MQ=255)
                                                                                                                                               aCTGCTAGGTCGGGCCGAGCATTTCCAGGGCGGCATAGGTGGTGAACAGCTTCATGGTCGAGGCCGGGTTCACCGATACGTCGGCGTTGTAGTAGGTGGGGTTGCCCGGGCCGTCGAGCGGAATCAGCACCAGCGACAGATCGTTGCCGGt   >  1:6214033/1‑151 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GGTCACGCCATTGGCGCGCAGGTCACGCATCAACAGCCAGAGCTTCTCCATGTTCAGTTTCGGATCGCCGCCGCCCTTGAGGTAGAGGTTGCCGTTGAGCACTCCGTTCTTGAGCTGGCCGTCGGTGTAGAACTCGGTCTTCCACTGGTAGGTCGGGCCGAGCATTTCCAGGGCGGCATAGGTGGTGAACAGCTTCATGGTCGAGGCCGGGTTCACCGATACGTCGGCGTTGTAGTAGGTGGGGTTGCCCGGGCCGTCGAGCGGAATCAGCACCAGCGACAGATCGTTGCCGGTC  >  NC_008463/2158676‑2158970

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