Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,556,501 C→T 5.0% G266D (GGC→GAC)  eutC ← ethanolamine ammonia‑lyase, small subunit (light chain)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,556,5010CT 5.0% 272.0 / 3.4 80G266D (GGC→GAC) eutCethanolamine ammonia‑lyase, small subunit (light chain)
Reads supporting (aligned to +/- strand):  ref base C (32/43);  new base T (2/2);  total (35/45)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

TCACCCGCATAGCGGTGACAGACGGTCGAATCAAATCTAAAGCTGGCATGATGCCTCCTTATCGGGTCATGTTGATGCCGGATGCTTTCTGCTCCAGCATACGTTTCGCCAAATCCACAATGACGGCTGCGGCTTCAACCGGCGGCGTGCCGCCCTGGTGAATGTTAGAGATACAGGTGCGATCGGCCTCAACGGTGGTCGCCATACGCGGCGAGTAAACGGCGTAGCAGGAGAGGCTTTCTGACTGACCCAACCCCGGACGTTCGCCCACCAGCAGGATCACTACTTTCGCGCCGAG  >  NC_000913/2556352‑2556649
                                                                                                                                                     |                                                                                                                                                    
tCACCCGCATAGCGGTGACAGACGGTCGAATCAAATCTAAAGCTGGCATGATGCCTCCTTATCGGGTCATGTTGATGCCGGATGCTTTCTGCTCCAGCATACGTTTCGCCAAATCCACAATGACGGCTGCGGCTTCAACCGGCGGCGTgc                                                                                                                                                      >  1:1397418/1‑150 (MQ=255)
  aCCCGCATAGCGGTGACAGACGGTCGAATCAAATCTACAGCTGGCATGATGCCTCCTTATCGGGTCATGTTGATGCCGGATGCTTTCTGCTCCAGCATACGTTTCGCCAAATCCACAATGACGGCTGCGGCTTCAACCGGCGGCGTgccg                                                                                                                                                    <  2:355359/150‑1 (MQ=255)
  aCCCGCATAGCGGTGACAGACGGTCGAATCAAATCTAAAGCTGGCATGATGCCTCCTTATCGGGTCATGTTGATGCCGGATGCTTTCTGCTCCAGCATACGTTTCGCCAAATCCACAATGACGGCTGCGGCTTCAACCGGCGGCGTgccg                                                                                                                                                    <  1:38063/150‑1 (MQ=255)
  aCCCGCATAGCGGTGACAGACGGTCGAATCAAATCTAAAGCTGGCATGATGCCTCCTTATCGGGTCATGTTGATGCCGGATGCTTTCTGCTCCAGCATACGTTTCGCCAAATCCACAATGACGGCTGCGGCTTCAACCGGCGGCGTgccg                                                                                                                                                    <  1:1291764/150‑1 (MQ=255)
    ccGCATAGCGGTGACAGACGGTCGAATCAAATCTAAAGCTGGCATGATGCCTCCTTATCGGGTCATGTTGATGCCGGATGCTTTCTGCTCCAGCATACGTTTCGCCAAATCCACAATGACGGCTGCGGCTTCAACCGGCGGCGTgccgcc                                                                                                                                                  <  1:1344856/150‑1 (MQ=255)
    ccGCATAGCGGTGACAGACGGTCGAATCAAATCTAAAGCTGGCATGATGCCTCCTTATCGGGTCATGTTGATGCCGGATGCTTTCTGCTCCAGCATACGTTTCGCCAAATCCACAATGACGGCTGCGGCTTCAACCGGCGGCGTgccgcc                                                                                                                                                  <  1:1344864/150‑1 (MQ=255)
    ccGCATAGCGGTGACAGACGGTCGAATCAAATCTAAAGCTGGCATGATGCCTCCTTATCGGGTCATGTTGATGCCGGATGCTTTCTGCTCCAGCATACGTTTCGCCAAATCCACAATGACGGCTGCGGCTTCAACCGGCGGCGTgccgcc                                                                                                                                                  >  2:150766/1‑150 (MQ=255)
    ccGCATAGCGGTGACAGACGGTCGAATCAAATCTAAAGCTGGCATGATGCCTCCTTATCGGGTCATGTTGATGCCGGATGCTTTCTGCTCCAGCATACGTTTCGCCAAATCCACAATGACGGCTGCGGCTTCAACCGGCGGCGTgccgcc                                                                                                                                                  >  2:291608/1‑150 (MQ=255)
      gCATAGCGGTGACAGACGGTCGAATCAAATCTAAAGCTGGCATGATGCCTCCTTATCGGGTCATGTTGATGCCGGATGCTTTCTGCTCCAGCATACGTTTCGCCAAATCCACAATGACGGCTGCGGCTTCAACCGGCGGCACgccgccgg                                                                                                                                                <  2:839020/150‑3 (MQ=255)
         tAGCGGTGACAGACGGTCGAATCAAATCTAAAGCTGGCATGATGCCTCCTTATCGGGTCATGTTGATGCCGGATGCTTTCTGCTCCAGCATACGTTTCGCCAAATCCACAATGACGGCTGCGGCTTCAACCGGCGGCGTGCCGCCCTGGt                                                                                                                                             <  1:579201/150‑1 (MQ=255)
         tAGCGGTGACAGACGGTCGAATCAAATCTAAAGCTGGCATGATGCCTCCTTATCGGGTCATGTTGATGCCGGATGCTTTCTGCTCCAGCATACGTTTCGCCAAATCCACAATGACGGCTGCGGCTTCAACCGGCGGCGTGCCGCCCTGGt                                                                                                                                             <  1:201559/150‑1 (MQ=255)
             ggTGACAGACGGTCGAATCAAATCTAAAGCTGGCATGATGCCTCCTTATCGGGTCATGTTGATGCCGGATGCTTTCTGCTCCAGCATACGTTTCGCCAAATCCACAATGAGGGCTGCGGCTTCAACCGGCGGCGTGCCGCCCTGGTGAAt                                                                                                                                         <  1:1307622/150‑1 (MQ=255)
                    gACGGTCGAATCAAATCTAAAGCTGGCATGATGCCTCCTTATCGGGTCATGTTGATGCCGGATGCTTTCTGCTCCAGCATACGTTTCGCCAAATCCACAATGACGGCTGCGGCTTCAACCGGCGGCGTGCCGCCCTGGTGAATGTTagag                                                                                                                                  <  2:1203415/150‑1 (MQ=255)
                     aCGGTCGAATCAAATCTAAAGCTGGCATGATGCCTCCTTATCGGGTCATGTTGATGCCGGATGCTTTCTGCTCCAGCATACGTTTCGCCAAATCCACAATGACGGCTGCGGCTTCAACCGGCGGCGTGCCGCCCTGGTGAATGTTAgaga                                                                                                                                 >  2:467373/1‑150 (MQ=255)
                        gTCGAATCAAATCTAAAGCTGGCATGATGCCTCCTTATCGGGTCATGTTGATGCCGGATGCTTTCTGCTCCAGCATACGTTTCGCCAAATCCACAATGACGGCTGCGGCTTCAACCGGCGGCGTGCCGCCCTGGTGAATGTTAGAGATAc                                                                                                                              <  1:218232/150‑1 (MQ=255)
                           gAATCAAATCTAAAGCTGGCATGATGCCTCCTTATCGGGTCATGTTGATGCCGGATGCTTTCTGCTCCAGCATACGGTTCGCCAAATCCACAATGACGGCTGCGGCTTCAACCGGCGGCGTGCCGCCCTGGTGCATGTTAGAGATACAgg                                                                                                                           >  1:447213/1‑150 (MQ=255)
                            aaTCAAATCTAAAGCTGGCATGATGCCTCCTTATCGGGTCATGTTGATGCCGGATGCTTTCTGCTCCAGCATACGTTTCGCCAAATCCACAATGACGGCTGCGGCTTCAACCGGCGGCGTGCCGCCCTGGTGAATGTTAGAGATACAGGt                                                                                                                          <  1:1004947/150‑1 (MQ=255)
                              tCAAATCTAAAGCTGGCATGATGCCTCCTTATCGGGTCATGTTGATGCCGGATGCTTTCTGCTCCAGCATACGTTTCGCCAAATCCACAATGACGGCTGCGGCTTCAACCGGCGGCGTGCCGCCCTGGTGACTGTTAGAGATACAGGTGc                                                                                                                        >  1:192072/1‑150 (MQ=255)
                              tCAAATCTAAAGCTGGCATGATGCCTCCTTATCGGGTCATGTTGATGCCGGATGCTTTCTGCTCCAGCATACGTTTCGCCAAATCCACAATGACGGCTGCGGCTTCAACCGGCGGCGTGCCGCCCTGGTGAATGTTAGAGATACAGGTGc                                                                                                                        >  1:269473/1‑150 (MQ=255)
                              tCAAATCTAAAGCTGGCATGATGCCTCCTTATCGGGTCATGTTGATGCCGGATGCTTTCTGCTCCAGCATACGTTTCGCCAAATCCACAATGACGGCTGCGGCTTCAACCGGCGGCGTGCCGCCCTGGTGAATGTTAGAGATACAGGTGc                                                                                                                        >  2:1030501/1‑150 (MQ=255)
                              tCAAATCTAAAGCTGGCATGATGCCTCCTTATCGGGTCATGTTGATGCCGGATGCTTTCTGCTCCAGCATACGTTACGCCAAATCCACAATGACGGCTGCGGCTTCAACCGGCGGCGTGCCGCCCTGGTGAATGTTAGAGATACAGGTGc                                                                                                                        >  2:685023/1‑150 (MQ=255)
                                 aaTCTAAAGCTGGCATGATGCCTCCTTATCGGGTCATGTTGATGCCGGATGCTTTCTGCTCCAGCATACGTTTCGCCAAATCCACAATGACGGCTGCGGCTTCAACCGGCGGCGTGCCGCCCTGGTGAATGTTAGAGATACAGGTGCGAt                                                                                                                     <  2:253112/150‑1 (MQ=255)
                                 aaTCTAAAGCTGGCATGATGCCTCCTTATCGGGTCATGTTGATGCCGGATGCTTTCTGCTCCAGCATACGTTTCGCCAAATCCACAATGACGGCTGCGGCTTCAACCGGCGGCGTGCCGCCCTGGTGAATGTTAGAGATACAGGTGCGAt                                                                                                                     <  2:253453/150‑1 (MQ=255)
                                      aaaGCTGGCATGATGCCTCCTTATCGGGTCATGTTGATGCCGGATGCTTTCTGCTCCAGCATACGTTTCGCCAAATCCACAATGACGGCTGCGGCTTCAACCGGCGGCGTGCCGCCCTGGTGAATGTTAGAGATACAGGTGCGATCGGcc                                                                                                                <  2:1296368/150‑1 (MQ=255)
                                      aaaGCTGGCATGATGCCTCCTTATCGGGTCATGTTGATGCCGGATGCTTTCTGCTCCAGCATACGTTTCGCCAAATCCACAATGACGGCTGCGGCTTCAACCGGCGGCGTGCCGCCCTGGTGAATGTTAGAGATACAGGTGCGATCGGcc                                                                                                                <  2:761742/150‑1 (MQ=255)
                                       aaGCTGGCATGATGCCTCCTTATCGGGTCATGTTGATGCCGGATGCTTTCTGCTCCAGCATACGTTTCGCCAAATCCACAATGACGGCTGCGGCTTCAACCGGCGGCGTGCCGCCCTGGTGAATGTTAGAGATACAGGTGCGATCGGCCt                                                                                                               >  1:649289/1‑150 (MQ=255)
                                         gCTGGCATGATGCCTCCTTATCGGGTCATGTTGATGCCGGATGCTTTCTGCTCCAGCATACGTTTCGCCAAATCCACAATGACGGCTGCGGCTTCAACCGGCGGCGTGCCGCCCTGGTGAATGTTAGAGATACAGGTGCGATCGGCCTCa                                                                                                             >  1:20850/1‑150 (MQ=255)
                                          cTGGCATGATGCCTCCTTATCGGGTCATGTTGATGCCGGATGCTTTCTGCTCCAGCATACGTTTCGCCAAATCCACAATGACGGCTGCGGCTTCAACCGGCGGCGTGCCGCCCTGGTGAATGTTAGAGATACAGGTGCGATCGGCCTCaa                                                                                                            <  2:1141500/150‑1 (MQ=255)
                                          cTGGCATGATGCCTCCTTATCGGGTCATGTTGATGCCGGATGCTTTCTGCTCCAGCATACGTTTCGCCAAATCCACAATGACGGCTGCGGCTTCAACCGGCGGCGTGCCGCCCTGGTGAATGTTAGAGATACAGGTGCGATCGGCCTCaa                                                                                                            >  2:467419/1‑150 (MQ=255)
                                           tGGCATGATGCCTCCTTATCGGGTCATGTTGATGCCGGATGCTTTCTGCTCCAGCATACGTTTCGCCAAATCCACAATGACGGCTGCGGCTTCAACCGGCGGCGTGCCGCCCTGGTGAATGTTAGAGATACAGGTGCGATCGGCCTCAAc                                                                                                           <  2:365316/150‑1 (MQ=255)
                                             gCATGATGCCTCCTTATCGGGTCATGTTGATGCCGGATGCTTTCTGCTCCAGCATACGTTTCGCCAAATCCACAATGACGGCTGCGGCTTCAACCGGCGGCGTGCCGCCCTGGTGAATGTTAGAGATACAGGTGCGATCGGCCTCAACgg                                                                                                         >  1:723052/1‑150 (MQ=255)
                                             gCATGATGCCTCCTTATCGGGTCATGTTGATGCCGGATGCTTTCTGCTCCAGCATACGTTTCGCCAAATCCACAATGACGGCTGCGGCTTCAACCGGCGGCGTGCCGCCCGGGTGAATGTTAGAGATACAGGTGCGATCGGCCTCAACgg                                                                                                         >  2:188238/1‑150 (MQ=255)
                                              cATGATGCCTCCTTATCGGGTCATGTTGATGCCGGATGCTTTCTGCTCCAGCATACGTTTCGCCAAATCCACAATGACGGCTGCGGCTTCAACCGGCGGCGTGCCGCCCTGGTGAATGTTAGAGATACAGGTGCGATCGGCCTCAACggt                                                                                                        <  2:474287/150‑1 (MQ=255)
                                              cATGATGCCTCCTTATCGGGTCATGTTGATGCCGGATGCTTTCTGCTCCAGCATACGTTTCGCCAAATCCACAATGACGGCTGCGGCTTCAACCGGCGGCGTGCCGCCCTGGTGAATGTTAGAGATACAGGTGCGATCGGCCTCAACggt                                                                                                        <  2:476038/150‑1 (MQ=255)
                                               atgatgCCTCCTTATCGGGTCATGTTGATGCCGGATGCTTTCTGCTCCAGCATACGTTTCGCCAAATCCACAATGACGGCTGCGGCTTCAACCGGCGGCGTGCCGCCCTGGTGAATGTTAGAGATACAGGTGCGATCGGCCTCAACggtg                                                                                                       <  2:1247846/150‑1 (MQ=255)
                                                tgatgCCTCCTTATCGGGTCATGTTGATGCCGGATGCTTTCTGCTCCAGCATACGTTTCGCCAAATCCACAATGACGGCTGCGGCTTCAACCGGCGGCGTGCCGCCCTGGTGAATGTTAGAGATACAGGTGCGATCGGCCTCAACggtgg                                                                                                      >  2:1244627/1‑150 (MQ=255)
                                                tgatgCCTCCTTATCGGGTCATGTTGATGCCGGATGCTTTCTGCTCCAGCATACGTTTCGCCAAATCCACAATGACGGCTGCGGCTTCAACCGGCGGCGTGCCGCCCTGGTGAATGTTAGAGATACAGGTGCGATCGGCCTCAACggtgg                                                                                                      <  2:984931/150‑1 (MQ=255)
                                                 gatgCCTCCTTATCGGGTCATGTTGATGCCGGATGCTTTCTGCTCCAGCATACGTTTCGCCAAATCCACAATGACGGCTGCGGCTTCAACCGGCGGCGTGTCGCCCTGGTGAATGTTAGAGATACAGGTGCGATCGGCCTCAACggtggt                                                                                                     >  2:726400/1‑150 (MQ=255)
                                                 gatgCCTCCTTATCGGGTCATGTTGATGCCGGATGCTTTCTGCTCCAGCATACGTTTCGCCAAATCCACAATGACGGCTGCGGCTTCAACCGGCGGCGTGTCGCCCTGGTGAATGTTAGAGATACAGGTGCGATCGGCCTCAACggtggt                                                                                                     >  2:726419/1‑150 (MQ=255)
                                                     cctcctTATCGGGTCATGTTGATGCCGGATGCTTTCTGCTCCAGCATACGTTTCGCCAAATCCACAATGACGGCTGCGGCTTCAACCGGCGGCGTGCCGCCCTGGTGAATGTTAGAGATACAGGTGCGATCGGCCTCAACGGTGGTCGcc                                                                                                 >  2:1236262/1‑150 (MQ=255)
                                                     cctcctTATCGGGTCATGTTGATGCCGGATGCTTTCTGCTCCAGCATACGTTTCGCCAAATCCACAATGACGGCTGCGGCTTCAACCGGCGGCGTGCCGCCCTGGTGAATGTTAGAGATACAGGTGCGATCGGCCTCAACGGTGGTCGcc                                                                                                 >  2:155909/1‑150 (MQ=255)
                                                     cctcctTATCGGGTCATGTTGATGCCGGATGCTTTCTGCTCCAGCATACGTTTCGCCAAATCCACAATGACGGCTGCGGCTTCAACCGGCGGCGTGCCGCCCTGGTGAATGTTAGAGATACAGGTGCGATCGGCCTCAACGGTGGTCGcc                                                                                                 >  2:469105/1‑150 (MQ=255)
                                                         ctTATCGGGTCATGTTGATGCCGGATGCTTTCTGCTCCAGCATACGTTTCGCCAAATCCACAATGACGGCTGCGGCTTCAACCGGCGGCGTGCCGCCCTGGTGAATGTTAGAGATACAGGTGCGATCGGCCTCAACGGTGGTCGCCATAc                                                                                             <  2:1061397/150‑1 (MQ=255)
                                                          ttATCGGGTCATGTTGATGCCGGATGCTTTCTGCTCCAGCATACGTTTCGCCAAATCCACAATGACGGCTGCGGCTTCAACCGGCGGCGTGCCGCCCTGGTGAATGGTAGAGATACAGGTGCGATCGGCCTCAACGGTGGTCGCCATAcg                                                                                            >  2:914261/1‑150 (MQ=255)
                                                           tATCGGGTCATGTTGATGCCGGATGCTTTCTGCTCCAGCATACGTTTCGCCAAATCCACAATGACGGCTGCGGCTTCAACCGGCGGCGTGCCGCCCTGGTGAATGTTAGAGATACAGGTGCGATCGGCCTCAACGGTGGTCGCCATAcgc                                                                                           >  2:546612/1‑150 (MQ=255)
                                                             tCGGGTCATGTTGATGCCGGATGCTTTCTGCTCCAGCATACGTTTCGCCAAATCCACAATGACGGCTGCGGCTTCAACCGGCGGCGTGCCGCCCTGGTGAATGTTAGAGATACAGGTGCGATCGGCCTCAACGGTGGTCGCCATACgcgg                                                                                         <  2:1117515/150‑1 (MQ=255)
                                                             tCGGGTCATGTTGATGCCGGATGCTTTCTGCTCCAGCATACGTTTCGCCAAATCCACAATGACGGCTGCGGCTTCAACCGGCGGCGTGCCGCCCTGGTGAATGTTAGAGATACAGGTGCGATCGGCCTCAACGGTGGTCGCCATACgcgg                                                                                         <  2:8676/150‑1 (MQ=255)
                                                                ggTCATGTTGATGCCGGATGCTTTCTGCTCCAGCATACGTTTCGCCAAATCCACAATGACGGCTGCGGCTTCAACCGGCGGCGTGCCGCCCTGGTGAATGTTAGAGATACAGGTGCGATCGGCCTCAACGGTGGTCGCCATACGCGGCGa                                                                                      <  2:1092955/150‑1 (MQ=255)
                                                                ggTCATGTTGATGCCGGATGCTTTCTGCTCCAGCATACGGTTCGCCAAATCCACAATGACGGCTGCGGCTTCAACCGGCGGCGTGCCGCCCTGGTGAATGTTAGAGATACAGGTGCGATCGGCCTCAACGGTGGTCGCCATACGCGGCGa                                                                                      <  1:829682/150‑1 (MQ=255)
                                                                   cATGTTGATGCCGGATGCTTTCTGCTCCAGCATACGTTTCGCCAAATCCACAATGACGGCTGCGGCTTCAACCGGCGGCGTGCCGCCCTGGTGAATGTTAGAGATACAGGTGCGATCGGCCTCAACGGTGGTCGCCATACGCGGCGAGTa                                                                                   <  1:331791/150‑1 (MQ=255)
                                                                     tGTTGATGCCGGATGCTTTCTGCTCCAGCATACGTTTCGCCAAATCCACAATGACGGCTGCGGCTTCAACCGGCGGCGTGCCGCCCTGGTGAATGTTAGAGATACAGGTGCGATCGGCCTCAACGGTGGTCGCCATACGCGGCGAGTaaa                                                                                 <  2:1076348/150‑1 (MQ=255)
                                                                      gTTGATGCCGGATGCTTTCTGCTCCAGCATACGTTTCGCCAAATCCACAATGACGGCTGCGGCTTCAACCGGCGGCGTGCCGCCCTGGTGAATGTTAGAGATACAGGTGCGATCGGCCTCAACGGTGGTCGCCATACGCGGCGAGTAAAc                                                                                >  2:1266988/1‑150 (MQ=255)
                                                                         gaagCCGGATGCTTTCTGCTCCAGCATACGTTTCGCCAAATCCACAATGACGGCAGCGGCTTCAACCGGCGGCGTGCCGCCCTGGTGAATGTTAGAGATACAGGTGCGATCGGCCTCAACGGTGGTCGCCATACGCGGCGAGTAAACGGc                                                                             >  2:32136/4‑150 (MQ=255)
                                                                          aTGCCGGATGCTTTCTGCTCCAGCATACGTTTCGCCAAATCCACAATGACGGCTGCGGCTTCAACCGGCGGCGTGCCGCCCTGGTGAATGGTAGAGATACAGGTGCGATCGGCCTCAACGGTGGTCGCCATACGCGGCGAGTAAACGGCg                                                                            >  2:1389399/1‑150 (MQ=255)
                                                                             ccGGATGCTTTCTGCTCCAGCATACGTTTCGCCAAATCCACAATGACGGCTGCGGCTTCAACCGGCGGCGTGCCGCCCTGGTGAATGTTAGAGATACAGGTGCGATCGGCCTCAACGGTGGTCGCCATACGCGGCGAGTAAACGGCGTAg                                                                         >  2:1140907/1‑150 (MQ=255)
                                                                              cGGATGCTTTCTGCTCCAGCATACGTTTCGCCAAATCCACAATGACGGCTGTGGCTTCAACCGGCGGCGTGCCGCCCTGGTGAATGTTAGAGATACAGGTGCGATCGGCCTCAACGGTGGTCGCCATACGCGGCGAGTAAACGGCGTAGc                                                                        >  1:153341/1‑150 (MQ=255)
                                                                              cGGATGCTTTCTGCTCCAGCATACGTTTCGCCAAATCAACAATGACGGCTGCGGCTTCAACCGGCGGCGTGCCGCCCTGGTGAATGTTAGAGATACAGGTGCGATCGGCCTCAACGGTGGTCGCCATACGCGGCGAGTAAACTGCGTAGc                                                                        >  2:715967/1‑150 (MQ=255)
                                                                                 aTGCTTTCTGCTCCAGCATACGTTTCGCCAAATCCACAATGACGGCTGCGGCTTCAACCGGCGGCGTGCCGCCCTGGTGAATGTTAGAGATACAGGTGCGATCGGCCTCAACGGTGGTCGCCATACGCGGCGAGTAAACGGCGTAGCAgg                                                                     <  2:526338/150‑1 (MQ=255)
                                                                                    cTTTCTGCTCCAGCATACGTTTCGCCAAATCCACAATGACGGCTGCGGCTTCAACCGGCGGCGTGCCGCCCTGGTGAATGTTAGAGATACAGGTGCGATCGGCCTCAACGGTGGTCGCCATACGCGGCGAGTAAACGGCGTAGCAGgaga                                                                  <  1:998147/150‑1 (MQ=255)
                                                                                    cTTTCTGCTCCAGCATACGTTTCGCCAAATCCACAATGACGGCTGCGGCTTCAACCGGCGGCGTGCCGCCCTGGTGAATGTTAGAGATACAGGTGCGATCGGCCTCAACGGTGGTCGCCATACGCGGCGAGTAAACGGCGTAGCAGgaga                                                                  <  1:979356/150‑1 (MQ=255)
                                                                                    cTTTCTGCTCCAGCATACGTTTCGCCAAATCCACAATGACGGCTGCGGCTTCAACCGGCGGCGTGCCGCCCTGGTGAATGTTAGAGATACAGGTGCGATCGGCCTCAACGGTGGTCGCCATACGCGGCGAGTAAACGGCGTAGCAGgaga                                                                  <  1:828813/150‑1 (MQ=255)
                                                                                    cTTTCTGCTCCAGCATACGTTTCGCCAAATCCACAATGACGGCTGCGGCTTCAACCGGCGGCGTGCCGCCCTGGTGAATGTTAGAGATACAGGTGCGATCGGCCTCAACGGTGGTCGCCATACGCGGCGAGTAAACGGCGTAGCAGgaga                                                                  >  1:362296/1‑150 (MQ=255)
                                                                                       tCTGCTCCAGCATACGTTTCGCCAAATCCACAATTACGGCTGAGGCTTCAAACGGCGGCGTGAAGCACTGGTGACTGTTAGAGATACAGGTGCGATCTGCCTCAAAGGTGGTAGCCATACGAGGAGAGTAAACGGCGTAGCAGGAGAGGc                                                               >  1:1098485/1‑150 (MQ=255)
                                                                                       tCTGCTCCAGCATACGTTTCGCCAAATCCACAATGACGGCTGCGGCTTCAACCGGCGGCGTGCCGCCCTGGTGAATGTTAGAGATACAGGTGCGATCGGCCTCAACGGTGGTCGCCATACGCGGCGAGTAAACGGCGTAGCAGGAGTGGc                                                               <  1:1258961/150‑1 (MQ=255)
                                                                                          gCTCCAGCATACGTTTCGCCAAATCCACAATGACGGCTGCGGCTTCAACCGGCGGCGTGCCGCCCTGGTGAATGTTAGAGATACAGGTGCGATCGGCCTCAACGGTGGTCGCCATACGCGGCGAGTAAACGGCGTAGCAGGAGAGGCttt                                                            <  2:998844/150‑1 (MQ=255)
                                                                                          gCTCCAGCATACGTTTCGCCAAATCCACAATGACGGCTGCGGCTTCAACCGGCGGCGTGCCGCCCTGGTGAATGTTAGAGATACAGGTGCGATCGGCCTCAACGGTGGTCGCCATACGCGGCGAGTAAACGGCGTAGCAGGAGAGGCttt                                                            <  2:925219/150‑1 (MQ=255)
                                                                                           cTCCAGCCTACGTTTCGCCAAATCCACAATGACGGCTGCGGCTTCAACCGGCGGCGTGCCGCCCTGGTGAATGTTAGAGATACAGGTGCGATCGGCCTCAACGGTGGTCGCCATACGCGGCGAGTAAACGGCGTAGCAGGAGAGGCTTTc                                                           <  1:474770/150‑1 (MQ=255)
                                                                                           cTCCAGCATACGTTTCGCCAAATCCACAATGACGGCTGCGGCTTCAACCGGCGGCGTGCCGCCCTGGTGAATGTTAGAGATACAGGTGCGATCGGCCTCAACGGTGGTCGCCATACGCGGCGAGTAAACGGCGTAGCAGGAGAGGCTTTc                                                           >  1:1123558/1‑150 (MQ=255)
                                                                                           cTCCAGCATACGTTTCGCCAAATCCACAATGACGGCTGCGGCTTCAACCGGCGGCGTGCCGCCCTGGTGAATGTTAGAGATACAGGTGCGATCGGCCTCAACGGTGGTCGCCATACGCGGCGAGTAAACGGCGTAGCAGGAGAGGCTTTc                                                           >  1:1123581/1‑150 (MQ=255)
                                                                                           cTCCAGCATACGTTTCGCCAAATCCACAAAGACGGCTGCGTCTTCAACCGGCGGCGTGCCGCACTGGTGAATGTTAGAGATACAGGTGCGATCGGCCTCAACGGTGGTCGCCATACGCGGCGAGTAAACGGCGTAGCAGGAGAGGCTTTc                                                           >  1:1123574/1‑150 (MQ=255)
                                                                                               aGCATACGTTTCGCCAAATCCACAATGACGGCTGCGGCTTCAACCGGCGGCGTGCCGCCCTGGTGAATGTTAGAGATACAGGTGCGATCGGCCTCAACGGTGGTCGCCATACGCGGCGAGTAAACGGCGTAGCAGGAGAGGCTTTctgac                                                       <  1:801266/150‑1 (MQ=255)
                                                                                                gCATACGTTTCGCCAAATCCACAATGACGGCTGCGGCTTCAACCGGCGGCGTGCCGCCCTGGTGAATGTTAGAGATACAGGTGCGATCGGCCTCAACGGTGGTCGCCATACGCGGCGAGTAAACGGCGTAGCAGGAGAGGCTTTctgact                                                      <  1:518089/150‑1 (MQ=255)
                                                                                                 cATACGTTTCGCCAAATCCACAATGACGGCTGCGGCTTCAACCGGCGGCGTGCCGCCCTGGTGAATGTTAGAGATACAGGTGCGATCGGCCTCAACGGTGGTCGCCATACGCGGCGAGTAAACGGCGTAGCAGGAGAGGCTTTctgactg                                                     >  2:1107445/1‑150 (MQ=255)
                                                                                                           gCCAAATCCACAATGACGGCTGCGGCTTCAACCGGCGGCGTGCCGCCCTGGTGAATGTTAGAGATACAGGTGCGATCGGCCTCAACGGTGGTCGCCATACGCGGCGAGTAAACGGCGTAGCAGGAGAGGCTTTCTGACTGACCCAAcccc                                           <  2:127781/150‑1 (MQ=255)
                                                                                                           gCCAAATCCAAAATGACGGCTTCGGCTTCAACCGGCGGCGTGCCGCCCTGGTGAATGTTAGAGATACAGGTGCGATCGGCCTCAACGGTGGTCGCCATACGCGGCGAGTAAACGGCGTAGCAGGAGAGGCTTTCTGACTGACCCAAcccc                                           <  2:837720/150‑1 (MQ=255)
                                                                                                            ccAAATCCACAATGACGGCTGCGGCTTCAACCGGCGGCGTGCCGCCCTGGTGAATGTTAGAGATACAGGTGCGATCGGCCTCAACGGTGGTCGCCATACGCGGCGAGTAAACGGCGTAGCAGGAGAGGCTTTCTGACTGACCCAACCCCg                                          >  2:1260049/1‑150 (MQ=255)
                                                                                                            ccAAATCCACAATGACGGCTGCGGCTTCAACCGGCGGCGTGCCGCCCTGGTGAATGTTAGAGATACAGGTGCGATCGGCCTCAACGGTGGTCGCCATACGCGGCGAGTAAACGGCGTAGCAGGAGAGGCTTTCTGACTGACCCAACCCCg                                          >  2:360099/1‑150 (MQ=255)
                                                                                                                   cacaATGACGGCTGCGGCTTCAACCGGCGGCGTGCCGCCCTGGTGAATGTTAGAGATACAGGTGCGATCGGCCTCAACGGTGGTCGCCATACGCGGCGAGTAAACGGCGTAGCAGGAGAGGCTTTCTGACTGACCCAACCCCGGACGTTc                                   <  2:986916/150‑1 (MQ=255)
                                                                                                                                  ggCTTCAACCGGCGGCGTGCCGCCCTGGTGAATGTTAGAGATACAGGTGCGATCGGCCTCAACGGTGGTCGCCATACGCGGCGAGTAAACGGCGTAGCAGGAGAGGCTTTCTGACTGACCCAACCCCGGACGTTCGCCCACCAGCAGGAt                    <  1:25642/150‑1 (MQ=255)
                                                                                                                                  ggCTTCAACCGGCGGCGTGCCGCCCTGGTGAATGTTAGAGATACAGGTGCGATCGGCCTCAACGGTGGTCGCCATACGCGGCGAGTAAACGGCGTAGCAGGAGAGGCTTTCTGACTGACCCAACCCCGGACGTTCGCCCACCAGCAGGAt                    <  1:735327/150‑1 (MQ=255)
                                                                                                                                    cTTCAACCGGCGGCGTGCCGCCCTGGTGAATGTTAGAGACACAGGTGCGATCGGCCTCAACGGTGGTCGCCATACGCGGCGAGTAAACGGCGTAGCAGGAGAGGCTTTCTGACTGACCCAACCCCGGACGTTCGCCCACCAGCAGGATCa                  <  1:570658/150‑1 (MQ=255)
                                                                                                                                     ttCAACCGGCGGCGTGCCGCCCTGGTGAATGTTAGAGATACAGGTGCGATCGGCCTCAACGGTGGTCGCCATACGCGGCGAGTAAACGGCGTAGCAGGAGAGGCTTTCTGACTGACCCAACCCCGGACGTTCGCCCACCAGCAGGATCac                 >  2:378751/1‑150 (MQ=255)
                                                                                                                                      tCAACCGGCGGCGTGTCGCCCTGGTGAATGTTAGAGATACAGGTGCGATCGGCCTCAACGGTGGTCGCCATACGCGGCGAGTAAACGGCGTAGCAGGAGAGGCTTTCTGACTGACCCAACCCCGGACGTTCGCCCACCAGCAGGATCact                <  1:726400/150‑1 (MQ=255)
                                                                                                                                      tCAACCGGCGGCGTGTCGCCCTGGTGAATGTTAGAGATACAGGTGCGATCGGCCTCAACGGTGGTCGCCATACGCGGCGAGTAAACGGCGTAGCAGGAGAGGCTTTCTGACTGACCCAACCCCGGACGTTCGCCCACCAGCAGGATCact                <  1:726419/150‑1 (MQ=255)
                                                                                                                                            ggcggcGTGCCGCCCTGGTGAATGTTAGAGATACAGGTGCGATCGGCCTCAACGGTGGTCGCCATACGCGGCGAGTAAACGGCGTAGCAGGAGAGGCTTTCTGACTGACCCAACCCCGGACGTTCGCCCACCAGCAGGATCACTACTTTc          <  1:799991/150‑1 (MQ=255)
                                                                                                                                                gcgTGCCGCCCTGGTGAATGTTAGAGATACAGGTGCGATCGGCCTCAACGGTGGTCGCCATACGCGGCGa                                                                                      <  2:669592/70‑1 (MQ=255)
                                                                                                                                                gcgTGCCGCCCTGGTGAATGTTAGAGATACAGGTGCGATCGGCCTCAACGGTGGTCGCCATACGCGGCGa                                                                                      >  1:669592/1‑70 (MQ=255)
                                                                                                                                                   tGCCGCCCTGGTGTATGTTAGAGATACAGGTGCGATCGGCCTCAACGGTGGTCGCCATACGCGGCGAGTAAACGGCGTAGCAGGAGAGGCTTTCTGACTGACCCAACCCCGGACGTTCGCCCACCAGCAGGATCACTACTTTCGCGccta   <  1:935259/150‑3 (MQ=255)
                                                                                                                                                   tGCCGCCCTGGTGAATGTTAGAGATACAGGTGCGATCGGCCTCAACGGTGGTCGCCATACGCGGCGAGTAAACGGCGTAGCAGGAGAGGCTTTCTGACTGACCCAACCCCGGACGTTCGCCCACCAGCAGGATCACTACTTTCGCGCCGa   <  1:538149/150‑1 (MQ=255)
                                                                                                                                                    gccgccCTGGTGAATGTTAGAGATACAGGTGCGATCGGCCTCAACGGTGGTCGCCATACGCGGCGAGTAAACGGCGTAGCAGGAGAGGCTTTCTGACTGACCCAACCCCGGACGTTCGCCCACCAGCAGGATCACTACTTTCGCGCCGAg  >  2:797054/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TCACCCGCATAGCGGTGACAGACGGTCGAATCAAATCTAAAGCTGGCATGATGCCTCCTTATCGGGTCATGTTGATGCCGGATGCTTTCTGCTCCAGCATACGTTTCGCCAAATCCACAATGACGGCTGCGGCTTCAACCGGCGGCGTGCCGCCCTGGTGAATGTTAGAGATACAGGTGCGATCGGCCTCAACGGTGGTCGCCATACGCGGCGAGTAAACGGCGTAGCAGGAGAGGCTTTCTGACTGACCCAACCCCGGACGTTCGCCCACCAGCAGGATCACTACTTTCGCGCCGAG  >  NC_000913/2556352‑2556649

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: