Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,152,555 C→T 5.1% T107I (ACC→ATC)  metC → cystathionine beta‑lyase/L‑cysteine desulfhydrase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,152,5550CT 5.1% 251.7 / 2.0 79T107I (ACC→ATC) metCcystathionine beta‑lyase/L‑cysteine desulfhydrase
Reads supporting (aligned to +/- strand):  ref base C (39/35);  new base T (2/2);  total (42/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.19e-01

ACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCTGCAGCCAAACACTAAAATCGTGTTTCTGGAATCGCCAG  >  NC_000913/3152406‑3152704
                                                                                                                                                     |                                                                                                                                                     
aCGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGAc                                                                                                                                                       <  2:598570/150‑1 (MQ=255)
 cGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGAcc                                                                                                                                                      <  1:25159/150‑1 (MQ=255)
 cGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGAcc                                                                                                                                                      <  1:1827774/150‑1 (MQ=255)
 cGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGAcc                                                                                                                                                      >  2:1925268/1‑150 (MQ=255)
    cgcgGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAac                                                                                                                                                   <  1:733059/150‑1 (MQ=255)
            gTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGGTCATCAacacatgatc                                                                                                                                           <  1:711456/150‑7 (MQ=255)
            gTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTa                                                                                                                                           <  1:1944804/150‑1 (MQ=255)
             ttAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGGTCATCAacacatgatcg                                                                                                                                          <  2:1818032/150‑8 (MQ=255)
                aCCCATTTCTCCTTACCACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCTTGCTATTTCCCTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACATGGCGATCATGTGTTGATGACCAACACCGCCTATGaa                                                                                                                                       <  2:1998445/150‑1 (MQ=255)
                 cccATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAAc                                                                                                                                      >  1:860755/1‑150 (MQ=255)
                 cccATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAAc                                                                                                                                      >  2:678751/1‑150 (MQ=255)
                  ccATTTCTCCTTACAACAAGCGATGGGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTGGCTTTTCTCGTACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAAcc                                                                                                                                     <  2:753148/150‑1 (MQ=255)
                        ctcCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCa                                                                                                                               >  1:1625155/1‑150 (MQ=255)
                         tcCTTACAACAAGCGATGTGTGAACTGGAAGGTGTCGCAGGCTGCGTGCTATTTCCCTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAg                                                                                                                              <  1:789843/150‑1 (MQ=255)
                         tcCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAg                                                                                                                              <  2:1657949/150‑1 (MQ=255)
                          ccTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAgg                                                                                                                             >  1:757697/1‑150 (MQ=255)
                          ccTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAgg                                                                                                                             >  1:1146220/1‑150 (MQ=255)
                           cTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGa                                                                                                                            >  2:330023/1‑150 (MQ=255)
                               caacaaGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTc                                                                                                                        >  2:1728592/1‑150 (MQ=255)
                                 acaaGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTg                                                                                                                      <  1:739623/150‑1 (MQ=255)
                                  caaGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGt                                                                                                                     <  1:208227/150‑1 (MQ=255)
                                   aaGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTa                                                                                                                    >  2:1071226/1‑150 (MQ=255)
                                   aaGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTa                                                                                                                    >  1:1662780/1‑150 (MQ=255)
                                   aaGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTa                                                                                                                    >  1:1655809/1‑150 (MQ=255)
                                    cgCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAg                                                                                                                   <  1:669968/149‑1 (MQ=255)
                                     gCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGc                                                                                                                  <  2:1625155/150‑1 (MQ=255)
                                             tgAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCt                                                                                                          >  1:657927/1‑150 (MQ=255)
                                              gAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTc                                                                                                         <  2:860755/150‑1 (MQ=255)
                                                aCTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAg                                                                                                       <  1:156628/150‑1 (MQ=255)
                                                aCTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAg                                                                                                       <  1:149574/150‑1 (MQ=255)
                                                 cTGGAAGGTGGAGCAGGCTGCGTGCTATTTCCCTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGc                                                                                                      >  2:870647/1‑150 (MQ=255)
                                                      aGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACACGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACt                                                                                                 <  2:986398/150‑1 (MQ=255)
                                                         tGGCGCAGGCTGCGTGCTATTTCCCTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATcc                                                                                                           >  2:896336/1‑137 (MQ=255)
                                                         tGGCGCAGGCTGCGTGCTATTTCCCTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATcc                                                                                                           <  1:896336/137‑1 (MQ=255)
                                                              cAGGCTGCGTGCTATTTCCCTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTaa                                                                                         <  2:657927/150‑1 (MQ=255)
                                                               aGGCTGCGTGCTATTTCCCTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAAc                                                                                        >  1:1816930/1‑150 (MQ=255)
                                                                ggCTGCGTGCTATTTCCCTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACg                                                                                       >  1:828075/1‑150 (MQ=255)
                                                                  cTGCGTGCTATTTCCCTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTa                                                                                          <  2:1019322/145‑1 (MQ=255)
                                                                  cTGCGTGCTATTTCCCTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTa                                                                                          >  1:1019322/1‑145 (MQ=255)
                                                                          tATTTCCCTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGt                                                                             <  2:1146220/150‑1 (MQ=255)
                                                                           aTTTCCCTGGGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAATATCCTCAGCAAACTGGCCGTAACGACATCATGGtt                                                                            <  2:876621/150‑1 (MQ=255)
                                                                           aTTTCCCTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGtt                                                                            >  2:311499/1‑150 (MQ=255)
                                                                              tCCCTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGa                                                                         >  1:1592409/1‑150 (MQ=255)
                                                                               cccTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGAt                                                                        >  2:993822/1‑150 (MQ=255)
                                                                                ccTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATc                                                                       <  2:959354/150‑1 (MQ=255)
                                                                                 cTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATcc                                                                      >  2:1951886/1‑150 (MQ=255)
                                                                                       ggcggcAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGAt                                                                <  1:1330014/150‑1 (MQ=255)
                                                                                              gCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGcc                                                         <  1:330023/150‑1 (MQ=255)
                                                                                               cGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCg                                                        >  2:1512961/1‑150 (MQ=255)
                                                                                                ggTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGa                                                       <  1:311499/150‑1 (MQ=255)
                                                                                                      tAATTCCATTCTTTCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGt                                                 >  2:1562258/1‑150 (MQ=255)
                                                                                                      tAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGt                                                 >  2:1010161/1‑150 (MQ=255)
                                                                                                       aaTTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGtt                                                <  1:870647/150‑1 (MQ=255)
                                                                                                        aTTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAATCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTa                                               <  2:355056/150‑1 (MQ=255)
                                                                                                          tCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAg                                             >  2:673545/1‑150 (MQ=255)
                                                                                                           ccATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGc                                            >  1:1195413/1‑150 (MQ=255)
                                                                                                            cATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCa                                           >  2:1565815/1‑150 (MQ=255)
                                                                                                                cTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCTg                                       >  2:1517536/1‑150 (MQ=255)
                                                                                                                attGCTTTTTTCTTACAGGTCTATCATGTGTTGATGACCAACACCGCCTATGAACCGATTCAGGATTTCTGTAGCAAAATGCTCAGCAAAGTTTGCGTAACGACATCATTGTTTGATCCGCTGATTGGTGGCGATATCGTTATGGATCTg                                       <  2:1592409/149‑1 (MQ=255)
                                                                                                                   gCTTTTATCGAACACGGCGATCATGTGTTTATGAACAACACCGCCTATGACCCGAGTCAAAATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCTACAg                                    >  2:900612/1‑150 (MQ=255)
                                                                                                                    cTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCTGCAGc                                   >  2:271626/1‑150 (MQ=255)
                                                                                                                      tttATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCTGCAGCCa                                 <  1:1010161/150‑1 (MQ=255)
                                                                                                                         aTCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCTGCAGCCAAac                              >  1:807386/1‑150 (MQ=255)
                                                                                                                         aTCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCTGCAGCCAAac                              <  1:1071226/150‑1 (MQ=255)
                                                                                                                            gAACAGGGCGATCATGTGTTGATTACCAACACAGCCTATGAACCGAGTCAGGATTTCTGTAGCAAACTCCTCACCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGCTGCCGATATCGTTAAGCAACTGCAGCCAAACACt                           >  2:1243881/1‑150 (MQ=255)
                                                                                                                            gAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCTGCAGCCAAACACt                           >  1:1971107/1‑150 (MQ=255)
                                                                                                                            gAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCTGCAGACAAACACt                           >  2:1358096/1‑150 (MQ=255)
                                                                                                                              aCAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCTGCAGCCAAACACTaa                         >  1:905996/1‑150 (MQ=255)
                                                                                                                               cAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCTGCAGCCAAACACTaaa                        >  2:175847/1‑150 (MQ=255)
                                                                                                                                 gggCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCTGCAGCCAAACACTAAAAt                      >  2:1115322/1‑150 (MQ=255)
                                                                                                                                    cGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCTGCAGCCAAACACTAAAATCgt                   <  1:813747/150‑1 (MQ=255)
                                                                                                                                    ataggcggtgtTGGTCATCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCTGCAGCCAAACACTAAAATCgt                   >  2:774550/8‑150 (MQ=255)
                                                                                                                                     gATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCTGCAGCCAAACACTAAAATCgtg                  <  1:673545/150‑1 (MQ=255)
                                                                                                                                      aggcggtgtTGGTCATCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCTGCAGCCAAACACTAAAATCgtgt                 >  2:1653891/6‑150 (MQ=255)
                                                                                                                                          tgtgTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCTGCAGCCAAACACTAAAATCGTGTTTCt             <  2:1816930/150‑1 (MQ=255)
                                                                                                                                            tgtTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCTGCAGCCAAACACTAAAATCGTGTTTCTgg           <  1:415337/150‑1 (MQ=255)
                                                                                                                                            tgtTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCTGCAGCCAAACACTAAAATCGTGTTTCTgg           <  1:1728592/150‑1 (MQ=255)
                                                                                                                                             gtTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGCCATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCTGCAGCCAAACACTAAAATCGTGTTTCTGGa          >  2:1686552/1‑150 (MQ=255)
                                                                                                                                             gtTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCTGCAGCCAAACACTAAAATCGTGTTTCTGGa          >  2:614106/1‑150 (MQ=255)
                                                                                                                                             gtTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCTGCAGCCAAACACTAAAATCGTGTTTCTGGa          <  2:1195413/150‑1 (MQ=255)
                                                                                                                                              ttGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCTGCAGCCAAACACTAAAATCGTGTTTCTGGaa         >  2:1859587/1‑150 (MQ=255)
                                                                                                                                                gatgaCCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCTGCAGCCAAACACTAAAATCGTGTTTCTGGAATc       <  1:175847/150‑1 (MQ=255)
                                                                                                                                                 atgaCCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCTGCAGCCAAACACTAAAATCGTGTTTCTGGAATCg      >  2:822644/1‑150 (MQ=255)
                                                                                                                                                     ccAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCTGCAGCCAAACACTAAAATCGTGTTTCTGGAATCGCCAg  >  1:1289217/1‑150 (MQ=255)
                                                                                                                                                     ccAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCTGCAGCCAAACACTAAAATCGTGTTTCTGGAATCGCCAg  <  2:1198955/150‑1 (MQ=255)
                                                                                                                                                     ccAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCTGCAGCCAAACACTAAAATCGTGTTTCTGGAATCGCCAg  >  2:434929/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ACGGCGCGGAACGTTAACCCATTTCTCCTTACAACAAGCGATGTGTGAACTGGAAGGTGGCGCAGGCTGCGTGCTATTTCCCTGCGGGGCGGCAGCGGTTGCTAATTCCATTCTTGCTTTTATCGAACAGGGCGATCATGTGTTGATGACCAACACCGCCTATGAACCGAGTCAGGATTTCTGTAGCAAAATCCTCAGCAAACTGGGCGTAACGACATCATGGTTTGATCCGCTGATTGGTGCCGATATCGTTAAGCATCTGCAGCCAAACACTAAAATCGTGTTTCTGGAATCGCCAG  >  NC_000913/3152406‑3152704

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