Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,046,050 T→G 12.3% intergenic (‑149/+118) 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,0500TG12.3% 296.5 / 9.9 143intergenic (‑149/+118)ygfF/gcvPputative NAD(P)‑dependent oxidoreductase/glycine decarboxylase, PLP‑dependent, subunit P of glycine cleavage complex
Reads supporting (aligned to +/- strand):  ref base T (37/85);  new base G (11/10);  total (48/95)
Fisher's exact test for biased strand distribution p-value = 7.75e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 146 total aligned reads displayed.

TTGTTTTTCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATCGGTACGCAGGAGC  >  NC_000913/3045902‑3046199
                                                                                                                                                    |                                                                                                                                                     
ttGTTTTTCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCCGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGtt                                                                                                                                                      <  1:2371931/150‑1 (MQ=255)
   tttttCCGATAAAATTTATTTCTTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTccc                                                                                                                                                   <  1:1022625/150‑1 (MQ=255)
   tttttCCGAGAAAAATTACTTCGTTGACGTATCTTTGAGATTGCTCACCAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCGCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTccc                                                                                                                                                   <  2:2181937/150‑1 (MQ=255)
   tttttCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTccc                                                                                                                                                   <  2:927411/150‑1 (MQ=255)
    ttttCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTcccc                                                                                                                                                  >  2:2099081/1‑150 (MQ=255)
    tttaCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCTCGCCAGTGCCGGCCTGGTgccc                                                                                                                                                  >  1:911400/1‑150 (MQ=255)
     tttCCGAGAAAAATTATTTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCt                                                                                                                                                 <  2:1819905/150‑1 (MQ=255)
     tttCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCt                                                                                                                                                 <  2:1090739/150‑1 (MQ=255)
       tCCGAGACAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCa                                                                                                                                               <  2:1008539/150‑1 (MQ=255)
       tCCGAGAAATATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCa                                                                                                                                               <  2:886795/150‑1 (MQ=255)
       tCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCa                                                                                                                                               <  1:372741/150‑1 (MQ=255)
       tCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGAGTGCTCACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCa                                                                                                                                               <  2:2469990/150‑1 (MQ=255)
         cGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCAcc                                                                                                                                             >  2:523311/1‑150 (MQ=255)
         cGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGGTCCCCTCAcc                                                                                                                                             >  1:527075/1‑150 (MQ=255)
             aaaaaTTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCACACCCaaa                                                                                                                                         >  1:1268112/1‑150 (MQ=255)
                  ttaattCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATCTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGGTCCCCACACCCAAAaccac                                                                                                                                    >  1:1337153/1‑148 (MQ=255)
                   taattCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCctct                                                                                                                                   <  2:2678733/150‑1 (MQ=255)
                   taattCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCctct                                                                                                                                   <  2:2654926/150‑1 (MQ=255)
                    aattCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCtctc                                                                                                                                  <  1:930280/150‑1 (MQ=255)
                          ttGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCTGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCaaa                                                                                                                            <  1:1461466/150‑1 (MQ=255)
                          ttGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCaaa                                                                                                                            <  1:2234065/150‑1 (MQ=255)
                               gTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCAAACCCACACCCCAAAGGGGc                                                                                                                       >  1:2392700/1‑150 (MQ=255)
                                 aTCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGa                                                                                                                     <  1:618925/150‑1 (MQ=255)
                                    tttGTGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGGGAGGGTAAATCTTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAggg                                                                                                                  <  1:2815299/150‑1 (MQ=255)
                                    tttGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAggg                                                                                                                  <  2:2458247/150‑1 (MQ=255)
                                     ttGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGGTACCCACACCCCAAACCAAACCACAAAGGGGCgggggg                                                                                                                 >  2:500169/1‑150 (MQ=255)
                                        cgaTTGCTACCCCACGGGCTCCATCTTCAAGGGTGGCACGCTTGTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATAGTCGCGCCAGTGCCGGCCGGGTCCCCTCACCCTAACCCTCTCCCCACAGGGGCGAGGGGACt                                                                                                              <  2:597425/149‑1 (MQ=21)
                                        agaTTGCTCACACACTGGCTCAATCTTCAATGTGGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGGTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACt                                                                                                              <  2:597429/150‑1 (MQ=255)
                                        agaTTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGGTCCCCTCACCCTAACCCCCGCCCCTTTGGGGAGAgggttag                                                                                                              <  1:1601399/150‑5 (MQ=255)
                                         gaTTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTg                                                                                                             >  1:1943588/1‑150 (MQ=255)
                                            tGCTACCCAACTGGCTCAATCTTCAATGTTGGCACGCTTTACCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTca                                                                                                          <  1:300816/150‑2 (MQ=255)
                                            tGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGGTGCCCACACCCAAACCCACACCCCAAAGGGGCGCGGGGACTGTCt                                                                                                          >  1:2004414/1‑150 (MQ=255)
                                              cTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCCGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGGTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGa                                                                                                        <  1:1848896/150‑1 (MQ=255)
                                              cTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGa                                                                                                        <  2:880386/150‑1 (MQ=255)
                                               tAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCAAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAg                                                                                                       >  2:171318/1‑150 (MQ=255)
                                                 aCAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCCAACCCTAAACCGCACCCCAAAAGGGCGAGGGGACGGTCTGAGcc                                                                                                     >  1:1050848/1‑149 (MQ=255)
                                                  cAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCAc                                                                                                    <  1:2365867/150‑1 (MQ=255)
                                                   ataCTGGCTCAATCTTCATTGTTGGCACGCTTTTCCCCCTAGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCAAAAAGGGTCGAGGGGACTGTCTTAGCACt                                                                                                   <  1:2761472/148‑1 (MQ=25)
                                                   aaaCTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACt                                                                                                   <  1:2546660/150‑1 (MQ=255)
                                                     aCTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCAAACCCACACCCCCAAGGGGCGAGGGGACTGTATGAGCACttt                                                                                                 >  2:2851638/1‑150 (MQ=255)
                                                         gCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGTGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTa                                                                                             <  2:1023706/150‑1 (MQ=255)
                                                         gCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTa                                                                                             <  2:2905379/150‑1 (MQ=255)
                                                         gCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTa                                                                                             <  2:220813/150‑1 (MQ=255)
                                                         gCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTa                                                                                             <  2:220703/150‑1 (MQ=255)
                                                          cTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTAc                                                                                            >  2:2156163/1‑150 (MQ=255)
                                                          cTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCAAACCCTCTCCCCAAAGGGGCGCGGGGACTGTCGGAGGACTTTTGTAc                                                                                            >  1:410905/1‑150 (MQ=255)
                                                          cTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCAAACCCACTCCCCAAAGGGGCGAGGGGACTGTGTGAGCACTTTTGTAc                                                                                            >  1:410910/1‑150 (MQ=255)
                                                           tCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCAAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTGTTGTACt                                                                                           >  2:1233306/1‑150 (MQ=255)
                                                             aaTCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCCGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCGGTTCCCCACACCCAAACCCACACCCCAAACGGGCGCGGGGACAGTCTGTGCACgtgtgttaatt                                                                                         >  2:2322226/1‑145 (MQ=255)
                                                             aaTCGTCAATGTTGGCACGCTTGTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATAGTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACttt                                                                                         <  2:1430787/150‑1 (MQ=255)
                                                               tCTTCAATGTTGGCACGCTTTTCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGATTGTCTGAGCACTTTTGTACTTTGTc                                                                                      <  2:1492930/150‑1 (MQ=255)
                                                                cTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGGTCCCCACACCCAAAAACACACCCCAAAAGGGCGAGGGGCCAGTATGAGCACAGTTGTAATTTTTc                                                                                      >  1:899115/1‑150 (MQ=255)
                                                                 ttCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGGTGCCCGCCCCCAAACCCACAACCCCAAGGGGCGAGGGGGCGGGCTGTGCACTTTTGTTCtttttaa                                                                                     >  1:878476/1‑146 (MQ=255)
                                                                   cAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATc                                                                                   <  1:2615159/150‑1 (MQ=255)
                                                                      tGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGa                                                                                <  1:2545513/150‑1 (MQ=255)
                                                                      tGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGa                                                                                <  1:1667973/150‑1 (MQ=255)
                                                                       gTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGAc                                                                               <  1:1713361/150‑1 (MQ=255)
                                                                                 ttttCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGgcg                                                                     <  1:1058458/150‑1 (MQ=255)
                                                                                  tttCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGgcgc                                                                    <  1:1794243/150‑1 (MQ=255)
                                                                                  tttCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGgcgc                                                                    <  1:2165552/150‑1 (MQ=255)
                                                                                        ccTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAgc                                                              <  1:1905864/150‑1 (MQ=255)
                                                                                        ccTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCGAACCCGAAACCACTCCCCCAAGGGGCGCGGGGCCGGTCTGGGCACTGTTTTACGTTGGCATCTGACTAAAAAGGCGCCGAAgc                                                              >  2:587430/1‑150 (MQ=255)
                                                                                         cTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTGCCCACACCCGAACCCTCACCCCCAAGGGGCGCGGGGGCGGTCTGTGGCCATTTTTACTTTGTGAACTGACTAAAAAGGCGCCgccgcg                                                             >  1:1404184/1‑145 (MQ=255)
                                                                                         cTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCCCAGGGGCGAGGGGACTGTCTGAGCACGTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGCAgcg                                                             >  2:2644279/1‑150 (MQ=255)
                                                                                          tCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAgcgc                                                            >  2:2922795/1‑150 (MQ=255)
                                                                                          tCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCACACCCCAAAGGGGCGGGGGGACTGTCGTAGCACGGTTGTACGTTGTCATCTGACTAAAAATGCGCCGAAgcgc                                                            >  1:718434/1‑150 (MQ=25)
                                                                                                 ttCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGGTCCCCGCACCCCAACCCGCTCCCCCAAGGGGCCAGGGGACTGTCTGAGCACTTTTGTACATTGTCATATGACTAAAAAAGCGCCGAAGCGCCGTTAGa                                                     >  2:918751/1‑150 (MQ=255)
                                                                                                  tCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGaa                                                    <  1:1598885/150‑1 (MQ=255)
                                                                                                      ggAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATa                                                <  1:1529762/150‑1 (MQ=255)
                                                                                                      ggAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGGTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATa                                                <  2:1911213/150‑1 (MQ=255)
                                                                                                         gagGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTc                                             >  2:1938621/1‑150 (MQ=255)
                                                                                                         gagGGCCGGGGTGAGGGTAAATAGTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTc                                             <  2:574021/150‑1 (MQ=255)
                                                                                                          agGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCg                                            <  2:1247550/150‑1 (MQ=255)
                                                                                                           gggCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGGTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGa                                           <  2:1684252/150‑1 (MQ=255)
                                                                                                             gCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAAt                                         <  2:1654990/150‑1 (MQ=255)
                                                                                                                 gggTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGt                                     <  2:2659915/150‑1 (MQ=255)
                                                                                                                 gggTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCACACCCAAACCCACACCCCCAACGGGCGAGGGGGCTGTCTGAGCACGTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCGTTAGAAAATAGTCGAAACact                                     >  2:325187/1‑148 (MQ=255)
                                                                                                                 gggTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCACACCCAAAACCACTACCCAAAGGGGCGAGGGGACTGTCTGAGCACGTTTGTACTTTTTCAACAGACTAAAAATGCGCCGAAGCGCCTTTAGAAAATAGTCGAAACAGt                                     >  2:2809560/1‑150 (MQ=255)
                                                                                                                   gTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGGTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGa                                   <  1:2289956/150‑1 (MQ=255)
                                                                                                                     gAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAAt                                 <  1:575967/150‑1 (MQ=255)
                                                                                                                       gggTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTa                               <  1:523311/150‑1 (MQ=255)
                                                                                                                         gTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACt                             <  2:2004414/150‑1 (MQ=255)
                                                                                                                         gTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACt                             <  1:1571141/150‑1 (MQ=255)
                                                                                                                         gTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACt                             <  1:1618100/150‑1 (MQ=255)
                                                                                                                         gTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACt                             <  1:2648657/150‑1 (MQ=255)
                                                                                                                           aaaTATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCGAACCCGCTCCCCAAAGGGGCGAGGGGACGGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAATGCGCCCAAGCGCCGTTAGAAAAAAGTCGAATCAGTGAATTACTgg                           >  2:123943/1‑150 (MQ=38)
                                                                                                                               atTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTAtt                       <  2:1222776/150‑1 (MQ=255)
                                                                                                                                ttCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTc                      <  1:215995/150‑1 (MQ=255)
                                                                                                                                 tCGCGCCAGTGCCGGCCTGTTCCCCGCACCCTAAACCTCACCCCAAAGGGGCGGGGGGCCTGTCAGAGCACTTTTGTACTTTGTCAAATTACTAAAAAGGCGCCCAAGCGCCGTTAGAAAAAACTCGAAACAGTGAAATACAGGTaatag                     >  2:140020/1‑146 (MQ=255)
                                                                                                                                  cgcgCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGc                    <  1:476151/150‑1 (MQ=255)
                                                                                                                                       cAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATc               <  2:2352517/150‑1 (MQ=255)
                                                                                                                                       cAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCGCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATc               <  2:1408565/150‑1 (MQ=255)
                                                                                                                                       cAGTGCCGGCCTGGTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATc               <  2:2858943/150‑1 (MQ=255)
                                                                                                                                           gccggccgGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATCGGTa           <  2:2146913/150‑1 (MQ=255)
                                                                                                                                             cGGCCTGTTCCCCACACCCAAACCCTCACCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATCGGTACg         >  2:1225822/1‑150 (MQ=255)
                                                                                                                                              ggCCGGTTCCCCTCACCCTAACCCTCTCCCCACAGGGGCGAGGGGACTGTCGGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAAGAGTCGAATCAGTGAATTACTGGTATTCGCTAATCGGTACGc        <  2:1859624/150‑1 (MQ=255)
                                                                                                                                               gCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATCGGTACGCa       <  2:2392700/150‑1 (MQ=255)
                                                                                                                                                ccTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATCGGTACGCAg      <  1:1768254/150‑1 (MQ=255)
                                                                                                                                                  tGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCTCTAATCGGTACGCAGGa    >  2:124451/1‑150 (MQ=255)
                                                                                                                                                    ttCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATCGGTACGCAGGAGc  >  2:1272650/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TTGTTTTTCCGAGAAAAATTAATTCGTTGACGTATCTTTGAGATTGCTAACAAACTGGCTCAATCTTCAATGTTGGCACGCTTTTCCCCCTCGCCCTTTCAGGGAGAGGGCCGGGGTGAGGGTAAATATTCGCGCCAGTGCCGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAAGGGGCGAGGGGACTGTCTGAGCACTTTTGTACTTTGTCATCTGACTAAAAAGGCGCCGAAGCGCCTTTAGAAAATAGTCGAATCAGTGAATTACTGGTATTCGCTAATCGGTACGCAGGAGC  >  NC_000913/3045902‑3046199

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: