Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,046,068 T→A 17.7% intergenic (‑167/+100) ygfF ← / ← gcvP putative NAD(P)‑dependent oxidoreductase/glycine decarboxylase, PLP‑dependent, subunit P of glycine cleavage complex

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,046,0680TA17.7% 1097.5 / 66.1 455intergenic (‑167/+100)ygfF/gcvPputative NAD(P)‑dependent oxidoreductase/glycine decarboxylase, PLP‑dependent, subunit P of glycine cleavage complex
Reads supporting (aligned to +/- strand):  ref base T (109/249);  new base A (88/3);  total (201/254)
Fisher's exact test for biased strand distribution p-value = 1.08e-33
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 475 total aligned reads displayed.

TAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATCGGTACGCAGGAGCAGAACAGGTTACGGTCGC  >  NC_000913/3045921‑3046217
                                                                                                                                                   |                                                                                                                                                     
taattCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCctct                                                                                                                                                     >  2:2824269/1‑150 (MQ=255)
actattttctaaaggcgcttcggcgcctttttagtcagatgacaaagtacaaaagtgctcagacagtCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCctct                                                                                                                                                     >  1:3333157‑M1/68‑150 (MQ=255)
 aattCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCtctc                                                                                                                                                    <  2:7421212/150‑1 (MQ=255)
   ttCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTccc                                                                                                                                                  <  2:6196927/150‑1 (MQ=255)
          aCGTATCTTTGTGATTGCTAACAAACTGGCTAAATCTTCAATGTTGGCACTCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCTCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAggg                                                                                                                                           <  2:1689295/150‑1 (MQ=255)
          aCGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCAAACCCACACCCCAAAggg                                                                                                                                           >  1:4240053/1‑150 (MQ=255)
          aCGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTACCCTCACCCGAACCCTCACCCCAAAggg                                                                                                                                           >  2:3919601/1‑150 (MQ=255)
             tATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAAACCTATACCCAAAGGGGCg                                                                                                                                        >  1:4713002/1‑150 (MQ=255)
                cTTTGAGATTGCTAACAAACTGGCTCAATCTTAAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAgg                                                                                                                                     <  1:2156571/150‑1 (MQ=255)
                 tttGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCTGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAggg                                                                                                                                    <  1:6204264/150‑1 (MQ=255)
                 tttGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAggg                                                                                                                                    >  2:1137206/1‑150 (MQ=255)
                  ttGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAgggg                                                                                                                                   >  2:6824809/1‑150 (MQ=255)
                       aTTGCTAACAAACTGGCTCAATCTCCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCCCCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGt                                                                                                                              <  1:2586236/150‑1 (MQ=255)
                            tAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAg                                                                                                                                      >  1:130966/1‑137 (MQ=255)
                            tAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAg                                                                                                                         >  2:5182874/1‑150 (MQ=255)
                            tAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAg                                                                                                                         >  2:7061457/1‑150 (MQ=255)
                            tAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCAAACCCACAACCCAAACGGGGGAGGGGCCGGTCTGAg                                                                                                                         >  2:8087747/1‑150 (MQ=255)
                            tAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGGTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAg                                                                                                                         <  1:1103600/150‑1 (MQ=255)
                             aaCAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCGAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGc                                                                                                                        >  2:183828/1‑150 (MQ=255)
                              aCAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCCAAGGGGCGAGGGGACTGTCTGAGCa                                                                                                                       >  2:1395331/1‑150 (MQ=255)
                              aCAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCa                                                                                                                       <  1:4929936/150‑1 (MQ=255)
                              aCAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCGAAACCACTCCCCAAAGGGGCGAGGGGACTGTCTGAGCa                                                                                                                       >  1:4653866/1‑150 (MQ=255)
                                          aaTCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCACACCCTAACCCACACCCCCAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACttt                                                                                                           >  1:5428472/1‑150 (MQ=255)
                                          aaTCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTGTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACttt                                                                                                           <  2:4188536/150‑1 (MQ=255)
                                           aTCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCGGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTg                                                                                                          >  2:5241552/1‑150 (MQ=255)
                                             cTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTc                                                                                                        <  1:6892265/150‑1 (MQ=255)
                                             cTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTc                                                                                                        <  2:5000416/150‑1 (MQ=255)
                                              ttCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCAAAACCACTCCCCAAAGGGGGGGGGGGACTGTCTGTGCACGTTTGTACGTTTTCa                                                                                                       >  1:7995785/1‑150 (MQ=255)
                                                cAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCACACCCAAACCCTATCCCCCAACGGGCGAGGGGACAGTCTGAGCACTTTTGTACTTTGTCATc                                                                                                     >  2:3225922/1‑150 (MQ=255)
                                                cAATGTTGGCACGCTTTTCCCCCTCGCCCTGTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGGTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATc                                                                                                     <  1:7502644/150‑1 (MQ=255)
                                                 aaTGTTGGCACTCATTTCCCCCTCCCCCAAACAGGGAGAGGCCCGGGGTGAGCGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTATGTACATTGTCATCt                                                                                                    <  2:2171262/150‑1 (MQ=25)
                                                   tGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGa                                                                                                  <  2:4425306/150‑1 (MQ=255)
                                                   tGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGa                                                                                                  <  2:7528494/150‑1 (MQ=255)
                                                   tGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGa                                                                                                  <  2:348770/150‑1 (MQ=255)
                                                    gTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGAc                                                                                                 <  1:6557854/150‑1 (MQ=255)
                                                    gTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCGCACCCCAAACCACACACCAAAACGGGGGGGGGACTGTGTGAGCACATTTGTACTTTGTCATCAGaa                                                                                                 >  1:7135403/1‑149 (MQ=255)
                                                       ggCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTaa                                                                                              <  2:3361567/150‑1 (MQ=255)
                                                       ggCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTaa                                                                                              <  2:4313895/150‑1 (MQ=255)
                                                         cACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTaaaa                                                                                            >  2:5576627/1‑150 (MQ=255)
                                                          aCGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGAATaaaaa                                                                                           <  2:6515391/150‑1 (MQ=255)
                                                            gCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGCCGGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAgg                                                                                         >  2:1038289/1‑150 (MQ=255)
                                                             cTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGgc                                                                                        <  2:582243/150‑1 (MQ=255)
                                                              ttttCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGGTCCCCTCACCCGAAACCACGCCCCACACGGGCGAGGGGCAGGTCTGAGCCCGTTTGGGATTTGTGATATTACTtaaaatggg                                                                                       >  1:7860465/1‑146 (MQ=255)
                                                               tttCCCCCTCGCCCTTTCTGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGgcgc                                                                                      <  1:3699010/150‑1 (MQ=255)
                                                               tttCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGgcgc                                                                                      <  2:1567381/150‑1 (MQ=255)
                                                               tttCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGgcgc                                                                                      <  1:3927889/150‑1 (MQ=255)
                                                               tttCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCACACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACGAAAAAGgcgc                                                                                      >  2:8390902/1‑150 (MQ=255)
                                                                    CCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAg                                                                                 >  2:7718459‑M1/1‑150 (MQ=255)
                                                                     CCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAgc                                                                                <  2:8170428‑M1/150‑1 (MQ=255)
                                                                       TCGCCCTTTCAGGGAGTGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAgcgc                                                                              >  1:164232‑M1/1‑150 (MQ=255)
                                                                         GCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCGAACCCTCTCCCCCAAGGGGCGAGGGGACGGGCTGAGCACGGTTGTACGTTGTCAACTGACTAAAAAGGCGCCGAAGCgccg                                                                            <  2:567878‑M1/150‑2 (MQ=255)
                                                                          CCCTTTCCGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCtt                                                                           >  1:8385651‑M1/1‑150 (MQ=255)
                                                                            CTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTa                                                                         >  2:8316078‑M1/1‑150 (MQ=255)
                                                                              ttCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGa                                                                       <  1:98824/150‑1 (MQ=255)
                                                                              ttCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGa                                                                       <  1:6735476/150‑1 (MQ=255)
                                                                                 aGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGaaaa                                                                    <  1:8251279/150‑1 (MQ=255)
                                                                                 aGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGaaaa                                                                    <  1:5267294/150‑1 (MQ=255)
                                                                                 aGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGGTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGaaaa                                                                    <  2:413771/150‑1 (MQ=255)
                                                                                 aGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGGTCCCCGCCCCCAAACCCACACCCCCACAGGGGGGGGGGGCTGTCGGAGCACATTTTTACTTTGTGATATGTCTAAAAAAGCCCCCCAACGCCGGTAGaaaa                                                                    >  1:7204051/1‑150 (MQ=255)
                                                                                  gggAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAAt                                                                   <  1:4127005/150‑1 (MQ=255)
                                                                                   ggAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATa                                                                  <  1:2824269/150‑1 (MQ=255)
                                                                                   ggAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATa                                                                  <  1:2280827/150‑1 (MQ=255)
                                                                                     agagGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGt                                                                <  2:6991211/150‑1 (MQ=255)
                                                                                      gagGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTc                                                               <  2:4009732/150‑1 (MQ=255)
                                                                                         ggCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCACTACCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAAAAGTCGaa                                                            >  2:210507/1‑150 (MQ=255)
                                                                                         ggCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGGTGCCCACACCCAAAACCACACCCCCAACGGGGGAGGGGCCTGTGTGAGCACATTTGTAATTTGTCATAAGACAAAAAATTCGCCGCAGCGCCGTTAGAAAAAAgaagaa                                                            >  1:5644301/1‑145 (MQ=255)
                                                                                          gCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTTACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAAt                                                           >  2:5965748/1‑150 (MQ=255)
                                                                                          gCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAAt                                                           <  2:3678185/150‑1 (MQ=255)
                                                                                          gCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAAt                                                           <  2:7873250/150‑1 (MQ=255)
                                                                                          gCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCACACCCAAACCCAATCCCCAAAGGGGCGACGGGACTGTCAGAGCACTTTTGTACAATTTCATCTGACTAAAAAGGCGCCCAAGCGCCGTTAGCAAATAGACGaaa                                                           >  2:5138762/1‑149 (MQ=255)
                                                                                              gggTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGt                                                       <  1:1997972/150‑1 (MQ=255)
                                                                                              gggTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGt                                                       <  1:3185666/150‑1 (MQ=255)
                                                                                              gggTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGt                                                       <  1:4075673/150‑1 (MQ=255)
                                                                                                gTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGa                                                     <  2:7739497/150‑1 (MQ=255)
                                                                                                gTGAGGGTAAATAGTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGa                                                     <  1:5656118/150‑1 (MQ=255)
                                                                                                 tGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCGCACCCTAACCCACACCCCAAAGGGGCGAGGGGACTGTCTGAGCACGTTTGTACTTTGTAATCTGAATAAAAAGGCGCCGAAGCGCCATTAGAAAATAGTCGAATCAGTGaa                                                    >  2:1528429/1‑150 (MQ=39)
                                                                                                      gTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACt                                               <  1:1153646/150‑1 (MQ=255)
                                                                                                      gTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACt                                               >  2:4791022/1‑150 (MQ=255)
                                                                                                      gTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAAACCTCACCCCAAAGGGGCGCGGGGACTGTCTGAGCACGTTTGTACTTTGTAATCTGACTAAAAAGGCGCCGAAGCGCCGTTAGAAAATAGTCGAATAAGTTAATTACt                                               >  2:2644868/1‑150 (MQ=255)
                                                                                                            cgtCGCGCCCGTGCCGGCCTGTTCCCCTCACCCTAACCCTCGCCCCCCAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATCGTCGAATCAGTGAATTACTGGTAtt                                         <  2:2954344/148‑1 (MQ=38)
                                                                                                              tCGCGCCAGTGCCGGCCTGTTCCCCACACCCAAACCCACTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCGTTAGAAAAAAGTCGAAACAGTGAATTAATGGTAttag                                       >  2:3289345/1‑148 (MQ=38)
                                                                                                               cgcgCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGc                                      <  2:2383425/150‑1 (MQ=255)
                                                                                                               cgcgCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGc                                      <  1:7393040/150‑1 (MQ=255)
                                                                                                                 cgcCAGTGCCGGCCTGTGCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTa                                    <  2:7931579/150‑1 (MQ=255)
                                                                                                                    cAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATc                                 <  1:4326957/150‑1 (MQ=255)
                                                                                                                       tGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATCGGt                              <  2:4713002/150‑1 (MQ=255)
                                                                                                                         ccGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATCGGTAc                            <  2:6307207/150‑1 (MQ=255)
                                                                                                                             ccTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATCGGTACGCAg                        <  1:2696964/150‑1 (MQ=255)
                                                                                                                             ccTGGTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATCGGTACGCAg                        <  1:5806762/150‑1 (MQ=255)
                                                                                                                                 ttCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATCGGTACGCAGGAGc                    >  2:714658/1‑150 (MQ=255)
                                                                                                                                    cccTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATCGGTACGCAGGAGCAGa                 >  2:6671701/1‑150 (MQ=255)
                                                                                                                                     ccTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATCGGTACGCAGGAGCAGaa                <  2:2068278/150‑1 (MQ=255)
                                                                                                                                      cTCACCCTAACCCTCTCCCCAACGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATCGGTACGCAGGAGCAGAAc               <  1:175364/150‑1 (MQ=255)
                                                                                                                                      cTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAATGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATCGGTACGCAGGAGCAGAAc               >  1:5993975/1‑150 (MQ=255)
                                                                                                                                      cTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATCGGTACGCAGGAGCAGAAc               >  1:6071755/1‑150 (MQ=255)
                                                                                                                                        cACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATCGGTACGCAGGAGCAGAACAg             <  1:2480386/150‑1 (MQ=255)
                                                                                                                                           ccTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATCGGTACGCAGGAGCAGAACAGGtt          >  1:8149242/1‑150 (MQ=255)
                                                                                                                                            cTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCCTCTGACTAAACAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATCGGTACGCAGGAGCAGAACAGGTTa         <  2:2722152/150‑1 (MQ=255)
                                                                                                                                            cTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATCGGTACGCAGGAGCAGAACAGGTTa         <  2:1278321/150‑1 (MQ=255)
                                                                                                                                                   tctcCCCCAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCGTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATCGGTACGCAGGAGCAGAACAGGTTACGGTCGc  >  2:261464/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
TAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATCGGTACGCAGGAGCAGAACAGGTTACGGTCGC  >  NC_000913/3045921‑3046217

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