Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 950,646 A→G 6.1% V145V (GTT→GTC) ‡ pflA ← pyruvate formate‑lyase 1‑activating enzyme; [formate‑C‑acetyltransferase 1]‑activating enzyme; PFL activase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913950,6460AG 6.1% 335.3 / 9.8 117V145V (GTT→GTC) ‡pflApyruvate formate‑lyase 1‑activating enzyme; [formate‑C‑acetyltransferase 1]‑activating enzyme; PFL activase
Reads supporting (aligned to +/- strand):  ref base A (56/52);  new base G (2/5);  total (59/58)
Fisher's exact test for biased strand distribution p-value = 2.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.50e-01

Only 100 of 119 total aligned reads displayed.

TCAACGTTGCCCATATCACGGGTAAATTCACCGAGGCGATGCGCTGAATCGTCATCGTCAGACCAGCCTGGGACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGTTCATCAATCACCGGATCGTAACGACGAACAAAACCGTTGGTGTCCAGACAGGTATGAATGCCTTCTTTT  >  NC_000913/950497‑950793
                                                                                                                                                     |                                                                                                                                                   
tcAACGTTGCCCATATCACGGGTAAATTCACCGAGGCGATGCGCTGAATCGTCATCGTCAGACCAGCCTGGGACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGaa                                                                                                                                                     >  2:849602/1‑150 (MQ=255)
 cAACGTTGCCCATATCACGGGTAAATTCACCGAGGCGATGCGCTGAATCGTCATCGTCAGACCAGCCTGGGACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGaaa                                                                                                                                                    >  2:1529887/1‑150 (MQ=255)
   aCGTAGCCCATATCCCGGGTAAATTCACCGAGGCGATGCGCTGAATCGTCATCGTCAGACCAGCCTGGGACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACt                                                                                                                                                  <  1:1512658/150‑1 (MQ=255)
      ttGCCCATATCACGGGTAAATTCACCGAGGCGATGCGCTGAATCGTCATCGTCAGACCAGCCTGGGACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCa                                                                                                                                               <  2:431553/150‑1 (MQ=255)
      ttGCCCATATCACGGGTAAATTCACCGAGGCGATGCGCTGAATCGTCATCGTCAGACCAGCCTGGGACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCa                                                                                                                                               <  1:1378943/150‑1 (MQ=255)
         cccATATCACGGGTAAATTCACCGAGGCGATGCGCTGAATCGTCATCGTCAGACCAGCCTGGGACAACAACGTAGCGGATCCACCCCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAAcc                                                                                                                                            >  1:57916/1‑150 (MQ=255)
         cccATATCACGGGTAAATTCACCGAGGCGATGCGCTGAATCGTCATCGTCAGACCAGCCTGGGACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAGTTTCCAAcc                                                                                                                                            >  1:1043188/1‑150 (MQ=255)
         cccATATCACGGGTAAATTCACCGAGGCGATGCGCTGAATCGTCATCGTCAGACCAGCCTGGGACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAAcc                                                                                                                                            >  1:945206/1‑150 (MQ=255)
          ccATATCACGGGTAAATTCACCGAGGCGATGCGCTGAATCGTCATCGTCAGACCAGCCTGGGACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCa                                                                                                                                           >  1:788598/1‑150 (MQ=255)
          ccATATCACGGGTAAATTCACCGAGGCGATGCGCTGAATCGTCATCGTCAGACCAGCCTGGGACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCa                                                                                                                                           >  2:713062/1‑150 (MQ=255)
           cATATCACGGGTAAATTCACCGAGGCGATGCGCTGAATCGTCATCGTCAGACCAGCCTGGGACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAGTTTCCAACCac                                                                                                                                          <  2:1043188/150‑2 (MQ=255)
           cATATCACGGGTAAATTCACCGAGGCGATGCGCTGAATCGTCATCGTCAGACCAGCCTGGGACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAGTTTCCAACCac                                                                                                                                          <  1:698567/150‑2 (MQ=255)
            atatCACGGGTAAATTCACCGAGGCGATGCGCTGAATCGTCATCGTCAGACCAGCCTGGGACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAGTTTCCAACCacc                                                                                                                                         <  1:1128787/150‑3 (MQ=255)
            atatCACGGGTAAATTCACCGAGGCGATGCGCTGAATCGTCATCGTCAGACCAGCCTGGGACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGa                                                                                                                                         >  2:1249697/1‑150 (MQ=255)
                cACGGGTAAATTCACCGAGGCGATGCGCTGAATCGTCATCGTCAGACCTGCCTGGGACAACAACGTAGCGGATCCAAAACTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGAtttt                                                                                                                                     <  1:847701/150‑1 (MQ=255)
                cACGGGTAAATTCACCGAGGCGATGCGCTGAATCGTCATCGTCAGACCAGCCTGGGACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGAtttt                                                                                                                                     <  1:1390512/150‑1 (MQ=255)
                   gggTAAATTCACCGAGGCGATGCGCTGAATCGTCATCGTCAGACCAGCCTGGGACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTtggt                                                                                                                                  >  2:700657/1‑150 (MQ=255)
                   gggTAAATTCACCGAGGCGATGCGCTGAATCGTCATCGTCAGACCAGCCTGGGACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTtggt                                                                                                                                  >  2:1677454/1‑150 (MQ=255)
                        aaTTCACCGAGGCGATGCGCTGAATCGTCATCGTCAGACCAGCCTGGGACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGAtc                                                                                                                             >  1:385313/1‑150 (MQ=255)
                        aaTTCACCGAGGCGATGCGCTGAATCGTCATCGTCAGACCAGCCTGGGACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGAtc                                                                                                                             >  1:391659/1‑150 (MQ=255)
                        aaTTCACCGAGGCGATGCGCTGAATCGTCATCGTCAGACCAGCCTGGGACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGAtc                                                                                                                             >  2:641889/1‑150 (MQ=255)
                         aTTCACCGAGGCGATGCGCTGAATCGTCATCGTCAGACCAGCCTGGGACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGAtct                                                                                                                            >  2:529251/1‑150 (MQ=255)
                            cACCGCGGCGATGCGCTGAATCGTCATCGTCAGACCAGCCTGGGACAACTCCGTAGCGGATCCACTCCTTCACATTTTTGTTCGCCAGATATTTAGCGATCGCCAGCGTGCGGTGGTTGGATACTCCAACCAGATTTTGGTGGATCtcgt                                                                                                                         <  1:234285/150‑1 (MQ=255)
                            cACCGAGGCGATGCGCTGAATCGTCATCGTCAGACCAGCCTGGGACAACAACGTAGCGGATCCACCCCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCtcgt                                                                                                                         >  1:1505704/1‑150 (MQ=255)
                            cACCGAGGCGATGCGCTGAATCGTCATCGTCAGACCAGCCTGGGACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTTGCGACCTCCAGCGAGCGGTGGTTGGAAACTACCAGCAGATTTTGGTGGATCtcgt                                                                                                                         >  1:126371/1‑150 (MQ=255)
                            cACCGAGGCGATGCGCTGAATCGTCATCGTCAGACCAGCCTGGGACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCtcgt                                                                                                                         <  1:1001377/150‑1 (MQ=255)
                            cACCGAGGCGATGCGCGGAATCGTCATCGTCAGACCAGCCTGGGACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCtcgt                                                                                                                         >  1:1699246/1‑150 (MQ=255)
                                  ggCGATGCGCTGAATCGTCATCGTCAGACCAGCCTGGGACAACAACGTAGCGGATCCACACCTTCACATTTTAGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCa                                                                                                                   <  1:1249697/150‑1 (MQ=255)
                                    cGATGCGCTGAATCGTCATCGTCAGACCAGCCTGGGACAACAACGTAGCGGCTCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCAtg                                                                                                                 >  2:1631797/1‑149 (MQ=255)
                                     gATGCGCTGAATCGTCATCGTCAGACCAGCCTGGGACAACAAAGTAGGGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCt                                                                                                                <  1:147496/150‑1 (MQ=255)
                                        gcgcTGAATCGTCATCGTCAGACCAGCCTGGGACTACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGGGGATCTCGTCGTTCATCTGtt                                                                                                             <  1:301014/150‑1 (MQ=255)
                                               atcgtcatcgtcAGACCAGCCTGTGACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATc                                                                                                      <  1:849602/150‑1 (MQ=255)
                                                 cgtcatcgtcaGACCAGCCTGGGACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGa                                                                                                    >  1:551988/1‑150 (MQ=255)
                                                    catcgtcaGACCAGCCTGGGACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCa                                                                                                 <  1:1371182/150‑1 (MQ=255)
                                                       cgtcaGACCAGCCTGGGACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTa                                                                                              >  2:1024174/1‑150 (MQ=255)
                                                        gtcaGACCAGCCTGGGACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTAc                                                                                             <  1:1613026/150‑1 (MQ=255)
                                                         tcaGACCAGCCTGGGACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTAcc                                                                                            >  1:1477798/1‑150 (MQ=255)
                                                         tcaGACCAGCCTGGGACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTAcc                                                                                            <  1:641889/150‑1 (MQ=255)
                                                           aGACCAGCCTGGGACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAg                                                                                          <  1:529251/150‑1 (MQ=255)
                                                              ccAGCCTGGGACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTc                                                                                       >  1:1696172/1‑150 (MQ=255)
                                                              ccAGCCTGGGACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTc                                                                                       <  1:713062/150‑1 (MQ=255)
                                                              ccAGCCTGGGACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTc                                                                                       <  2:461496/150‑1 (MQ=255)
                                                                 gCCTGGGACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGt                                                                                    <  1:1024174/150‑1 (MQ=255)
                                                                 gCCTGGGACAACAACGTAGCGGAGCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGt                                                                                    <  1:927326/150‑1 (MQ=255)
                                                                   cTGGGACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTa                                                                                  >  2:12180/1‑150 (MQ=255)
                                                                     gggACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACt                                                                                <  2:788598/150‑1 (MQ=255)
                                                                     gggACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACt                                                                                <  1:1653758/150‑1 (MQ=255)
                                                                     gggACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATATGTTTGAGATCGAGCATTACCAGGTCGCTTACt                                                                                >  2:400552/1‑150 (MQ=255)
                                                                      gaacaacaACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACtt                                                                               >  2:871292/3‑150 (MQ=255)
                                                                         caacaaCGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGCTACTTCca                                                                            <  1:605368/150‑1 (MQ=255)
                                                                         caacaaCGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATAGAGCATTACCAGGTCGGTTACTTCca                                                                            >  2:337258/1‑150 (MQ=255)
                                                                          aacaacGTAGCGGATCCACACCTTCATATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCcag                                                                           >  2:1680462/1‑150 (MQ=255)
                                                                          aacaacGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCcag                                                                           <  2:1542530/150‑1 (MQ=255)
                                                                             aacGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGAATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCcagcag                                                                        <  2:1426885/150‑1 (MQ=255)
                                                                               cGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGtt                                                                      >  2:446867/1‑150 (MQ=255)
                                                                               cGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGtt                                                                      >  1:1424230/1‑150 (MQ=255)
                                                                               cGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGACGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGtt                                                                      >  2:612034/1‑150 (MQ=255)
                                                                                  aGCGGATCGACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGTtcat                                                                   <  1:340877/150‑1 (MQ=255)
                                                                                  aGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGATCATTACCAGGTCGGTTACTTCCAGCAGTtcat                                                                   <  1:831896/150‑1 (MQ=255)
                                                                                  aGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGTtcat                                                                   <  1:1677454/150‑1 (MQ=255)
                                                                                   gCGGCTCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGTtcatc                                                                  <  1:644560/150‑1 (MQ=255)
                                                                                   gCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGTtcatc                                                                  <  2:1659306/150‑1 (MQ=255)
                                                                                   gCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGTtcatc                                                                  <  2:1477798/150‑1 (MQ=255)
                                                                                      gATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCAACTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGTTAatcaat                                                               >  2:1134101/1‑150 (MQ=255)
                                                                                        tCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATAGAGCATTACCAGGTCGGTTACTTCCAGCAGTTCatcaatca                                                             >  1:1372002/1‑150 (MQ=255)
                                                                                          gacacCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTGGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGTTCATCAATCAcc                                                           <  1:46646/149‑1 (MQ=255)
                                                                                            cacCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTAGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTAATTCAAGCAGTTCATCAATCACCgg                                                         >  2:877191/1‑150 (MQ=255)
                                                                                              ccTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGTTCATCAATCACCGGAt                                                       >  1:574711/1‑150 (MQ=255)
                                                                                              ccTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGTTCATCAATCACCGGAt                                                       >  1:1210859/1‑150 (MQ=255)
                                                                                               cTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGATCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGTTCATCAATCACCGGATc                                                      <  2:1138065/150‑1 (MQ=255)
                                                                                                    caTTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGTTCATCAATCACCGGATCGTAac                                                 <  1:900879/150‑1 (MQ=255)
                                                                                                     gctggagTTCGCTAAATATCTGGCGAACTCCAGCGGGCGGTGGTTGGAAACTCCCACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTAAAAGGTCGGTTACTTACAGCAGTTCTTCAATCAGCGGATCGTAtcg                                                >  2:1075627/7‑150 (MQ=255)
                                                                                                     aTTTTTGTTCGCCAGATATTTAGCGAACTCCATCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGTTCATCAATCACCGGATCGTAacg                                                <  1:81077/150‑1 (MQ=255)
                                                                                                        tttGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGTTCATCAATCACCGGATCGTAacgacg                                             <  1:700657/150‑1 (MQ=255)
                                                                                                            ttCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGTTCATCAATCACCGGATCGTAACGACGAACa                                         <  2:416276/150‑1 (MQ=255)
                                                                                                               gCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTAATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGTTCATCAATCACCGGATCGTAACGACTAACaaaa                                      >  1:1624552/1‑150 (MQ=255)
                                                                                                                ccAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGTTCATCAATCACCGGATCGTAACGACGAACAAAAc                                     <  1:612034/150‑1 (MQ=255)
                                                                                                                   gATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGTTCATCAATCACCGGATCGTAACGACGAACAAAACCGt                                  <  2:385313/150‑1 (MQ=255)
                                                                                                                   gATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGTTCATCAATCACCGGATCGTAACGACGAACAAAACCGt                                  <  2:391659/150‑1 (MQ=255)
                                                                                                                     tatTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTAACAGGTCGGTTACTTCCAGCAGGTCATCAATCACCGGATCGTAACGACGAACAAAACCGTTg                                >  1:1408032/1‑150 (MQ=255)
                                                                                                                     taattAGCGAACTCCAGCGTGCGGTGCTTGTAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGTTCATCAATAACCGGATCGTAACGACGCACAAAACCGTTg                                >  2:315079/4‑150 (MQ=255)
                                                                                                                         tAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTGGGTGGATCTCGTCGTTCATCTGTGTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGTTCATCAATCACCGGATCGTAACGACGAACAAAACCGTTGgtgt                            >  1:1060984/1‑150 (MQ=255)
                                                                                                                          aGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGTTCATCAATCACCGGATCGTAACGACGAACAAAACCGTTGGTGTc                           >  2:34657/1‑150 (MQ=255)
                                                                                                                          aGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGTTCATCAATCACCGGATCGTAACGACGAACAAAACCGTTGGTGTc                           >  1:290724/1‑150 (MQ=255)
                                                                                                                          aGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGACGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGTTCATCAATCACCGGATCGTAACGACGAACAAAACCGTTGGTGTc                           <  1:1475652/150‑1 (MQ=255)
                                                                                                                           gCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGTTCATCAATCACCGGATCGTAacgacg                                             >  2:1201231/1‑131 (MQ=255)
                                                                                                                           gCGAACTCCAGCGTGCGGTGGTGGGAAACTCCAACCAGATTTTGGTGGATCTCGTGGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGTTCATCAATCACCGGATCGTAACGACGAACAAAACCGTTGGTGTcc                          <  1:871292/150‑1 (MQ=255)
                                                                                                                             gAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGTTCATCAATCACCGGATCGTAACGACGAACAAAACCGTTGGTGTCCAg                        <  1:631934/150‑1 (MQ=255)
                                                                                                                                cTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGTTCATCAATCACCGGATCGTAACGACGAACAAAACCGTTGGTGTCCAGACa                     >  1:539075/1‑150 (MQ=255)
                                                                                                                                 tCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGTTCATCAATCACCGGATCGTAACGACGAACAAAACCGTTGGTGTCCAGACAg                    <  1:1529887/150‑1 (MQ=255)
                                                                                                                                  ccAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGTTCATCAATCACCGGATCGTAACGACGAACAAAACCGTTGGTGTCCAGACAgg                   >  2:362080/1‑150 (MQ=255)
                                                                                                                                  ccAGCGTGCGGTGGTTGGAAAATCCAACCAGATTTTGGTGGATCTCGTCGGTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGTTCATCAATCACCGTATCGTAACGACGAACAAAACCGTTGGTGTCCAGACAgt                   <  2:1064742/150‑2 (MQ=255)
                                                                                                                                      cGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGTTCATCAATCACCGGATCGTAACGACGAACAAAACCGTTGGTGTCCAGACAGGTATg               <  2:57916/150‑1 (MQ=255)
                                                                                                                                      cGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGTTCATCAATCACCGGATCGTAACGACGAACAAAACCGTTGGTGTCCAGACAGGTATg               >  2:376802/1‑150 (MQ=255)
                                                                                                                                        tGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGTTCATCAATCACCGGATCGTAACGACGAACAAAACCGTTGGTGTCCAGACAGGTATGaa             <  2:1424230/150‑1 (MQ=255)
                                                                                                                                          cGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGTTCATCAATCACCGGATCGTAACGACGAACAAAACCGTTGGTGTCCAGACAGGTATGAATg           <  1:1547422/150‑1 (MQ=255)
                                                                                                                                            ctggtTGGAGTTTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGTTCATCAATCACCGGATCGTAACGACGAACAAAACCGTTGGTGTCCAGACAGGTATGAATGcc         >  1:1187399/2‑150 (MQ=255)
                                                                                                                                                ttGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGTTCATCAATCACCGGATCGTAACGACGAACAAAACCGGTGGTGTCCAGACAGGTATGAATGCcttct     >  1:1240030/1‑150 (MQ=255)
                                                                                                                                                  ggAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGTTCATCAATCACCGGATCGTAACGACGAACAAAACCGTTGGTGTCCAGACAGGTATGAATGCCTTCttt   <  1:12180/150‑1 (MQ=255)
                                                                                                                                                   gAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGTTCATCAATCACCGGATCGTTACGACGAACAAAACCGTTGGTGTCCAGACAGGTATGAATGCCTTCtttt  >  1:232049/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
TCAACGTTGCCCATATCACGGGTAAATTCACCGAGGCGATGCGCTGAATCGTCATCGTCAGACCAGCCTGGGACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGTTCATCAATCACCGGATCGTAACGACGAACAAAACCGTTGGTGTCCAGACAGGTATGAATGCCTTCTTTT  >  NC_000913/950497‑950793

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