Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 620,891 C→T 8.3% P184L (CCG→CTG)  entC → isochorismate synthase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273620,8910CT 8.3% 656.0 / 31.6 192P184L (CCG→CTG) entCisochorismate synthase 1
Reads supporting (aligned to +/- strand):  ref base C (82/94);  new base T (7/9);  total (89/103)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

Only 100 of 193 total aligned reads displayed.

CGCCGCCGCACTTACCGCCACGCCGCAGGTCGACAAAGTGGTGTTGTCACGGTTGATTGATATCACCACTGACGCCGCCATTGATAGTGGCGTATTGCTGGAACGGTTGATTGCGCAAAACCCGGTTAGTTACAACTTCCATGTTCCGCTGGCTGATGGTGGCGTCCTGCTGGGGGCCAGCCCGGAACTGCTGCTACGTAAAGACGGCGAGCGTTTTAGCTCCATTCCGTTAGCCGGTTCCGCGCGTCGTCAGCCGGATGAAGTGCTCGATCGCGAAGCAGGTAATCGTCTGC  >  CP009273/620745‑621037
                                                                                                                                                  |                                                                                                                                                  
cgccgccgcACTTACCGCCACGCCGCAGGTCGACAAAGTGGTGTTGTCACGGTTGATTGATATCACCACTGACGCCGCCATTGATAGTGGCGTATTGCTGGAACGGTTGATTGCGCAAAACCCGGTTAGTTACAACTTCCATGTTCCg                                                                                                                                                   <  1:1169652/148‑1 (MQ=255)
 gccgccgcACTTACCGCCACGCCGCAGGTCGACAAAGTGGTGTTGTCACGGTTGATTGATATCACCACTGACGCCGCCATTGATAGTGGCGTATTGCTGGAACGGTTGATTGCGCAAAACCCGGTTAGTTACAACTTCCATGTTCCgc                                                                                                                                                  <  1:3422899/148‑1 (MQ=255)
 gccgccgcACTTACCGCCACGCCGCAGGTCGACAAAGTGGTGTTGTCACGGTTGATTGATATCACCACTGACGCCGCCATTGATAGTGGCGTATTGCTGGAACGGTTGATTGCGCAAAACCCGGTTAGTTACAACTTCCATGTTCCgc                                                                                                                                                  <  2:1340073/148‑1 (MQ=255)
     ccgcACTTACCGCCACGCCGCAGGTCGACAAAGTGGTGTTGTCACGGTTGATTGATATCACCACTGACGCCGCCATTGATAGTGGCGTATTGCTGGAACGGTTGATTGCGCAAAACCCGGTTAGTTACAACTTCCATGTTCCgctggc                                                                                                                                              <  1:2448095/148‑1 (MQ=255)
       gcACTTACCGCCACGCCGCAGGTCGACAAAGTGGTGTTGTCACGGTTGATTGATATCACCACTGACGCCGCCATTGATAGTGGCGTATTGCTGGAACGGTTGATTGCGCAAAACCCGGTTAGTTACAACTTCCATGTTCCgctggctg                                                                                                                                            <  1:534437/148‑1 (MQ=255)
          cTTACCGCCACGCCGCAGGTCGACAAAGTGGTGTTGTCACGGTTGATTGATATCACCACTGACGCCGCCATTGATAGTGGCGTATTGCTGGAACGGTTGATTGCGCAAAACCCGGTTAGTTACAACTTCCATGTTCCGCTGGCTGAtg                                                                                                                                         <  2:437201/148‑1 (MQ=255)
          cTTACCGCCACGCCGCAGGTCGACAAAGTGGTGTTGTCACGGTTGATTGATATCACCACTGACGCCGCCATTGATAGTGGCGTATTGCTGGAACGGTTGATTGCGCAAAACCCGGTTAGTTACAACTTCCATGTTCCGCTGGCTGAtg                                                                                                                                         <  2:1802131/148‑1 (MQ=255)
           ttACCGCCACGCCGCAGGTCGACAAAGTGGTGTTGTCACGGTTGATTGATATCACCACTGACGCCGCCATTGATAGTGGCGTATTGCTGGAACGGTTGATTGCGCAAAACCCGGTTAGTTACAACTTCCATGTTCCGCTGGCTGAt                                                                                                                                          >  2:2437191/1‑146 (MQ=255)
            tACCGCCACGCCGCAGGTCGACAAAGTGGTGTTGTCACGGTTGATTGATATCACCACTGACGCCGCCATTGATAGTGGCGTATTGCTGGAACGGTTGATTGCGCAAAACCCGGTTAGTTACAACTTCCATGTTCCGCTGGCTGAtggt                                                                                                                                       <  2:378035/148‑1 (MQ=255)
             aCCGCCACGCCGCAGGTCGACAAAGTGGTGTTGTCACGGTTGATTGATATCACCACTGACGCCGCCATTGATAGTGGCGTATTGCTGGAACGGTTGATTGCGCAAAACCCGGTTAGTTACAACTTCCATGTTCCGCTGGCTGAtggtg                                                                                                                                      >  1:1647139/1‑148 (MQ=255)
                  cACGCCGCAGGTCGACAAAGTGGTGTTGTCACGGTTGATTGATATCACCACTGACGCCGCCATTGATAGTGGCGTATTGCTGGAACGGTTGATTGCGCAAAACCCGGTTAGTTACAACTTCCATGTTCCGCTGGCTGATGGTGGCGTc                                                                                                                                 >  2:332160/1‑148 (MQ=255)
                   aCGCCGCAGGTCGACAAAGTGGTGTTGTCACGGTTGATTGATATCACCACTGACGCCGCCATTGATAGTGGCGTATTGCTGGAACGGTTGATTGCGCAAAACCCGGTTAGTTACAACTTCCATGTTCCGCTGGCTGATGGTGGCGTcc                                                                                                                                <  1:703025/148‑1 (MQ=255)
                    cgccgcAGGTCGACAAAGTGGTGTTGTCACGGTTGATTGATATCACCACTGACGCCGCCATTGATAGTGGCGTATTGCTGGAACGGTTGATTGCGCAAAACCCGGTTAGTTACAACTTCCATGTTCCGCTGGCTGATGGTGGCGTCct                                                                                                                               >  2:3256232/1‑148 (MQ=255)
                           ggTCGACAAAGTGGTGTTGTCACGGTTGATTGATATCACCACTGACGCCGCCATTGATAGTGGCGTATTGCTGGAACGGTTGATTGCGCAAAACCCGGTTAGTTACAACTTCCATGTTCCGCTGGCTGATGGTGGCGTCCTGCTgggg                                                                                                                        >  1:489480/1‑148 (MQ=255)
                            gTCGACAAAGTGGTGTTGTCACGGTTGATTGATATCACCACTGACGCCGCCATTGATAGTGGCGTATTGCTGGAACGGTTGATTGCGCAAAACCCGGTTAGTTACAACTTCCATGTTCCGCTGGCTGATGGTGGCGTCCTGCTggggg                                                                                                                       >  1:3113954/1‑148 (MQ=255)
                            gTCGACAAAGTGGTGTTGTCACGGTTGATTGATATCACCACTGACGCCGCCATTGATAGTGGCGTATTGCTGGAACGGTTGATTGCGCAAAACCCGGTTAGTTACAACTTCCATGTTCCGCTGGCTGATGGTGGCGTCCTGCTggggg                                                                                                                       >  1:2505762/1‑148 (MQ=255)
                                aCAAAGTGGTGTTGTCACGGTTGATTGATATCACCACTGACGCCGCCATTGATAGTGGCGTATTGCTGGAACGGTTGATTGCGCAAAACCCGGTTAGTTACAACTTCCATGTTCTGCTGGCTGATGGTGGCGTCCTGCTGGGGGCCAg                                                                                                                   >  1:1535074/1‑148 (MQ=255)
                                aCAAAGTGGTGTTGTCACGGTTGATTGATATCACCACTGACGCCGCCATTGATAGTGGCGTATTGCTGGAACGGTTGATTGCGCAAAACCCGGTTAGTTACAACTTCCATGTTCCGCTGGCTGATGGTGGCGTCCTGCTGGGGGCCAg                                                                                                                   >  2:395669/1‑148 (MQ=255)
                                aCAAAGTGGTGTTGTCACGGTTGATTGATATCACCACTGACGCCGCCATTGATAGTGGCGTATTGCTGGAACGGTTGATTGCGCAAAACCCGGTTAGTTACAACTTCCATGTTCCGCTGGCTGATGGTGGCGTCCTGCTGGGGGCCAg                                                                                                                   >  1:3221141/1‑148 (MQ=255)
                                 cAAAGTGGTGTTGTCACGGTTGATTGATATCACCACTGACGCCGCCATTGATAGTGGCGTATTGCTGGAACGGTTGATTGCGCAAAACCCGGTTAGTTACAACTTCCATGTTCCGCTGGCTGATGGTGGCGTCCTGCTGGGGGCCAGc                                                                                                                  <  2:1369207/148‑1 (MQ=255)
                                  aaaGTGGTGTTGTCACGGTTGATTGATATCACCACTGACGCCGCCATTGATAGTGGCGTATTGCTGGAACGGTTGATTGCGCAAAACCCGGTTAGTTACAACTTCCATGTTCCGCTGGCTGATGGTGGCGTCCTGCTGGGGGCCAGcc                                                                                                                 <  2:2909329/148‑1 (MQ=255)
                                  aaaGTGGTGTTGTCACGGTTGATTGATATCACCACTGACGCCGCCATTGATAGTGGCGTATTGCTGGAACGGTTGATTGCGCAAAACCCGGTTAGTTACAACTTCCATGTTCCGCTGGCTGATGGTGGCGTCCTGCTGGGGGCCAGcc                                                                                                                 <  2:3001703/148‑1 (MQ=255)
                                     gtggtgTTGTCACGGTTGATTGATATCACCACTGACGCCGCCATTGATAGTGGCGTATTGCTGGAACGGTTGATTGCGCAAAACCCGGTTAGTTACAACTTCCATGTTCCGCTGGCTGATGGTGGCGTCCTGCTGGGGGCCAGCCCgg                                                                                                              >  1:1990912/1‑148 (MQ=255)
                                     gtggtgTTGTCACGGTTGATTGATATCACCACTGACGCCGCCATTGATAGTGGCGTATTGCTGGAACGGTTGATTGCGCAAAACCCGGTTAGTTACAACTTCCATGTTCCGCTGGCTGATGGTGGCGTCCTGCTGGGGGCCAGCCCgg                                                                                                              >  1:510200/1‑148 (MQ=255)
                                     gtggtgTTGTCACGGTTGATTGATATCACCACTGACGCCGCCATTGATAGTGGCGTATTGCTGGAACGGTTGATTGCGCAAAACCCGGTTAGTTACAACTTCCATGTTCCGCTGGCTGATGGTGGCGTCCTGCTGGGGGCCAGCCCgg                                                                                                              >  1:508251/1‑148 (MQ=255)
                                      tggtgTTGTCACGGTTGATTGATATCACCACTGACGCCGCCATTGATAGTGGCGTATTGCTGGAACGGTTGATTGCGCAAAACCCGGTTAGTTACAACTTCCATGTTCTGCTGGCTGATGGTGGCGTCCTGCTGGGGGCCAGCCCGGa                                                                                                             >  1:2842794/1‑148 (MQ=255)
                                       ggtgTTGTCACGGTTGATTGATATCACCACTGACGCCGCCATTGATAGTGGCGTATTGCTGGAACGGTTGATTGCGCAAAACCCGGTTAGTTACAACTTCCATGTTCCGCTGGCTGATGGTGGCGTCCTGCTGGGGGCCAGCCCGGaa                                                                                                            >  2:1055063/1‑148 (MQ=255)
                                       ggtgTTGTCACGGTTGATTGATATCACCACTGACGCCGCCATTGATAGTGGCGTATTGCTGGAACGGTTGATTGCGCAAAACCCGGTTAGTTACAACTTCCATGTTCCGCTGGCTGATGGTGGCGTCCTGCTGGGGGCCAGCCCGGaa                                                                                                            <  1:869505/148‑1 (MQ=255)
                                       ggtgTTGTCACGGTTGATTGATATCACCACTGACGCCGCCATTGATAGTGGCGTATTGCTGGAACGGTTGATTGCGCAAAACCCGGTTAGTTACAACTTCCATGTTCCGCTGGCTGATGGTGGCGTCCTGCTGGGGGCCAGCCCGGaa                                                                                                            <  2:88502/148‑1 (MQ=255)
                                          gttgtCACGGTTGATTGATATCACCACTGACGCCGCCATTGATAGTGGCGTATTGCTGGAACGGTTGATTGCGCAAAACCCGGTTAGTTACAACTTCCATGTTCCGCTGGCTGATGGTGGCGTCCTGCTGGGGGCCAGCCCGGAActg                                                                                                         <  1:2677495/148‑1 (MQ=255)
                                          gttgtCACGGTTGATTGATATCACCACTGACGCCGCCATTGATAGTGGCGTATTGCTGGAACGGTTGATTGCGCAAAACCCGGTTAGTTACAACTTCCATGTTCCGCTGGCTGATGGTGGCGTCCTGCTGGGGGCCAGCCCGGAActg                                                                                                         <  1:39402/148‑1 (MQ=255)
                                             gtCACGGTTGATTGATATCACCACTGACGCCGCCATTGATAGTGGCGTATTGCTGGAACGGTTGATTGCGCAAAACCCGGTTAGTTACAACTTCCATGTTCCGCTGGCTGATGGTGGCGTCCTGCTGGGGGCCAGCCCGGAActgctg                                                                                                      <  1:2619515/148‑1 (MQ=255)
                                             gtCACGGTTGATTGATATCACCACTGAAGCCGCCATTGATAGTGGCGTATTGCTGGAACGGTTGATTGCGCAAAACCCGGTTAGTTACAACTTCCATGTTCCGCTGGCTGATGGTGGCGTCCTGCTGGGGGCCAGCCCGGAActgctg                                                                                                      <  1:598441/148‑1 (MQ=255)
                                              tCACGGTTGATTGATATCACCACTGACGCCGCCATTGATAGTGGCGTATTGCTGGAACGGTTGATTGCGCAAAACCCGGTTAGTTACAACTTCCATGTTCCGCTGGCTGATGGTGGCGTCCTGCTGGGGGCCAGCCCGGAACtgctgc                                                                                                     <  2:2069029/148‑1 (MQ=255)
                                                aCGGTTGATTGATATCACCACTGACGCCGCCATTGATAGTGGCGTATTGCTGGAACGGTTGATTGCGCAAAACCCGGTTAGTTACAACTTCCATGTTCCGCTGGCTGATGGTGGCGTCCTGCTGGGGGCCAGCCCGGAACTGCTGCTa                                                                                                   <  2:489485/148‑1 (MQ=255)
                                                 cGGTTGATTGATATCACCACTGACGCCGCCATTGATAGTGGCGTATTGCTGGAACGGTTGATTGCGCAAAACCCGGTTAGTTACAACTTCCATGTTCCGCTGGCTGATGGTGGCGTCCTGCTGGGGGCCAGCCCGGAACTGCTGCTAc                                                                                                  >  1:815412/1‑148 (MQ=255)
                                                  ggTTGATTGATATCACCACTGACGCCGCCATTGATAGTGGCGTATTGCTGGAACGGTTGATTGCGCAAAACCCGGTTAGTTACAACTTCCATGTTCCGCTGGCTGATGGTGGCGTCCTGCTGGGGGCCAGCCCGGAACTGCTGCTACg                                                                                                 <  2:1320560/148‑1 (MQ=255)
                                                    ttgattgaTATCACCACTGACGCCGCCATTGATAGTGGCGTATTGCTGGAACGGTTGATTGCGCAAAACCCGGTTAGTTACAACTTCCATGTTCCGCTGGCTGATGGTGGCGTCCTGCTGGGGGCCAGCCCGGAACTGCTGCTACGTa                                                                                               >  1:2961457/1‑148 (MQ=255)
                                                    ttgattgaTATCACCACTGACGCCGCCATTGATAGTGGCGTATTGCTGGAACGGTTGATTGCGCAAAACCCGGTTAGTTACAACTTCCATGTTCCGCTGGCTGATGGTGGCGTCCTGCTGGGGGCCAGCCCGGAACTGCTGCTACGTa                                                                                               >  1:1498650/1‑148 (MQ=255)
                                                       attgatATCACCACTGACGCCGCCATTGATAGTGGCGTATTGCTGGAACGGTTGATTGCGCAAAACCCGGTTAGTTACAACTTCCATGTTCCGCTGGCTGATGGTGGCGTCCTGCTGGGGGCCAGCCCGGAACTGCTGCTACGTAAAg                                                                                            >  1:3550896/1‑148 (MQ=255)
                                                           atatCACCACTGACGCCGCCATTGATAGTGGCGTATTGCTGGAACGGTTGATTGCGCAAAACCCGGTTAGTTACAACTTCCATGTTCCGCTGGCTGATGGTGGCGTCCTGCTGTGGGCCAGCCCGGAACTGCTGCTACGTAAAGACgg                                                                                        <  2:2634599/148‑1 (MQ=255)
                                                           atatCACCACTGACGCCGCCATTGATAGTGGCGTATTGCTGGAACGGTTGATTGCGCAAAACCCGGTTAGTTACAACTTCCATGTTCCGCTGGCTGATGGTGGCGTCCTGCTGGGGGCCAGCCCGGAACTGCTGCTACGTAAAGACgg                                                                                        <  2:2988196/148‑1 (MQ=255)
                                                           atatCACCACTGACGCCGCCATTGATAGTGGCGTATTGCTGGAACGGTTGATTGCGCAAAACCCGGTTAGTTACAACTTCCATGTTCCGCTGGCTGATGGTGGCGTCCTGCTGGGGGCCAGCCCGGAACTGCTGCTACGTAAAGACgg                                                                                        <  2:1407086/148‑1 (MQ=255)
                                                           atatCACCACTGACGCCGCCATTGATAGTGGCGTATTGCTGGAACGGTTGATTGCGCAAAACCCGGTTAGTTACAACTTCCATGTTCCGCTGGCTGATGGTGGCGTCCTGCTGGGGGCCAGCCCGGAACTGCTGCTACGTAAAGACgg                                                                                        <  2:844636/148‑1 (MQ=255)
                                                              tCACCACTGACGCCGCCATTGATAGTGGCGTATTGCTGGAACGGTTGATTGCGCAAAACCCGGTTAGTTACAACTTCCATGTTCCGCTGGCTGATGGTGGCGTCCTGCTGGGGGCCAGCCCGGAACTGCTGCTACGTAAAGACGGCGa                                                                                     >  2:170646/1‑148 (MQ=255)
                                                              tCACCACTGACGCCGCCATTGATAGTGGCGTATTGCTGGAACGGTTGATTGCGCAAAACCCGGTTAGTTACAACTTCCATGTTCCGCTGGCTGATGGTGGCGTCCTGCTGGGGGCCAGCCCGGAACTGCTGCTACGTAAAGACGGCGa                                                                                     >  2:3007720/1‑148 (MQ=255)
                                                               caccacTGACGCCGCCATTGATAGTGGCGTATTGCTGGAACGGTTGATTGCGCAAAACCCGGTTAGTTACAACTTCCATGTTCCGCTGGCTGATGGTGGCGTCCTGCTGGGGGCCAGCCCGGAACTGCTGCTACGTAAAGACGGCGAg                                                                                    >  1:1233360/1‑148 (MQ=255)
                                                               caccacTGACGCCGCCATTGATAGTGGCGTATTGCTGGAACGGTTGATTGCGCAAAACCCGGTTAGTTACAACTTCCATGTTCCGCTGGCTGATGGTGGCGTCCTGCTGGGGGCCAGCCCGGAACTGCTGCTACGTAAAGACGGCGAg                                                                                    >  1:3112464/1‑148 (MQ=255)
                                                                accacTGACGCCGCCATTGATAGTGGCGTATTGCTGGAACGGTTGATTGCGCAAAACCCGGTTAGTTACAACTTCCATGTTCCGCTGGCTGATGGTGGCGTCCTGCTGGGGGCCAGCCCGGAACTGCTGCTACGTAAAGACGGCGAGc                                                                                   >  1:2117545/1‑148 (MQ=255)
                                                                 ccacTGACGCCGCCATTGATAGTGGCGTATTGCTGGAACGGTTGATTGCGCAAAACCCGGTTAGTTACAACTTCCATGTTCCGCTGGCTGATGGTGGCGTCCTGCTGGGGGCCAGCCCGGAACTGCTGCTACGTAAAGACGGCGAGCg                                                                                  <  2:510200/148‑1 (MQ=255)
                                                                 ccacTGACGCCGCCATTGATAGTGGCGTATTGCTGGAACGGTTGATTGCGCAAAACCCGGTTAGTTACAACTTCCATGTTCCGCTGGCTGATGGTGGCGTCCTGCTGGGGGCCAGCCCGGAACTGCTGCTACGTAAAGACGGCGAGCg                                                                                  >  1:2975970/1‑148 (MQ=255)
                                                                  cacTGACGCCGCCATTGATAGTGGCGTATTGCTGGAACGGTTGATTGCGCAAAACCCGGTTAGTTACAACTTCCATGTTCCGCTGGCTGATGGTGGCGTCCTGCTGGGGGCCAGCCCGGAACTGCTGCTACGTAAAGACGGCGAGCGt                                                                                 <  2:2679594/148‑1 (MQ=255)
                                                                      gACGCCGCCATTGATAGTGGCGTATTGCTGGAACGGTTGATTGCGCAAAACCCGGTTAGTTACAACTTCCATGTTCCGCTGGCTGATGGTGGCGTCCTGCTGGGGGCCAGCCCGGAACTGCTGCTACGTAAAGACGGCGAGCGTTTTa                                                                             >  1:2319170/1‑148 (MQ=255)
                                                                       aCGCCGCCATTGATAGTGGCGTATTGCTGGAACGGTTGATTGCGCAAAACCCGGTTAGTTACAACTTCCATGTTCCGCTGGCTGATGGTGGCGTCCTGCTGGGGGCCAGCCCGGAACTGCTGCTACGTAAAGACGGCGAGCGTTTTAg                                                                            <  1:2879954/148‑1 (MQ=255)
                                                                       aCGCCGCCATTGATAGTGGCGTATTGCTGGAACGGTTGATTGCGCAAAACCCGGTTAGTTACAACTTCCATGTTCCGCTGGCTGATGGTGGCGTCCTGCTGGGGGCCAGCCCGGAACTGCTGCTACGTAAAGACGGCGAGCGTTTTAg                                                                            <  1:459674/148‑1 (MQ=255)
                                                                         gccgccATTGATAGTGGCGTATTGCTGGAACGGTTGATTGCGCAAAACCCGGTTAGTTACAACTTCCATGTTCCGCTGGCTGATGGTGGCGTCCTGCTGGGGGCCAGCCCGGAACTGCTGCTACGTAAAGACGGCGAGCGTTTTAGCt                                                                          >  2:276843/1‑148 (MQ=255)
                                                                           cgccATTGATAGTGGCGTATTGCTGGAACGGTTGATTGCGCAAAACCCGGTTAGTTACAACTTCCATGTTCCGCTGGCTGATGGTGGCGTCCTGCTGGGGGCCAGCCCGGAACTGCTGCTACGTAAAGACGGCGAGCGTTTTAGCTcc                                                                        <  2:3313488/148‑1 (MQ=255)
                                                                            gccATTGATAGTGGCGTATTGCTGGAACGGTTGATTGCGCAAAACCCGGTTAGTTACAACTTCCATGTTCCGCTGGCTGATGGTGGCGTCCTGCTGGGGGCCAGCCCGGAACTGCTGCTACGTAAAGACGGCGAGCGTTTTAGCTCCa                                                                       <  2:2095139/148‑1 (MQ=255)
                                                                              cATTGATAGTGGCGTATTGCTGGAACGGTTGATTGCGCAAAACCCGGTTAGTTACAACTTCCATGTTCCGCTGGCTGATGGTGGCGTCCTGCTGGGGGCCAGCCCGGAACTGCTGCTACGTAAAGACGGCGAGCGTTTTAGCTCCAtt                                                                     <  1:2379049/148‑1 (MQ=255)
                                                                                ttGATAGTGGCGTATTGCTGGAACGGTTGATTGCGCAAAACCCGGTTAGTTACAACTTCCATGTTCCGCTGGCTGATGGTGGCGTCCTGCTGGGGGCCAGCCCGGAACTGCTGCTACGTAAAGACGGCGAGCGTTTTAGCTCCATTcc                                                                   >  2:1855667/1‑148 (MQ=255)
                                                                                  gATAGTGGCGTATTGCTGGAACGGTTGATTGCGCAAAACCCGGTTAGTTACAACTTCCATGTTCCGCTGGCTGATGGTGGCGTCCTGCTGGGGGCCAGCCCGGAACTGCTGCTACGTAAAGACGGCGAGCGTTTTAGCTCCATTCCGt                                                                 <  1:606033/148‑1 (MQ=255)
                                                                                  gATAGTGGCGTATTGCTGGAACGGTTGATTGCGCAAAACCCGGTTAGTTACAACTTCCATGTTCCGCTGGCTGATGGTGGCGTCCTGCTGGGGGCCAGCCCGGAACTGCTGCTACGTAAAGACGGCGAGCGTTTTAGCTCCATTCCGt                                                                 <  1:2746458/148‑1 (MQ=255)
                                                                                   aTAGTGGCGTATTGCTGGAACGGTTGATTGCGCAAAACCCGGTTAGTTACAACTTCCATGTTCTGCTGGCTGATGGTGGCGTCCTGCTGGGGGCCAGCCCGGAACTGCTGCTACGTAAAGACGGCGAGCGTTTTAGCTCCATTCCGtt                                                                <  2:3347689/148‑1 (MQ=255)
                                                                                   aTAGTGGCGTATTGCTGGAACGGTTGATTGCGCAAAACCCGGTTAGTTACAACTTCCATGTTCCGCTGGCTGATGGTGGCGTCCTGCTGGGGGCCAGCCCGGAACTGCTGCTACGTAAAGACGGCGAGCGTTTTAGCTCCATTCCGtt                                                                >  2:247090/1‑148 (MQ=255)
                                                                                    tAGTGGCGTATTGCTGGAACGGTTGATTGCGCAAAACCCGGTTAGTTACAACTTCCATGTTCCGCTGGCTGATGGTGGCGTCCTGCTGGGGGCCAGCCCGGAACTGCTGCTACGTAAAGACGGCGAGCGTTTTAGCTCCATTCCGTTa                                                               >  1:987445/1‑148 (MQ=255)
                                                                                       tGGCGTATTGCTGGAACGGTTGATTGCGCAAAACCCGGTTAGTTACAACTTCCATGTTCCGCTGGCTGATGGTGGCGTCCTGCTGGGGGCCAGCCCGGAACTGCTGCTACGTAAAGACGGCGAGCGTTTTAGCTCCATTCCGTTAGcc                                                            >  2:1338516/1‑148 (MQ=255)
                                                                                       tGGCGTATTGCTGGAACGGTTGATTGCGCAAAACCCGGTTAGTTACAACTTCCATGTTCCGCTGGCTGATGGTGGCGTCCTGCTGGGGGCCAGCCCGGAACTGCTGCTACGTAAAGACGGCGAGCGTTTTAGCTCCATTCCGTTAGcc                                                            >  2:1338479/1‑148 (MQ=255)
                                                                                           gTATTGCTGGAACGGTTGATTGCGCAAAACCCGGTTAGTTACAACTTCCATGTTCCGCTGGCTGATGGTGGCGTCCTGCTGGGGGCCAGCCCGGAACTGCTGCTACGTAAAGACGGCGAGCGTTTTAGCTCCATTCCGTTAGCCGGtt                                                        <  2:1116290/148‑1 (MQ=255)
                                                                                           gTATTGCTGGAACGGTTGATTGCGCAAAACCCGGTTAGTTACAACTTCCATGTTCCGCTGGCTGAGGGTGGCGTCCTGCTGGGGGCCAGCCCGGAACTGCTGCTACGTAAAGACGGCGAGCGTTTTAGCTCCATTCCGTTAGCCGGtt                                                        <  2:1116151/148‑1 (MQ=255)
                                                                                            tATTGCTGGAACGGTTGATTGCGCAAAACCCGGTTAGTTACAACTTCCATGTTCCGCTGGCTGATGGTGGCGTCCTGCTGGGGGCCAGCCCGGAACTGCTGCTACGTAAAGACGGCGAGCGTTTTAGCTCCATTCCGTTAGCCGGTTc                                                       <  1:1380454/148‑1 (MQ=255)
                                                                                                  tGGAACGGTTGATTGCGCAAAACCCGGTTAGTTACAACTTCCATGTTCTGCTGGCTGATGGTGGCGTCCTGCTGGGGGCCAGCCCGGAACTGCTGCTACGTAAAGACGGCGAGCGTTTTAGCTCCATTCCGTTAGCCGGTTCcgcgcg                                                 <  1:594340/148‑1 (MQ=255)
                                                                                                   ggAACGGTTGATTGCGCAAAACCCGGTTAGTTACAACTTCCATGTTCCGCTGGCTGATGGTGGCGTCCTGCTGGGGGCCAGCCCGGAACTGCTGCTACGTAAAGACGGCGAGCGTTTTAGCTCCATTCCGTTAGCCGGTTCCGCGcgt                                                <  1:2355199/148‑1 (MQ=255)
                                                                                                    gAACGGTTGATTGCGCAAAACCCGGTTAGTTACAACTTCCATGTTCCGCTGGCTGATGGTGGCGTCCTGCTGGGGGCCAGCCCGGAACTGCTGCTACGTAAAGACGGCGAGCGTTTTAGCTCCATTCCGTTAGCCGGTTCCGCGcgtc                                               <  1:1462599/148‑1 (MQ=255)
                                                                                                     aaCGGTTGATTGCGCAAAACCCGGTTAGTTACAACTTCCATGTTCCGCTGGCTGATGGTGGCGTCCTGCTGGGGGCCAGCCCGGAACTGCTGCTACGTAAAGACGGCGAGCGTTTTAGCTCCATTCCGTTAGCCGGTTCCGCGcgtcg                                              <  1:407024/148‑1 (MQ=255)
                                                                                                     aaCGGTTGATTGCGCAAAACCCGGTTAGTTACAACTTCCATGTTCCGCTGGCTGATGGTGGCGTCCTGCTGGGGGCCAGCCCGGAACTGCTGCTACGTAAAGACGGCGAGCGTTTTAGCTCCATTCCGTTAGCCGGTTCCGCGcgtcg                                              <  2:1111183/148‑1 (MQ=255)
                                                                                                      ccGGTTGATTGCGCAACACCCGGTTAGTTACAACTTCCATGTTCCGCTGGCTGATGGTGGCGTCCTGCTGGGGGCCAGCCCGGAACTGCTGCTACGTAAAGACGGCGAGCGTTTTAGCTCCATTCCGTTAGCCGGTTCCGCGcgtcgt                                             <  2:2056369/147‑1 (MQ=255)
                                                                                                          ttGATTGCGCAAAACCCGGTTAGTTACAACTTCCATGTTCCGCTGGCTGATGGTGGCGTCCTGCTGGGGGCCAGCCCGGAACTGCTGCTACGTAAAGACGGCGAGCGTTTTAGCTCCATTCCGTTAGCCGGTTCCGCGCGTCGTCAGc                                         >  2:2231717/1‑148 (MQ=255)
                                                                                                               tGCGCAAACCCCGGTTAGTTACAACTTCCATGTTCCGCTGGCTGATGGTGGCGTCCTGCTGGGGGCCAGCCCGGAACTGCTGCTACGTAAAGACGGCGAGCGTTTTAGCTCCATTCCGTTAGCCGGTTCCGCGCGTCGTCAGCCGGAt                                    <  2:413069/148‑1 (MQ=255)
                                                                                                               tGCGCAAAACCCGGTTAGTTACAACTTCCATGTTCCGCTGGCTGATGGTGGCGTCCTGCTGGGGGCCAGCCCGGAACTGCTGCTACGTAAAGACGGCGAGCGTTTTAGCTCCATTCCGTTAGCCGGTTCCGCGCGTCGTCAGCCGGAt                                    <  2:2210724/148‑1 (MQ=255)
                                                                                                                     aaaCCCGGTTAGTTACAACTTCCATGTTCCGCTGGCTGATGGTGGCGTCCTGCTGGGGGCCAGCCCGGAACTGCTGCTACGTAAAGACGGCGAGCGTTTTAGCTCCATTCCGTTAGCCGGTTCCGCGCGTCGTCAGCCGGATGAAGTg                              >  2:2353315/1‑148 (MQ=255)
                                                                                                                     aaaCCCGGTTAGTTACAACTTCCATGTTCCGCTGGCTGATGGTGGCGTCCTGCTGGGGGCCAGCCCGGAACTGCTGCTACGTAAAGACGGCGAGCGTTTTAGCTCCATTCCGTTAGCCGGTTCCGCGCGTCGTCAGCCGGATGAAGTg                              <  2:2599480/148‑1 (MQ=255)
                                                                                                                     aaaCCCGGTTAGTTACAACTTCCATGTTCCGCTGGCTGATGGTGGCGTCCTGCTGGGGGCCAGCCCGGAACTGCTGCTACGTAAAGACGGCGAGCGTTTTAGCTCCATTCCGTTAGCCGGTTCCGCGCGTCGTCAGCCGGATGAAGTg                              >  1:2607319/1‑148 (MQ=255)
                                                                                                                       aCCCGGTTAGTTACAACTTCCATGTTCTGCTGGCTGATGGTGGCGTCCTGCTGGGGGCCAGCCCGGAACTGCTGCTACGTAAAGACGGCGAGCGTTTTAGCTCCATTCCGTTAGCCGGTTCCGCGCGTCGTCAGCCGGATGAAGTGCt                            <  2:3076248/148‑1 (MQ=255)
                                                                                                                       aCCCGGTTAGTTACAACTTCCATGTTCCGCTGGCTGATGGTGGCGTCCTGCTGGGGGCCAGCCCGGAACTGCTGCTACGTAAAGACGGCGAGCGTTTTAGCTCCATTCCGTTAGCCGGTTCCGCGCGTCGTCAGCCGGATGAAGTGCt                            >  1:36846/1‑148 (MQ=255)
                                                                                                                        cccGGTTAGTTACAACTTCCATGTTCCGCTGGCTGATGGTGGCGTCCTGCTGGGGGCCAGCCCGGAACTGCTGCTACGTAAAGACGGCGAGCGTTTTAGCTCCATTCCGTTAGCCGGTTCCGCGCGTCGTCAGCCGGATGAAGTGCTc                           >  1:25413/1‑148 (MQ=255)
                                                                                                                           ggTTAGTTACAACTTCCATGTTCCGCTGGCTGATGGTGGCGTCCTGCTGGGGGCCAGCCCGGAACTGCTGCTACGTAAAGACGGCGAGCGTTTTAGCTCCATTCCGTTAGCCGGTTCCGCGCGTCGTCAGCCGGATGAAGTGCTCGAt                        <  1:3256232/148‑1 (MQ=255)
                                                                                                                            gttagttaCAACTTCCATGTTCCGCTGGCTGATGGTGGCGTCCTGCTGGGGGCCAGCCCGGAACTGCTGCTACGTAAAGACGGCGAGCGTTTTAGCTCCATTCCGTTAGCCGGTTCCGCGCGTCGTCAGCCGGATGAAGTGCTCGATc                       >  2:2544375/1‑148 (MQ=255)
                                                                                                                              tagttaCAACTTCCATGTTCCGCTGGCTGATGGTGGCGTCCTGCTGGGGGCCAGCCCGGAACTGCTGCTACGTAAAGACGGCGAGCGTTTTAGCTCCATTCCGTTAGCCGGTTCCGCGCGTCGTCAGCCGGATGAAGTGCTCGATcgc                     >  2:1761827/1‑148 (MQ=255)
                                                                                                                                 ttaCAACTTCCATGTTCCGCTGGCTGATGGTGGCGTCCTGCTGGGGGCCAGCCCGGAACTGCTGCTACGTAAAGACGGCGAGCGTTTTAGCTCCATTCCGTTAGCCGGTTCCGCGCGTCGTCAGCCGGATGAAGTGCTCGATCGCGaa                  >  1:3183647/1‑148 (MQ=255)
                                                                                                                                    cAACTTCCATGTTCCGCTGGCTGATGGTGGCGTCCTGCTGGGGGCCAGCCCGGAACTGCTGCTACGTAAAGACGGCGAGCGTTTTAGCTCCATTCCGTTAGCCGGTTCCGCGCGTCGTCAGCCGGATGAAGTGCTCGATCGCGAAGCa               <  2:2185437/148‑1 (MQ=255)
                                                                                                                                      aCTTCCATGTTCCGCTGGCTGATGGTGGCGTCCTGCTGGGGGCCAGCCCGGAACTGCTGCTACGTAAAGACGGCGAGCGTTTTAGCTCCATTCCGTTAGCCGGTTCCGCGCGTCGTCAGCCGGATGAAGTGCTCGATCGCGAAGCAgg             >  1:2433888/1‑148 (MQ=255)
                                                                                                                                        ttCCATGTTCCGCTGGCTGATGGTGGCGTCCTGCTGGGGGCCAGCCCGGAACTGCTGCTACGTAAAGACGGCGAGCGTTTTAGCTCCATTCCGTTAGCCGGTTCCGCGCGTCGTCAGCCGGATGAAGTGCTCGATCGCGAAGCAGGTa           >  2:1975078/1‑148 (MQ=255)
                                                                                                                                            aTGTTCCGCTGGCTGATGGTGGCGTCCTGCTGGGGGCCAGCCCGGAACTGCTGCTACGTAAAGACGGCGAGCGTTTTAGCTCCATTCCGTTAGCCGGTTCCGCGCGTCGTCAGCCGGATGAAGTGCTCGATCGCGAAGCAGGTAATCg       >  1:2196696/1‑148 (MQ=255)
                                                                                                                                             tGTTCCGCTGGCTGATGGTGGCGTCCTGCTGGGGGCCAGCCCGGAACTGCTGCTACGTAAAGACGGCGAGCGTTTTAGCTCCATTCCGTTAGCCGGTTCCGCGCGTCGTCAGCCGGATGAAGTGCTCGATCGCGAAGCAGGTAATCGt      >  2:2559457/1‑148 (MQ=255)
                                                                                                                                             tGTTCCGCTGGCTGATGGTGGCGTCCTGCTGGGGGCCAGCCCGGAACTGCTGCTACGTAAAGACGGCGAGCGTTTTAGCTCCATTCCGTTAGCCGGTTCCGCGCGTCGTCAGCCGGATGAAGTGCTCGATCGCGAAGCAGGTAATCGt      <  1:2270269/148‑1 (MQ=255)
                                                                                                                                              gTTCCGCTGGCTGATGGTGGCGTCCTGCTGGGGGCCAGCCCGGAACTGCTGCTACGTAAAGACGGCGAGCGTTTTAGCTCCATTCCGTTAGCCGGTTCCGCGCGTCGTCAGCCGGATGAAGTGCTCGATCGCGAAGCAGGTAATCGTc     >  1:1882372/1‑148 (MQ=255)
                                                                                                                                              gTTCCGCTGGCTGATGGTGGCGTCCTGCTGGGGGCCAGCCCGGAACTGCTGCTACGTAAAGACGGCGAGCGTTTTAGCTCCATTCCGTTAGCCGGTTCCGCGCGTCGTCAGCCGGATGAAGTGCTCGATCGCGAAGCAGGTAATCGTc     >  1:1016589/1‑148 (MQ=255)
                                                                                                                                                tCCGCTGGCTGATGGTGGCGTCCTGCTGGGGGCCAGCCCGGAACTGCTGCTACGTAAAGACGGCGAGCGTTTTAGCTCCATTCCGTTAGCCGGTTCCGCGCGTCGTCAGCCGGATGAAGTGCTCGATCGCGAAGCAGGTAATCGTctg   >  1:1570426/1‑148 (MQ=255)
                                                                                                                                                tCCGCTGGCTGATGGTGGCGTCCTGCTGGGGGCCAGCCCGGAACTGCTGCTACGTAAAGACGGCGAGCGTTTTAGCTCCATTCCGTTAGCCGGTTCCGCGCGTCGTCAGCCGGATGAAGTGCTCGATCGCGAAGCAGGTAATCGTctg   <  1:1855667/148‑1 (MQ=255)
                                                                                                                                                tCCGCTGGCTGATGGTGGCGTCCTGCTGGGGGCCAGCCCGGAACTGCTGCTACGTAAAGACGGCGAGCGTTTTAGCTCCATTCCGTTAGCCGGTTCCGCGCGTCGTCAGCCGGATGAAGTGCTCGATCGCGAAGCAGGTAATCGTctg   >  2:192301/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
CGCCGCCGCACTTACCGCCACGCCGCAGGTCGACAAAGTGGTGTTGTCACGGTTGATTGATATCACCACTGACGCCGCCATTGATAGTGGCGTATTGCTGGAACGGTTGATTGCGCAAAACCCGGTTAGTTACAACTTCCATGTTCCGCTGGCTGATGGTGGCGTCCTGCTGGGGGCCAGCCCGGAACTGCTGCTACGTAAAGACGGCGAGCGTTTTAGCTCCATTCCGTTAGCCGGTTCCGCGCGTCGTCAGCCGGATGAAGTGCTCGATCGCGAAGCAGGTAATCGTCTGC  >  CP009273/620745‑621037

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