Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 950,647 A→T 5.9% V145D (GTT→GAT) ‡ 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,6470AT 5.9% 345.5 / 7.5 117V145D (GTT→GAT) ‡pflApyruvate formate‑lyase 1‑activating enzyme; [formate‑C‑acetyltransferase 1]‑activating enzyme; PFL activase
Reads supporting (aligned to +/- strand):  ref base A (57/53);  new base T (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.49e-01

Only 100 of 119 total aligned reads displayed.

CAACGTTGCCCATATCACGGGTAAATTCACCGAGGCGATGCGCTGAATCGTCATCGTCAGACCAGCCTGGGACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGTTCATCAATCACCGGATCGTAACGACGAACAAAACCGTTGGTGTCCAGACAGGTATGAATGCCTTCTTTTTTG  >  NC_000913/950498‑950796
                                                                                                                                                     |                                                                                                                                                     
cAACGTTGCCCATATCACGGGTAAATTCACCGAGGCGATGCGCTGAATCGTCATCGTCAGACCAGCCTGGGACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGaaa                                                                                                                                                       >  2:1529887/1‑150 (MQ=255)
  aCGTAGCCCATATCCCGGGTAAATTCACCGAGGCGATGCGCTGAATCGTCATCGTCAGACCAGCCTGGGACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACt                                                                                                                                                     <  1:1512658/150‑1 (MQ=255)
     ttGCCCATATCACGGGTAAATTCACCGAGGCGATGCGCTGAATCGTCATCGTCAGACCAGCCTGGGACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCa                                                                                                                                                  <  1:1378943/150‑1 (MQ=255)
     ttGCCCATATCACGGGTAAATTCACCGAGGCGATGCGCTGAATCGTCATCGTCAGACCAGCCTGGGACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCa                                                                                                                                                  <  2:431553/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:498128/150‑2 (MQ=255)
          cATATCACGGGTAAATTCACCGAGGCGATGCGCTGAATCGTCATCGTCAGACCAGCCTGGGACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAGTTTCCAACCac                                                                                                                                             <  1:698567/150‑2 (MQ=255)
          cATATCACGGGTAAATTCACCGAGGCGATGCGCTGAATCGTCATCGTCAGACCAGCCTGGGACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAGTTTCCAACCac                                                                                                                                             <  2:1043188/150‑2 (MQ=255)
           atatCACGGGTAAATTCACCGAGGCGATGCGCTGAATCGTCATCGTCAGACCAGCCTGGGACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAGTTTCCAACCacc                                                                                                                                            <  1:1128787/150‑3 (MQ=255)
               cACGGGTAAATTCACCGAGGCGATGCGCTGAATCGTCATCGTCAGACCTGCCTGGGACAACAACGTAGCGGATCCAAAACTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGAtttt                                                                                                                                        <  1:847701/150‑1 (MQ=255)
               cACGGGTAAATTCACCGAGGCGATGCGCTGAATCGTCATCGTCAGACCAGCCTGGGACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGAtttt                                                                                                                                        <  1:1390512/150‑1 (MQ=255)
                  gggTAAATTCACCGAGGCGATGCGCTGAATCGTCATCGTCAGACCAGCCTGGGACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTtggt                                                                                                                                     >  2:1677454/1‑150 (MQ=255)
                       aaTTCACCGAGGCGATGCGCTGAATCGTCATCGTCAGACCAGCCTGGGACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGAtc                                                                                                                                >  2:641889/1‑150 (MQ=255)
                       aaTTCACCGAGGCGATGCGCTGAATCGTCATCGTCAGACCAGCCTGGGACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGAtc                                                                                                                                >  1:391659/1‑150 (MQ=255)
                       aaTTCACCGAGGCGATGCGCTGAATCGTCATCGTCAGACCAGCCTGGGACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGAtc                                                                                                                                >  1:385313/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)
                                   cGATGCGCTGAATCGTCATCGTCAGACCAGCCTGGGACAACAACGTAGCGGATCCACACCGTCACATGTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTCCGGTGGTTGGACACTCCAAAGAGATTTTGGTGGATCTCGTCGTTGATc                                                                                                                    >  2:326317/1‑150 (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)
                                                   catcgtcaGACCAGCCTGGGACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCa                                                                                                    >  2:1547422/1‑150 (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:713062/150‑1 (MQ=255)
                                                             ccAGCCTGGGACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTc                                                                                          >  1:1696172/1‑150 (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)
                                                                            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:1210859/1‑150 (MQ=255)
                                                                                             ccTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGTTCATCAATCACCGGAt                                                          >  1:574711/1‑150 (MQ=255)
                                                                                              cTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGATCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGTTCATCAATCACCGGATc                                                         <  2:1138065/150‑1 (MQ=255)
                                                                                                   caTTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGTTCATCAATCACCGGATCGTAac                                                    <  1:900879/150‑1 (MQ=255)
                                                                                                    aTTTTTGTTCGCCAGATATTTAGCGAACTCCATCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGTTCATCAATCACCGGATCGTAacg                                                   <  1:81077/150‑1 (MQ=255)
                                                                                                       tttGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGTTCATCAATCACCGGATCGTAacgacg                                                <  1:700657/150‑1 (MQ=255)
                                                                                                          gTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGTTCATCAATCACCGGATCGTAACGACGAAc                                             >  2:1124899/1‑150 (MQ=255)
                                                                                                           ttCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGTTCATCAATCACCGGATCGTAACGACGAACa                                            <  2:416276/150‑1 (MQ=255)
                                                                                                            tCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGTTCATCAATCACCGGATCGTAACGCCGAACaa                                           >  2:400472/1‑150 (MQ=255)
                                                                                                              gCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTAATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGTTCATCAATCACCGGATCGTAACGACTAACaaaa                                         >  1:1624552/1‑150 (MQ=255)
                                                                                                               ccAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGTTCATCAATCACCGGATCGTAACGACGAACAAAAc                                        <  1:612034/150‑1 (MQ=255)
                                                                                                                  gATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGTTCATCAATCACCGGATCGTAACGACGAACAAAACCGt                                     <  2:385313/150‑1 (MQ=255)
                                                                                                                    tatTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTAACAGGTCGGTTACTTCCAGCAGGTCATCAATCACCGGATCGTAACGACGAACAAAACCGTTg                                   >  1:1408032/1‑150 (MQ=255)
                                                                                                                     atTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCCGGGCGGTTACTTCCAGCAGTTCATCAATCACCGGATCGTAACGACGAACAAAACCGTTgg                                  >  2:709446/1‑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)
                                                                                                                                  cAGCGTGCGGTGGTTGGAAACTCCAACCAGTTTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGTTCATCAATCACCGTATCGTAACGACGAACAAAACCGTTGGTGTCCAGACAGGt                     <  2:126371/150‑1 (MQ=255)
                                                                                                                                     cGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGTTCATCAATCACCGGATCGTAACGACGAACAAAACCGTTGGTGTCCAGACAGGTATg                  >  2:376802/1‑150 (MQ=255)
                                                                                                                                     cGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGTTCATCAATCACCGGATCGTAACGACGAACAAAACCGTTGGTGTCCAGACAGGTATg                  <  2:57916/150‑1 (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)
                                                                                                                                                     aCTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGTTCATCAATCACCGGATCGTAACGACGAACAAAACCGTTGGTGTCCAGACAGGTATGAATGCCTTCTTTTTTg  <  1:981435/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CAACGTTGCCCATATCACGGGTAAATTCACCGAGGCGATGCGCTGAATCGTCATCGTCAGACCAGCCTGGGACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGTTCATCAATCACCGGATCGTAACGACGAACAAAACCGTTGGTGTCCAGACAGGTATGAATGCCTTCTTTTTTG  >  NC_000913/950498‑950796

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: