Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 950,648 C→T 5.9% V145I (GTT→ATT) ‡ 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,6480CT 5.9% 378.2 / 6.3 118V145I (GTT→ATT) ‡pflApyruvate formate‑lyase 1‑activating enzyme; [formate‑C‑acetyltransferase 1]‑activating enzyme; PFL activase
Reads supporting (aligned to +/- strand):  ref base C (57/53);  new base T (2/5);  total (59/59)
Fisher's exact test for biased strand distribution p-value = 2.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 120 total aligned reads displayed.

ACGTTGCCCATATCACGGGTAAATTCACCGAGGCGATGCGCTGAATCGTCATCGTCAGACCAGCCTGGGACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGTTCATCAATCACCGGATCGTAACGACGAACAAAACCGTTGGTGTCCAGACAGGTATGAATGCCTTCTTTTTTGC  >  NC_000913/950500‑950797
                                                                                                                                                    |                                                                                                                                                     
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)
       ccATATCACGGGTAAATTCACCGAGGCGATGCGCTGAATCGTCATCGTCAGACCAGCCTGGGACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCa                                                                                                                                               >  2:713062/1‑150 (MQ=255)
       ccATATCACGGGTAAATTCACCGAGGCGATGCGCTGAATCGTCATCGTCAGACCAGCCTGGGACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCa                                                                                                                                               >  1:788598/1‑150 (MQ=255)
        cATATCACGGGTAAATTCACCGAGGCGATGCGCTGAATCGTCATCGTCAGACCAGCCTGGGACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAGTTTCCAACCac                                                                                                                                              <  2:498128/150‑2 (MQ=255)
        cATATCACGGGTAAATTCACCGAGGCGATGCGCTGAATCGTCATCGTCAGACCAGCCTGGGACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAGTTTCCAACCac                                                                                                                                              <  2:1043188/150‑2 (MQ=255)
        cATATCACGGGTAAATTCACCGAGGCGATGCGCTGAATCGTCATCGTCAGACCAGCCTGGGACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAGTTTCCAACCac                                                                                                                                              <  1:698567/150‑2 (MQ=255)
         atatCACGGGTAAATTCACCGATGCGATGCGCTGAATCGTCATCGTCAGACCAGCCTGGGACAACAACGTAGCGGATCCACACCTTCACATTTTTGATCGCCAGATATTTAGCGAACTCCAGCGGGCGGTGGTGGGAGTTTCCAACCacc                                                                                                                                             <  1:151769/150‑3 (MQ=255)
         atatCACGGGTAAATTCACCGAGGCGATGCGCTGAATCGTCATCGTCAGACCAGCCTGGGACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAGTTTCCAACCacc                                                                                                                                             <  1:1128787/150‑3 (MQ=255)
         atatCACGGGTAAATTCACCGAGGCGATGCGCTGAATCGTCATCGTCAGACCAGCCTGGGACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGa                                                                                                                                             >  2:1249697/1‑150 (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:385313/1‑150 (MQ=255)
                     aaTTCACCGAGGCGATGCGCTGAATCGTCATCGTCAGACCAGCCTGGGACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGAtc                                                                                                                                 >  1:391659/1‑150 (MQ=255)
                      aTTCACCGAGGCGATGCGCTGAATCGTCATCGTCAGACCAGCCTGGGACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGAtct                                                                                                                                >  2:529251/1‑150 (MQ=255)
                         cACCGCGGCGATGCGCTGAATCGTCATCGTCAGACCAGCCTGGGACAACTCCGTAGCGGATCCACTCCTTCACATTTTTGTTCGCCAGATATTTAGCGATCGCCAGCGTGCGGTGGTTGGATACTCCAACCAGATTTTGGTGGATCtcgt                                                                                                                             <  1:234285/150‑1 (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)
                                        cTGAATCGTCATCGTCAGACCAGCCTGGGACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATATTGGTGGATCTCGTCGTTCATCTGTTTga                                                                                                              >  1:653741/1‑150 (MQ=255)
                                            atcgtcatcgtcAGACCAGCCTGTGACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATc                                                                                                          <  1:849602/150‑1 (MQ=255)
                                              cgtcatcgtcaGACCAGCCTGGGACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGa                                                                                                        >  1:551988/1‑150 (MQ=255)
                                                 catcgtcaGACCAGCCTGGGACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCa                                                                                                     <  1:1371182/150‑1 (MQ=255)
                                                  atcgtcaGACCAGCCTGGGACAACAACGTAGCGGATCCACACCGTCACATTTTTGTTCGCCAGATATTTAGCGACCTCCAGCGTGCGGTGGTTGGAAACTCAAACCAGATTTAGGTCGATCTCGACGTTCATCTGATTGAGAACGCGCAt                                                                                                    >  1:92192/1‑150 (MQ=255)
                                                    cgtcaGACCAGCCTGGGACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTa                                                                                                  >  2:1024174/1‑150 (MQ=255)
                                                     gtcaGACCAGCCTGGGACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTAc                                                                                                 <  1:1613026/150‑1 (MQ=255)
                                                      tcaGACCAGCCTGGGACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTAcc                                                                                                >  1:1659306/1‑150 (MQ=255)
                                                      tcaGACCAGCCTGGGACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTAcc                                                                                                <  1:641889/150‑1 (MQ=255)
                                                      tcaGACCAGCCTGGGACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTAcc                                                                                                >  1:1477798/1‑150 (MQ=255)
                                                        aGACCAGCCTGGGACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAg                                                                                              <  1:529251/150‑1 (MQ=255)
                                                           ccAGCCTGGGACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTc                                                                                           >  1:1696172/1‑150 (MQ=255)
                                                           ccAGCCTGGGACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTc                                                                                           <  1:713062/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)
                                                                cTGGGACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTa                                                                                      <  1:268317/150‑1 (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                                                                          >  1:1424230/1‑150 (MQ=255)
                                                                            cGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGtt                                                                          >  2:446867/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)
                                                                                       cacaCCTTCACATTTTTGTTCGCCAGATATTTATCGAACGCCAGCGTGAGGATGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGGCGTGCATCTGTGTGAGATCGAGCATTACCAGGTCGGTTCCTTCCAGCAGTTCATCAATCAcc                                                               <  1:41219/150‑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)
                                                                                             ttCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGTTCATCAATCACCGGATCg                                                         >  1:1309659/1‑150 (MQ=255)
                                                                                                 caTTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGTTCATCAATCACCGGATCGTAac                                                     <  1:900879/150‑1 (MQ=255)
                                                                                                  gctggagTTCGCTAAATATCTGGCGAACTCCAGCGGGCGGTGGTTGGAAACTCCCACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTAAAAGGTCGGTTACTTACAGCAGTTCTTCAATCAGCGGATCGTAtcg                                                    >  2:1075627/7‑150 (MQ=255)
                                                                                                  aTTTTTGTTCGCCAGATATTTAGCGAACTCCATCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGTTCATCAATCACCGGATCGTAacg                                                    <  1:81077/150‑1 (MQ=255)
                                                                                                         ttCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGTTCATCAATCACCGGATCGTAACGACGAACa                                             <  2:416276/150‑1 (MQ=255)
                                                                                                          tCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGTTCATCAATCACCGGATCGTAACGCCGAACaa                                            >  2:400472/1‑150 (MQ=255)
                                                                                                            gCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGTTCATCAATCCCCGGATCGTAACGACGAACaaaa                                          >  1:820226/1‑150 (MQ=255)
                                                                                                                gATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGTTCATCAATCACCGGATCGTAACGACGAACAAAACCGt                                      <  2:385313/150‑1 (MQ=255)
                                                                                                                  tatTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGTTCATCAATCACCGGATCGTAACGACGAACAAAACCGTTg                                    >  1:1414576/1‑150 (MQ=255)
                                                                                                                  taattAGCGAACTCCAGCGTGCGGTGCTTGTAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGTTCATCAATAACCGGATCGTAACGACGCACAAAACCGTTg                                    >  2:315079/4‑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                                                 <  1:1201231/131‑1 (MQ=255)
                                                                                                                        gCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGTTCATCAATCACCGGATCGTAacgacg                                                 >  2:1201231/1‑131 (MQ=255)
                                                                                                                        gCGAACTCCAGCGTGCGGTGGTGGGAAACTCCAACCAGATTTTGGTGGATCTCGTGGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGTTCATCAATCACCGGATCGTAACGACGAACAAAACCGTTGGTGTcc                              <  1:871292/150‑1 (MQ=255)
                                                                                                                          gAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGTTCATCAATCACCGGATCGTAACGACGAACAAAACCGTTGGTGTCCAg                            <  1:631934/150‑1 (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: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)
                                                                                                                                             ttGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGTTCATCAATCACCGGATCGTAACGACGAACAAAACCGGTGGTGTCCAGACAGGTATGAATGCcttct         >  1:1240030/1‑150 (MQ=255)
                                                                                                                                               ggAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGTTCATCAATCACCGGATCGTAACGACGAACAAAACCGTTGGTGTCCAGACAGGTATGAATGCCTTCttt       <  1:12180/150‑1 (MQ=255)
                                                                                                                                                gAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCAACTGTTTGAGATAGAGCATTACCAGGTCTGTTACTACCAGCAGTTCAACAATAACCTGATCGTAACGCAGAACAAAACCGTACGTGTGCAGACAGGAAAGAATGCattcttat      >  1:594317/1‑148 (MQ=255)
                                                                                                                                                   aCTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGTTCATCAATCACCGGATCGTAACGACGAACAAAACCGTTGGTGTCCAGACAGGTATGAATGCCTTCTTTTTTg   <  1:981435/150‑1 (MQ=255)
                                                                                                                                                    cTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGTTCATCAATCACCGGATCGTAACGACGAACAAAACCGTTGGTGTCCAGACAGGTATGAATGCCTTCTTTTTTGc  >  2:1427731/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
ACGTTGCCCATATCACGGGTAAATTCACCGAGGCGATGCGCTGAATCGTCATCGTCAGACCAGCCTGGGACAACAACGTAGCGGATCCACACCTTCACATTTTTGTTCGCCAGATATTTAGCGAACTCCAGCGTGCGGTGGTTGGAAACTCCAACCAGATTTTGGTGGATCTCGTCGTTCATCTGTTTGAGATCGAGCATTACCAGGTCGGTTACTTCCAGCAGTTCATCAATCACCGGATCGTAACGACGAACAAAACCGTTGGTGTCCAGACAGGTATGAATGCCTTCTTTTTTGC  >  NC_000913/950500‑950797

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: