Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 488,177 G→C 5.3% A43G (GCC→GGC)  pdxA ← 4‑hydroxythreonine‑4‑phosphate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947488,1770GC 5.3% 1082.7 / 36.1 318A43G (GCC→GGC) pdxA4‑hydroxythreonine‑4‑phosphate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (162/139);  new base C (8/9);  total (170/148)
Fisher's exact test for biased strand distribution p-value = 6.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.50e-01

Only 100 of 322 total aligned reads displayed.

GGCCTTGTCGAGCTGGCCTGGCACCACCGGGGCCGCCAAGGGGGTATCCCAGACGTACAGGCTACCTGCGGGTGCCGGGAGGTCAGGCCATTGGCCCGGGGCTACCGGCAGCAGGCTGACGGCCAGCCCCAGCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCTGGTGATGGCGATCAGGGGGTGGGGCTGGGCGTCTGCGGCGAGCAGCAGGCACAGGTCGGGGCCTATGCCGGCCGGCTCGCCGGGGGTGACGGCGAAGCGCAGGGGCTTCACTGGGCGGCCTGGTCAGCGCCAG  >  NC_002947/488029‑488326
                                                                                                                                                    |                                                                                                                                                     
  ccTTGTCGAGCTGGCCTGGCACCACCGGGGCCGCCAAGGGGGTATCCCAGACGTACAGGCTACCTGCGGGTGCCGGGAGGTCAGGCCATTGGCCCGGGGCTACCGGCAGCAGGCTGACGGCCAGCCCCAGCTGCGTGGCCCGCTCGGCGa                                                                                                                                                    <  2:4340552/150‑1 (MQ=255)
  ccTTGTCGAGCTGGCCTGGCACCACCGGGGCCGCCAAGGGGGTATCCCAGACGTACAGGCTACCTGCGGGTGCCGGGAGGTCAGGCCATTGGCCCGGGGCTACCGGCAGCAGGCTGACGGCCAGCCCCAGCTGCGTGGCCCGCTCGGCGa                                                                                                                                                    >  2:732479/1‑150 (MQ=255)
  ccTTGTCGAGCTGGCCTGGCACCACCGGGGCCGCCAAGGGGGTATCCCAGACGTACAGGCTACCTGCGGGTGCCGGGAGGTCAGGCCATTGGCCCGGGGCTACCGGCAGCAGGCTGACGGCCAGCCCCAGCTGCGTGGCCCGCTCGGCGa                                                                                                                                                    <  1:338925/150‑1 (MQ=255)
      gTCGAGCTGGCCTGGCACCACCGGGGCCGCCAAGGGGGTATCCCAGACGTACAGGCTACCTGCGGGTGCCGGGAGGTCAGGCCATTGGCCCGGGGCTACCGGCAGCAGGCTGACGGCCAGCCCCAGCTGCGTGGCCCGCTCGGCGAGCAg                                                                                                                                                <  1:3551076/150‑1 (MQ=255)
          aGCTGGCCTGGCACCACCGGGGCCGCCAAGGGGGTATCCCAGACGTACAGGCTACCTGCGGGTGCCGGGAGGTCAGGCCATTGGCCCGGGGCTACCGGCAGCAGGCTGACGGCCAGCCCCAGCTGCGTGACCTGCTCGCCGAGCGGGCCa                                                                                                                                            <  2:2813228/150‑1 (MQ=255)
           gCTGGCCTGGCACCACCGGGGCCGCCAAGGGGGTATCCCAGACGTACAGGCTACCTGCGGGTGCCGGGAGGTCAGGCCATTGGCCCGGGGCTACCGGCAGCAGGCTGACGGCCAGCCCCAGCTGCGTGGCCTGCTCGCCGAGCGGGCCAc                                                                                                                                           <  1:2267344/150‑1 (MQ=255)
            ctggcctggcACCACCGGGGCCGCCAAGGGGGTATCCCAGACGTACAGGCTACCTGCGGGTGCCGGGAGGTCAGGCCATTGGCCCGGGGCTACCGGCAGCAGGCTGACGGCCAGCCCCAGCTGCGTGACCTGCTCGCCGAGCGGGCCACg                                                                                                                                          <  2:4703609/150‑1 (MQ=255)
              ggcctggcACCACCGGGGCCGCCAAGGGGGTATCCCAGACGTACAGGCTACCTGCGGGTGCCGGGAGGTCAGGCCATTGGCCCGGGGCTACCGGCAGCAGGCTGACGGCCAGCCCCAGCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGc                                                                                                                                        <  1:3089254/150‑1 (MQ=255)
              ggcctggcACCACCGGGGCCGCCAAGGGGGTATCCCAGACGTACAGGCTACCTGCGGGTGCCGGGAGGTCAGGCCATTGGCCCGGGGCTACCGGCAGCAGGCTGACGGCCAGCCCCAGCTGCGTGACCTGCTCGCCGAGCGGGCCACgca                                                                                                                                        <  1:1294037/150‑3 (MQ=255)
               gcctggcACCACCGGGGCCGCCAAGGGGGTATCCCAGACGTACAGGCTACCTGCGGGTGCCGGGAGGTCAGGCCATTGGCCCGGGGCTACCGGCAGCAGGCTGACGGCCAGCCCCAGCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCt                                                                                                                                       <  1:2663374/150‑1 (MQ=255)
                   ggcACCACCGGGGCCGCCAAGGGGGTATCCCAGACGTACAGGCTACCTGCGGGTGCCGGGAGGTCAGGCCATTGGCCCGGGGCTACCGGCAGCAGGCTGACGGCCAGCCCCAGCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCTGGTg                                                                                                                                   <  1:1430857/150‑1 (MQ=255)
                   ggcACCACCGGGGCCGCCAAGGGGGTATCCCAGACGTACAGGCTACCTGCGGGTGCCGGGAGGTCAGGCCATTGGCCCGGGGCTACCGGCAGCAGGCTGACGGCCAGCCCCAGCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCTGGTg                                                                                                                                   <  1:6291021/150‑1 (MQ=255)
                      accaccGGGGCCGCCAAGGGGGTATCCCAGACGTACAGGCTACCTGCGGGTGCCGGGAGGTCAGGCCATTGGCCCGGGGCTACCGGCAGCAGGCTGACGGCCAGCCCCAGCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCTGGTGATg                                                                                                                                >  2:3712931/1‑150 (MQ=255)
                       ccaccGGGGCCGCCCAGGGGGTATCCCAGACGTACAGGCTACCTGCGGGTGCCGGGAGGTCAGGCCATTGGCCCGGGGCTACCGGCAGCAGGCTGACGGCCAGCCCCAGCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCTGGTGATgg                                                                                                                               >  2:5335856/1‑150 (MQ=255)
                       ccaccGGGGCCGCCAAGGGGGTATCCCAGACGTACAGGCTACCTGCGGGTGCCGGGAGGTCAGGCCATTGGCCCGGGGCTACCGGCAGCAGGCTGACGGCCAGCCCCAGCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCTGGTGATgg                                                                                                                               <  2:1972750/150‑1 (MQ=255)
                       ccaccGGGGCCGCCAAGGGGGTATCCCAGACGTACAGGCTACCTGCGGGTGCCGGGAGGTCAGGCCATTGGCCCGGGGCTACCGGCAGCAGGCTGACGGCCAGCCCCAGCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCTGGTGATgg                                                                                                                               <  1:1187166/150‑1 (MQ=255)
                         accGGGGCCGCCAAGGGGGTATCCCAGACGTACAGGCTACCTGCGGGTGCCGGGAGGTCAGGCCATTGGCCCGGGGCTACCGGCAGCAGGCTGACGGCCAGCCCCAGCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCTGGTGATGGCg                                                                                                                             <  2:1694980/150‑1 (MQ=255)
                         accGGGGCCGCCAAGGGGGTATCCCAGACGTACAGGCTACCTGCGGGTGCCGGGAGGTCAGGCCATTGGCCCGGGGCTACCGGCAGCAGGCTGACGGCCAGCCCCAGCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCTGGTGATGGCg                                                                                                                             >  2:6338072/1‑150 (MQ=255)
                           cGGGGCCGCCAAGGGGGTATCCCAGACGTACAGGCTACCTGCGGGTGCCGGGAGGTCAGGCCCTTGGCCCGGGGCTACCGGCAGCAGGCTGACGGCCAGCCCCAGCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCTGGTGATGGCGAt                                                                                                                           >  2:3828476/1‑150 (MQ=255)
                             gggCCGCCAAGGGGGTATCCCAGACGTACAGGCTACCTGCGGGTGCCGGGAGGTCAGGCCATTGGCCCGGGGCTACCGGCAGCAGGCTGACGGCCAGCCCCAGCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCTGGTGATGGCGATCa                                                                                                                         >  1:4695120/1‑150 (MQ=255)
                             gggCCGCCAAGGGGGTATCCCAGACGTACAGGCTACCTGCGGGTGCCGGGAGGTCAGGCCATTGGCCCGGGGCTACCGGCAGCAGGCTGACGGCCAGCCCCAGCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCTGGTGATGGCGATCa                                                                                                                         >  2:621268/1‑150 (MQ=255)
                               gccgccAAGGGGGTATCCCAGACGTACAGGCTACCTGCGGGTGCCGGGAGGTCAGGCCATTGGCCCGGGGCTACCGGCAGCAGGCTGACGGCCAGCCCCAGCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCTGGTGATGGCGATCAgg                                                                                                                       <  2:900560/150‑1 (MQ=255)
                                   ccAAGGGGGTATCCCAGACGTACAGGCTACCTGCGGGTGCCGGGAGGTCAGGCCATTGGCCCGGGGCTACCGGCAGCAGGCTGACGGCCAGCCCCAGCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCTGGTGATGGCGATCAGGGGGt                                                                                                                   <  1:199877/150‑1 (MQ=255)
                                     aaGGGGGTATCCCAGACGTACAGGCTACCTGCGGGTGCCGGGAGGTCAGGCCATTGGCCCGGGGCTACCGGCAGCAGGCTGACGGCCAGCCCCAGCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCTGGTGATGGCGATCAGGGGGTgg                                                                                                                 <  2:3277938/150‑1 (MQ=255)
                                      aGGGGGTATCCCAGACGTACAGGCTACCTGCGGGTGCCGGGAGGTCAGGCCATTGGCCCGGGGCTACCGGCAGCAGGCTGACGGCCAGCCCCAGCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCTGGTGATGGCGATCAGGGGGTggg                                                                                                                <  2:5627447/150‑1 (MQ=255)
                                        ggggTATCCCAGACGTACAGGCTACCTGCGGGTGCCGGGAGGTCAGGCCATTGGCCCGGGGCTACCGGCAGCAGGCTGACGGCCAGCCCCAGCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCTGGTGATGGCGATCAGGGGGTGGGGc                                                                                                              >  2:5779706/1‑150 (MQ=255)
                                        ggggTATCCCAGACGTACAGGCTACCTGCGGGTGCCGGGAGGTCAGGCCATTGGCCCGGGGCTACCGGCAGCAGGCTGACGGCCAGCCCCAGCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCTGGTGATGGCGATCAGGGGGTGGGGc                                                                                                              >  2:1843571/1‑150 (MQ=255)
                                          ggTATCCCAGACGTACAGGCTACCTGCGGGTGCCGGGAGGTCAGGCCATTGGCCCGGGGCTACCGGCAGCAGGCTGACGGCCAGCCCCAGCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCTGGTGATGGCGATCAGGGGGTGGGGCTg                                                                                                            >  2:2623806/1‑150 (MQ=255)
                                           gTATCCCAGACGTACAGGCTACCTGCGGGTGCCGGGAGGTCAGGCCATTGGCCCGGGGCTACCGGCAGCAGGCTGACGGCCAGCCCCAGCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCTGGTGATGGCGATCAGGGGGTGGGGCTgg                                                                                                           <  2:1464326/150‑1 (MQ=255)
                                           gTATCCCAGACGTACAGGCTACCTGCGGGTGCCGGGAGGTCAGGCCATTGGCCCGGGGCTACCGGCAGCAGGCTGACGGCCAGCCCCAGCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCTGGTGATGGCGATCAGGGGGTGGGGCTgg                                                                                                           <  1:1126341/150‑1 (MQ=255)
                                            tATCCCAGACGTACAGGCTACCTGCGGGTGCCGGGAGGTCAGGCCATTGGCCCGGGGCTACCGGCAGCAGGCTGACGGCCAGCCCCAGCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCTGGTGATGGCGATCAGGGGGTGGGGCTggg                                                                                                          >  1:4015611/1‑150 (MQ=255)
                                             aTCCCAGACGTACAGGCTACCTGCGGGTGCCGGGAGGTCAGGCCATTGGCCCGGGGCTACCGGCAGCAGGCTGACGGCCAGCCCCAGCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCGGGTGATGGCGATCAGGGGGGGGGGCTGGGc                                                                                                         >  2:2999623/1‑150 (MQ=255)
                                              tCCCAGACGTACAGGCTACCTGCGGGTGCCGGGAGGTCAGGCCATTGGCCCGGGGCTACCGGCAGCAGGCTGACGGCCAGCCCCAGCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCTGGTGATGGCGATCAGGGGGTGGGGCTGGGCg                                                                                                        >  2:5560852/1‑150 (MQ=255)
                                               cccAGACGTACAGGCTACCTGCGGGTGCCGGGAGGTCAGGCCATTGGCCCGGGACTACCGGCAGCAGGCTGACGGCCAGCCCCAGCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCTGGTGATGGCGATCAGGGGGTGGGGCTGGGCGt                                                                                                       >  1:4204436/1‑150 (MQ=255)
                                                ccAGACGTACAGGCTACCTGCGGGTGCCGGGAGGTCAGGCCATTGGCCCGGGGCTACCGGCAGCAGGCTGACGGCCAGCCCCAGCTGCGTGGCCTGCTCGCCGAGcgg                                                                                                                                                <  2:4001480/108‑3 (MQ=255)
                                                ccAGACGTACAGGCTACCTGCGGGTGCCGGGAGGTCAGGCCATTGGCCCGGGGCTACCGGCAGCAGGCTGACGGCCAGCCCCAGCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCTGGTg                                                                                                                                   <  2:2034362/121‑1 (MQ=255)
                                                ccAGACGTACAGGCTACCTGCGGGTGCCGGGAGGTCAGGCCATTGGCCCGGGGCTACCGGCAGCAGGCTGACGGCCAGCCCCAGCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCTGGTGATGGCGATCAGGGGGTGGGGCTGGGCGTc                                                                                                      >  1:4293367/1‑150 (MQ=255)
                                                ccAGACGTACAGGCTACCTGCGGGTGCCGGGAGGTCAGGCCATTGGCCCGGGGCTACCGGCAGCAGGCTGACGGCCAGCCCCAGCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCTGGTGATGGCGATCAGGGGGTGGGGCTGGGCGTc                                                                                                      >  2:3205193/1‑150 (MQ=255)
                                                ccAGACGTACAGGCTACCTGCGGGTGCCGGGAGGTCAGGCCATTGGCCCGGGGCTACCGGCAGCAGGCTGACGGCCAGCCCCAGCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCTGGTGATGGCGATCAGGGGGTGGGGCTGGGCGTc                                                                                                      >  1:5904695/1‑150 (MQ=255)
                                                 cAGACGTACAGGCTACCTGCGGGTGCCGGGAGGTCAGGCCATTGGCCCGGGGCTACCGGCAGCAGGCTGACGGCCAGCCCCAGCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCTGGTGATGGCGATCAGGGGGTGGGGCTGGGCGTCt                                                                                                     <  1:4907601/150‑1 (MQ=255)
                                                  aGACGTACAGGCTACCTGCGGGTGCCGGGAGGTCAGGCCATTGGCCCGGGGCTACCGGCAGCAGGCTGACGGCCAGCCCCAGCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCTGGTGATGGCGATCAGGGGGGGGGGCTGGGCGTCTg                                                                                                    >  2:3282501/1‑150 (MQ=255)
                                                        aCAGGCTACCTGCGGGTGCCGGGAGGTCAGGCCATTGGCCCGGGGCTACCGGCAGCAGGCTGACGGCCAGCCCCAGCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCTGGTGATGGCGATCAGGGGGTGGGGCTGGGCGTCTGCGGCGa                                                                                              <  2:4095422/150‑1 (MQ=255)
                                                        aCAGGCTACCTGCGGGTGCCGGGAGGTCAGGCCATTGGCCCGGGGCTACCGGCAGCAGGCTGACGGCCAGCCCCAGCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCTGGTGATGGCGATCAGGGGGTGGGGCTGGGCGTCTGCGGCGa                                                                                              <  1:2876941/150‑1 (MQ=255)
                                                         cAGGCTACCTGCGGGTGCCGGGAGGTCAGGCCATTGGCCCGGGGCTACCGGCAGCAGGCTGACGGCCAGCCCCAGCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCTGGTGATGGCGATCAGGGGGTGGGGCTGGGCGTCTGCGGCGag                                                                                             >  1:737579/1‑150 (MQ=255)
                                                          aGGCTACCTGCGGGTGCCGGGAGGTCAGGCCATTGGCCCGGGGCTACCGGCAGCAGGCTGACGGCCAGCCCCAGCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCTGGTGATGGCGATCAGGGGGTGGGGCTGGGCGTCTGCGGCGagc                                                                                            <  1:327782/150‑1 (MQ=255)
                                                          aGGCTACCTGCGGGTGCCGGGAGGTCAGGCCATTGGCCCGGGGCTACCGGCAGCAGGCTGACGGCCAGCCCCAGCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCTGGTGATGGCGATCAGGGGGTGGGGCTGGGCGTCTGCGGCGagc                                                                                            >  1:6026472/1‑150 (MQ=255)
                                                          aGGCTACCTGCGGGTGCCGGGAGGTCAGGCCATTGGCCCGGGGCTACCGGCAGCAGGCTGACGGCCAGCCCCAGCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCTGGTGATGGCGATCAGGGGGTGGGGCTGGGCGTCTGCGGCGagc                                                                                            <  1:624241/150‑1 (MQ=255)
                                                           ggCTACCTGCGGGTGCCGGGAGGTCAGGCCATTGGCCCGGGGCTACCGGCAGCAGGCTGACGGCCAGCCCCAGCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCTGGTGATGGCGATCAGGGGGTGGGGCTGGGCGTCTGCGGCGagca                                                                                           <  2:4587736/150‑1 (MQ=255)
                                                            gCTACCTGCGGGTGCCGGGAGGTCAGGCCATTGGCCCGGGGCTACCGGCAGCAGGCTGACGGCCAGCCCCAGCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCTGGTGATGGCGAt                                                                                                                           >  2:5153630/1‑117 (MQ=255)
                                                              tACCTGCGGGTGCCGGGAGGTCAGGCCATTGGCCCGGGGCTACCGGCAGCAGGCTGACGGCCAGCCCCAGCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCTGGTGATGGCGATCAGGGGGTGGGGCTGGGCGTCTGCGGCGagccgcc                                                                                        >  2:4835846/1‑149 (MQ=255)
                                                                ccTGCGGGTGCCGGGAGGTCAGGCCATTGGCCCGGGGCTACCGGCAGCAGGCTGACGGCCAGCCCCAGCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCTGGTGATGGCGATCAGGGGGTGGGGCTGGGCGTCTGCGGCGAGCAGCAgg                                                                                      <  1:4827931/150‑1 (MQ=255)
                                                                 cTGCGGGTGCCGGGAGGTCAGGCCATTGGCCCGGGGCTACCGGCAGCAGGCTGACGGCCAGCCCCAGCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCTGGTGATGGCGATCAGGGGGTGGGGCTGGGCGTCTGCGGCGAGCAGCAGGc                                                                                     >  1:4420290/1‑150 (MQ=255)
                                                                   gCGGGTGCCGGGAGGTCAGGCCATTGGCCCGGGGCTACCGGCAGCAGGCTGACGGCAAGCCCCAGCTGCGTGGCCCGCTCGGCGAGCAGGTCACGACTGGTGATGGCGATCAGGGGGTGGGGCTGGGCGACTGCGGCGAGCAGCAGGcac                                                                                   >  1:90858/1‑150 (MQ=255)
                                                                     gggTGCCGGGAGGTCAGGCCATTGGCCCGGGGCTACCGGCAGCAGGCTGACGGCCAGCCCCAGCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCTGGTGATGGCGATCAGGGGGTGGGGCTGGGCGTCTGCGGCGAGCAGCAGGCACAg                                                                                 <  2:3425989/150‑1 (MQ=255)
                                                                       gTGCCGGGAGGTCAGGCCATTGGCCCGGGGCTACCGGCAGCAGGCTGACGGCCAGCCCCAGCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCTGGTGATGGCGATCAGGGGGTGGGGCTGGGCGTCTGCGGCGAGCAGCAGGCACAGGt                                                                               >  2:522095/1‑150 (MQ=255)
                                                                        tGCCGGGAGGTCAGGCCATTGGCCCGGGGCTACCGGCAGCAGGCTGACGGCCAGCCCCAGCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCTGGTGATGGCGATCAGGGGGTGGGGCTGGGCGTCTGCGGCGAGCAGCAGGCACAGGTc                                                                              <  1:5240282/150‑1 (MQ=255)
                                                                        tGCCGGGAGGTCAGGCCATTGGCCCGGGGCTACCGGCAGCAGGCTGACGGCCAGCCCCAGCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCTGGTGATGGCGATCAGGGGGTGGGGCTGGGCGTCTGCGGCGAGCAGCAGGCACAGGTc                                                                              <  2:2216149/150‑1 (MQ=255)
                                                                           cGGGAGGTCAGGCCATTGGCCCGGGGCTACCGGCAGCAGGCTGACGGCCAGCCCCAGCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCTGGTGATGGCGATCAGGGGGTGGGGCTGGGCGTCTGCGGCGAGCAGCAGGCACAGGTCggg                                                                           >  2:2415712/1‑150 (MQ=255)
                                                                           cGGGAGGTCAGGCCATTGGCCCGGGGCTACCGGCAGAAGGCTGACGGCCAGCCCCAGCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCTGGTGATGGCGATCAGGGGGTGGGGCTGGGCGTCTGCGGCGAGCAGCAGGCACAGGTCggg                                                                           <  1:6074699/150‑1 (MQ=255)
                                                                            gggAGGTCAGGCCATTGGCCCGGGGCTACCGGCAGCAGGCTGACGGCCAGCCCCAGCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCTGGTGATGGCGATCAGGGGGTGGGGCTGGGCGTCTGCGGCGAGCAGCAGGCACAGGTCgggg                                                                          <  1:5506047/150‑1 (MQ=255)
                                                                             ggAGGTCAGGCCATTGGCCCGGGGCTACCGGCAGCAGGCTGACGGCCAGCCCCAGCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCTGGTGATGGCGATCAGGGGGTGGGGCTGGGCGTCTGCGGCGAGCAGCAGGCACAGGTCGGGGc                                                                         >  1:2975929/1‑150 (MQ=255)
                                                                                ggTCAGGCCATTGGCCCGGGGCTACCGGCAGCAGGCTGACGGCCAGCCCCAGCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCTGGTGATGGCGATCAGGGGGTGGGGCTGGGCGTCTGCGGCGAGCAGCAGGCACAGGTCGGGGCCTa                                                                      <  1:2546555/150‑1 (MQ=255)
                                                                                  tCAGGCCATTGGCCCGGGGCTACCGGCAGCAGGCTGACGGCCAGCCCCAGCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCTGGTGATGGCGATCAGGGGGTGGGGCTGGGCGTCTGCGGCGAGCAGCAGGCACAGGTCGGGGCCTATg                                                                    <  1:5959466/150‑1 (MQ=255)
                                                                                  tCAGGCCATTGGCCCGGGGCTACCGGCAGCAGGCTGACGGCCAGCCCCAGCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCTGGTGATGGCGATCAGGGGGTGGGGCTGGGCGTCTGCGGCGAGCAGCAGGCACAGGTCGGGGCCTATg                                                                    >  1:2709263/1‑150 (MQ=255)
                                                                                   cAGGCCATTGGCCCGGGGCTACCGGCAGCAGGCTGACGGCCAGCCCCAGCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCTGGTGATGGCGATCAGGGGGTGGGGCTGGGCGTCTGCGGCGAGCAGCAGGCCCAGGTCGGGGCCTATgc                                                                   >  1:5047131/1‑150 (MQ=255)
                                                                                   cAGGCCATTGGCCCGGGGCTACCGGCAGCAGGCTGACGGCCAGCCCCAGCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCTGGTGATGGCGATCAGGGGGTGGGGCTGGGCGTCTGCGGCGAGCAGCAGGCACAGGTCGGGGCCTATgc                                                                   <  1:3839898/150‑1 (MQ=255)
                                                                                   cAGGCCATTGGCCCGGGGCTACCGGCAGCAGGCTGACGGCCAGCCCCAGCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCTGGTGATGGCGATCAGGGGGTGGGGCTGGGCGTCTGCGGCGAGCAGCAGGCACAGGTCGGGGCCTATgc                                                                   >  2:1359162/1‑150 (MQ=255)
                                                                                      gCCATTGGCCCGGGGCTACCGGCAGCAGGCTGACGGCCAGCCCCAGCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCTGGTGATGGCGATCAGGGGGTGGGGCTGGGCGTCTGCGGCGAGCAGCAGGCACAGGTCGGGGCCTATgccgg                                                                >  2:867768/1‑150 (MQ=255)
                                                                                       ccATTGGCCCGGGGCTACCGGCAGCAGGCTGACGGCCAGCCCCAGCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCTGGTGATGGCGATCAGGGGGTGGGGCTGGGCGTCTGCGGCGAGCAGCAGGCACAGGTCGGGGCCTATgccggc                                                               >  1:2394990/1‑150 (MQ=255)
                                                                                        cATTGACCCGGGGCTACCGGCAGCAGGCTGACGGCCAGCCCCAGCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCTGGTGATGGCGATCAGGGGGTGGGGCTGGGCGTCTGCGGCGAGCAGCAGGCACAGGTCGGGGCCTATgccggcc                                                              >  1:1661614/1‑150 (MQ=255)
                                                                                              cccGGGGCTACCGGCAGCAGGCTGACGGCCAGCCCCAGCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCTGGTGATGGCGATCAGGGGGTGGGGCTGGGCGTCTGCGGCGAGCAGCAGGCACAGGTCGGGGCCTATGCCGGCCGGCTCg                                                        >  2:1569741/1‑150 (MQ=255)
                                                                                              cccGGGGCTACCGGCAGCAGGCTGACGGCCAGCCCCAGCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCTGGTGATGGCGATCAGGGGGTGGGGCTGGGCGTCTGCGGCGAGCAGCAGGCACAGGTCGGGGCCTATGCCGGCCGGCTCg                                                        >  2:5705994/1‑150 (MQ=255)
                                                                                               ccGGGGCTACCGGCAGCAGGCTGACGGCCAGCCCCAGCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCTGGTGATGGCGATCAGGGGGTGGGGCTGGGCGTCTGCGGCGAGCAGCAGGCACAGGTCGGGGCCTATGCCGGCCGGCTCGc                                                       >  2:5904703/1‑150 (MQ=255)
                                                                                               ccGGGGCTACCGGCAGCAGGCTGACGGCCAGCCCCAGCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCTGGTGATGGCGATCAGGGGGTGGGGCTGGGCGTCTGCGGCGAGCAGCAGGCACAGGTCGGGGCCTATGCCGGCCGGCTCGc                                                       >  1:1002315/1‑150 (MQ=255)
                                                                                               ccGGGGCTACCGGCAGCAGGCTGACGGCCAGCCCCAGCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCTGGTGATGGCGATCAGGGGGTGGGGCTGGGCGTCTGCGGCGAGCAGCAGGCACAGGTCGGGGCCTATGCCGGCCGGCTCGc                                                       >  1:3643265/1‑150 (MQ=255)
                                                                                                    gCTACCGGCAGCAGGCTGACGGCCAGCCCCAGCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCTGGTGCTGGCGATCAGGGGGTGGGGCTGGGCGTCTGCGGCGAGCAGCAGGCACAGGTCGGGGCCTATGCCGGCCGGCTCGCCgggg                                                  <  2:2542494/150‑1 (MQ=255)
                                                                                                    gCTACCGGCAGCAGGCTGACGGCCAGCCCCAGCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCTGGTGATGGCGATCAGGGGGTGGGGCTGGGCGTCTGCGGCGAGCAGCAGGCACAGGTCGGGGCCTATGCCGGCCGGCTCGCCgggg                                                  >  1:5053300/1‑150 (MQ=255)
                                                                                                           gcagcaGGCTGACGGCCAGCCCCAGCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCTGGTGATGGCGATCAGGGGGTGGGGCTGGGCGTCTGCGGCGAGCAGCAGGCACAGGTCGGGGCCTATGCCGGCCGGCTCGCCGGGGGTGACgg                                           >  1:184901/1‑150 (MQ=255)
                                                                                                           gcagcaGGCTGACGGCCAGCCCCAGCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCGGGTGATGGCGATCCGGGGGTGGGGCTGGGCGTCTGCGGCGAGCAGCAGGCACAGGTCGGGGCCTATGCCGGCCGGCTCGCCGGGGGTGACgg                                           >  2:6612542/1‑150 (MQ=255)
                                                                                                             agcagGCTGACGGCCAGCCCCAGCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCTGGTGATGGCGATCAGGGGGTGGGGCTGGGCGTCTGCGGCGAGCAGCAGGCACAGGTCGGGGCCTATGCCGGCCGGCTCGCCGGGGGTGACGGCg                                         >  1:1627446/1‑150 (MQ=255)
                                                                                                                  gCTGACGGCCAGCCCCAGCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCTGGTGATGGCGATCAGGGGGTGGGGCTGGGCGTCTGCGGCGAGCAGCAGGCACAGGTCGGGGCCTATGCCGGCCGGCTCGCCGGGGGTGACGGCGAAgcg                                    >  1:3753806/1‑150 (MQ=255)
                                                                                                                  gCTGACGGCCAGCCCCAGCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCTGGTGATGGCGATCAGGGGGTGGGGCGGGGCGTCTGCGGCGAGCAGCAGGCACAGGTCGGGGCCTCTGCCGGCCGGCTCGCCGGGGGTGACGGCGAAgcg                                    >  2:356872/1‑150 (MQ=255)
                                                                                                                    tGACGGCCAGCCCCAGCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCTGGTGATGGCGATCAGGGGGTGGGGCTGGGCGTCTGCGGCGAGCAGCAGGCACAGGTCGGGGCCTATGCCGGCCGGCTCGCCGGGGGTGACGGCGAAGCGCa                                  >  1:2202771/1‑150 (MQ=255)
                                                                                                                    tGACGGCCAGCCCCAGCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCTGGTGATGGCGATCAGGGGGTGGGGCTGGGCGTCTGCGGCGAGCAGCAGGCACAGGTCGGGGCCTATGCCGGCCGGCTCGCCGGGGGTGACGGCGAAGCGCa                                  >  2:6429385/1‑150 (MQ=255)
                                                                                                                    tGACGGCCAGCCCCAGCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCTGGTGATGGCGATCAGGGGGTGGGGCTGGGCGTCTGCGGCGAGCAGCAGGCACAGGTCGGGGCCTATGCCGGCCGGCTCGCCGGGGGGGACGGCGAAgcgcc                                  >  1:6448634/1‑149 (MQ=255)
                                                                                                                      aCGGCCAGCCCCAGCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCTGGTGATGGCGATCAGGGGGTGGGGCTGGGCGTCTGCGGCGAGCAGCAGGCACAGGTCGGGGCCTATGCCGGCCGGCTCGCCGGGGGTGACGGCGAAGCGCAgg                                >  1:5597169/1‑150 (MQ=255)
                                                                                                                       cGGCCAGCCCCAGCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCTGGTGATGGCGATCAGGGGGTGGGGCTGGGCGTCTGCGGCGAGCAGCAGGCACAGGTCGGGGCCTATGCCGGCCGGCTCGCCGGGGGTGACGGCGAAGCGCAggg                               <  1:3514480/150‑1 (MQ=255)
                                                                                                                            aGCCCCAGCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCTGGTGATGGCGATCAGGGGGTGGGGATGGGCGTCTGCTGCGAGCAGCAGGCACAGGTCGGTGCCTATGCCGGCCGGCTCGCCGGGGGTGACGGCGAAGCGCAGGGGCTTc                          <  1:3939845/150‑1 (MQ=255)
                                                                                                                                ccAGCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCTGGTGATGGCGATCAGGGGGTGGGGCTGGGCGTCTGCGGCGAGCAGCAGGCACAGGTCGGGGCCTATGCCGGCCGGCTCGCCGGGGGTGACGGCGAAGCGCAGGGGCTTCACTg                      >  1:5873758/1‑150 (MQ=255)
                                                                                                                                   gCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCTGGTGATGGCGATCAGGGGGTGGGGCTGGGCGTCTGCGGCGAGCAGCAGGCACAGGTCGGGGCCTATGCCGGCCGGCTCGCCGGGGGTGACGGCGAAGCGCAGGGGCTTCACTGggc                   <  2:4886143/150‑1 (MQ=255)
                                                                                                                                   gCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCTGGTGATGGCGATCAGGGGGTGGGGCTGGGCGTCTGCGGCGAGCAGCAGGCACAGGTCGGGGCCTATGCCGGCCGGCTCGCCGGGGGTGACGGCGAAGCGCAGGGGCTTCACTGggc                   <  1:6338072/150‑1 (MQ=255)
                                                                                                                                      gCGTGGCCCGCTCGGCGAGCAGGTCACGGCTGGTGATGGCGATCAGGGGGTGGGGCTGGGCGTCTGCGGCGAGCAGCAGGCACAGGTCGGGGCCTATGCCGGCCGGCTCGCCGGGGGTGACGGCGAAGCGCAGGGGCTTCACTGggcggc                >  1:6751097/1‑150 (MQ=255)
                                                                                                                                        gTGGCCCGCTCGGCGAGCAGGTCACGGCTGGTGATGGCGATCAGGGGGTGGGGCTGGGCGTCTGCGGCGAGCAGCAGGCACAGGTCGGGGCCTATGCCGGCCGGCTCGCCGGGGGTGACGGCGAAGCGCAGGGGCTTCACTGGGCGGCCt              <  2:1401230/150‑1 (MQ=255)
                                                                                                                                             ccGCTCGGCGAGCAGGTCACGGCTGGTGATGGCGATCAGGGGGTGGGGCTGGGCGTCTGCGGCGAGCAGCAGGCACAGGTCGGGGCCTATGCCGGCCGGCTCGCCGGGGGTGACGGCGAAGCGCAGGGGCTTCACTGGGCGGCCTGGTCa         >  1:6595470/1‑150 (MQ=255)
                                                                                                                                             ccGCTCGGCGAGCAGGTCACGGCTGGTGATGGCGATCAGGGGGTGGGGCTGGGCGTCTGCGGCGAGCAGCAGGCACAGGTCGGGGCCTATGCCGGCCGGCTCGCCGGGGGTGACGGCGAAGCGCAGGGGCTTCACTGGGCGGCCTGGTCa         >  1:6222789/1‑150 (MQ=255)
                                                                                                                                             ccGCTCGGCGAGCAGGTCACGGCTGGTGATGGCGATCAGGGGGTGGGGCTGGGCGTCTGCGGCGAGCAGCAGGCACAGGTCGGGGCCTATGCCGGCCGGCTCGCCGGGGGTGACGGCGAAGCGCAGGGGCTTCACTGGGCGGCCTGGTCa         >  1:5669495/1‑150 (MQ=255)
                                                                                                                                                cTCGGCGAGCAGGTCACGGCTGGTGATGGCGATCAGGGGGTGGGGCTGGGCGTCTGCGGCGAGCAGCAGGCACAGGTCGGGGCCTATGCCGGCCGGCTCGCCGGGGGTGACGGCGAAGCGCAGGGGCTTCACTGGGCGGCCTGGt           <  2:5100984/145‑1 (MQ=255)
                                                                                                                                                cTCGGCGAGCAGGTCACGGCTGGTGATGGCGATCAGGGGGTGGGGCTGGGCGTCTGCGGCGAGCAGCAGGCACAGGTCGGGGCCTATGCCGGCCGGCTCGCCGGGGGTGACGGCGAAGCGCAGGGGCTTCACTGGGCGGCCTGGTCAgcg      >  1:1724752/1‑150 (MQ=255)
                                                                                                                                                  cGGCGAGCAGGTCACGGCTGGTGATGGCGATCAGGGGGTGGGGCTGGGCGTCTGCGGCGAGCAGCAGGCACAGGTCGGGGCCTATGCCGGCCGGTTCGCCGGGGGTGACGGCGAAGCGCAGGGGCTTCACTGGGCGGCCTGGTCAGCGcc    <  2:4428561/150‑1 (MQ=255)
                                                                                                                                                    gCGAGCAGGTCACGGCTGGTGATGGCGATCAGGGGGTGGGGCTGGGCGTCTGCGGCGAGCAGCAGGCACAGGTCGGGGCCTATGCCGGCCGGCTCGCCGGGGGTGACGGCGAAGCGCAGGGGCTTCACTGGGCGGCCTGGTCAGCGCCAg  <  1:2095599/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GGCCTTGTCGAGCTGGCCTGGCACCACCGGGGCCGCCAAGGGGGTATCCCAGACGTACAGGCTACCTGCGGGTGCCGGGAGGTCAGGCCATTGGCCCGGGGCTACCGGCAGCAGGCTGACGGCCAGCCCCAGCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCTGGTGATGGCGATCAGGGGGTGGGGCTGGGCGTCTGCGGCGAGCAGCAGGCACAGGTCGGGGCCTATGCCGGCCGGCTCGCCGGGGGTGACGGCGAAGCGCAGGGGCTTCACTGGGCGGCCTGGTCAGCGCCAG  >  NC_002947/488029‑488326

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: