Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 488,171 C→T 5.1% R45Q (CGG→CAG)  pdxA ← 4‑hydroxythreonine‑4‑phosphate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947488,1710CT 5.1% 1005.3 / 25.4 311R45Q (CGG→CAG) pdxA4‑hydroxythreonine‑4‑phosphate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (159/136);  new base T (7/9);  total (166/145)
Fisher's exact test for biased strand distribution p-value = 4.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 318 total aligned reads displayed.

CTGCGTTGGCCTTGTCGAGCTGGCCTGGCACCACCGGGGCCGCCAAGGGGGTATCCCAGACGTACAGGCTACCTGCGGGTGCCGGGAGGTCAGGCCATTGGCCCGGGGCTACCGGCAGCAGGCTGACGGCCAGCCCCAGCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCTGGTGATGGCGATCAGGGGGTGGGGCTGGGCGTCTGCGGCGAGCAGCAGGCACAGGTCGGGGCCTATGCCGGCCGGCTCGCCGGGGGTGACGGCGAAGCGCAGGGGCTTCACTGGGCGGCCTGGTCA  >  NC_002947/488022‑488319
                                                                                                                                                     |                                                                                                                                                    
    gTTGGCCTTGTCGAGCTGGCCTGGCACCACCGGGGCCGCCAAGGGGGTATCCCAGACGTACAGGCTACCTGCGGGTGCCGGGAGGTCAGGCCATTGGCCCGGGGCTACCGGCAGCAGGCTGACGGCCAGCCCCAGCTGCGTGGCCCGCTc                                                                                                                                                  >  2:5211153/1‑150 (MQ=255)
     ttGGCCTTGTCGAGCTGGCCTGGCACCACCGGGGCCGCCAAGGGGGTATCCCAGACGTACAGGCTACCTGCGGGTGCCGGGAGGTCAGGCCATTGGCCCGGGGCTACCGGCAGCAGGCTGACGGCCAGCCCCAGCTGCGTGGCCCGCTCg                                                                                                                                                 >  1:1687933/1‑150 (MQ=255)
       ggCCTTGTCGAGCTGGCCTGGCACCACCGGGGCCGCCAAGGGGGTATCCCAGACGTACAGGCTACCTGCGGGTGCCGGGAGGTCAGGCCATTGGCCCGGGGCTACCGGCAGCAGGCTGACGGCCAGCCCCAGCTGCGTGGCCCGCTCGGc                                                                                                                                               <  2:5105361/150‑1 (MQ=255)
          cTTGTCGAGCTGGCCTGGCACCACCGGGGCCGCCAAGGGGGTATCCCAGTCGTACAGGCTACCTGCGGGTGCCGGGAGGTCAGGCCATTGGCCCGGGGCTACCGGCAGCAGGCTGACGGCCAGCCCCAGCTGCGTGGCCCGCTCGGCGat                                                                                                                                            <  2:55049/150‑2 (MQ=255)
            tGTCGAGCTGGCCTGGCACCACCGGGGCCGCCAAGGGGGTATCCCAGACGTACAGGCTACCTGCGGGTGCCGGGAGGTCAGGCCATTGGCCCGGGGCTACCGGCAGCAGGCTGACGGCCAGCCCCAGCTGCGTGGCCCGCTCGGCGAGCa                                                                                                                                          <  1:2924406/150‑1 (MQ=255)
              tCGAGCTGGCCTGGCACCACCGGGGCCGCCAAGGGGGTATCCCAGACGTACAGGCTACCTGCGGGTGCCGGGAGGTCAGGCCATTGGCCCGGGGCTACCGGCAGCAGGCTGACGGCCAGCCCCAGCTGCGTGGCCCGCTCGGCGAGCAgg                                                                                                                                        <  2:783191/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)
                      gcctggcACCACCGGGGCCGCCAAGGGGGTATCCCAGACGTACAGGCTACCTGCGGGTGCCGGGAGGTCAGGCCATTGGCCCGGGGCTACCGGCAGCAGGCTGACGGCCAGCCCCAGCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCt                                                                                                                                >  2:5334785/1‑150 (MQ=255)
                         tggcACCACCGGGGCCGCCAAGGGGGTATCCCAGACGTACAGGCTACCTGCGGGTGCCGGGAGGTCAGGCCATTGGCCCGGGGCTACCGGCAGCAGGCTGACGGCCAGCCCCAGCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCTGGt                                                                                                                             >  1:5876105/1‑150 (MQ=255)
                             accaccGGGGCCGCCAAGGGGGTATCCCAGACGTACAGGCTACCTGCGGGTGCCGGGAGGTCAGGCCATTGGCCCGGGGCTACCGGCAGCAGGCTGACGGCCAGCCCCAGCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCTGGTGATg                                                                                                                         >  1:6624162/1‑150 (MQ=255)
                              ccaccGGGGCCGCCAAGGGGGTATCCCAGACGTACAGGCTACCTGCGGGTGCCGGGAGGTCAGGCCATTGGCCCGGGGCTACCGGCAGCAGGCTGACGGCCAGCCCCAGCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCTGGTGATgg                                                                                                                        <  2:1972750/150‑1 (MQ=255)
                              ccaccGGGGCCGCCAAGGGGGTATCCCAGACGTACAGGCTACCTGCGGGTGCCGGGAGGTCAGGCCATTGGCCCGGGGCTACCGGCAGCAGGCTGACGGCCAGCCCCAGCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCTGGTGATgg                                                                                                                        <  1:1187166/150‑1 (MQ=255)
                               caccGGGGCCGCCAAGGGGGTATCCCAGACGTACAGGCTACCTGCGGGTGCCGGGAGGTCAGGCCATTGGCCCGGGGCTACCGGCAGCAGGCTGACGGCCAGCCCCAGCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCTGGTGATGGc                                                                                                                       >  2:3537200/1‑150 (MQ=255)
                                accGGGGCCGCCAAGGGGGTATCCCAGACGTACAGGCTACCTGCGGGTGCCGGGAGGTCAGGCCATTGGCCCGGGGCTACCGGCAGCAGGCTGACGGCCAGCCCCAGCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCTGGTGATGGCg                                                                                                                      >  2:6338072/1‑150 (MQ=255)
                                 ccGGGGCCGCCAAGGGGGTATCCCAGACGTACAGGCTACCTGCGGGTGCCGGGAGGTCAGGCCATTGGCCCGGGGCTACCGGCAGCAGGCTGACGGCCAGCCCCAGCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCTGGTGATGGCGa                                                                                                                     >  1:717169/1‑150 (MQ=255)
                                  cGGGGCCGCCAAGGGGGTATCCCAGACGTACAGGCTACCTGCGGGTGCCGGGAGGTCAGGCCATTGGCCCGGGGCTACCGGCAGCAGGCTGACGGCCAGCCCCAGCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCTGGTGATGGCGAt                                                                                                                    >  2:408475/1‑150 (MQ=255)
                                    gggCCGCCAAGGGGGTATCCCAGACGTACAGGCTACCTGCGGGTGCCGGGAGGTCAGGCCATTGGCCCGGGGCTACCGGCAGCAGGCTGACGGCCAGCCCCAGCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCTGGTGATGGCGATCa                                                                                                                  >  2:621268/1‑150 (MQ=255)
                                      gccgccAAGGGGGTATCCCAGACGTACAGGCTACCTGCGGGTGCCGGGAGGTCAGGCCATTGGCCCGGGGCTACCGGCAGCAGGCTGACGGCCAGCCCCAGCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCTGGTGATGGCGATCAgg                                                                                                                <  1:6256950/150‑1 (MQ=255)
                                          ccAAGGGGGTATCCCAGACGTACAGGCTACCTGCGGGTGCCGGGAGGTCAGGCCATTGGCCCGGGGCTACCGGCAGCAGGCTGACGGCCAGCCCCAGCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCTGGTGATGGCGATCAGGGGGt                                                                                                            <  1:199877/150‑1 (MQ=255)
                                           cAAGGGGGTATCCCAGACGTACAGGCTACCTGCGGGTGCCGGGAGGTCAGGCCATTGGCCCGGGGCTACCGGCAGCAGGCTGACGGCCAGCCCCAGCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCTGGTGATGGCGATCAGGGGGTg                                                                                                           >  1:6227972/1‑150 (MQ=255)
                                            aaGGGGGTATCCCAGACGTACAGGCTACCTGCGGGTGCCGGGAGGTCAGGCCATTGGCCCGGGGCTACCGGCAGCAGGCTGACGGCCAGCCCCAGCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCTGGTGATGGCGATCAGGGGGTgg                                                                                                          <  2:3277938/150‑1 (MQ=255)
                                            aaGGGGGTATCCCAGACGTACAGGCTACCTGCGGGTGCCGGGAGGTCAGGCCATTGGCCCGGGGCTACCGGCAGCAGGCTGACGGCCAGCCCCAGCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCTGGTGATGGCGATCAGGGGGTgg                                                                                                          <  1:5750547/150‑1 (MQ=255)
                                            aaGGGGGTATCCCAGACGTACAGGCTACCTGCGGGTGCCGGGAGGTCAGGCCATTGGCCCGGGGCTACCGGCAGCAGGCTGACGGCCAGCCCCAGCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCTGGTGATGGCGATCAGGGGGTgg                                                                                                          <  1:4877297/150‑1 (MQ=255)
                                               ggggTATCCCAGACGTACAGGCTACCTGCGGGTGCCGGGAGGTCAGGCCATTGGCCCGGGGCTACCGGCAGCAGGCTGACGGCCAGCCCCAGCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCTGGTGATGGCGATCAGGGGGTGGGGc                                                                                                       >  2:1843571/1‑150 (MQ=255)
                                               ggggTATCCCAGACGTACAGGCTACCTGCGGGTGCCGGGAGGTCAGGCCATTGGCCCGGGGCTACCGGCAGCAGGCTGACGGCCAGCCCCAGCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCTGGTGATGGCGATCAGGGGGTGGGGc                                                                                                       >  2:5779631/1‑150 (MQ=255)
                                                 ggTATCCCAGACGTACAGGCTACCTGCGGGTGCCGGGAGGTCAGGCCATTGGCCCGGGGCTACCGGCAGCAGGCTGACGGCCAGCCCCAGCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCTGGTGATGGCGATCAGGGGGTGGGGCTg                                                                                                     >  2:2623806/1‑150 (MQ=255)
                                                  gTATCCCAGACGTACAGGCTACCTGCGGGTGCCGGGAGGTCAGGCCATTGGCCCGGGGCTACCGGCAGCAGGCTGACGGCCAGCCCCAGCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCTGGTGATGGCGATCAGGGGGTGGGGCTgg                                                                                                    <  1:1126341/150‑1 (MQ=255)
                                                    aTCCCAGACGTACAGGCTACCTGCGGGTGCCGGGAGGTCAGGCCATTGGCCCGGGGCTACCGGCAGCAGGCTGACGGCCAGCCCCAGCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCTGGTGATGGCGATCAGGGGGTGGGGCTGGGc                                                                                                  <  1:135597/150‑1 (MQ=255)
                                                    aTCCCAGACGTACAGGCTACCTGCGGGTGCCGGGAGGTCAGGCCATTGGCCCGGGGCTACCGGCAGCAGGCTGACGGCCAGCCCCAGCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCTGGTGATGGCGATCAGGGGGTGGGGCTGGGc                                                                                                  <  1:3525199/150‑1 (MQ=255)
                                                     tCCCAGACGTACAGGCTACCTGCGGGTGCCGGGAGGTCAGGCCATTGGCCCGGGGCTACCGGCAGCAGGCTGACGGCCAGCCCCAGCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCTGGTGATGGCGATCAGGGGGTGGGGCTGGGCg                                                                                                 >  1:4280942/1‑150 (MQ=255)
                                                       ccAGACGTACAGGCTACCTGCGGGTGCCGGGAGGTCAGGCCATTGGCCCGGGGCTACCGGCAGCAGGCTGACGGCCAGCCCCAGCTGCGTGGCCTGCTCGCCGAGcgg                                                                                                                                         <  2:4001216/108‑3 (MQ=255)
                                                       ccAGACGTACAGGCTACCTGCGGGTGCCGGGAGGTCAGGCCATTGGCCCGGGGCTACCGGCAGCAGGCTGACGGCCAGCCCCAGCTGCGTGGCCTGCTCGCCGAGcgg                                                                                                                                         >  1:4001480/1‑106 (MQ=255)
                                                       ccAGACGTACAGGCTACCTGCGGGTGCCGGGAGGTCAGGCCATTGGCCCGGGGCTACCGGCAGCAGGCTGACGGCCAGCCCCAGCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCTGGTg                                                                                                                            <  2:2034362/121‑1 (MQ=255)
                                                       ccAGACGTACAGGCTACCTGCGGGTGCCGGGAGGTCAGGCCATTGGCCCGGGGCTACCGGCAGCAGGCTGACGGCCAGCCCCAGCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCTGGTGATGGCGATCAGGGGGTGGGGCTGGGCGTc                                                                                               >  2:5233943/1‑150 (MQ=255)
                                                       ccAGACGTACAGGCTACCTGCGGGTGCCGGGAGGTCAGGCCATTGGCCCGGGGCTACCGGCAGCAGGCTGACGGCCAGCCCCAGCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCTGGTGATGGCGATCAGGGGGTGGGGCTGGGCGTc                                                                                               >  2:3205193/1‑150 (MQ=255)
                                                         aGACGTACAGGCTACCTGCGGGTGCCGGGAGGTCAGGCCATTGGCCCGGGGCTACCGGCAGCAGGCTGACGGCCAGCCCCAGCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCTGGTGATGGCGATCAGGGGGGGGGGCTGGGCGTCTg                                                                                             >  2:3282501/1‑150 (MQ=255)
                                                           aCGTACAGGCTACCTGCGGGTGCCGGGAGGTCAGGCCATTGGCCCGGGGCTACCGGCAGCAGGCTGACGGCCAGCCCCAGCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCTGGTGATGGCGATCAGGGGGTGGGGCTGGGCGTCTgcg                                                                                           <  2:1047574/150‑1 (MQ=255)
                                                             gTACAGGCTACCTGCGGGTGCCGGGAGGTCAGGCCATTGGCCCGGGGCTACCGGCAGCAGGCTGACGGCCAGCCCCAGCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCTGGTGATGGCGATCAGGGGGTGGGGCTGGGCGTCTgcggc                                                                                         <  1:5666847/150‑1 (MQ=255)
                                                               aCAGGCTACCTGCGGGTGCCGGGAGGTCAGGCCATTGGCCCGGGGCTACCGGCAGCAGGCTGACGGCCAGCCCCAGCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCTGGTGATGGCGATCAGGGGGTGGGGCTGGGCGTCTGCGGCGa                                                                                       <  1:2479845/150‑1 (MQ=255)
                                                               aCAGGCTACCTGCGGGTGCCGGGAGGTCAGGCCATTGGCCCGGGGCTACCGGCAGCAGGCTGACGGCCAGCCCCAGCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCTGGTGATGGCGATCAGGGGGTGGGGCTGGGCGTCTGCGGCGa                                                                                       <  1:3702406/150‑1 (MQ=255)
                                                               aCAGGCTACCTGCGGGTGCCGGGAGGTCAGGCCATTGGCCCGGGGCTACCGGCAGCAGGCTGACGGCCAGCCCCAGCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCTGGTGATGGCGATCAGGGGGTGGGGCTGGGCGTCTGCGGCGa                                                                                       <  1:6310083/150‑1 (MQ=255)
                                                                 aGGCTACCTGCGGGTGCCGGGAGGTCAGGCCATTGGCCCGGGGCTACCGGCAGCAGGCTGACGGCCAGCCCCAGCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCTGGTGATGGCGATCAGGGGGTGGGGCTGGGCGTCTGCGGCGagc                                                                                     >  1:6027140/1‑150 (MQ=255)
                                                                  ggCTACCTGCGGGTGCCGGGAGGTCAGGCCATTGGCCCGGGGCTACCGGCAGCAGGCTGACGGCCAGCCCCAGCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCTGGTGATGGCGATCAGGGGGTGGGGCTGGGCGTCTGCGGCGagca                                                                                    >  1:880113/1‑150 (MQ=255)
                                                                  ggCTACCTGCGGGTGCCGGGAGGTCAGGCCATTGGCCCGGGGCTACCGGCAGCAGGCTGACGGCCAGCCCCAGCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCTGGTGATGGCGATCAGGGGGTGGGGCTGGGCGTCTGCGGCGagca                                                                                    <  2:4587736/150‑1 (MQ=255)
                                                                   gCTACCTGCGGGTGCCGGGAGGTCAGGCCATTGGCCCGGGGCTACCGGCAGCAGGCTGACGGCCAGCCCCAGCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCTGGTGATGGCGAt                                                                                                                    <  1:5153630/117‑1 (MQ=255)
                                                                        cTGCGGGTGCCGGGAGGTCAGGCCATTGGCCCGGGGCTACCGGCAGCAGGCTGACGGCCAGCCCCAGCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCTGGTGATGGCGATCAGGGGGTGGGGCTGGGCGTCTGCGGCGAGCAGCAGGc                                                                              >  1:4420290/1‑150 (MQ=255)
                                                                        cTGCGGGTGCCGGGAGGTCAGGCCATTGGCCCGGGGCTACCGGCAGCAGGCTGACGGCCAGCCCCAGCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCTGGTGATGGCGATCAGGGGGTGGGGCTGGGCGTCTGCGGCGAGCAGCAGGc                                                                              >  2:2340726/1‑150 (MQ=255)
                                                                            gggTGCCGGGAGGTCAGGCCATTGGCCCGGGGCTACCGGCAGCAGGCTGACGGCCAGCCCCAGCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCTGGTGATGGCGATCAGGGGGTGGGGCTGGGCGTCTGCGGCGAGCAGCAGGCACAg                                                                          <  2:3425989/150‑1 (MQ=255)
                                                                            gggTGCCGGGAGGTCAGGCCATTGGCCCGGGGCTACCGGCAGCAGGCTGACGGCCAGCCCCAGCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCTGGTGATGGCGATCAGGGGGTGGGGCTGGGCGTCTGCGGCGAGCAGCAGGCACAg                                                                          >  1:4961356/1‑150 (MQ=255)
                                                                              gTGCCGGGAGGTCAGGCCATTGGCCCGGGGCTACCGGCAGCAGGCTGACGGCCAGCCCCAGCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCTGGTGATGGCGATCAGGGGGTGGGGCTGGGCGTCTGCGGCGAGCAGCAGGCACAGGt                                                                        <  1:5917901/150‑1 (MQ=255)
                                                                                 ccGGGAGGTCAGGCCATTGGCCCGGGGCTACCGGCAGCAGGCTGACGGCCAGCCCCAGCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCTGGTGATGGCGATCAGGGGGTGGGGCTGGGCGTCTGCGGCGAGCAGCAGGCACAGGTCgg                                                                     <  1:1497648/150‑1 (MQ=255)
                                                                                  cGGGAGGTCAGGCCATTGGCCCGGGGCTACCGGCAGCAGGCTGACGGCCAGCCCCAGCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCTGGTGATGGCGATCAGGGGGTGGGGCTGGGCGTCTGCGGCGAGCAGCAGGCACAGGTCggg                                                                    >  1:3584775/1‑150 (MQ=255)
                                                                                  cGGGAGGTCAGGCCATTGGCCCGGGGCTACCGGCAGCAGGCTGACGGCCAGCCCCAGCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCTGGTGATGGCGATCAGGGGGTGGGGCTGGGCGTCTGCGGCGAGCAGCAGGCACAGGTCggg                                                                    >  2:2415712/1‑150 (MQ=255)
                                                                                   gggAGGTCAGGCCATTGGCCCGGGGCTACCGGCAGCAGGCTGACGGCCAGCCCCAGCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCTGGTGATGGCGATCAGGGGGTgg                                                                                                          >  1:3280455/1‑111 (MQ=255)
                                                                                   gggAGGTCAGGCCATTGGCCCGGGGCTACCGGCAGCAGGCTGACGGCCAGCCCCAGCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCTGGTGATGGCGATCAGGGGGTGGGGCTGGGCGTCTGCGGCGAGCAGCAGGCACAGGTCgggg                                                                   <  1:5505799/150‑1 (MQ=255)
                                                                                    ggAGGTCAGGCCATTGGCCCGGGGCTACCGGCAGCAGGCTGACGGCCAGCCCCAGCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCTGGTGATGGCGATCAGGGGGTGGGGCTGGGCGTCTGCGGCGAGCAGCAGGCACAGGTCGGGGc                                                                  <  2:4388214/150‑1 (MQ=255)
                                                                                     gAGGTCAGGCCATTGGCCCGGGGCTACCGGCAGCAGGCTGACGGCCAGCCCCAGCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCTGGTGATGGCGATCAGGGGGTGGGGCTGGGCGTCTGCGGCGAGCAGCAGGCACAGGTCGGGGcc                                                                 >  2:306373/1‑150 (MQ=255)
                                                                                      aGGTCAGGCCATTGGCCCGGGGCTACCGGCAGCAGGCTGACGGCCAGCCCCAGCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCTGGTGATGGCGATCAGGGGGTGGGGCTGGGCGTCTGCGGCGAGCAGCAGGCACAGGTCGGGGCCt                                                                <  2:5627261/150‑1 (MQ=255)
                                                                                       ggTCAGGCCATTGGCCCGGGGCTACCGGCAGCAGGCTGACGGCCAGCCCCAGCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCTGGTGATGGCGATCAGGGGGTGGGGCTGGGCGTCTGCGGCGAGCAGCAGGCACAGGTCGGGGCCTa                                                               <  1:6114258/150‑1 (MQ=255)
                                                                                          cAGGCCATTGGCCCGGGGCTACCGGCAGCAGGCTGACGGCCAGCCCCAGCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCTGGTGATGGCGATCAGGGGGTGGGGCTGGGCGTCTGCGGCGAGCAGCAGGCACAGGTCGGGGCCTATgc                                                            >  2:1359162/1‑150 (MQ=255)
                                                                                           aGGCCATTGGCCCGGGGCTACCGGCAGCAGGCTGACGGCCAGCCCCAGCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCTGGTGATGGCGATCAgg                                                                                                                >  1:2554819/1‑97 (MQ=255)
                                                                                              ccATTGGCCCGGGGCTACCGGCAGCAGGCTGACGGCCAGCCCCAGCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCTGGTGATGGCGATCAGGGGGTGGGGCTGGGCGTCTGCGGCGAGCAGCAGGCACAGGTCGGGGCCTATgccggc                                                        >  1:504944/1‑150 (MQ=255)
                                                                                               cATTGGCCCGGGGCTACCGGCAGCAGGCTGACGGCCAGCCCCAGCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCTGGTGATGGCGATCAGGGGGTGGGGCTGGGCGTCTGCGGCGAGCAGCAGGCACAGGTCGGGGCCTATgccggcc                                                       >  1:2685156/1‑150 (MQ=255)
                                                                                                 ttGGCCCGGGGCTACCGGCAGCAGGCTGACGGCCAGCCCCAGCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCTGGTGATGGCGATCAGGGGGTGGGGCTGGGCGTCTGCGGCGAGCAGCAGGCACAGGTCGGGGCCTATGccggccgg                                                     <  1:312270/150‑1 (MQ=255)
                                                                                                    gCCCGGGGCTACCGGCAGCAGGCTGACGGCCAGCCCCAGCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCGGGTGATGGCGATCAGGGGGTGGGGCTGGGCGTCTGCGGCGAGCAGCAGGCACAGGTCGGGGCCTATGCCGGCCGGCTc                                                  >  1:2060376/1‑150 (MQ=255)
                                                                                                    gCCCGGGGCTACCGGCAGCAGGCTGACGGCCAGCCCCAGCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCGGGTGATGGCGATCAGGGGGTGGGGCTGGGCGTCTGCGGCGAGCAGCAGGCACAGGTCGGGGCCTAGGCCGGCCGGCTc                                                  >  2:1699655/1‑150 (MQ=255)
                                                                                                     cccGGGGCTACCGGCAGCAGGCTGACGGCCAGCCCCAGCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCTGGTGATGGCGATCAGGGGGTGGGGCTGGGCGTCTGCGGCGAGCAGCAGGCACAGGTCGGGGCCTATGCCGGCCGGCTCg                                                 >  2:5705994/1‑150 (MQ=255)
                                                                                                      ccGGGGCTACCGGCAGCAGGCTGACGGCCAGCCCCAGCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCTGGTGATGGCGATCAGGGGGTGGGGCTGGGCGTCTGCGGCGAGCAGCAGGCACAGGTCGGGGCCTATGCCGGCCGGCTCGc                                                >  1:1002315/1‑150 (MQ=255)
                                                                                                         gggCTACCGGCAGCAGGCTGACGGCCAGCCCCAGCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCTGGTGATGGCGATCAGGGGGTGGGGCTGGGCGTCTGCGGCGAGCAGCAGGCACAGGTCGGGGCCTATGCCGGCCGGCTCGCCgg                                             >  1:3504825/1‑150 (MQ=255)
                                                                                                           gCTACCGGCAGCAGGCTGACGGCCAGCCCCAGCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCTGGTGCTGGCGATCAGGGGGTGGGGCTGGGCGTCTGCGGCGAGCAGCAGGCACAGGTCGGGGCCTATGCCGGCCGGCTCGCCgggg                                           <  2:2542494/150‑1 (MQ=255)
                                                                                                           gCTACCGGCAGCAGGCTGACGGCCAGCCCCAGCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCTGGTGATGGCGATCAGGGGGTGGGGCTGGGCGTCTGCGGCGAGCAGCAGGCACAGGTCGGGGCCTATGCCGGCCGGCTCGCCgggg                                           <  2:4169388/150‑1 (MQ=255)
                                                                                                                cGGCAGCAGGCTGACGGCCAGCCCCAGCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCTGGTGATGGCGATAAGGGGGTGGGGCTGGGCGTCTGCGGCGAGCAGCAGGCACAGGTCGGGGCCTATGCCGGCCGGCTCGCCGGGGGTGAc                                      >  1:6371902/1‑150 (MQ=255)
                                                                                                                   cagcagGCTGACGGCCAGCCCCAGCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCTGGTGATGGCGATCAGGGGGTGGGGCTGGGCGTCTGCGGCGAGCAGCAGGCACAGGTCGGGGCCTATGCCGGCCGGCTCGCCGGGGGTGACGGc                                   <  1:5563871/150‑1 (MQ=255)
                                                                                                                    agcagGCTGACGGCCAGCCCCAGCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCTGGTGATGGCGATCAGGGGGTGGGGCTGGGCGTCTGCGGCGAGCAGCAGGCACAGGTCGGGGCCTATGCCGGCCGGCTCGCCGGGGGTGACGGCg                                  >  1:1627446/1‑150 (MQ=255)
                                                                                                                      cagGCTGACGGCCAGCCCCAGCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCTGGTGATGGCGATCAGGGGGTGGGGCTGGGCGTCTGCGGCGAGCAGCAGGCACAGGTCGGGGCCTATGCCGGCCGGCTCGCCGGGGGTGACGGCGaa                                <  1:3929175/150‑1 (MQ=255)
                                                                                                                      cagGCTGACGGCCAGCCCCAGCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCTGGTGATGGCGATCAGGGGGTGGGGCTGGGCGTCTGCGGCGAGCAGCAGGCACAGGTCGGGGCCTATGCCGGCCGGCTCGCCGGGGGTGACGGCGaa                                >  1:3032056/1‑150 (MQ=255)
                                                                                                                        ggCTGACGGCCAGCCCCAGCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCTGGTGATGGCGATCAGGGGGTGGGGCTGGGCGTCTGCGGCGAGCAGCAGGCACAGGTCGGGGCCTATGCCGGCCGGCTCGCCGGGGGTGACGGCGAAgc                              <  1:3828476/150‑1 (MQ=255)
                                                                                                                        ggCTGACGGCCAGCCCCAGCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCTGGTGATGGCGATCAGGGGGTGGGGCTGGGCGTCTGCGGCGAGCAGCAGGCACAGGTCGGGGCCTATGCCGGCCGGCTCGCCGGGGGTGACGGCGAAgc                              <  1:408475/150‑1 (MQ=255)
                                                                                                                         gCTGACGGCCAGCCCCAGCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCTGGTGATGGCGATCAGGGGGTGGGGCTGGGCGTCGGCGGCGAGCAGCAGGCACAGGTCGGGGCCTATGCCGGCCGGCTCGCCGGGGGTGACGGCGAAgcg                             >  2:149601/1‑150 (MQ=255)
                                                                                                                           tGACGGCCAGCCCCAGCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCTGGTGATGGCGATCAGGGGGTGGGGCTGGGCGTCTGCGGCGAGCAGCAGGCACAGGTCGGGGCCTATGCCGGCCGGCTCGCCGGGGGTGACGGCGAAGCGCa                           >  2:6429385/1‑150 (MQ=255)
                                                                                                                           tGACGGCCAGCCCCAGCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCTGGTGATGGCGATCAGGGGGTGGGGCTGGGCGTCTGCGGCGAGCAGCAGGCACAGGTCGGGGCCTATGCCGGCCGGCTCGCCGGGGGTGACGGCGAAGCGCa                           <  2:4812020/150‑1 (MQ=255)
                                                                                                                           tGACGGCCAGCCCCAGCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCTGGTGATGGCGATCAGGGGGTGGGGCTGGGCGTCTGCGGCGAGCAGCAGGCACAGGTCGGGGCCTATGCCGGCCGGCTCGCCGGGGGTGACGGCGAAGCGCa                           >  1:2202771/1‑150 (MQ=255)
                                                                                                                           tGACGGCCAGCCCCAGCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCTGGTGATGGCGATCAGGGGGTGGGGCTGGGCGTCTGCGGCGAGCAGCAGGCACAGGTCGGGGCCTATGCCGGCCGGCTCGCCGGGGGGGACGGCGAAgcgcc                           >  1:6448634/1‑149 (MQ=255)
                                                                                                                           tGACGGCCAGCCCCAGCTCCGTGGCCCGCTCGGCGAGCAGGTCACGGCTGGTGATGGCGATCAGGGGGTGGGGCTGGGCGTCTGCGGCGAGCAGCAGGCACAGGTCGGGGCCTATGCCGGCCGGCTCGCCGGGGGTGACGGCGAAGCGCa                           >  2:5925362/1‑150 (MQ=255)
                                                                                                                             aCGGCCAGCCCCAGCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCTGGTGAGGGCGATCAGGGGGTGGGGCTGGGCGTCTGCGGCGAGCAGCAGGCACAGGTCGGGGCCTATGCCGGCCGGCTCGCCGGGGGTGACGGCGAAGCGCAgg                         <  2:887388/150‑1 (MQ=255)
                                                                                                                                    gCCCCAGCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCTGGTGATGGCGATCAGGGGGTGGGGCTGGGCGTCTGCGGCGAGCAGCAGGCACAGGTCGGGGCCTATGCCGGCCGGCTCGCCGGGGGTGACGGCGAAGCGCAGGGGCTTCa                  <  1:4788345/150‑1 (MQ=255)
                                                                                                                                      cccAGCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCTGGTGATGGCGATCAGGGGGTGGGGCTGGGCGTCTGCGGCGAGCAGCAGGCACAGGTCGGGGCCTATGCCGGCCGGCTCGCCGGGGGTGACGGCGAAGCGCAGGGGCTTCACt                <  2:5800748/150‑1 (MQ=255)
                                                                                                                                        cAGCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCTGGTGATGGCGATCAGGGGGTGGGGCTGGGCGTCTGCGGCGAGCAGCAGGCACAGGTCGGGGCCTATGCCGGCCGGCTCGCCGGGGGTGACGGCGAAGCGCAGGGGCTTCACTgg              <  2:3807860/150‑1 (MQ=255)
                                                                                                                                        cAGCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCTGGTGATGGCGATCAGGGGGTGGGGCTGGGCGTCTGCGGCGAGCAGCAGGCACAGGTCGGGGCCTATGCCGGCCGGCTCGCCGGGGGTGACGGCGAAGCGCAGGGGCTTCACTgg              >  2:6867175/1‑150 (MQ=255)
                                                                                                                                         aGCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCTGGTGATGGCGATCAGGGGGTGGGGCTGGGCGTCTGCGGCGAGCAGCAGGCACAGGTCGGGGCCTATGCCGGCCGGCTCGCCGGGGGTGACGGCGAAGCGCAGGGGCTTCACTggg             <  1:5971713/150‑1 (MQ=255)
                                                                                                                                          gCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCTGGTGATGGCGATCAGGGGGTGGGGCTGGGCGTCTGCGGCGAGCAGCAGGCACAGGTCGGGGCCTATGCCGGCCGGCTCGCCGGGGGTGACGGCGAAGCGCAGGGGCTTCACTGggc            <  2:4886143/150‑1 (MQ=255)
                                                                                                                                          gCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCTGGTGATGGCGATCAGGGGGTGGGGCTGGGCGTCTGCGGCGAGCAGCAGGCACAGGTCGGGGCCTATGCCGGCCGGCTCGCCGGGGGTGACGGCGAAGCGCAGGGGCTTCACTGggc            <  2:1687933/150‑1 (MQ=255)
                                                                                                                                          cagccGTGACCTGCTCGCCGAGCGGGCCACGGCTGGTGATGGCGATCAGGGGGTGGGGCTGGGCGTCTGCGGCGAGCAGCAGGCACAGGTCGGGGCCTATGCCGGCCGGCTCGCCGGGGGTGACGGCGAAGCGCAGGGGCTTCACTGggc            >  1:27713/5‑150 (MQ=255)
                                                                                                                                           agccGTGACCTGCTCGCCGAGCGGGCCACGGCTGGTGATGGCGATCAGGGGGTGGGGCTGGGCGTCTGCGGCGAGCAGCAGGCACAGGTCGGGGCCTATGCCGGCCGGCTCGCCGGGGGTGACGGCGAAGCGCAGGGGCTTCACTGggcg           >  1:1731303/4‑150 (MQ=255)
                                                                                                                                               gTGGCCCGCTCGGCGAGCAGGTCACGGCTGGTGATGGCGATCAGGGGGTGGGGCTGGGCGTCTGCGGCGAGCAGCAGGCACAGGTCGGGGCCTATGCCGGCCGGCTCGCCGGGGGTGACGGCGAAGCGCAGGGGCTTCACTGGGCGGCCt       <  2:1401230/150‑1 (MQ=255)
                                                                                                                                               gTGGCCCGCTCGGCGAGCAGGTCACGGCTGGTGATGGCGATCAGGGGGTGGGGCTGGGCGTCTGCGGCGAGCAGCAGGCACAGGTCGGGGCCTATGCCGGCCGGCTCGCCGGGGGTGACGGCGAAGCGCAGGGGCTTCACTGGGCGGCCt       >  1:4473868/1‑150 (MQ=255)
                                                                                                                                                tGGCCCGCTCGGCGAGCAGGTCACGGCTGGTGATGGCGATCAGGGGGTGGGGCTGGGCGTCTGCGGCGAGCAGCAGGCACAGGTCGGGGCCTATGCCGGCCGGCTCGCCGGGGGTGACGGCGAAGCGCAGGGGCTTCACTGGGCGGCCTg      >  2:370370/1‑150 (MQ=255)
                                                                                                                                                    ccGCTCGGCGAGCAGGTCACGGCTGGTGATGGCGATCAGGGGGTGGGGCTGGGCGTCTGCGGCGAGCAGCAGGCACAGGTCGGGGCCTATGCCGGCCGGCTCGCCGGGGGTGACGGCGAAGCGCAGGGGCTTCACTGGGCGGCCTGGTCa  >  2:6561475/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CTGCGTTGGCCTTGTCGAGCTGGCCTGGCACCACCGGGGCCGCCAAGGGGGTATCCCAGACGTACAGGCTACCTGCGGGTGCCGGGAGGTCAGGCCATTGGCCCGGGGCTACCGGCAGCAGGCTGACGGCCAGCCCCAGCTGCGTGGCCCGCTCGGCGAGCAGGTCACGGCTGGTGATGGCGATCAGGGGGTGGGGCTGGGCGTCTGCGGCGAGCAGCAGGCACAGGTCGGGGCCTATGCCGGCCGGCTCGCCGGGGGTGACGGCGAAGCGCAGGGGCTTCACTGGGCGGCCTGGTCA  >  NC_002947/488022‑488319

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: